Improving the WRF model's simulation over sea ice surface   through coupling with a complex thermodynamic sea ice model

Abstract. Sea ice plays an important role in the air–ice–ocean interaction, but it is often represented simply in many regional atmospheric models. The Noah sea ice model, which has been widely used in the Weather Research and Forecasting (WRF) model, exhibits cold bias in simulating the Arctic sea ice temperature when validated against the Surface Heat Budget of the Arctic Ocean (SHEBA) in situ observations. According to sensitivity tests, this bias is attributed not only to the simulation of snow depth and turbulent fluxes but also to the heat conduction within snow and ice. Compared with the Noah sea ice model, the high-resolution thermodynamic snow and ice model (HIGHTSI) has smaller bias in simulating the sea ice temperature. HIGHTSI is further coupled with the WRF model to evaluate the possible added value from better resolving the heat transport and solar penetration in sea ice from a complex thermodynamic sea ice model. The cold bias in simulating the surface temperature over sea ice in winter by the original Polar WRF is reduced when HIGHTSI rather than Noah is coupled with the WRF model, and this also leads to a better representation of surface upward longwave radiation and 2 m air temperature. A discussion on the impact of specifying sea ice thickness in the WRF model is presented. Consistent with previous research, prescribing the sea ice thickness with observational information would result in the best simulation among the available methods. If no observational information is available, using an empirical method based on the relationship between sea ice concentration and sea ice thickness could mimic the large-scale spatial feature of sea ice thickness. The potential application of a thermodynamic sea ice model in predicting the change in sea ice thickness in a RCM is limited by the lack of sea ice dynamic processes in the model and the coarse assumption on the initial value of sea ice thickness.

Download Full-text

Improving the WRF model's (version 3.6.1) simulation over sea ice surface through coupling with a complex thermodynamic sea ice model (HIGHTSI)

Geoscientific Model Development ◽

10.5194/gmd-9-2239-2016 ◽

2016 ◽

Vol 9 (6) ◽

pp. 2239-2254 ◽

Cited By ~ 3

Author(s):

Yao Yao ◽

Jianbin Huang ◽

Yong Luo ◽

Zongci Zhao

Keyword(s):

Sea Ice ◽

Heat Budget ◽

Wrf Model ◽

The Arctic ◽

Wave Radiation ◽

Ice Thickness ◽

Sea Ice Thickness ◽

The Impact ◽

Ice Model ◽

Snow And Ice

Abstract. Sea ice plays an important role in the air–ice–ocean interaction, but it is often represented simply in many regional atmospheric models. The Noah sea ice scheme, which is the only option in the current Weather Research and Forecasting (WRF) model (version 3.6.1), has a problem of energy imbalance due to its simplification in snow processes and lack of ablation and accretion processes in ice. Validated against the Surface Heat Budget of the Arctic Ocean (SHEBA) in situ observations, Noah underestimates the sea ice temperature which can reach −10 °C in winter. Sensitivity tests show that this bias is mainly attributed to the simulation within the ice when a time-dependent ice thickness is specified. Compared with the Noah sea ice model, the high-resolution thermodynamic snow and ice model (HIGHTSI) uses more realistic thermodynamics for snow and ice. Most importantly, HIGHTSI includes the ablation and accretion processes of sea ice and uses an interpolation method which can ensure the heat conservation during its integration. These allow the HIGHTSI to better resolve the energy balance in the sea ice, and the bias in sea ice temperature is reduced considerably. When HIGHTSI is coupled with the WRF model, the simulation of sea ice temperature by the original Polar WRF is greatly improved. Considering the bias with reference to SHEBA observations, WRF-HIGHTSI improves the simulation of surface temperature, 2 m air temperature and surface upward long-wave radiation flux in winter by 6, 5 °C and 20 W m−2, respectively. A discussion on the impact of specifying sea ice thickness in the WRF model is presented. Consistent with previous research, prescribing the sea ice thickness with observational information results in the best simulation among the available methods. If no observational information is available, we present a new method in which the sea ice thickness is initialized from empirical estimation and its further change is predicted by a complex thermodynamic sea ice model. The ice thickness simulated by this method depends much on the quality of the initial guess of the ice thickness and the role of the ice dynamic processes.

Download Full-text

Suppression of Arctic sea ice growth by winter clouds and snowfall

