Engineering topics
Kwok, Ron
Publications and source records attributed to Kwok, Ron.
Evaluation of Arctic Sea Ice Thickness Simulated by Arctic Ocean Model Intercomparison Project Models
Six Arctic Ocean Model Intercomparison Project model simulations are compared with estimates of sea ice thickness derived from pan-Arctic satellite freeboard measurements (2004-2008); airborne electromagnetic measurements (2001-2009); ice draft data from moored instruments in Fram Strait, the Greenland Sea, and the Beaufort Sea (1992-2008) and from submarines (1975-2000); and drill hole data from the Arctic basin, Laptev, and East Siberian marginal seas (1982-1986) and coastal stations (1998-2009). Despite an assessment of six models that differ in numerical methods, resolution, domain, forcing, and boundary conditions, the models generally overestimate the thickness of measured ice thinner than approximately 2 mand underestimate the thickness of ice measured thicker than about approximately 2m. In the regions of flat immobile landfast ice (shallow Siberian Seas with depths less than 25-30 m), the models generally overestimate both the total observed sea ice thickness and rates of September and October ice growth from observations by more than 4 times and more than one standard deviation, respectively. The models do not reproduce conditions of fast ice formation and growth. Instead, the modeled fast ice is replaced with pack ice which drifts, generating ridges of increasing ice thickness, in addition to thermodynamic ice growth. Considering all observational data sets, the better correlations and smaller differences from observations are from the Estimating the Circulation and Climate of the Ocean, Phase II and Pan-Arctic Ice Ocean Modeling and Assimilation System models.
Laser Altimetry Sampling Strategies over Sea Ice
With the conclusion of the science phase of the Ice, Cloud and land Elevation Satellite (ICESat) mission in late 2009, and the planned launch of ICESat-2 in late 2015, NASA has recently established the IceBridge program to provide continuity between missions. A major goal of IceBridge is to obtain a sea-ice thickness time series via airborne surveys over the Arctic and Southern Oceans. Typically two laser altimeters, the Airborne Topographic Mapper (ATM) and the Land, Vegetation and Ice Sensor (LVIS), are utilized during IceBridge flights. Using laser altimetry simulations of conventional analogue systems such as ICESat, LVIS and ATM, with the multi-beam system proposed for ICESat-2, we investigate differences in measurements gathered at varying spatial resolutions and the impact on sea-ice freeboard. We assess the ability of each system to reproduce the elevation distributions of two seaice models and discuss potential biases in lead detection and sea-surface elevation, arising from variable footprint size and spacing. The conventional systems accurately reproduce mean freeboard over 25km length scales, while ICESat-2 offers considerable improvements over its predecessor ICESat. In particular, its dense along-track sampling of the surface will allow flexibility in the algorithmic approaches taken to optimize the signal-to-noise ratio for accurate and precise freeboard retrieval.
Evaluation of Arctic Sea Ice Thickness Simulated by AOMIP Models
We compare results from six AOMIP model simulations with estimates of sea ice thickness obtained from ICESat, moored and submarine-based upward looking sensors, airborne electromagnetic measurements and drill holes. Our goal is to find patterns of model performance to guide model improvement. The satellite data is pan-arctic from 2004-2008, ice-draft data is from moored instruments in Fram Strait, the Greenland Sea and the Beaufort Sea from 1992-2008 and from submarines from 1975-2000. The drill hole data are from the Laptev and East Siberian marginal seas from 1982-1986 and from coastal stations from 1998-2009. While there are important caveats when comparing modeled results with measurements from different platforms and time periods such as these, the models agree well with moored ULS data. In general, the AOMIP models underestimate the thickness of measured ice thicker than about 2 m and overestimate thickness of ice thinner than 2 m. The simulated results are poor over the fast ice and marginal seas of the Siberian shelves. Averaging over all observational data sets, the better correlations and smaller differences from observed thickness are from the ECCO2 and UW models.
Baffin Bay Ice Drift and Export: 2002-2007
Multiyear estimates of sea ice drift in Baffin Bay and Davis Strait are derived for the first time from the 89 GHz channel of the AMSR-E instrument. Uncertainties in the drift estimates, assessed with Envisat ice motion, are approximately 2-3 km/day. A persistent atmospheric trough, between the coast of Greenland and Baffin Island, drives the prevailing southward drift pattern with average daily displacements in excess of 18-20 km during winter. Over the 5-year record, the ice export ranges between 360 and 675 x 10(exp 3) km(exp 2), with an average of 530 x 10(exp 3) km(exp 2). Sea ice area inflow from the Nares Strait, Lancaster Sound and Jones Sound potentially contribute up to a third of the net area outflow while ice production at the North Water Polynya contributes the balance. Rough estimates of annual volume export give approximately 500-800 km(exp 3). Comparatively, these are approximately 70% and approximately 30% of the annual area and Strait.
On large outflows of Arctic sea ice into the Barents Sea
Winter outflows of Arctic sea ice into the Barents Sea are estimated using a 10-year record of satellite ice motion and thickness. The mean winter volume export through the Svalbard/Franz Josef Land passage is 40 km3, and ranges from -280 km3 to 340 km3. A large outflow in 2003 is preconditioned by an unusually high concentration of thick perennial ice over the Nansen Basin at the end of the 2002 summer. With a deep atmospheric low situated over the eastern Barents Sea in winter, the result is an increased export of Arctic ice. The Oct-Mar ice area flux, at 110 x 10 to the third power km3, is not only unusual in magnitude but also remarkable in that >70% of the area is multiyear ice; the ice volume flux at~340 km3 is almost one-fifth of the ice flux through the Fram Strait. Another large outflow of Arctic sea ice through this passage, comparable to that in 2003, is found in 1996. This southward flux of sea ice represents one of two major sources of freshwater in the Barents Sea; the other is the eastward flux of water via the Norwegian Coastal Current. The possible consequences of variable freshwater input on the Barents Sea hydrography and its impact on transformation of Atlantic Water en route to the Arctic Ocean are examined with a 25-year coupled ice-ocean model.
ICESat Observations of Arctic Sea Ice: A First Look
Analysis of near-coincident ICESat and RADARSAT imagery shows that the retrieved elevations from the laser altimeter are sensitive to new openings (containing thin ice or open water) in the sea ice cover as well as to surface relief of old and first-year ice. The precision of the elevation estimates, measured over relatively flat sea ice, is approx. 2 cm. Using the thickness of thin-ice in recent openings to estimate sea level references, we obtain the sea-ice freeboard along the altimeter tracks. This step is necessitated by the large uncertainties in the sea surface topography compared to that required for accurate determination of freeboard. Unknown snow depth introduces the largest uncertainty in the conversion of freeboard to ice thickness. Surface roughness is also derived, for the first time, from the variability of successive elevation estimates along the altimeter track. Overall, these ICESat measurements provide an unprecedented view of the Arctic Ocean ice cover at length scales at and above the spatial dimension of the altimeter footprint of approx. 70 m.
ICESat Observations of Arctic Sea Ice: A First Look
Analysis of near-coincident ICESat and RADARSAT imagery shows that the retrieved elevations from the laser altimeter are sensitive to new openings (containing thin ice or open water) in the sea ice cover as well as to surface relief of old and first-year ice. The precision of the elevation estimates, measured over relatively flat sea ice, is approx. 2 cm Using the thickness of thin-ice in recent openings to estimate sea level references, we obtain the sea-ice free-board along the altimeter tracks. This step is necessitated by the large uncertainties in the time-varying sea surface topography compared to that required for accurate determination of free-board. Unknown snow depth introduces the largest uncertainty in the conversion of free-board to ice thickness. Surface roughness is also derived, for the first time, from the variability of successive elevation estimates along the altimeter track Overall, these ICESat measurements provide an unprecedented view of the Arctic Ocean ice cover at length scales at and above the spatial dimension of the altimeter footprint.
Sea ice investigations from Seasat to the present
In this paper, we provide a short review of sea ice investigations starting from Seasat. We focus particularly on the detailed and quantitative measurements of the sea ice motion field, which were some of the earliest results from Seasat and have sugsequently been shown to be of critical value to the derivation of several key climatically important variables. Other key investigations discussed include examination of the seasonal melt cycle, ice extent and concentration, and estimates of thickness from the proxy measurements of ice type and age and more directly from freeboard. We end with a brief discussion on how thse measurements might be improved in the future.
Spatial Patterns of Variability in Antarctic Surface Temperature: Connections to the Southern Hemisphere Annular Mode and the Southern Oscillation
The 17-year (1982-1998) trend in surface temperature shows a general cooling over the Antarctic continent, warming of the sea ice zone, with moderate changes over the oceans. Warming of the peripheral seas is associated with negative trends in the regional sea ice extent. Effects of the Southern Hemisphere Annular Mode (SAM) and the extrapolar Southern Oscillation (SO) on surface temperature are quantified through regression analysis. Positive polarities of the SAM are associated with cold anomalies over most of Antarctica, with the most notable exception of the Antarctic Peninsula. Positive temperature anomalies and ice edge retreat in the Pacific sector are associated with El Nino episodes. Over the past two decades, the drift towards high polarity in the SAM and negative polarity in the SO indices couple to produce a spatial pattern with warmer temperatures in the Antarctic Peninsula and peripheral seas, and cooler temperatures over much of East Antarctica.
Ice Flow in the North East Greenland Ice Stream
Early observations with ERS-1 SAR image data revealed a large ice stream in North East Greenland (Fahnestock 1993). The ice stream has a number of the characteristics of the more closely studied ice streams in Antarctica, including its large size and gross geometry. The onset of rapid flow close to the ice divide and the evolution of its flow pattern, however, make this ice stream unique. These features can be seen in the balance velocities for the ice stream (Joughin 1997) and its outlets. The ice stream is identifiable for more than 700 km, making it much longer than any other flow feature in Greenland. Our research goals are to gain a greater understanding of the ice flow in the northeast Greenland ice stream and its outlet glaciers in order to assess their impact on the past, present, and future mass balance of the ice sheet. We will accomplish these goals using a combination of remotely sensed data and ice sheet models. We are using satellite radar interferometry data to produce a complete maps of velocity and topography over the entire ice stream. We are in the process of developing methods to use these data in conjunction with existing ice sheet models similar to those that have been used to improve understanding of the mechanics of flow in Antarctic ice streams.
Variability of Fram Strait Ice Flux and North Atlantic Oscillation
An important term in the mass balance of the Arctic Ocean sea ice is the ice export. We estimated the winter sea ice export through the Fram Strait using ice motion from satellite passive microwave data and ice thickness data from moored upward looking sonars. The average winter area flux over the 18-year record (1978-1996) is 670,000 square km, approximately 7% of the area of the Arctic Ocean. The winter area flux ranges from a minimum of 450,000 sq. km in 1984 to a maximum of 906,000 sq km in 1995. The daily, monthly and interannual variabilities of the ice area flux are high. There is an upward trend in the ice area flux over the 18-year record. The average winter volume flux over the winters of October 1990 through May 1995 is 1745 cubic km ranging from a low of 1375 cubic km in 1990 to a high of 2791 cubic km in 1994. The sea-level pressure gradient across the Fram Strait explains more than 80% of the variance in the ice flux over the 18-year record. We use the coefficients from the regression of the time-series of area flux versus pressure gradient across the Fram Strait and ice thickness data to estimate the summer area and volume flux. The average 12-month area flux and volume flux are 919,000 sq km and 2366 cubic km. We find a significant correlation (R =0.86) between the area flux and positive phases of the North Atlantic Oscillation (NAO) index over the months of December through March. Correlation between our six years of volume flux estimates and the NAO index gives R =0.56. During the high NAO years, a more intense Icelandic low increases the gradient in the sea-level pressure by almost 1 mbar across the Fram Strait thus increasing the atmospheric forcing on ice transport. Correlation is reduced during the negative NAO years because of decreased dominance of this large-scale atmospheric pattern on the sea-level pressure gradient across the Fram Strait. Additional information is contained in the original.
Interferometric estimation of ice sheet motion and topography
With ERS-1/2 satellite radar interferometry, it is possible to make measurements of glacier motion with high accuracy and fine spatial resolution. Interferometric techniques were applied to map velocity and topography for several outlet glaciers in Greenland. For the Humboldt and Petermann glaciers, data from several adjacent tracks were combined to make a wide-area map that includes the enhanced flow regions of both glaciers. The discharge flux of the Petermann glacier upstream of the grounding line was estimated, thereby establishing the potential use of ERS-1/2 interferometric data for monitoring ice-sheet discharge. Interferograms collected along a single track are sensitive to only one component of motion. By utilizing data from ascending and descending passes and by making a surface-parallel flow assumption, it is possible to measure the full three-dimensional vector flow field. The application of this technique for an area on the Ryder glacier is demonstrated. Finally, ERS-1/2 interferograms were used to observe a mini-surge on the Ryder glacier that occurred in autumn of 1995.
Interferometric Study of the Ice Stream in Interior Northeast Greenland
Initial work in Greenland with ERS-1 synthetic aperture radar (SAR) imagery led to the discovery of a large ice stream which carries much of the ice discharge from the northeast quadrant of the ice sheet. The ice stream was identified based on the morphology of the surface, which indicated the influence of localized rapid ice flow. Continued investigations using repeat-pass interferometry from ERS-1 and ERS-2 allowed the character of the rapid motion to be determined, mapping the onset region in the deep interior of the ice sheet as well as more rapid motion downstream. Markers were placed on the ice sheet and are being tracked with the global positioning system (GPS) survey to provide tie-points for the interferometry, which will allow the absolute velocities to be determined in this very remote area.
Balance Velocities for the Greenland Ice Sheet
We present a map of balance velocities for the Greenland ice sheet. The resolution of the DEM, which was derived from the radar altimetery data, yields far greater detail than in earlier balance velocity estimates for Greenland.
Balance Velocities of the Greenland Ice Sheet
We present a map of balance velocities for the Greenland ice sheet. The resolution of the underlying DEM, which was derived primarily from radar altimetry data, yields far greater detail than earlier balance velocity estimates for Greenland. The velocity contours reveal in striking detail the location of an ice stream in northeastern Greenland, which was only recently discovered using satellite imagery. Enhanced flow associated with all of the major outlets is clearly visible, although small errors in the source data result in less accurate estimates of the absolute flow speeds. Nevertheless, the balance map is useful for ice-sheet modelling, mass balance studies, and field planning.
Estimation of Ice Sheet Motion Using Satellite Radar Interferometry
We examine the sources of error in interferometrically-derived ice-motion measurements, including those errors due to inaccurate estimates of the interferometric baseline.