Engineering Papers⌕ Search

Engineering topics

Kwok, R.

Publications and source records attributed to Kwok, R..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

The Role of Cyclone Activity in Snow Accumulation on Arctic Sea Ice

Identifying the mechanisms controlling the timing and magnitude of snow accumulation on sea ice is crucial for understanding snow’s net effect on the surface energy budget and sea-ice mass balance. Here, we analyze the role of cyclone activity on the seasonal buildup of snow on Arctic sea ice using model, satellite, and in situ data over 1979–2016. On average, 44% of the variability in monthly snow accumulation was controlled by cyclone snowfall and 29% by sea-ice freeze-up. However, there were strong spatio-temporal differences. Cyclone snowfall comprised ~50% of total snowfall in the Pacific compared to 83% in the Atlantic. While cyclones are stronger in the Atlantic, Pacific snow accumulation is more sensitive to cyclone strength. These findings highlight the heterogeneity in atmosphere-snow-ice interactions across the Arctic, and emphasize the need to scrutinize mechanisms governing cyclone activity to better understand their effects on the Arctic snow-ice system with anthropogenic warming.

Webster, M. A.↗
NASA NTRS2017

Potential Basin-Scale Estimates of Arctic Snow Depth with Sea Ice Freeboards from CryoSat-2 and ICESat-2: An Exploratory Analysis

The potential of deriving snow depth estimates using differences in freeboard heights from CryoSat-2 (CS-2) and ICESat-2 (IS-2) is examined. In our analysis, we use lidar freeboard from the Airborne Topographic Mapper (ATM) on Operation IceBridge (OIB) as proxy of IS-2 total (snow plus ice) freeboard. Snow depths are estimates from the OIB snow radar. Differences in height between the total (ATM) and ice (CS-2) freeboards are related to snow depth by the refractive index of the snow layer, which is dependent on snow density. For two years (2014 and 2015), regression of the ATM and CS-2 freeboard differences against OIB snow depth gives correlations of approximately 0.80, estimated eta (sub s) of 1.21, and standard errors of 8 centimeters. The resulting refractive index, eta (sub s), can be compared to that expected of the Arctic snow cover in early spring (1.25 plus or minus 0.05). The expected biases and variability in the regression analysis are discussed. Results suggest that snow depth can be estimated from the freeboard differences. The benefits of adjusting the orbit of CS-2 for providing more optimized overlaps between IS-2 and CS-2 are considered.

Sea ice↗
NASA NTRS2014

Profiling Sea Ice with a Multiple Altimeter Beam Experimental Lidar (MABEL)

The sole instrument on the upcoming ICESat-2 altimetry mission is a micropulse lidar that measures the time-of-flight of individual photons from laser pulses transmitted at 532 nm. Prior to launch, MABEL serves as an airborne implementation for testing and development. In this paper, we provide a first examination of MABEL data acquired on two flights over sea ice in April 2012: one north of the Arctic coast of Greenland, and the other in the East Greenland Sea.We investigate the phenomenology of photon distributions in the sea ice returns. An approach to locate the surface and estimate its elevation in the distributions is described, and its achievable precision assessed. Retrieved surface elevations over relatively flat leads in the ice cover suggest that precisions of several centimeters are attainable. Restricting the width of the elevation window used in the surface analysis can mitigate potential biases in the elevation estimates due to subsurface returns at 532 nm. Comparisons of nearly coincident elevation profiles from MABEL with those acquired by an analog lidar show good agreement.Discrimination of ice and open water, a crucial step in the determination of sea ice free board and the estimation of ice thickness, is facilitated by contrasts in the observed signal background photon statistics. Future flight lines will sample a broader range of seasonal ice conditions for further evaluation of the year-round profiling capabilities and limitations of the MABEL instrument.

Lidar↗
NASA NTRS2012

Arctic Sea Ice Freeboard from Icebridge Acquisitions in 2009: Estimates and Comparisons with ICEsat

