skip to main content


Title: Improving model-satellite comparisons of sea ice melt onset with a satellite simulator
Abstract. Seasonal transitions in Arctic sea ice, such as the melt onset, have been found to be useful metrics for evaluating sea ice in climate models against observations. However, comparisons of melt onset dates between climate models and satellite observations are indirect. Satellite data products of melt onset rely on observed brightness temperatures, while climate models do not currently simulate brightness temperatures, and must therefore define melt onset with other modeled variables. Here we adapt a passive microwave sea ice satellite simulator, the Arctic Ocean Observation Operator (ARC3O), to produce simulated brightness temperatures that can be used to diagnose the timing of the earliest snowmelt in climate models, as we show here using Community Earth System Model version 2 (CESM2) ocean-ice hindcasts. By producing simulated brightness temperatures and earliest snowmelt estimation dates using CESM2 and ARC3O, we facilitate new and previously impossible comparisons between the model and satellite observations by removing the uncertainty that arises due to definition differences. Direct comparisons between the model and satellite data allow us to identify an early bias across large areas of the Arctic at the beginning of the CESM2 ocean-ice hindcast melt season, as well as improve our understanding of the physical processes underlying seasonal changes in brightness temperatures. In particular, the ARC3O allows us to show that satellite algorithm-based melt onset dates likely occur after significant snowmelt has already taken place.  more » « less
Award ID(s):
1847398
NSF-PAR ID:
10357455
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
16
Issue:
8
ISSN:
1994-0424
Page Range / eLocation ID:
3235 to 3248
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract. Arctic sea ice experiences a dramatic annual cycle, and seasonal ice loss and growth can be characterized by various metrics: melt onset, breakup, opening, freeze onset, freeze-up, and closing. By evaluating a range of seasonal sea ice metrics, CMIP6 sea ice simulations can be evaluated in more detail than by using traditional metrics alone, such as sea ice area. We show that models capture the observed asymmetry in seasonal sea ice transitions, with spring ice loss taking about 1–2 months longer than fall ice growth. The largest impacts of internal variability are seen in the inflow regions for melt and freeze onset dates, but all metrics show pan-Arctic model spreads exceeding the internal variability range, indicating the contribution of model differences. Through climate model evaluation in the context of both observations and internal variability, we show that biases in seasonal transition dates can compensate for other unrealistic aspects of simulated sea ice. In some models, this leads to September sea ice areas in agreement with observations for the wrong reasons. 
    more » « less
  3. Abstract

    Observations show predictive skill of the minimum sea ice extent (Min SIE) from late winter anomalous offshore ice drift along the Eurasian coastline, leading to local ice thickness anomalies at the onset of the melt season—a signal then amplified by the ice–albedo feedback. We assess whether the observed seasonal predictability of September sea ice extent (Sept SIE) from Fram Strait Ice Area Export (FSIAE; a proxy for Eurasian coastal divergence) is present in global climate model (GCM) large ensembles, namely the CESM2-LE, GISS-E2.1-G, FLOR-LE, CNRM-CM6-1, and CanESM5. All models show distinct periods where winter FSIAE anomalies are negatively correlated with the May sea ice thickness (May SIT) anomalies along the Eurasian coastline, and the following Sept Arctic SIE, as in observations. Counterintuitively, several models show occasional periods where winter FSIAE anomalies are positively correlated with the following Sept SIE anomalies when the mean ice thickness is large, or late in the simulation when the sea ice is thin, and/or when internal variability increases. More important, periods with weak correlation between winter FSIAE and the following Sept SIE dominate, suggesting that summer melt processes generally dominate over late-winter preconditioning and May SIT anomalies. In general, we find that the coupling between the winter FSIAE and ice thickness anomalies along the Eurasian coastline at the onset of the melt season is a ubiquitous feature of GCMs and that the relationship with the following Sept SIE is dependent on the mean Arctic sea ice thickness.

     
    more » « less
  4. Abstract

    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
  5. The Yukon River basin encompasses over 832,000 km2 of boreal Arctic Alaska and northwest Canada, providing a major transportation corridor and multiple natural resources to regional communities. The river seasonal hydrology is defined by a long winter frozen season and a snowmelt-driven spring flood pulse. Capabilities for accurate monitoring and forecasting of the annual spring freshet and river ice breakup (RIB) in the Yukon and other northern rivers is limited, but critical for understanding hydrologic processes related to snow, and for assessing flood-related risks to regional communities. We developed a regional snow phenology record using satellite passive microwave remote sensing to elucidate interactions between the timing of upland snowmelt and the downstream spring flood pulse and RIB in the Yukon. The seasonal snow metrics included annual Main Melt Onset Date (MMOD), Snowoff (SO) and Snowmelt Duration (SMD) derived from multifrequency (18.7 and 36.5 GHz) daily brightness temperatures and a physically-based Gradient Ratio Polarization (GRP) retrieval algorithm. The resulting snow phenology record extends over a 29-year period (1988–2016) with 6.25 km grid resolution. The MMOD retrievals showed good agreement with similar snow metrics derived from in situ weather station measurements of snowpack water equivalence (r = 0.48, bias = −3.63 days) and surface air temperatures (r = 0.69, bias = 1 day). The MMOD and SO impact on the spring freshet was investigated by comparing areal quantiles of the remotely sensed snow metrics with measured streamflow quantiles over selected sub-basins. The SO 50% quantile showed the strongest (p < 0.1) correspondence with the measured spring flood pulse at Stevens Village (r = 0.71) and Pilot (r = 0.63) river gaging stations, representing two major Yukon sub-basins. MMOD quantiles indicating 20% and 50% of a catchment under active snowmelt corresponded favorably with downstream RIB (r = 0.61) from 19 river observation stations spanning a range of Yukon sub-basins; these results also revealed a 14–27 day lag between MMOD and subsequent RIB. Together, the satellite based MMOD and SO metrics show potential value for regional monitoring and forecasting of the spring flood pulse and RIB timing in the Yukon and other boreal Arctic basins. 
    more » « less