Prediction of monthly Arctic sea ice concentrations using satellite and reanalysis data based on convolutional neural networks

Changes in Arctic sea ice affect atmospheric circulation, ocean current, and polar ecosystems. There have been unprecedented decreases in the amount of Arctic sea ice due to global warming. In this study, a novel 1-month sea ice concentration (SIC) prediction model is proposed, with eight predictors using a deep-learning approach, convolutional neural networks (CNNs). This monthly SIC prediction model based on CNNs is shown to perform better predictions (mean absolute error - MAE - of 2.28 %, anomaly correlation coefficient - ACC - of 0.98, root-mean-square error - RMSE - of 5.76 %, normalized RMSE - nRMSE - of 16.15 %, and NSE - Nash-Sutcliffe efficiency - of 0.97) than a random-forest-based (RF-based) model (MAE of 2.45 %, ACC of 0.98, RMSE of 6.61 %, nRMSE of 18.64 %, and NSE of 0.96) and the persistence model based on the monthly trend (MAE of 4.31 %, ACC of 0.95, RMSE of 10.54 %, nRMSE of 29.17 %, and NSE of 0.89) through hindcast validations. The spatio-temporal analysis also confirmed the superiority of the CNN model. The CNN model showed good SIC prediction results in extreme cases that recorded unforeseen sea ice plummets in 2007 and 2012 with RMSEs of less than 5.0 %. This study also examined the importance of the input variables through a sensitivity analysis. In both the CNN and RF models, the variables of past SICs were identified as the most sensitive factor in predicting SICs. For both models, the SIC-related variables generally contributed more to predict SICs over ice-covered areas, while other meteorological and oceanographic variables were more sensitive to the prediction of SICs in marginal ice zones. The proposed 1-month SIC prediction model provides valuable information which can be used in various applications, such as Arctic shipping-route planning, management of the fishing industry, and long-term sea ice forecasting and dynamics

An assessment of key model parametric uncertainties in projections of Greenland Ice Sheet behavior

Lack of knowledge about the values of ice sheet model input parameters introduces substantial uncertainty into projections of Greenland Ice Sheet contributions to future sea level rise. Computer models of ice sheet behavior provide one of several means of estimating future sea level rise due to mass loss from ice sheets. Such models have many input parameters whose values are not well known. Recent studies have investigated the effects of these parameters on model output, but the range of potential future sea level increases due to model parametric uncertainty has not been characterized. Here, we demonstrate that this range is large, using a 100-member perturbed-physics ensemble with the SICOPOLIS ice sheet model. Each model run is spun up over 125 000 yr using geological forcings and subsequently driven into the future using an asymptotically increasing air temperature anomaly curve. All modeled ice sheets lose mass after 2005 AD. Parameters controlling surface melt dominate the model response to temperature change. After culling the ensemble to include only members that give reasonable ice volumes in 2005 AD, the range of projected sea level rise values in 2100 AD is ~40 % or more of the median. Data on past ice sheet behavior can help reduce this uncertainty, but none of our ensemble members produces a reasonable ice volume change during the mid-Holocene, relative to the present. This problem suggests that the model's exponential relation between temperature and precipitation does not hold during the Holocene, or that the central-Greenland temperature forcing curve used to drive the model is not representative of conditions around the ice margin at this time (among other possibilities). Our simulations also lack certain observed physical processes that may tend to enhance the real ice sheet's response. Regardless, this work has implications for other studies that use ice sheet models to project or hindcast the behavior of the Greenland Ice Sheet

Multi-physics ensemble snow modelling in the western Himalaya

