skip to main content

Title: Effects of Increasing the Category Resolution of the Sea Ice Thickness Distribution in a Coupled Climate Model on Arctic and Antarctic Sea Ice Mean State

Many modern sea ice models used in global climate models represent the subgrid‐scale heterogeneity in sea ice thickness with an ice thickness distribution (ITD), which improves model realism by representing the significant impact of the high spatial heterogeneity of sea ice thickness on thermodynamic and dynamic processes. Most models default to five thickness categories. However, little has been done to explore the effects of the resolution of this distribution (number of categories) on sea‐ice feedbacks in a coupled model framework and resulting representation of the sea ice mean state. Here, we explore this using sensitivity experiments in CESM2 with the standard 5 ice thickness categories and 15 ice thickness categories. Increasing the resolution of the ITD in a run with preindustrial climate forcing results in substantially thicker Arctic sea ice year‐round. Analyses show that this is a result of the ITD influence on ice strength. With 15 ITD categories, weaker ice occurs for the same average thickness, resulting in a higher fraction of ridged sea ice. In contrast, the higher resolution of thin ice categories results in enhanced heat conduction and bottom growth and leads to only somewhat increased winter Antarctic sea ice volume. The spatial resolution of the ICESat‐2 satellite mission provides a new opportunity to compare model outputs with observations of seasonal evolution of the ITD in the Arctic (ICESat‐2; 2018–2021). Comparisons highlight significant differences from the ITD modeled with both runs over this period, likely pointing to underlying issues contributing to the representation of average thickness.