10.21203/rs.3.rs-84779/v1 ◽

2021 ◽

Author(s):

Won-il Lim ◽

Hyo-Seok Park ◽

Andrew Stewart ◽

Kyong-Hwan Seo

Keyword(s):

Sea Ice ◽

Longwave Radiation ◽

Arctic Sea Ice ◽

The Arctic ◽

Ice Thickness ◽

Sea Ice Thickness ◽

Ice Growth ◽

Downward Longwave Radiation ◽

Arctic Sea ◽

The Impact

Abstract The ongoing Arctic warming has been pronounced in winter and has been associated with an increase in downward longwave radiation. While previous studies have demonstrated that poleward moisture flux into the Arctic strengthens downward longwave radiation, less attention has been given to the impact of the accompanying increase in snowfall. Here, utilizing state-of-the art sea ice models, we show that typical winter snowfall anomalies of 1.0 cm, accompanied by positive downward longwave radiation anomalies of ~5 W m-2 can decrease sea ice thickness by around 5 cm in the following spring over the Eurasian Seas. This basin-wide ice thinning is followed by a shrinking of summer ice extent in extreme cases. In the winter of 2016–17, anomalously strong warm/moist air transport combined with ~2.5 cm increase in snowfall decreased spring ice thickness by ~10 cm and decreased the following summer sea ice extent by 5–30%. Projected future reductions in the thickness of Arctic sea ice and snow will amplify the impact of anomalous winter snowfall events on winter sea ice growth and seasonal sea ice thickness.

Download Full-text

Using CloudSat snowfall rate observations to constrain and characterize the uncertainties of Arctic snow-on-sea-ice

10.5194/egusphere-egu2020-11436 ◽

2020 ◽

Author(s):

Alex Cabaj ◽

Paul Kushner ◽

Alek Petty ◽

Stephen Howell ◽

Christopher Fletcher

Keyword(s):

Sea Ice ◽

Arctic Ocean ◽

Snow Depth ◽

Arctic Sea Ice ◽

The Arctic ◽

Ice Thickness ◽

Sea Ice Thickness ◽

The Arctic Ocean ◽

Ice Model

Snow on Arctic sea ice plays multiple&#8212;and sometimes contrasting&#8212;roles in several feedbacks between sea ice and the global climate system. For example, the presence of snow on sea ice may mitigate sea ice melt by increasing the sea ice albedo and enhancing the ice-albedo feedback. Conversely, snow can inhibit sea ice growth by insulating the ice from the atmosphere during the sea ice growth season. In addition to its contribution to sea ice feedbacks, snow on sea ice also poses a challenge for sea ice observations. In particular, snow contributes to uncertainties in retrievals of sea ice thickness from satellite altimetry measurements, such as those from ICESat-2. Snow-on-sea-ice models can produce basin-wide snow depth estimates, but these models require snowfall input from reanalysis products. In-situ snowfall measurements are absent over most of the Arctic Ocean, so it can be difficult to determine which reanalysis snowfall product is best suited to be used as input for a snow-on-sea-ice model.In the absence of in-situ snowfall rate measurements, measurements from satellite instruments can be used to quantify snowfall over the Arctic Ocean. The CloudSat satellite, which is equipped with a 94 GHz Cloud Profiling Radar instrument, measures vertical radar reflectivity profiles from which snowfall rates can be retrieved. This instrument provides the most extensive high-latitude snowfall rate observation dataset currently available. CloudSat&#8217;s near-polar orbit enables it to make measurements at latitudes up to 82&#176;N, with a 16-day repeat cycle, over the time period from 2006-2016.We present a calibration of reanalysis snowfall to CloudSat observations over the Arctic Ocean, which we then apply to reanalysis snowfall input for the NASA Eulerian Snow On Sea Ice Model (NESOSIM). This calibration reduces the spread in snow depths produced by NESOSIM when different reanalysis inputs are used. In light of this calibration, we revise the NESOSIM parametrizations of wind-driven snow processes, and we characterize the uncertainties in NESOSIM-generated snow depths resulting from uncertainties in snowfall input. We then extend this analysis further to estimate the resulting uncertainties in sea ice thickness retrieved from ICESat-2 when snow depth estimates from NESOSIM are used as input for the retrieval.

Download Full-text