Combining multiple data sources with multi-physics simulation frameworks offers new potential to extend snow model inter-comparison efforts to the Himalaya. As such, this study evaluates the sensitivity of simulated regional snow cover and runoff dynamics to different snowpack process representations. The evaluation is based on a spatially distributed version of the Factorial Snowpack Model (FSM) set up for the Astore catchment in the upper Indus basin. The FSM multi-physics model was driven by climate fields from the High Asia Refined Analysis (HAR) dynamical downscaling product. Ensemble performance was evaluated primarily using MODIS remote sensing of snow-covered area, albedo and land surface temperature. In line with previous snow model inter-comparisons, no single FSM configuration performs best in all of the years simulated. However, the results demonstrate that performance variation in this case is at least partly related to inaccuracies in the sequencing of inter-annual variation in HAR climate inputs, not just FSM model limitations. Ensemble spread is dominated by interactions between parameterisations of albedo, snowpack hydrology and atmospheric stability effects on turbulent heat fluxes. The resulting ensemble structure is similar in different years, which leads to systematic divergence in ablation and mass balance at high elevations. While ensemble spread and errors are notably lower when viewed as anomalies, FSM configurations show important differences in their absolute sensitivity to climate variation. Comparison with observations suggests that a subset of the ensemble should be retained for climate change projections, namely those members including prognostic albedo and liquid water retention, refreezing and drainage processes

Cloud and ice in the planetary scale circulation and in climate

The roles of the cryosphere, and of cloud-radiative interactions are investigated. The effects clouds and ice have in the climate system are examined. The cloud radiation research attempts explain the modes of interaction (feedback) between raditive transfer, cloud formation, and atmospheric dynamics. The role of sea ice in weather and climate is also discussed. Models are used to describe the ice and atmospheric dynamics under study

Representing moisture fluxes and phase changes in glacier debris cover using a reservoir approach

Due to the complexity of treating moisture in supraglacial debris, surface energy balance models to date have neglected moisture infiltration and phase changes in the debris layer. The latent heat flux (QL) is also often excluded due to the uncertainty in determining the surface vapour pressure. To quantify the importance of moisture on the surface energy and climatic mass balance (CMB) of debris-covered glaciers, we developed a simple reservoir parameterization for the debris ice and water content, as well as an estimation of the latent heat flux. The parameterization was incorporated into a CMB model adapted for debris-covered glaciers. We present the results of two point simulations, using both our new “moist” and the conventional “dry” approaches, on the Miage Glacier, Italy, during summer 2008 and fall 2011. The former year coincides with available in situ glaciological and meteorological measurements, including the first eddy-covariance measurements of the turbulent fluxes over supraglacial debris, while the latter contains two refreeze events that permit evaluation of the influence of phase changes. The simulations demonstrate a clear influence of moisture on the glacier energy and mass-balance dynamics. When water and ice are considered, heat transmission to the underlying glacier ice is lower, as the effective thermal diffusivity of the saturated debris layers is reduced by increases in both the density and the specific heat capacity of the layers. In combination with surface heat extraction by QL, subdebris ice melt is reduced by 3.1% in 2008 and by 7.0% in 2011 when moisture effects are included. However, the influence of the parameterization on the total accumulated mass balance varies seasonally. In summer 2008, mass loss due to surface vapour fluxes more than compensates for the reduction in ice melt, such that the total ablation increases by 4.0 %. Conversely, in fall 2011, the modulation of basal debris temperature by debris ice results in a decrease in total ablation of 2.1 %. Although the parameterization is a simplified representation of the moist physics of glacier debris, it is a novel attempt at including moisture in a numerical model of debris-covered glaciers and one that opens up additional avenues for future research

Estimation of subsurface porosities and thermal conductivities of polygonal tundra by coupled inversion of electrical resistivity, temperature, and moisture content data

