More Like this

1. Improvements in one-dimensional grounding-line parameterizations in an ice-sheet model with lateral variations (PSUICE3D v2.1)

   Pollard, D. ; DeConto, R. ( January 2020 , Geoscientific model development)

   null (Ed.)

   The use of a boundary-layer parameterization of buttressing and ice flux across grounding lines in a two dimensional ice-sheet model is improved by allowing general orientations of the grounding line. This and another modification to the model’s grounding-line parameterization are assessed in three settings: rectangular fjord-like domains – the third Marine Ice Sheet Model Intercomparison Project (MISMIPC) and Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d) – and future simulations of West Antarctic ice retreat under Representative Concentration Pathway (RCP) 8.5-based climates. The new modifications are found to have significant effects on the fjord-like results, which are now within the envelopes of other models in the MISMIP+ and MISMIP3d intercomparisons. In contrast, the modifications have little effect on West Antarctic retreat, presumably because dynamics in the wider major Antarctic basins are adequately represented by the model’s previous simpler one-dimensional formulation. As future grounding lines retreat across very deep bedrock topography in the West Antarctic simulations, buttressing is weak and deviatoric stress measures exceed the ice yield stress, implying that structural failure at these grounding lines would occur. We suggest that these grounding-line quantities should be examined in similar projections by other ice models to better assess the potential for future structural failure. 

   more » « less
2. A Generalized Interpolation Material Point Method for Shallow Ice Shelves. 1: Shallow Shelf Approximation and Ice Thickness Evolution

   https://doi.org/10.1029/2020MS002277  

   Huth, Alex ; Duddu, Ravindra ; Smith, Ben ( August 2021 , Journal of Advances in Modeling Earth Systems)

   We develop a generalized interpolation material point method (GIMPM) for the shallow shelf approximation (SSA) of ice flow. The GIMPM, which can be viewed as a particle version of the finite element method, is used here to solve the shallow shelf approximations of the momentum balance and ice thickness evolution equations. We introduce novel numerical schemes for particle splitting and integration at domain boundaries to accurately simulate the spreading of an ice shelf. The advantages of the proposed GIMPM‐SSA framework include efficient advection of history or internal state variables without diffusion errors, automated tracking of the ice front and grounding line at sub‐element scales, and a weak formulation based on well‐established conventions of the finite element method with minimal additional computational cost. We demonstrate the numerical accuracy and stability of the GIMPM using 1‐D and 2‐D benchmark examples. We also compare the accuracy of the GIMPM with the standard material point method (sMPM) and a reweighted form of the sMPM. We find that the grid‐crossing error is very severe for SSA simulations with the sMPM, whereas the GIMPM successfully mitigates this error. While the grid‐crossing error can be reasonably reduced in the sMPM by implementing a simple material point reweighting scheme, this approach it not as accurate as the GIMPM. Thus, we illustrate that the GIMPM‐SSA framework is viable for the simulation of ice sheet‐shelf evolution and enables boundary tracking and error‐free advection of history or state variables, such as ice thickness or damage.

    

   more » « less