Parameter calibration and uncertainty analysis for snow depths from the NASA Eulerian Snow On Sea Ice Model and derived sea ice thickness from ICESat-2

10.5194/egusphere-egu21-13791 ◽

2021 ◽

Author(s):

Alex Cabaj ◽

Paul Kushner ◽

Alek Petty

Keyword(s):

Sea Ice ◽

Snow Depth ◽

Arctic Sea Ice ◽

The Arctic ◽

Ice Thickness ◽

Blowing Snow ◽

Sea Ice Thickness ◽

Ice Growth ◽

Arctic Sea ◽

Ice Model

Snow on Arctic sea ice plays many, sometimes contrasting roles in Arctic climate feedbacks. During the sea ice growth season, the presence of snow on sea ice can enhance ice growth by increasing the sea ice albedo, or conversely, inhibit sea ice growth by insulating the ice from the cold atmosphere. Furthermore, estimates of snow depth on Arctic sea ice are also a key input for deriving sea ice thickness from altimetry measurements, such as satellite lidar altimetry measurements from ICESat-2. Due to the logistical challenges of making measurements in as remote a region as the Arctic, snow depth on Arctic sea ice is difficult to observationally constrain. The NASA Eulerian Snow On Sea Ice Model (NESOSIM) can be used to provide snow depth and density estimates over Arctic sea ice with pan-Arctic coverage within a relatively simple framework. The latest version of NESOSIM, version 1.1, is a 2-layer model with simple representations of the processes of accumulation, wind packing, loss due to blowing snow, and redistribution due to sea ice motion. Relative to version 1.0, NESOSIM 1.1 features an extended model domain, and reanalysis snowfall input scaled to observed snowfall retrieved from CloudSat satellite radar reflectivity measurements. In this work, we present a systematic calibration, and an accompanying estimate in the uncertainty of the free parameters in NESOSIM, targeting airborne snow radar measurements from Operation IceBridge. We further investigate uncertainties in snow depth and the resulting uncertainties in derived sea ice thickness from ICESat-2 altimetry measurements using NESOSIM snow depths.

Download Full-text

The implications of selected processing methods on satellite altimetry derived sea ice thickness state and trends in the seasonal ice zone

10.5194/egusphere-egu21-12240 ◽

2021 ◽

Author(s):

Isolde Glissenaar ◽

Jack Landy ◽

Alek Petty ◽

Nathan Kurtz ◽

Julienne Stroeve

Keyword(s):

Sea Ice ◽

Snow Depth ◽

Satellite Altimetry ◽

Region Of Interest ◽

The Arctic ◽

Ice Thickness ◽

Sea Ice Thickness ◽

Baffin Bay ◽

The Impact

The ice cover of the Arctic Ocean is increasingly becoming dominated by seasonal sea ice. It is important to focus on the processing of altimetry ice thickness data in thinner seasonal ice regions to understand seasonal sea ice behaviour better. This study focusses on Baffin Bay as a region of interest to study seasonal ice behaviour.We aim to reconcile the spring sea ice thickness derived from multiple satellite altimetry sensors and sea ice charts in Baffin Bay and produce a robust long-term record (2003-2020) for analysing trends in sea ice thickness. We investigate the impact of choosing different snow depth products (the Warren climatology, a passive microwave snow depth product and modelled snow depth from reanalysis data) and snow redistribution methods (a sigmoidal function and an empirical piecewise function) to retrieve sea ice thickness from satellite altimetry sea ice freeboard data.The choice of snow depth product and redistribution method results in an uncertainty envelope around the March mean sea ice thickness in Baffin Bay of 10%. Moreover, the sea ice thickness trend ranges from -15 cm/dec to 20 cm/dec depending on the applied snow depth product and redistribution method. Previous studies have shown a possible long-term asymmetrical trend in sea ice thinning in Baffin Bay. The present study shows that whether a significant long-term asymmetrical trend was found depends on the choice of snow depth product and redistribution method. The satellite altimetry sea ice thickness results with different snow depth products and snow redistribution methods show that different processing techniques can lead to different results and can influence conclusions on total and spatial sea ice thickness trends. Further processing work on the historic radar altimetry record is needed to create reliable sea ice thickness products in the marginal ice zone.

Download Full-text

Near Real Time Arctic sea ice thickness and volume from CryoSat-2

10.5194/tc-2016-21 ◽