Studies indicate greenhouse gas emissions following permafrost thaw will amplify current rates of atmospheric warming, a process referred to as the permafrost carbon feedback. However, large uncertainties exist regarding the timing and magnitude of the permafrost carbon feedback, in part due to uncertainties associated with subsurface permafrost parameterization and structure. Development of robust parameter estimation methods for permafrost-rich soils is becoming urgent under accelerated warming of the Arctic. Improved parameterization of the subsurface properties in land system models would lead to improved predictions and a reduction of modeling uncertainty. In this work we set the groundwork for future parameter estimation (PE) studies by developing and evaluating a joint PE algorithm that estimates soil porosities and thermal conductivities from time series of soil temperature and moisture measurements and discrete in-time electrical resistivity measurements. The algorithm utilizes the Model-Independent Parameter Estimation and Uncertainty Analysis toolbox and coupled hydrological-thermal-geophysical modeling. We test the PE algorithm against synthetic data, providing a proof of concept for the approach. We use specified subsurface porosities and thermal conductivities and coupled models to set up a synthetic state, perturb the parameters, and then verify that our PE method is able to recover the parameters and synthetic state. To evaluate the accuracy and robustness of the approach we perform multiple tests for a perturbed set of initial starting parameter combinations. In addition, we varied types and quantities of data to better understand the optimal dataset needed to improve the PE method. The results of the PE tests suggest that using multiple types of data improve the overall robustness of the method. Our numerical experiments indicate that special care needs to be taken during the field experiment setup so that (1) the vertical distance between adjacent measurement sensors allows the signal variability in space to be resolved and (2) the longer time interval between resistivity snapshots allows signal variability in time to be resolved

Fluctuations of a Greenlandic tidewater glacier driven by changes in atmospheric forcing : observations and modelling of Kangiata Nunaata Sermia, 1859–present

Acknowledgements. The authors wish to thank Stephen Price, Mauri Pelto, and the anonymous reviewer for their reviews and comments that helped to improve the manuscript. RACMO2.1 data were provided by Jan van Angelen and Michiel van den Broeke, IMAU, Utrecht University. MAR v3.2 data used for runoff calculations were provided by Xavier Fettweis, Department of Geography, University of Liège. The photogrammetric DEM used in Figs. 1 and 3 was provided by Kurt H. Kjær, Centre for GeoGenetics, University of Copenhagen. This research was financially supported by J. M. Lea’s PhD funding, NERC grant number NE/I528742/1. Support for F. M. Nick was provided through the Conoco-Phillips/Lundin Northern Area Program CRIOS project (Calving Rates and Impact on Sea Level).Peer reviewedPublisher PD

Impact of rheology on probabilistic forecasts of sea ice trajectories: application for search and rescue operations in the Arctic

We present a sensitivity analysis and discuss the probabilistic forecast capabilities of the novel sea ice model neXtSIM used in hindcast mode. The study pertains to the response of the model to the uncertainty on winds using probabilistic forecasts of ice trajectories. neXtSIM is a continuous Lagrangian numerical model that uses an elasto-brittle rheology to simulate the ice response to external forces. The sensitivity analysis is based on a Monte Carlo sampling of 12 members. The response of the model to the uncertainties is evaluated in terms of simulated ice drift distances from their initial positions, and from the mean position of the ensemble, over the mid-term forecast horizon of 10 days. The simulated ice drift is decomposed into advective and diffusive parts that are characterised separately both spatially and temporally and compared to what is obtained with a free-drift model, that is, when the ice rheology does not play any role in the modelled physics of the ice. The seasonal variability of the model sensitivity is presented and shows the role of the ice compactness and rheology in the ice drift response at both local and regional scales in the Arctic. Indeed, the ice drift simulated by neXtSIM in summer is close to the one obtained with the free-drift model, while the more compact and solid ice pack shows a significantly different mechanical and drift behaviour in winter. For the winter period analysed in this study, we also show that, in contrast to the free-drift model, neXtSIM reproduces the sea ice Lagrangian diffusion regimes as found from observed trajectories. The forecast capability of neXtSIM is also evaluated using a large set of real buoy's trajectories and compared to the capability of the free-drift model. We found that neXtSIM performs significantly better in simulating sea ice drift, both in terms of forecast error and as a tool to assist search and rescue operations, although the sources of uncertainties assumed for the present experiment are not sufficient for complete coverage of the observed IABP positions

A double continuum hydrological model for glacier applications

Peer reviewedPublisher PD