scholarly journals Adjoint sensitivities of sub-ice-shelf melt rates to ocean circulation under the Pine Island Ice Shelf, West Antarctica

2012 ◽  
Vol 53 (60) ◽  
pp. 59-69 ◽  
Author(s):  
Patrick Heimbach ◽  
Martin Losch

AbstractWe investigate the sensitivity of sub-ice-shelf melt rates under Pine Island Ice Shelf, West Antarctica, to changes in the oceanic state using an adjoint ocean model that is capable of representing the flow in sub-ice-shelf cavities. The adjoint code is based on algorithmic differentiation (AD) of the Massachusetts Institute of Technology’s ocean general circulation model (MITgcm). The adjoint model was extended by adding into the AD process the corresponding sub-ice-shelf cavity code, which implements a three-equation thermodynamic melt-rate parameterization to infer heat and freshwater fluxes at the ice-shelf/ocean boundary. The inferred sensitivities reveal dominant timescales of 30–60 days over which the shelf exit is connected to the deep interior via advective processes. They exhibit rich three-dimensional time-evolving patterns that can be understood in terms of a combination of the buoyancy forcing by inflowing water masses, the cavity geometry and the effect of rotation and topography in steering the flow in the presence of prominent features in the bedrock bathymetry. Dominant sensitivity pathways are found over a sill, as well as ‘shadow regions’ of very low sensitivities. To the extent that these transient patterns are robust they carry important information for decision-making in observation deployment and monitoring.

2003 ◽  
Vol 15 (1) ◽  
pp. 13-23 ◽  
Author(s):  
DAVID M. HOLLAND ◽  
STANLEY S. JACOBS ◽  
ADRIAN JENKINS

We applied a modified version of the Miami isopycnic coordinate ocean general circulation model (MICOM) to the ocean cavity beneath the Ross Ice Shelf to investigate the circulation of ocean waters in the sub-ice shelf cavity, along with the melting and freezing regimes at the base of the ice shelf. Model passive tracers are utilized to highlight the pathways of waters entering and exiting the cavity, and output is compared with data taken in the cavity and along the ice shelf front. High Salinity Shelf Water on the western Ross Sea continental shelf flows into the cavity along the sea floor and is transformed into Ice Shelf Water upon contact with the ice shelf base. Ice Shelf Water flows out of the cavity mainly around 180°, but also further east and on the western side of McMurdo Sound, as observed. Active ventilation of the region near the ice shelf front is forced by seasonal variations in the density structure of the water column to the north, driving rapid melting. Circulation in the more isolated interior is weaker, leading to melting at deeper ice and refreezing beneath shallower ice. Net melting over the whole ice shelf base is lower than other estimates, but is likely to increase as additional forcings are added to the model.


2008 ◽  
Vol 21 (11) ◽  
pp. 2558-2572 ◽  
Author(s):  
Paul R. Holland ◽  
Adrian Jenkins ◽  
David M. Holland

Abstract A three-dimensional ocean general circulation model is used to study the response of idealized ice shelves to a series of ocean-warming scenarios. The model predicts that the total ice shelf basal melt increases quadratically as the ocean offshore of the ice front warms. This occurs because the melt rate is proportional to the product of ocean flow speed and temperature in the mixed layer directly beneath the ice shelf, both of which are found to increase linearly with ocean warming. The behavior of this complex primitive equation model can be described surprisingly well with recourse to an idealized reduced system of equations, and it is shown that this system supports a melt rate response to warming that is generally quadratic in nature. This study confirms and unifies several previous examinations of the relation between melt rate and ocean temperature but disagrees with other results, particularly the claim that a single melt rate sensitivity to warming is universally valid. The hypothesized warming does not necessarily require a heat input to the ocean, as warmer waters (or larger volumes of “warm” water) may reach ice shelves purely through a shift in ocean circulation. Since ice shelves link the Antarctic Ice Sheet to the climate of the Southern Ocean, this finding of an above-linear rise in ice shelf mass loss as the ocean steadily warms is of significant importance to understanding ice sheet evolution and sea level rise.


Author(s):  
V. B. Zalesny ◽  
S. N. Moshonkin