2016 ◽

Cited By ~ 3

Author(s):

R. L. Tilling ◽

A. Ridout ◽

A. Shepherd

Keyword(s):

Sea Ice ◽

Real Time ◽

Arctic Sea Ice ◽

Arctic Climate ◽

The Arctic ◽

Ice Thickness ◽

Satellite Orbits ◽

Sea Ice Thickness ◽

Arctic Sea ◽

Arctic Sea Ice Thickness

Abstract. Timely observations of sea ice thickness help us to understand Arctic climate, and can support maritime activities in the Polar Regions. Although it is possible to calculate Arctic sea ice thickness using measurements acquired by CryoSat-2, the latency of the final release dataset is typically one month, due to the time required to determine precise satellite orbits. We use a new fast delivery CryoSat-2 dataset based on preliminary orbits to compute Arctic sea ice thickness in near real time (NRT), and analyse this data for one sea ice growth season from October 2014 to April 2015. We show that this NRT sea ice thickness product is of comparable accuracy to that produced using the final release CryoSat-2 data, with an average thickness difference of 5 cm, demonstrating that the satellite orbit is not a critical factor in determining sea ice freeboard. In addition, the CryoSat-2 fast delivery product also provides measurements of Arctic sea ice thickness within three days of acquisition by the satellite, and a measurement is delivered, on average, within 10, 7 and 6 km of each location in the Arctic every 2, 14 and 28 days respectively. The CryoSat-2 NRT sea ice thickness dataset provides an additional constraint for seasonal predictions of Arctic climate change, and will allow industries such as tourism and transport to navigate the polar oceans with safety and care.

Download Full-text

New insight from CryoSat-2 sea ice thickness for sea ice modelling

10.5194/tc-2018-159 ◽

2018 ◽

Author(s):

David Schröder ◽

Danny L. Feltham ◽

Michel Tsamados ◽

Andy Ridout ◽

Rachel Tilling

Keyword(s):

Sea Ice ◽

Thickness Distribution ◽

Arctic Sea Ice ◽

Ice Thickness ◽

Atmospheric Conditions ◽

Sea Ice Concentration ◽

Sea Ice Thickness ◽

Ice Growth ◽

Melt Pond ◽

The Impact

Abstract. Estimates of Arctic sea ice thickness are available from the CryoSat-2 (CS2) radar altimetry mission during ice growth seasons since 2010. We derive the sub-grid scale ice thickness distribution (ITD) with respect to 5 ice thickness categories used in a sea ice component (CICE) of climate simulations. This allows us to initialize the ITD in stand-alone simulations with CICE and to verify the simulated cycle of ice thickness. We find that a default CICE simulation strongly underestimates ice thickness, despite reproducing the inter-annual variability of summer sea ice extent. We can identify the underestimation of winter ice growth as being responsible and show that increasing the ice conductive flux for lower temperatures (bubbly brine scheme) and accounting for the loss of drifting snow results in the simulated sea ice growth being more realistic. Sensitivity studies provide insight into the impact of initial and atmospheric conditions and, thus, on the role of positive and negative feedback processes. During summer, atmospheric conditions are responsible for 50 % of September sea ice thickness variability through the positive sea ice and melt pond albedo feedback. However, atmospheric winter conditions have little impact on winter ice growth due to the dominating negative conductive feedback process: the thinner the ice and snow in autumn, the stronger the ice growth in winter. We conclude that the fate of Arctic summer sea ice is largely controlled by atmospheric conditions during the melting season rather than by winter temperature. Our optimal model configuration does not only improve the simulated sea ice thickness, but also summer sea ice concentration, melt pond fraction, and length of the melt season. It is the first time CS2 sea ice thickness data have been applied successfully to improve sea ice model physics.

Download Full-text

Investigating the variability of the Arctic sea ice thickness in response to a stochastic thermodynamic atmospheric forcing

Climate Dynamics ◽

10.1007/s003820000096 ◽

2001 ◽

Vol 17 (2-3) ◽

pp. 107-125 ◽

Cited By ~ 14

Author(s):

B. L'Hévéder ◽

M.-N. Houssais

Keyword(s):

Sea Ice ◽

Arctic Sea Ice ◽

The Arctic ◽

Ice Thickness ◽

Atmospheric Forcing ◽