During the spring of 2009, the Airborne Topographic Mapper (ATM) system on the IceBridge mission acquired cross-basin surveys of surface elevations of Arctic sea ice. In this paper, the total freeboard derived from four 2000 km transects are examined and compared with those from the 2009 ICESat campaign. Total freeboard, the sum of the snow and ice freeboards, is the elevation of the air-snow interface above the local sea surface. Prior to freeboard retrieval, signal dependent range biases are corrected. With data from a near co-incident outbound and return track on 21 April, we show that our estimates of the freeboard are repeatable to within 4 cm but dependent locally on the density and quality of sea surface references. Overall difference between the ATM and ICESat freeboards for the four transects is 0.7 (8.5) cm (quantity in bracket is standard deviation), with a correlation of 0.78 between the data sets of one hundred seventy-eight 50 km averages. This establishes a level of confidence in the use of ATM freeboards to provide regional samplings that are consistent with ICESat. In early April, mean freeboards are 41 cm and 55 cm over first year and multiyear sea ice (MYI), respectively. Regionally, the lowest mean ice freeboard (28 cm) is seen on 5 April where the flight track sampled the large expanse of seasonal ice in the western Arctic. The highest mean freeboard (71 cm) is seen in the multiyear ice just west of Ellesmere Island from 21 April. The relatively large unmodeled variability of the residual sea surface resolved by ATM elevations is discussed.

Airborne Topographic Mapper (ATM)↗
NASA NTRS2011

Airborne Surveys of Snow Depth over Arctic Sea Ice

During the spring of 2009, an ultrawideband microwave radar was deployed as part of Operation IceBridge to provide the first cross-basin surveys of snow thickness over Arctic sea ice. In this paper, we analyze data from three approx 2000 km transects to examine detection issues, the limitations of the current instrument, and the regional variability of the retrieved snow depth. Snow depth is the vertical distance between the air \snow and snow-ice interfaces detected in the radar echograms. Under ideal conditions, the per echogram uncertainty in snow depth retrieval is approx 4 - 5 cm. The finite range resolution of the radar (approx 5 cm) and the relative amplitude of backscatter from the two interfaces limit the direct retrieval of snow depths much below approx 8 cm. Well-defined interfaces are observed over only relatively smooth surfaces within the radar footprint of approx 6.5 m. Sampling is thus restricted to undeformed, level ice. In early April, mean snow depths are 28.5 +/- 16.6 cm and 41.0 +/- 22.2 cm over first-year and multiyear sea ice (MYI), respectively. Regionally, snow thickness is thinner and quite uniform over the large expanse of seasonal ice in the Beaufort Sea, and gets progressively thicker toward the MYI cover north of Ellesmere Island, Greenland, and the Fram Strait. Snow depth over MYI is comparable to that reported in the climatology by Warren et al. Ongoing improvements to the radar system and the utility of these snow depth measurements are discussed.

Kwok, R.↗
NASA NTRS2007

Ross Sea Polynyas: Response of Ice Concentration Retrievals to Large Areas of Thin Ice

For a 3-month period between May and July of 2005, we examine the response of the Advanced Microwave Scanning Radiometer (AMSR-E) Enhanced NASA Team 2 (NT2) and AMSR-E Bootstrap (ABA) ice concentration algorithms to large areas of thin ice of the Ross Sea polynyas. Coincident Envisat Synthetic Aperture Radar (SAR) coverage of the region during this period offers a detailed look at the development of the polynyas within several hundred kilometers of the ice front. The high-resolution imagery and derived ice motion fields show bands of polynya ice, covering up to approximately 105 km(sup 2) of the Ross Sea, that are associated with wind-forced advection. In this study, ice thickness from AMSR-E 36 GHz polarization information serves as the basis for examination of the response. The quality of the thickness of newly formed sea ice (<10 cm) from AMSR-E is first assessed with thickness estimates derived from ice surface temperatures from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. The effect of large areas of thin ice in lowering the ice concentration estimates from both NT2/ABA approaches is clearly demonstrated. Results show relatively robust relationships between retrieved ice concentrations and thin ice thickness estimates that differ between the two algorithms. These relationships define the approximate spatial coincidence of ice concentration and thickness isopleths. Using the 83% (ABA) and 91% (NT2) isopleths as polynya boundaries, we show that the computed coverage compares well with that using the estimated 10-cm thickness contour. The thin ice response characterized here suggests that in regions with polynyas, the retrieval results could be used to provide useful geophysical information, namely thickness and coverage.