Ocean general circulation model (OGCM) of the INM RAS with embedded k turbulent model is developed. The solution of the k model equations depends on the frequencies of buoyancy and velocity shift which are generated by the OGCM. The coefficients of vertical turbulence in OGCM depend on k and omega. The numerical algorithms of both models are based on the splitting method for physical processes. The model equations are split into two stages, describing the three-dimensional transport-diffusion of the kinetic energy of turbulence and frequency and their local generation-dissipation. The system of ordinary differential equations arising at the second stage is solved analytically, which ensures the efficiency of the algorithm. Analytical solution also written for the vertical turbulence coefficient equation. The model is used to study the sensitivity of the model circulation of the North AtlanticArctic Ocean to variations in the parameters of vertical turbulence. Experiments show that varying the coefficients of the analytical model solution can improve the adequacy of the simulation.


2021 ◽  
Author(s):  
Lennart Ramme ◽  
Jochem Marotzke

Abstract. When a snowball Earth deglaciates through a very high atmospheric CO2 concentration, the resulting inflow of freshwater leads to a stably stratified ocean, and the strong greenhouse conditions drive the climate into a very warm state. Here, we use a coupled atmosphere-ocean general circulation model, applying different scenarios for the evolution of atmospheric CO2, to conduct the first simulation of the climate and the three-dimensional ocean circulation in the aftermath of the Marinoan snowball Earth. The simulations show that the strong freshwater stratification breaks up on a timescale in the order of 103 years, mostly independent of the applied CO2 scenario. This is driven by the upwelling of salty waters in high latitudes, mainly the northern hemisphere, where a strong circumpolar current dominates the circulation. In the warmest CO2 scenario, the simulated Marinoan supergreenhouse climate reaches a global mean surface temperature of about 30 °C under an atmospheric CO2 concentration of 15 × 103 parts per million by volume, which is a moderate temperature compared to previous estimates. Consequently, the thermal expansion of seawater causes a sea-level rise of only 8 m, with most of it occurring during the first 3000 years. Our results imply that the surface temperatures of that time were potentially not as threatening for early metazoa as previously assumed. Furthermore, the short destratification timescale found in this study implies a very rapid accumulation of Marinoan cap dolostones, given that they were deposited in a freshwater environment.


2021 ◽  
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Huiskamp ◽  
Stefan Petri ◽  
Torsten Albrecht ◽  
...  

<p>The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high resolution configurations, limiting these studies to individual glaciers or regions over short time scales of decades to a few centuries. To study global and long term interactions, we developed a framework to couple the dynamic ice sheet model PISM with the global ocean general circulation model MOM5 via the ice-shelf cavity module PICO. Since ice-shelf cavities are not resolved by MOM5, but parameterized with the box model PICO, the framework allows the ice sheet and ocean model to be run at resolution of 16 km and 3 degrees, respectively. We present first results from our coupled setup and discuss stability, feedbacks, and interactions of the Antarctic Ice Sheet and the global ocean system on millennial time scales.</p>


2019 ◽  
Vol 49 (5) ◽  
pp. 1141-1157 ◽  
Author(s):  
Patrick Wagner ◽  
Siren Rühs ◽  
Franziska U. Schwarzkopf ◽  
Inga Monika Koszalka ◽  
Arne Biastoch

AbstractTo model tracer spreading in the ocean, Lagrangian simulations in an offline framework are a practical and efficient alternative to solving the advective–diffusive tracer equations online. Differences in both approaches raise the question of whether both methods are comparable. Lagrangian simulations usually use model output averaged in time, and trajectories are not subject to parameterized subgrid diffusion, which is included in the advection–diffusion equations of ocean models. Previous studies focused on diffusivity estimates in idealized models but could show that both methods yield similar results as long as the deformations-scale dynamics are resolved and a sufficient amount of Lagrangian particles is used. This study compares spreading of an Eulerian tracer simulated online and a cloud of Lagrangian particles simulated offline with velocities from the same ocean model. We use a global, eddy-resolving ocean model featuring 1/20° horizontal resolution in the Agulhas region around South Africa. Tracer and particles were released at one time step in the Cape Basin and below the mixed layer and integrated for 3 years. Large-scale diagnostics, like mean pathways of floats and tracer, are almost identical and 1D horizontal distributions show no significant differences. Differences in vertical distributions, seen in a reduced vertical spreading and downward displacement of particles, are due to the combined effect of unresolved subdaily variability of the vertical velocities and the spatial variation of vertical diffusivity. This, in turn, has a small impact on the horizontal spreading behavior. The estimates of eddy diffusivity from particles and tracer yield comparable results of about 4000 m2 s−1 in the Cape Basin.


