Search Results
You are looking at 11 - 20 of 22 items for
- Author or Editor: Paul R. Holland x
- Refine by Access: All Content x
Abstract
The transport of heat and salt through turbulent ice shelf–ocean boundary layers is a large source of uncertainty within ocean models of ice shelf cavities. This study uses small-scale, high-resolution, 3D numerical simulations to model an idealized boundary layer beneath a melting ice shelf to investigate the influence of ambient turbulence on double-diffusive convection (i.e., convection driven by the difference in diffusivities between salinity and temperature). Isotropic turbulence is forced throughout the simulations and the temperature and salinity are initialized with homogeneous values similar to observations. The initial temperature and the strength of forced turbulence are varied as controlling parameters within an oceanographically relevant parameter space. Two contrasting regimes are identified. In one regime double-diffusive convection dominates, and in the other convection is inhibited by the forced turbulence. The convective regime occurs for high temperatures and low turbulence levels, where it is long lived and affects the flow, melt rate, and melt pattern. A criterion for identifying convection in terms of the temperature and salinity profiles, and the turbulent dissipation rate, is proposed. This criterion may be applied to observations and theoretical models to quantify the effect of double-diffusive convection on ice shelf melt rates.
Abstract
The transport of heat and salt through turbulent ice shelf–ocean boundary layers is a large source of uncertainty within ocean models of ice shelf cavities. This study uses small-scale, high-resolution, 3D numerical simulations to model an idealized boundary layer beneath a melting ice shelf to investigate the influence of ambient turbulence on double-diffusive convection (i.e., convection driven by the difference in diffusivities between salinity and temperature). Isotropic turbulence is forced throughout the simulations and the temperature and salinity are initialized with homogeneous values similar to observations. The initial temperature and the strength of forced turbulence are varied as controlling parameters within an oceanographically relevant parameter space. Two contrasting regimes are identified. In one regime double-diffusive convection dominates, and in the other convection is inhibited by the forced turbulence. The convective regime occurs for high temperatures and low turbulence levels, where it is long lived and affects the flow, melt rate, and melt pattern. A criterion for identifying convection in terms of the temperature and salinity profiles, and the turbulent dissipation rate, is proposed. This criterion may be applied to observations and theoretical models to quantify the effect of double-diffusive convection on ice shelf melt rates.
Abstract
The ice shelf–ocean boundary current has an important control on heat delivery to the base of an ice shelf. Climate and regional models that include a representation of ice shelf cavities often use a coarse grid, and results have a strong dependence on resolution near the ice shelf–ocean interface. This study models the ice shelf–ocean boundary current with a nonhydrostatic z-level configuration at turbulence-permitting resolution (1 m). The z-level model performs well when compared against state-of-the-art large-eddy simulations, showing its capability in representing the correct physics. We show that theoretical results from a one-dimensional model with parameterized turbulence reproduce the z-level model results to a good degree, indicating possible utility as a turbulence closure. The one-dimensional model evolves to a state of marginal instability, and we use the z-level model to demonstrate how this is represented in three dimensions. Instabilities emerge that regulate the strength of the pycnocline and coexist with persistent Ekman rolls, which are identified prior to the flow becoming intermittently unstable. When resolution of the z-level model is degraded to understand the gridscale dependencies, the degradation is dominated by the established problem of excessive numerical diffusion. We show that at intermediate resolutions (2–4 m), the boundary layer structure can be partially recovered by tuning diffusivities. Last, we compare replacing prescribed melting with interactive melting that is dependent on the local ocean conditions. Interactive melting results in a feedback such that the system evolves more slowly, which is exaggerated at lower resolution.
Abstract
The ice shelf–ocean boundary current has an important control on heat delivery to the base of an ice shelf. Climate and regional models that include a representation of ice shelf cavities often use a coarse grid, and results have a strong dependence on resolution near the ice shelf–ocean interface. This study models the ice shelf–ocean boundary current with a nonhydrostatic z-level configuration at turbulence-permitting resolution (1 m). The z-level model performs well when compared against state-of-the-art large-eddy simulations, showing its capability in representing the correct physics. We show that theoretical results from a one-dimensional model with parameterized turbulence reproduce the z-level model results to a good degree, indicating possible utility as a turbulence closure. The one-dimensional model evolves to a state of marginal instability, and we use the z-level model to demonstrate how this is represented in three dimensions. Instabilities emerge that regulate the strength of the pycnocline and coexist with persistent Ekman rolls, which are identified prior to the flow becoming intermittently unstable. When resolution of the z-level model is degraded to understand the gridscale dependencies, the degradation is dominated by the established problem of excessive numerical diffusion. We show that at intermediate resolutions (2–4 m), the boundary layer structure can be partially recovered by tuning diffusivities. Last, we compare replacing prescribed melting with interactive melting that is dependent on the local ocean conditions. Interactive melting results in a feedback such that the system evolves more slowly, which is exaggerated at lower resolution.
Abstract
Leads are cracks in sea ice that often form because of deformation. During winter months, leads expose the ocean to the cold atmosphere, resulting in supercooling and the formation of frazil ice crystals within the mixed layer. Here the authors investigate the role of frazil ice formation in leads on the mass balance of the sea ice pack through the incorporation of a new module into the Los Alamos sea ice model (CICE). The frazil ice module considers an initial cooling of leads followed by a steady-state formation of uniformly distributed single size frazil ice crystals that precipitate to the ocean surface as grease ice. The grease ice is pushed against one of the lead edges by wind and water drag that the authors represent through a variable collection thickness for new sea ice. Simulations of the sea ice cover in the Arctic and Antarctic are performed and compared to a model that treats leads the same as the open ocean. The processes of ice formation in the new module slow down the refreezing of leads, resulting in a longer period of frazil ice production. The fraction of frazil-derived sea ice increases from 10% to 50%, corresponding better to observations. The new module has higher ice formation rates in areas of high ice concentration and thus has a greater impact within multiyear ice than it does in the marginal seas. The thickness of sea ice in the central Arctic increases by over 0.5 m, whereas within the Antarctic it remains unchanged.
Abstract
Leads are cracks in sea ice that often form because of deformation. During winter months, leads expose the ocean to the cold atmosphere, resulting in supercooling and the formation of frazil ice crystals within the mixed layer. Here the authors investigate the role of frazil ice formation in leads on the mass balance of the sea ice pack through the incorporation of a new module into the Los Alamos sea ice model (CICE). The frazil ice module considers an initial cooling of leads followed by a steady-state formation of uniformly distributed single size frazil ice crystals that precipitate to the ocean surface as grease ice. The grease ice is pushed against one of the lead edges by wind and water drag that the authors represent through a variable collection thickness for new sea ice. Simulations of the sea ice cover in the Arctic and Antarctic are performed and compared to a model that treats leads the same as the open ocean. The processes of ice formation in the new module slow down the refreezing of leads, resulting in a longer period of frazil ice production. The fraction of frazil-derived sea ice increases from 10% to 50%, corresponding better to observations. The new module has higher ice formation rates in areas of high ice concentration and thus has a greater impact within multiyear ice than it does in the marginal seas. The thickness of sea ice in the central Arctic increases by over 0.5 m, whereas within the Antarctic it remains unchanged.
Abstract
The exchange between the open ocean and sub–ice shelf cavities is important to both water mass transformations and ice shelf melting. Here, the authors use a high-resolution (500 m) numerical model to investigate to which degree eddies produced by frontal instability at the edge of a polynya are capable of transporting dense high-salinity shelf water (HSSW) underneath an ice shelf. The applied surface buoyancy flux and ice shelf geometry is based on Ronne Ice Shelf in the southern Weddell Sea, an area of intense wintertime sea ice production where a flow of HSSW into the cavity has been observed. Results show that eddies are able to enter the cavity at the southwestern corner of the polynya where an anticyclonic rim current intersects the ice shelf front. The size and time scale of simulated eddies are in agreement with observations close to the Ronne Ice Front. The properties and strength of the inflow are sensitive to the prescribed total ice production, flushing the ice shelf cavity at a rate of 0.2–0.4 × 106 m3 s−1 depending on polynya size and magnitude of surface buoyancy flux. Eddy-driven HSSW transport into the cavity is reduced by about 50% if the model grid resolution is decreased to 2–5 km and eddies are not properly resolved.
Abstract
The exchange between the open ocean and sub–ice shelf cavities is important to both water mass transformations and ice shelf melting. Here, the authors use a high-resolution (500 m) numerical model to investigate to which degree eddies produced by frontal instability at the edge of a polynya are capable of transporting dense high-salinity shelf water (HSSW) underneath an ice shelf. The applied surface buoyancy flux and ice shelf geometry is based on Ronne Ice Shelf in the southern Weddell Sea, an area of intense wintertime sea ice production where a flow of HSSW into the cavity has been observed. Results show that eddies are able to enter the cavity at the southwestern corner of the polynya where an anticyclonic rim current intersects the ice shelf front. The size and time scale of simulated eddies are in agreement with observations close to the Ronne Ice Front. The properties and strength of the inflow are sensitive to the prescribed total ice production, flushing the ice shelf cavity at a rate of 0.2–0.4 × 106 m3 s−1 depending on polynya size and magnitude of surface buoyancy flux. Eddy-driven HSSW transport into the cavity is reduced by about 50% if the model grid resolution is decreased to 2–5 km and eddies are not properly resolved.
Abstract
Unlike the rapid sea ice losses reported in the Arctic, satellite observations show an overall increase in Antarctic sea ice concentration over recent decades. However, observations of decadal trends in Antarctic ice thickness, and hence ice volume, do not currently exist. In this study a model of the Southern Ocean and its sea ice, forced by atmospheric reanalyses, is used to assess 1992–2010 trends in ice thickness and volume. The model successfully reproduces observations of mean ice concentration, thickness, and drift, and decadal trends in ice concentration and drift, imparting some confidence in the hindcasted trends in ice thickness. The model suggests that overall Antarctic sea ice volume has increased by approximately 30 km3 yr−1 (0.4% yr−1) as an equal result of areal expansion (20 × 103 km2 yr−1 or 0.2% yr−1) and thickening (1.5 mm yr−1 or 0.2% yr−1). This ice volume increase is an order of magnitude smaller than the Arctic decrease, and about half the size of the increased freshwater supply from the Antarctic Ice Sheet. Similarly to the observed ice concentration trends, the small overall increase in modeled ice volume is actually the residual of much larger opposing regional trends. Thickness changes near the ice edge follow observed concentration changes, with increasing concentration corresponding to increased thickness. Ice thickness increases are also found in the inner pack in the Amundsen and Weddell Seas, where the model suggests that observed ice-drift trends directed toward the coast have caused dynamical thickening in autumn and winter. Modeled changes are predominantly dynamic in origin in the Pacific sector and thermodynamic elsewhere.
Abstract
Unlike the rapid sea ice losses reported in the Arctic, satellite observations show an overall increase in Antarctic sea ice concentration over recent decades. However, observations of decadal trends in Antarctic ice thickness, and hence ice volume, do not currently exist. In this study a model of the Southern Ocean and its sea ice, forced by atmospheric reanalyses, is used to assess 1992–2010 trends in ice thickness and volume. The model successfully reproduces observations of mean ice concentration, thickness, and drift, and decadal trends in ice concentration and drift, imparting some confidence in the hindcasted trends in ice thickness. The model suggests that overall Antarctic sea ice volume has increased by approximately 30 km3 yr−1 (0.4% yr−1) as an equal result of areal expansion (20 × 103 km2 yr−1 or 0.2% yr−1) and thickening (1.5 mm yr−1 or 0.2% yr−1). This ice volume increase is an order of magnitude smaller than the Arctic decrease, and about half the size of the increased freshwater supply from the Antarctic Ice Sheet. Similarly to the observed ice concentration trends, the small overall increase in modeled ice volume is actually the residual of much larger opposing regional trends. Thickness changes near the ice edge follow observed concentration changes, with increasing concentration corresponding to increased thickness. Ice thickness increases are also found in the inner pack in the Amundsen and Weddell Seas, where the model suggests that observed ice-drift trends directed toward the coast have caused dynamical thickening in autumn and winter. Modeled changes are predominantly dynamic in origin in the Pacific sector and thermodynamic elsewhere.
Abstract
Both the Southern Annular Mode (SAM) and El Niño–Southern Oscillation (ENSO) are critical factors contributing to Antarctic sea ice variability on interannual time scales. However, their joint effects on sea ice are complex and remain unclear for each austral season. In this study, satellite sea ice concentration (SIC) observations and atmospheric reanalysis data are utilized to assess the impacts of combined SAM and ENSO on seasonal Antarctic sea ice changes. The joint SAM–ENSO impacts on southern high latitudes are principally controlled by the strength and position of the wave activity and associated atmospheric circulation anomalies affected by their interactions. In-phase events (La Niña/positive SAM and El Niño/negative SAM) are characterized with an SIC dipole located in the Weddell/Bellingshausen Seas and Amundsen/Ross Seas, while out-of-phase events (El Niño/positive SAM and La Niña/negative SAM) experience significant SIC anomalies in the Indian Ocean and western Pacific Ocean. Sea ice budget analyses are conducted to separate the dynamic and thermodynamic contributions inducing the sea ice intensification anomalies. The results show that in-phase intensification anomalies also display a pattern similar to the SIC dipole and are mainly driven by the direct thermodynamic forcing at the ice edge and thermodynamic responses to meridional sea ice drift in the inner pack, especially in autumn and winter. Dynamic processes caused by zonal sea ice drift also play an important role during out-of-phase conditions in addition to the same mechanisms during in-phase conditions.
Abstract
Both the Southern Annular Mode (SAM) and El Niño–Southern Oscillation (ENSO) are critical factors contributing to Antarctic sea ice variability on interannual time scales. However, their joint effects on sea ice are complex and remain unclear for each austral season. In this study, satellite sea ice concentration (SIC) observations and atmospheric reanalysis data are utilized to assess the impacts of combined SAM and ENSO on seasonal Antarctic sea ice changes. The joint SAM–ENSO impacts on southern high latitudes are principally controlled by the strength and position of the wave activity and associated atmospheric circulation anomalies affected by their interactions. In-phase events (La Niña/positive SAM and El Niño/negative SAM) are characterized with an SIC dipole located in the Weddell/Bellingshausen Seas and Amundsen/Ross Seas, while out-of-phase events (El Niño/positive SAM and La Niña/negative SAM) experience significant SIC anomalies in the Indian Ocean and western Pacific Ocean. Sea ice budget analyses are conducted to separate the dynamic and thermodynamic contributions inducing the sea ice intensification anomalies. The results show that in-phase intensification anomalies also display a pattern similar to the SIC dipole and are mainly driven by the direct thermodynamic forcing at the ice edge and thermodynamic responses to meridional sea ice drift in the inner pack, especially in autumn and winter. Dynamic processes caused by zonal sea ice drift also play an important role during out-of-phase conditions in addition to the same mechanisms during in-phase conditions.
Abstract
In the Southern Ocean the Antarctic Circumpolar Current is significantly steered by large topographic features, and subpolar gyres form in their lee. The geometry of topographic features in the Southern Ocean is highly variable, but the influence of this variation on the large-scale flow is poorly understood. Using idealized barotropic simulations of a zonal channel with a meridional ridge, it is found that the ridge geometry is important for determining the net zonal volume transport. A relationship is observed between ridge width and volume transport that is determined by the form stress generated by the ridge. Gyre formation is also highly reliant on the ridge geometry. A steep ridge allows gyres to form within regions of unblocked geostrophic (f/H) contours, with an increase in gyre strength as the ridge width is reduced. These relationships among ridge width, gyre strength, and net zonal volume transport emerge to simultaneously satisfy the conservation of momentum and vorticity.
Abstract
In the Southern Ocean the Antarctic Circumpolar Current is significantly steered by large topographic features, and subpolar gyres form in their lee. The geometry of topographic features in the Southern Ocean is highly variable, but the influence of this variation on the large-scale flow is poorly understood. Using idealized barotropic simulations of a zonal channel with a meridional ridge, it is found that the ridge geometry is important for determining the net zonal volume transport. A relationship is observed between ridge width and volume transport that is determined by the form stress generated by the ridge. Gyre formation is also highly reliant on the ridge geometry. A steep ridge allows gyres to form within regions of unblocked geostrophic (f/H) contours, with an increase in gyre strength as the ridge width is reduced. These relationships among ridge width, gyre strength, and net zonal volume transport emerge to simultaneously satisfy the conservation of momentum and vorticity.
Abstract
The melt rate of Antarctic ice shelves is of key importance for rising sea levels and future climate scenarios. Recent observations beneath Larsen C Ice Shelf revealed an ocean boundary layer that was highly turbulent and raised questions on the effect of these rich flow dynamics on the ocean heat transfer and the ice shelf melt rate. Directly motivated by the field observations, we have conducted large-eddy simulations (LES) to further examine the ocean boundary layer beneath Larsen C Ice Shelf. The LES was initialized with uniform temperature and salinity (T–S) and included a realistic tidal cycle and a small basal slope. A new parameterization based on previous work was applied at the top boundary to model near-wall turbulence and basal melting. The resulting vertical T–S profiles, melt rate, and friction velocity matched well with the Larsen C Ice Shelf observations. The instantaneous melt rate varied strongly with the tidal cycle, with faster flow increasing the turbulence and mixing of heat toward the ice base. An Ekman layer formed beneath the ice base and, due to the strong vertical shear of the current, Ekman rolls appeared in the mixed layer and stratified region (depth ≈ 20–60 m). In an additional high-resolution simulation (conducted with a smaller domain) the Ekman rolls were associated with increased turbulent kinetic energy, but a relatively small vertical heat flux. Our results will help with interpreting field observations and parameterizing the ocean-driven basal melting of ice shelves.
Abstract
The melt rate of Antarctic ice shelves is of key importance for rising sea levels and future climate scenarios. Recent observations beneath Larsen C Ice Shelf revealed an ocean boundary layer that was highly turbulent and raised questions on the effect of these rich flow dynamics on the ocean heat transfer and the ice shelf melt rate. Directly motivated by the field observations, we have conducted large-eddy simulations (LES) to further examine the ocean boundary layer beneath Larsen C Ice Shelf. The LES was initialized with uniform temperature and salinity (T–S) and included a realistic tidal cycle and a small basal slope. A new parameterization based on previous work was applied at the top boundary to model near-wall turbulence and basal melting. The resulting vertical T–S profiles, melt rate, and friction velocity matched well with the Larsen C Ice Shelf observations. The instantaneous melt rate varied strongly with the tidal cycle, with faster flow increasing the turbulence and mixing of heat toward the ice base. An Ekman layer formed beneath the ice base and, due to the strong vertical shear of the current, Ekman rolls appeared in the mixed layer and stratified region (depth ≈ 20–60 m). In an additional high-resolution simulation (conducted with a smaller domain) the Ekman rolls were associated with increased turbulent kinetic energy, but a relatively small vertical heat flux. Our results will help with interpreting field observations and parameterizing the ocean-driven basal melting of ice shelves.
ABSTRACT
Open-ocean polynyas in the Weddell Sea of Antarctica are the product of deep convection, which transports Warm Deep Water (WDW) to the surface and melts sea ice or prevents its formation. These polynyas occur only rarely in the observational record but are a near-permanent feature of many climate and ocean simulations. A question not previously considered is the degree to which the Weddell polynya affects the nearby Filchner–Ronne Ice Shelf (FRIS) cavity. Here we assess these effects using regional ocean model simulations of the Weddell Sea and FRIS, where deep convection is imposed with varying area, location, and duration. In these simulations, the idealized Weddell polynyas consistently cause an increase in WDW transport onto the continental shelf as a result of density changes above the shelf break. This leads to saltier, denser source waters for the FRIS cavity, which then experiences stronger circulation and increased ice shelf basal melting. It takes approximately 14 years for melt rates to return to normal after the deep convection ceases. Weddell polynyas similar to those seen in observations have a modest impact on FRIS melt rates, which is within the range of simulated interannual variability. However, polynyas that are larger or closer to the shelf break, such as those seen in many ocean models, trigger a stronger response. These results suggest that ocean models with excessive Weddell Sea convection may not be suitable boundary conditions for regional models of the Antarctic continental shelf and ice shelf cavities.
ABSTRACT
Open-ocean polynyas in the Weddell Sea of Antarctica are the product of deep convection, which transports Warm Deep Water (WDW) to the surface and melts sea ice or prevents its formation. These polynyas occur only rarely in the observational record but are a near-permanent feature of many climate and ocean simulations. A question not previously considered is the degree to which the Weddell polynya affects the nearby Filchner–Ronne Ice Shelf (FRIS) cavity. Here we assess these effects using regional ocean model simulations of the Weddell Sea and FRIS, where deep convection is imposed with varying area, location, and duration. In these simulations, the idealized Weddell polynyas consistently cause an increase in WDW transport onto the continental shelf as a result of density changes above the shelf break. This leads to saltier, denser source waters for the FRIS cavity, which then experiences stronger circulation and increased ice shelf basal melting. It takes approximately 14 years for melt rates to return to normal after the deep convection ceases. Weddell polynyas similar to those seen in observations have a modest impact on FRIS melt rates, which is within the range of simulated interannual variability. However, polynyas that are larger or closer to the shelf break, such as those seen in many ocean models, trigger a stronger response. These results suggest that ocean models with excessive Weddell Sea convection may not be suitable boundary conditions for regional models of the Antarctic continental shelf and ice shelf cavities.
Abstract
The top 2000 m of the Southern Ocean has freshened and warmed over recent decades. However, the high-latitude (south of 50°S) southeast Pacific was observed to be cooler and fresher in the years 2008–10 compared to 2005–07 over a wide depth range including surface, mode, and intermediate waters. The causes and impacts of this event are analyzed using the ocean–sea ice data-assimilating Southern Ocean State Estimate (SOSE) and observationally based products. In 2008–10, a strong positive southern annular mode coincided with a negative El Niño–Southern Oscillation and a deep Amundsen Sea low. Enhanced meridional winds drove strong sea ice export from the eastern Ross Sea, bringing large amounts of ice to the Amundsen Sea ice edge. In 2008, together with increased precipitation, this introduced a strong freshwater anomaly that was advected eastward by the Antarctic Circumpolar Current (ACC), mixing along the way. This anomaly entered the ocean interior not only as Antarctic Intermediate Water, but also as lighter Southeast Pacific Subantarctic Mode Water (SEPSAMW). A numerical particle release experiment carried out in SOSE showed that the Ross Sea sector was the dominant source of particles reaching the SEPSAMW formation region. This suggests that large-scale climate fluctuations can induce strong interannual variability of volume and properties of SEPSAMW. These fluctuations act at different time scales: instantaneously via direct forcing and also lagged over advective time scales of several years from upstream regions.
Abstract
The top 2000 m of the Southern Ocean has freshened and warmed over recent decades. However, the high-latitude (south of 50°S) southeast Pacific was observed to be cooler and fresher in the years 2008–10 compared to 2005–07 over a wide depth range including surface, mode, and intermediate waters. The causes and impacts of this event are analyzed using the ocean–sea ice data-assimilating Southern Ocean State Estimate (SOSE) and observationally based products. In 2008–10, a strong positive southern annular mode coincided with a negative El Niño–Southern Oscillation and a deep Amundsen Sea low. Enhanced meridional winds drove strong sea ice export from the eastern Ross Sea, bringing large amounts of ice to the Amundsen Sea ice edge. In 2008, together with increased precipitation, this introduced a strong freshwater anomaly that was advected eastward by the Antarctic Circumpolar Current (ACC), mixing along the way. This anomaly entered the ocean interior not only as Antarctic Intermediate Water, but also as lighter Southeast Pacific Subantarctic Mode Water (SEPSAMW). A numerical particle release experiment carried out in SOSE showed that the Ross Sea sector was the dominant source of particles reaching the SEPSAMW formation region. This suggests that large-scale climate fluctuations can induce strong interannual variability of volume and properties of SEPSAMW. These fluctuations act at different time scales: instantaneously via direct forcing and also lagged over advective time scales of several years from upstream regions.