algorithms↗
NASA NTRS2007

A New Sea Ice Model: Pancakes to Metrics

We are developing, testing, and validating a new sea ice dynamics model that treats the ice cover as an elastic/decohesive material in the permanent pack and includes the correct frazil/pancake behavior in the marginal zone. Two salient features of present ice dynamics models are that they do not: 1) reproduce the oriented fracture patterns of openings and closings in the pack ice, and 2) accurately model the effects of frazil/pancake ice formation in the ice margin. These poorly modeled areas account for a substantial portion of the ice growth, turbulent heat flux to the atmosphere, salt flux to the ocean, and energy dissipation due to slippage, ridging, and rafting, in the Arctic. Existing sea ice models have shown limited success in predicting the degree to which a lead will open for prescribed or observed forcing conditions. An important aspect of the new model we are developing is that the existence of cracks, along with their orientation, opening, and closing, is predicted. To put this effort in perspective a short history of ice dynamics modeling and data collection is presented. The RGPS data set is used to validate the model. As part of the testing and validation of the model, we are working on a new metric for comparing linear features (leads and ridges) in the data and model to be used in data assimilation for this model. The model framework is presented as well as some results showing the creation and development of leads in a simulation of ice dynamics in the Beaufort Sea. Other presentations by the authors will show other results from this effort.

Coon, M.↗
NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗
NASA NTRS2006

Contrasts in Sea Ice Formation and Production in the Arctic Seasonal and Perennial Ice Zones

Four years (1997-2000) of RADARSAT Geophysical Processor System (RGPS) data are used to contrast the sea ice deformation and production regionally, and in the seasonal (SIZ) and perennial (PIZ) ice zones. Ice production is of seasonal ice in openings during the winter. 3-day estimates of these quantities are provided within Lagrangian elements initially 10 km on a side. A distinct seasonal cycle is seen in both zones with these estimates highest in the late fall and with seasonal minimums in the mid-winter. Regional divergence over the winter could be up to 30%. Spatially, the highest deformation is in the SIZ north of coastal Alaska. Both ice deformation and production are higher in the SIZ: deformation-related ice production in the SIZ (approx.0.5 m) is 1.5-2.3 times that of the PIZ (approx.0.3 m) - this is connected to ice strength and thickness. Atmospheric forcing and boundary layer structure contribute to only the seasonal and interannual variability. Seasonal ice growth in ice fractures accounts for approx.25-40% of the total ice production of the Arctic Ocean. By itself, this deformation-ice production relationship could be considered a negative feedback when thickness is perturbed. However, the overall effect on ice production in the face of increasing seasonal and thinner/weaker ice coverage could be modified by: local destabilization of the water column promoting overturning of warmer water due to increased brine rejection; and, the upwelling of the pynocline associated with increased occurrence of large shear motion in sea ice.

Kwok, R.↗
NASA NTRS2006

Upwelling of Arctic pycnocline associated with shear motion of sea ice