2018 ◽  
Vol 12 (2) ◽  
pp. 453-476 ◽  
Author(s):  
Rachael D. Mueller ◽  
Tore Hattermann ◽  
Susan L. Howard ◽  
Laurie Padman

Abstract. Recent modeling studies of ocean circulation in the southern Weddell Sea, Antarctica, project an increase over this century of ocean heat into the cavity beneath Filchner–Ronne Ice Shelf (FRIS). This increase in ocean heat would lead to more basal melting and a modification of the FRIS ice draft. The corresponding change in cavity shape will affect advective pathways and the spatial distribution of tidal currents, which play important roles in basal melting under FRIS. These feedbacks between heat flux, basal melting, and tides will affect the evolution of FRIS under the influence of a changing climate. We explore these feedbacks with a three-dimensional ocean model of the southern Weddell Sea that is forced by thermodynamic exchange beneath the ice shelf and tides along the open boundaries. Our results show regionally dependent feedbacks that, in some areas, substantially modify the melt rates near the grounding lines of buttressed ice streams that flow into FRIS. These feedbacks are introduced by variations in meltwater production as well as the circulation of this meltwater within the FRIS cavity; they are influenced locally by sensitivity of tidal currents to water column thickness (wct) and non-locally by changes in circulation pathways that transport an integrated history of mixing and meltwater entrainment along flow paths. Our results highlight the importance of including explicit tidal forcing in models of future mass loss from FRIS and from the adjacent grounded ice sheet as individual ice-stream grounding zones experience different responses to warming of the ocean inflow.


2010 ◽  
Vol 40 (7) ◽  
pp. 1441-1457 ◽  
Author(s):  
Zhu Min Lu ◽  
Rui Xin Huang

Abstract Based on the classical Ekman layer theory, a simple analytical solution of the steady flow induced by a stationary hurricane in a homogenous ocean is discussed. The model consists of flow converging in an inward spiral in the deeper layer and diverging in the upper layer. The simple analytical model indicates that both the upwelling flux and the horizontal transport increase linearly with increasing radius of maximum winds. Furthermore, they both have a parabolic relationship with the maximum wind speed. The Coriolis parameter also affects the upwelling flux: the response to a hurricane is stronger at low latitudes than that at middle latitudes. Numerical solutions based on a regional version of an ocean general circulation model are similar to the primary results obtained through the analytical solution. Thus, the simplifications made in formulating the analytical solution are reasonable. Although the analytical solution in this paper is sought for a rather idealized ocean, it can help to make results from the more complicated numerical model understandable. These conceptual models provide a theoretical limit structure of the oceanic response to a moving hurricane over a stratified ocean.


2017 ◽  
Vol 47 (8) ◽  
pp. 1941-1959 ◽  
Author(s):  
David S. Trossman ◽  
Brian K. Arbic ◽  
David N. Straub ◽  
James G. Richman ◽  
Eric P. Chassignet ◽  
...  

AbstractMotivated by the substantial sensitivity of eddies in two-layer quasigeostrophic (QG) turbulence models to the strength of bottom drag, this study explores the sensitivity of eddies in more realistic ocean general circulation model (OGCM) simulations to bottom drag strength. The OGCM results are interpreted using previous results from horizontally homogeneous, two-layer, flat-bottom, f-plane, doubly periodic QG turbulence simulations and new results from two-layer, β-plane QG turbulence simulations run in a basin geometry with both flat and rough bottoms. Baroclinicity in all of the simulations varies greatly with drag strength, with weak drag corresponding to more barotropic flow and strong drag corresponding to more baroclinic flow. The sensitivity of the baroclinicity in the QG basin simulations to bottom drag is considerably reduced, however, when rough topography is used in lieu of a flat bottom. Rough topography reduces the sensitivity of the eddy kinetic energy amplitude and horizontal length scales in the QG basin simulations to bottom drag to an even greater degree. The OGCM simulation behavior is qualitatively similar to that in the QG rough-bottom basin simulations, in that baroclinicity is more sensitive to bottom drag strength than are eddy amplitudes or horizontal length scales. Rough topography therefore appears to mediate the sensitivity of eddies in models to the strength of bottom drag. The sensitivity of eddies to parameterized topographic internal lee wave drag, which has recently been introduced into some OGCMs, is also briefly discussed. Wave drag acts like a strong bottom drag in that it increases the baroclinicity of the flow, without strongly affecting eddy horizontal length scales.


Sign in / Sign up

Export Citation Format

Share Document