skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Homogenization for convection-enhanced thermal transport in sea ice
Sea ice regulates heat exchange between the ocean and atmosphere in Earth’s polar regions. The thermal conductivity of sea ice governs this exchange, and is a key parameter in climate modelling. However, it is challenging to measure and predict due to its sensitive dependence on temperature, salinity and brine microstructure. Moreover, as temperature increases, sea ice becomes permeable, and fluid can flow through the porous microstructure. While models for thermal diffusion through sea ice have been obtained, advective contributions to transport have not been considered theoretically. Here, we homogenize a multiscale advection–diffusion equation that models thermal transport through porous sea ice when fluid flow is present. We consider two-dimensional models of convective flow and use an integral representation to derive bounds on the thermal conductivity as a function of the Péclet number. These bounds guarantee enhancement in the thermal conductivity due to the added flow. Further, we relate the Péclet number to temperature, making these bounds useful for global climate models. Our analytic approach offers a mathematical theory which can not only improve predictions of atmosphere–ice–ocean heat exchanges in climate models, but can provide a theoretical framework for a range of problems involving advection–diffusion processes in various fields of application.  more » « less
Award ID(s):
2206171 2136198
PAR ID:
10592779
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
The Royal Society Publishing
Date Published:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume:
480
Issue:
2296
ISSN:
1364-5021
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The mixed layer of polynyas is vital for local climate as it determines the exchange of properties and energy between ocean, sea ice, and atmosphere. However, its evolution is poorly understood, as it is controlled by complex interactions among these components, yet highly undersampled, especially outside summer. Here, we present a 2-month, high vertical-resolution, full-depth hydrographic dataset from the southeastern Amundsen Sea polynya in austral autumn (from mid-February to mid-April 2014) collected by a recovered seal tag. This novel dataset quantifies the changes in upper-ocean temperature and salinity stratification in this previously unobserved season. Our seal-tag measurements reveal that the mixed layer experiences deepening, salinification, and intense heat loss through surface fluxes. Heat and salt budgets suggest a sea ice formation rate of ∼3 cm per day. We use a one-dimensional model to reproduce the mixed layer evolution and further identify key controls on its characteristics. Our experiments with a range of reduced or amplified air–sea fluxes show that heat loss to the atmosphere and related sea ice formation are the principal determinants of stratification evolution. Additionally, our modeling demonstrates that horizontal advection is required to fully explain the mixed layer evolution, underlining the importance of the ice-covered neighboring region for determining sea ice formation rates in the Amundsen Sea polynya. Our findings suggest that the potential overestimation of sea ice production by satellite-based methods, due to the absence of oceanic heat flux, could be offset by horizontal advection inhibiting mixed layer deepening and sustaining sea ice formation. 
    more » « less
  2. Abstract We investigate the role of ocean heat transport (OHT) in driving the decadal variability of the Arctic climate by analyzing the pre‐industrial control simulation of a high‐resolution climate model. While the OHT variability at 65°N is greater in the Atlantic, we find that the decadal variability of Arctic‐wide surface temperature and sea ice area is much better correlated with Bering Strait OHT than Atlantic OHT. In particular, decadal Bering Strait OHT variability causes significant changes in local sea ice cover and air‐sea heat fluxes, which are amplified by shortwave feedbacks. These heat flux anomalies are regionally balanced by longwave radiation at the top of the atmosphere, without compensation by atmospheric heat transport (Bjerknes compensation). The sensitivity of the Arctic to changes in OHT may thus rely on an accurate representation of the heat transport through the Bering Strait, which is difficult to resolve in coarse‐resolution ocean models. 
    more » « less
  3. Abstract A number of feedbacks regulate the response of Arctic sea ice to local atmospheric warming. Using a realistic coupled ocean‐sea ice model and its adjoint, we isolate a mechanism by which significant ice growth at the end of the melt season may occur as a lagged response to Arctic atmospheric warming. A series of perturbation simulations informed by adjoint model‐derived sensitivity patterns reveal the enhanced ice growth to be accompanied by a reduction of snow thickness on the ice pack. Detailed analysis of ocean‐ice‐snow heat budgets confirms the essential role of the reduced snow thickness for persistence and delayed overshoot of ice growth. The underlying mechanism is a snow‐melt‐conductivity feedback, wherein atmosphere‐driven snow melt leads to a larger conductive ocean heat loss through the overlying ice layer. Our results highlight the need for accurate observations of snow thickness to constrain climate models and to initialize sea ice forecasts. 
    more » « less
  4. Abstract High latitudes, including the Bering Sea, are experiencing unprecedented rates of change. Long-term Bering Sea warming trends have been identified, and marine heatwaves (MHWs), event-scale elevated sea surface temperature (SST) extremes, have also increased in frequency and longevity in recent years. Recent work has shown that variability in air–sea coupling plays a dominant role in driving Bering Sea upper-ocean thermal variability and that surface forcing has driven an increase in the occurrence of positive ocean temperature anomalies since 2010. In this work, we characterize the drivers of the anomalous surface air–sea heat fluxes in the Bering Sea over the period 2010–22 using ERA5 fields. We show that the surface turbulent heat flux dominates the net surface heat flux variability from September to April and is primarily a result of near-surface air temperature and specific humidity anomalies. The airmass anomalies that account for the majority of the turbulent heat flux variability are a function of wind direction, with southerly (northerly) wind advecting anomalously warm (cool), moist (dry) air over the Bering Sea, resulting in positive (negative) surface turbulent flux anomalies. During the remaining months of the year, anomalies in the surface radiative fluxes account for the majority of the net surface heat flux variability and are a result of anomalous cloud coverage, anomalous lower-tropospheric virtual temperature, and sea ice coverage variability. Our results indicate that atmospheric variability drives much of the Bering Sea upper-ocean temperature variability through the mediation of the surface heat fluxes during the analysis period. Significance StatementA long-term ocean warming trend and a recent increase in marine heatwaves in the Bering Sea have been identified. Previous work showed that anomalies in the exchange of heat between the ocean and the atmosphere were the primary driver of Bering Sea temperature variability, but the processes responsible for the heat exchange anomalies were unknown. In this work, we show that the atmosphere is the primary driver of anomalies in the Bering Sea air–sea heat exchange and therefore plays an important role in altering the thermal state of the Bering Sea. Our results highlight the importance of understanding more about how the ocean and the atmosphere interact at high latitudes and how this relationship will be affected by future climate change. 
    more » « less
  5. Abstract. The Laurentian Great Lakes significantly influence the climate of the Midwest and Northeast United States due to their vast thermal inertia, moisture source potential, and complex heat and moisture flux dynamics. This study presents a newly developed coupled lake–ice–atmosphere (CLIAv1) modeling system for the Great Lakes by coupling the National Aeronautics and Space Administration (NASA) Unified Weather Research and Forecasting (NU-WRF) regional climate model (RCM) with the three-dimensional (3D) Finite Volume Community Ocean Model (FVCOM) and investigates the impact of coupled dynamics on simulations of the Great Lakes' winter climate. By integrating 3D lake hydrodynamics, CLIAv1 demonstrates superior performance in reproducing observed lake surface temperatures (LSTs), ice cover distribution, and the vertical thermal structure of the Great Lakes compared to the NU-WRF model coupled with the default 1D Lake Ice Snow and Sediment Simulator (LISSS). CLIAv1 also enhances the simulation of over-lake atmospheric conditions, including air temperature, wind speed, and sensible and latent heat fluxes, underscoring the importance of resolving complex lake dynamics for reliable regional Earth system projections. More importantly, the key contribution of this study is the identification of critical physical processes that influence lake thermal structure and ice cover – processes that are missed by 1D lake models but effectively resolved by 3D lake models. Through process-oriented numerical experiments, we identify key 3D hydrodynamic processes – ice transport, heat advection, and shear production in turbulence – that explain the superiority of 3D lake models to 1D lake models, particularly in cold season performance and lake–atmosphere interactions. Critically, all three of these processes are dynamically linked to water currents – spatially and temporally evolving flow fields that are structurally absent in 1D models. This study aims to advance our understanding of the physical mechanisms that underlie the fundamental differences between 3D and 1D lake models in simulating key hydrodynamic processes during the winter season, and it offers generalized insights that are not constrained by specific model configurations. 
    more » « less