more » « less
Award ID(s):
1724467 2138787
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Oceans
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    We explore the response of wintertime Arctic sea ice growth to strong cyclones and to large-scale circulation patterns on the daily scale using Earth system model output in phase 5 of the Coupled Model Intercomparison Project (CMIP5). A combined metrics ranking method selects three CMIP5 models that are successful in reproducing the wintertime Arctic dipole (AD) pattern. A cyclone identification method is applied to select strong cyclones in two subregions in the North Atlantic to examine their different impacts on sea ice growth. The total change of sea ice growth rate (SGR) is split into those respectively driven by the dynamic and thermodynamic atmospheric forcing. Three models reproduce the downward longwave radiation anomalies that generally match thermodynamic SGR anomalies in response to both strong cyclones and large-scale circulation patterns. For large-scale circulation patterns, the negative AD outweighs the positive Arctic Oscillation in thermodynamically inhibiting SGR in both impact area and magnitude. Despite the disagreement on the spatial distribution, the three CMIP5 models agree on the weaker response of dynamic SGR than thermodynamic SGR. As the Arctic warms, the thinner sea ice results in more ice production and smaller spatial heterogeneity of thickness, dampening the SGR response to the dynamic forcing. The higher temperature increases the specific heat of sea ice, thus dampening the SGR response to the thermodynamic forcing. In this way, the atmospheric forcing is projected to contribute less to change daily SGR in the future climate. 
    more » « less
  2. Abstract

    An Arctic sea ice‐ocean model is run with three uniform horizontal resolutions (6, 4, and 2 km) and identical sea ice and ocean model parameterizations, including an isotropic viscous‐plastic sea ice rheology, a mechanical ice strength parameterization, and an ice ridging parameterization. Driven by the same atmospheric forcing, the three model versions all produce similar spatial patterns and temporal variations of ice thickness and motion fields, resulting in almost identical magnitude and seasonal evolution of total ice volume and mean ice concentration, ice speed, and fractions of ice of various thickness categories over the Arctic Ocean. Increasing model resolution from 6 to 2 km does not significantly improve model performance when compared to NASA IceBridge ice thickness observations. This suggests that the large‐scale sea ice properties of the model are insensitive to varying high resolutions within the multifloe scale (2–10 km), and it may be unnecessary to adjust model parameters constantly with increasingly high resolutions. This is also true with models within the aggregate scale (10–75 km), indicating that model parameters used at coarse resolution may be used at high or multiscale resolution. However, even though the three versions all yield similar mean state of sea ice, they differ in representing anisotropic properties of sea ice. While they produce a basic pattern of major sea ice leads similar to satellite observations, their leads are distributed differently in space and time. Without changing model parameters and sea ice spatiotemporal variability, the 2‐km resolution model tends to capture more leads than the other two models.

    more » « less
  3. During the Arctic melt season, relatively fresh meltwater layers can accumulate under sea ice as a result of snow and ice melt, far from terrestrial freshwater inputs. Such under-ice meltwater layers, sometimes referred to as under-ice melt ponds, have been suggested to play a role in the summer sea ice mass balance both by isolating the sea ice from saltier water below, and by driving formation of ‘false bottoms’ below the sea ice. Such layers form at the interface of the fresher under-ice layer and the colder, saltier seawater below. During the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) expedition in the Central Arctic, we observed the presence of under-ice meltwater layers and false bottoms throughout July 2020 at primarily first-year ice locations. Here, we examine the distribution, prevalence, and drivers of under-ice ponds and the resulting false bottoms during this period. The average thickness of observed false bottoms and freshwater equivalent of under-ice meltwater layers was 0.08 m, with false bottom ice comprised of 74–87% FYI melt and 13–26% snow melt. Additionally, we explore these results using a 1D model to understand the role of dynamic influences on decoupling the ice from the seawater below. The model comparison suggests that the ice-ocean friction velocity was likely exceptionally low, with implications for air-ice-ocean momentum transfer. Overall, the prevalence of false bottoms was similar to or higher than noted during other observational campaigns, indicating that these features may in fact be common in the Arctic during the melt season. These results have implications for the broader ice-ocean system, as under-ice meltwater layers and false bottoms provide a source of ice growth during the melt season, potentially reduce fluxes between the ice and the ocean, isolate sea ice primary producers from pelagic nutrient sources, and may alter light transmission to the ocean below. 
    more » « less
  4. null (Ed.)
    Abstract. Ocean–sea-ice coupled models constrained by various observations provide different ice thickness estimates in the Antarctic. We evaluatecontemporary monthly ice thickness from four reanalyses in the Weddell Sea: the German contribution of the project Estimating the Circulation and Climate ofthe Ocean Version 2 (GECCO2), the Southern Ocean State Estimate (SOSE), the Ensemble Kalman Filter system based on the Nucleus for European Modelling of the Ocean (NEMO-EnKF) and the Global Ice–Ocean Modeling and Assimilation System (GIOMAS). The evaluation is performed againstreference satellite and in situ observations from ICESat-1, Envisat, upward-looking sonars and visual ship-based sea-ice observations. Compared withICESat-1, NEMO-EnKF has the highest correlation coefficient (CC) of 0.54 and lowest root mean square error (RMSE) of 0.44 m. Compared within situ observations, SOSE has the highest CC of 0.77 and lowest RMSE of 0.72 m. All reanalyses underestimate ice thickness near the coast ofthe western Weddell Sea with respect to ICESat-1 and in situ observations even though these observational estimates may be biased low. GECCO2 andNEMO-EnKF reproduce the seasonal variation in first-year ice thickness reasonably well in the eastern Weddell Sea. In contrast, GIOMAS ice thicknessperforms best in the central Weddell Sea, while SOSE ice thickness agrees most with the observations from the southern coast of the Weddell Sea. Inaddition, only NEMO-EnKF can reproduce the seasonal evolution of the large-scale spatial distribution of ice thickness, characterized by the thickice shifting from the southwestern and western Weddell Sea in summer to the western and northwestern Weddell Sea in spring. We infer that the thickice distribution is correlated with its better simulation of northward ice motion in the western Weddell Sea. These results demonstrate thepossibilities and limitations of using current sea-ice reanalysis for understanding the recent variability of sea-ice volume in the Antarctic. 
    more » « less
  5. Melt ponds on sea ice play an important role in the Arctic climate system. Their presence alters the partitioning of solar radiation: decreasing reflection, increasing absorption and transmission to the ice and ocean, and enhancing melt. The spatiotemporal properties of melt ponds thus modify ice albedo feedbacks and the mass balance of Arctic sea ice. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition presented a valuable opportunity to investigate the seasonal evolution of melt ponds through a rich array of atmosphere-ice-ocean measurements across spatial and temporal scales. In this study, we characterize the seasonal behavior and variability in the snow, surface scattering layer, and melt ponds from spring melt to autumn freeze-up using in situ surveys and auxiliary observations. We compare the results to satellite retrievals and output from two models: the Community Earth System Model (CESM2) and the Marginal Ice Zone Modeling and Assimilation System (MIZMAS). During the melt season, the maximum pond coverage and depth were 21% and 22 ± 13 cm, respectively, with distribution and depth corresponding to surface roughness and ice thickness. Compared to observations, both models overestimate melt pond coverage in summer, with maximum values of approximately 41% (MIZMAS) and 51% (CESM2). This overestimation has important implications for accurately simulating albedo feedbacks. During the observed freeze-up, weather events, including rain on snow, caused high-frequency variability in snow depth, while pond coverage and depth remained relatively constant until continuous freezing ensued. Both models accurately simulate the abrupt cessation of melt ponds during freeze-up, but the dates of freeze-up differ. MIZMAS accurately simulates the observed date of freeze-up, while CESM2 simulates freeze-up one-to-two weeks earlier. This work demonstrates areas that warrant future observation-model synthesis for improving the representation of sea-ice processes and properties, which can aid accurate simulations of albedo feedbacks in a warming climate. 
    more » « less