3. Wind- and Wave-Driven Reynolds Stress and Velocity Shear in the Upper Ocean for Idealized Misaligned Wind-Wave Conditions

   https://doi.org/10.1175/JPO-D-20-0157.1  

   Wang, Dong ; Kukulka, Tobias ( February 2021 , Journal of Physical Oceanography)

   null (Ed.)

   Abstract This study investigates the dynamics of velocity shear and Reynolds stress in the ocean surface boundary layer for idealized misaligned wind and wave fields using a large-eddy simulation (LES) model based on the Craik–Leibovich equations, which captures Langmuir turbulence (LT). To focus on the role of LT, the LES experiments omit the Coriolis force, which obscures a stress–current-relation analysis. Furthermore, a vertically uniform body force is imposed so that the volume-averaged Eulerian flow does not accelerate but is steady. All simulations are first spun-up without wind-wave misalignment to reach a fully developed stationary turbulent state. Then, a crosswind Stokes drift profile is abruptly imposed, which drives crosswind stresses and associated crosswind currents without generating volume-averaged crosswind currents. The flow evolves to a new stationary state, in which the crosswind Reynolds stress vanishes while the crosswind Eulerian shear and Stokes drift shear are still present, yielding a misalignment between Reynolds stress and Lagrangian shear (sum of Eulerian current and Stokes drift). A Reynolds stress budgets analysis reveals a balance between stress production and velocity–pressure gradient terms (VPG) that encloses crosswind Eulerian shear, demonstrating a complex relation between shear and stress. In addition, the misalignment between Reynolds stress and Eulerian shear generates a horizontal turbulent momentum flux (due to correlations of along-wind and crosswind turbulent velocities) that can be important in producing Reynolds stress (due to correlations of horizontal and vertical turbulent velocities). Thus, details of the Reynolds stress production by Eulerian and Stokes drift shear may be critical for driving upper-ocean currents and for accurate turbulence parameterizations in misaligned wind-wave conditions. 

   more » « less
4. Treatment of ice-shelf evolution combining flow and flexure

   https://doi.org/10.1017/jog.2021.39  

   MacAyeal, Douglas R. ; Sergienko, Olga V. ; Banwell, Alison F. ; Macdonald, Grant J. ; Willis, Ian C. ; Stevens, Laura A. ( October 2021 , Journal of Glaciology)

   Abstract We develop a two-dimensional, plan-view formulation of ice-shelf flow and viscoelastic ice-shelf flexure. This formulation combines, for the first time, the shallow-shelf approximation for horizontal ice-shelf flow (and shallow-stream approximation for flow on lubricated beds such as where ice rises and rumples form), with the treatment of a thin-plate flexure. We demonstrate the treatment by performing two finite-element simulations: one of the relict pedestalled lake features that exist on some debris-covered ice shelves due to strong heterogeneity in surface ablation, and the other of ice rumpling in the grounding zone of an ice rise. The proposed treatment opens new venues to investigate physical processes that require coupling between the longitudinal deformation and vertical flexure, for instance, the effects of surface melting and supraglacial lakes on ice shelves, interactions with the sea swell, and many others. 

   more » « less
5. Wind–Wave Misalignment Effects on Langmuir Turbulence in Tropical Cyclone Conditions

   https://doi.org/10.1175/JPO-D-19-0093.1  

   Wang, Dong ; Kukulka, Tobias ; Reichl, Brandon G. ; Hara, Tetsu ; Ginis, Isaac ( December 2019 , Journal of Physical Oceanography)

   This study utilizes a large-eddy simulation (LES) approach to systematically assess the directional variability of wave-driven Langmuir turbulence (LT) in the ocean surface boundary layer (OSBL) under tropical cyclones (TCs). The Stokes drift vector, which drives LT through the Craik–Leibovich vortex force, is obtained through spectral wave simulations. LT’s direction is identified by horizontally elongated turbulent structures and objectively determined from horizontal autocorrelations of vertical velocities. In spite of a TC’s complex forcing with great wind and wave misalignments, this study finds that LT is approximately aligned with the wind. This is because the Reynolds stress and the depth-averaged Lagrangian shear (Eulerian plus Stokes drift shear) that are key in determining the LT intensity (determined by normalized depth-averaged vertical velocity variances) and direction are also approximately aligned with the wind relatively close to the surface. A scaling analysis of the momentum budget suggests that the Reynolds stress is approximately constant over a near-surface layer with predominant production of turbulent kinetic energy by Stokes drift shear, which is confirmed from the LES results. In this layer, Stokes drift shear, which dominates the Lagrangian shear, is aligned with the wind because of relatively short, wind-driven waves. On the contrary, Stokes drift exhibits considerable amount of misalignments with the wind. This wind–wave misalignment reduces LT intensity, consistent with a simple turbulent kinetic energy model. Our analysis shows that both the Reynolds stress and LT are aligned with the wind for different reasons: the former is dictated by the momentum budget, while the latter is controlled by wind-forced waves.

    

   more » « less