High-resolution radar imagery shows that the dynamic response of winter sea ice to gradients in large-scale surface wind stress is often localized along quasi-linear fractures hundreds of kilometers long. Relative shearing motion across these narrow fractures can exceed 10 cm/s. In one event recorded during the drift of the SHEBA ice camp, we observed an intense zone of pycnocline upwelling (approx.14 m) associated with significant shear motion near the camp, while upward turbulent heat flux in the ocean boundary layer reached nearly 400 W/sq m, an order of magnitude greater than at any other time during the year-long drift. We attribute the upwelling to Ekman pumping associated with concentrated ice shear. Over the entire Arctic Ocean sea ice cover, this process could be responsible for significant heat exchange between the cold surface layer and warmer subsurface water at the ubiquitous fractures resulting from large-scale atmosphere-ice interactions.

McPhee, M. G.↗
NASA NTRS2004

Annual Cycles of Multiyear Sea Ice Coverage of the Arctic Ocean: 1999-2003

For the years 1999-2003, we estimate the time-varying perennial ice zone (PIZ) coverage and construct the annual cycles of multiyear (MY, including second year) ice coverage of the Arctic Ocean using QuikSCAT backscatter, MY fractions from RADARSAT, and the record of ice export from satellite passive microwave observations. An area balance approach extends the winter MY coverage from QuikSCAT to the remainder of the year. From these estimates, the coverage of MY ice at the beginning of each year is 3774 x 10(exp 3) sq km (2000), 3896 x 10(exp 3) sq km (2001), 4475 x 10(exp 3) sq km (2002), and 4122 x 10(exp 3) sq km (2003). Uncertainties in coverage are approx.150 x 10(exp 3) sq km. In the mean, on 1 January, MY ice covers approx.60% of the Arctic Ocean. Ice export reduces this coverage to approx.55% by 1 May. From the multiple annual cycles, the area of first-year (FY) ice that survives the intervening summers are 1192 x 10(exp 3) sq km (2000), 1509 x 10(exp 3) sq km (2001), and 582 x 10(exp 3) sq km (2002). In order for the MY coverage to remain constant from year to year, these replenishment areas must balance the overall area export and melt during the summer. The effect of the record minimum in Arctic sea ice area during the summer of 2002 is seen in the lowest area of surviving FY ice of the three summers. In addition to the spatial coverage, the location of the PIZ is important. One consequence of the unusual location of the PIZ at the end of the summer of 2002 is the preconditioning for enhanced export of MY ice into the Barents and Kara seas. Differences between the minimums in summer sea ice coverage from our estimates and passive microwave observations are discussed.

Arctic Ocean↗
NASA NTRS2004

Fram Strait sea ice outflow

We summarize 24 years (1978??2) of ice export estimates and examine, over a 9-year record, the associated variability in the time-varying upward-looking sonar (ULS) thickness distributions of the Fram Strait.

Fram strait ice export ice thickness↗
NASA NTRS2004

Fram Strait sea ice outflow

We summarize 24 years of ice export estimates and examine, over a 9-year record, the associated variability in the time-varying upward-looking sonar (ULS) thickness distributions of the Fram Strait. A more thorough assessment of the PMW (passive microwave) ice motion with 5 years of synthetic aperture radar (SAR)observations shows the uncertainties to be consistent with that found by Kwok and Rothrock [1999], giving greater confidence to the record of ice flux calculations.

Fram Strait ice export ice thickness↗
NASA NTRS2003

Sub-daily sea ice motion and deformation from RADARSAT observations

We find a persistent level of oscillatory sea ice motion and deformation, superimposed on the large-scale wind-driven field, in May 2002 (spring) and February 2003 (mid-winter), in the high Arctic over a region centered at approx.(85degreeN, 135degreeW). At this latitude, the RADARSAT wide-swath SAR coverage provides 4??equential observations every day, for ice motion retrieval, with a sampling interval at the orbital period of approx. 101 minutes.

ice motion ice mechanics oceanography↗
NASA NTRS2003

Sub-daily sea ice and deformation from RADARSAT observations

Periodic correlations in ice motion and deformation can be seen in lenght scales from 10 km and above, and suggest a 12 hour oscillation that is more likely associated with inertial rather than tidal frequencies.

oceanography ice mechanicscoriolis effects↗