Sea Ice Thickness ◽

Arctic Sea ◽

Arctic Sea Ice Thickness

Download Full-text

Winter Arctic sea ice freeboard, snow depth and thickness variability from ICESat-2 and NESOSIM

10.5194/egusphere-egu21-13779 ◽

2021 ◽

Author(s):

Alek Petty ◽

Nicole Keeney ◽

Alex Cabaj ◽

Paul Kushner ◽

Nathan Kurtz ◽

...

Keyword(s):

Sea Ice ◽

Snow Depth ◽

Arctic Sea Ice ◽

The Arctic ◽

Ice Thickness ◽

Sea Ice Thickness ◽

Ice Cloud ◽

Ice Drift ◽

Thickness Variability ◽

Arctic Ice

<div> <div> <div> <div> National Aeronautics and Space Administration's (NASA's) Ice, Cloud, and land Elevation Satellite&#8208; 2 (ICESat&#8208;2) mission was launched in September 2018 and is now providing routine, very high&#8208;resolution estimates of surface height/type (the ATL07 product) and freeboard (the ATL10 product) across the Arctic and Southern Oceans. In recent work we used snow depth and density estimates from the NASA Eulerian Snow on Sea Ice Model (NESOSIM) together with ATL10 freeboard data to estimate sea ice thickness across the entire Arctic Ocean. Here we provide an overview of updates made to both the underlying ATL10 freeboard product and the NESOSIM model, and the subsequent impacts on our estimates of sea ice thickness including updated comparisons to the original ICESat mission and ESA&#8217;s CryoSat-2. Finally we compare our Arctic ice thickness estimates from the 2018-2019 and 2019-2020 winters and discuss possible causes of these differences based on an analysis of atmospheric data (ERA5), ice drift (NSIDC) and ice type (OSI SAF). </div> </div> </div> </div>

Download Full-text

Estimating Arctic Sea Ice Thickness with CryoSat-2 Altimetry Data Using the Least Squares Adjustment Method

Sensors ◽

10.3390/s20247011 ◽

2020 ◽

Vol 20 (24) ◽

pp. 7011

Author(s):

Feng Xiao ◽

Fei Li ◽

Shengkai Zhang ◽

Jiaxing Li ◽

Tong Geng ◽

...

Keyword(s):

Sea Ice ◽

Least Squares ◽

Arctic Sea Ice ◽

The Arctic ◽

Ice Thickness ◽

Altimetry Data ◽

Sea Ice Thickness ◽

Arctic Sea ◽

Least Squares Adjustment ◽

Arctic Sea Ice Thickness

Satellite altimeters can be used to derive long-term and large-scale sea ice thickness changes. Sea ice thickness retrieval is based on measurements of freeboard, and the conversion of freeboard to thickness requires knowledge of the snow depth and snow, sea ice, and sea water densities. However, these parameters are difficult to be observed concurrently with altimeter measurements. The uncertainties in these parameters inevitably cause uncertainties in sea ice thickness estimations. This paper introduces a new method based on least squares adjustment (LSA) to estimate Arctic sea ice thickness with CryoSat-2 measurements. A model between the sea ice freeboard and thickness is established within a 5 km × 5 km grid, and the model coefficients and sea ice thickness are calculated using the LSA method. Based on the newly developed method, we are able to derive estimates of the Arctic sea ice thickness for 2010 through 2019 using CryoSat-2 altimetry data. Spatial and temporal variations of the Arctic sea ice thickness are analyzed, and comparisons between sea ice thickness estimates using the LSA method and three CryoSat-2 sea ice thickness products (Alfred Wegener Institute (AWI), Centre for Polar Observation and Modelling (CPOM), and NASA Goddard Space Flight Centre (GSFC)) are performed for the 2018–2019 Arctic sea ice growth season. The overall differences of sea ice thickness estimated in this study between AWI, CPOM, and GSFC are 0.025 ± 0.640 m, 0.143 ± 0.640 m, and −0.274 ± 0.628 m, respectively. Large differences between the LSA and three products tend to appear in areas covered with thin ice due to the limited accuracy of CryoSat-2 over thin ice. Spatiotemporally coincident Operation IceBridge (OIB) thickness values are also used for validation. Good agreement with a difference of 0.065 ± 0.187 m is found between our estimates and the OIB results.

Download Full-text