Optimization of a sea ice model using basinwide observations of Arctic sea ice thickness, extent, and velocity

A stand-alone sea ice model is tuned and validated using satellite-derived, basinwide observations of sea ice thickness, extent, and velocity from the years 1993 to 2001. This is the first time that basin-scale measurements of sea ice thickness have been used for this purpose. The model is based on the CICE sea ice model code developed at the Los Alamos National Laboratory, with some minor modifications, and forcing consists of 40-yr ECMWF Re-Analysis (ERA-40) and Polar Exchange at the Sea Surface (POLES) data. Three parameters are varied in the tuning process: Ca, the air–ice drag coefficient; P*, the ice strength parameter; and α, the broadband albedo of cold bare ice, with the aim being to determine the subset of this three-dimensional parameter space that gives the best simultaneous agreement with observations with this forcing set. It is found that observations of sea ice extent and velocity alone are not sufficient to unambiguously tune the model, and that sea ice thickness measurements are necessary to locate a unique subset of parameter space in which simultaneous agreement is achieved with all three observational datasets

Recent loss of floating ice and the consequent sea level contribution

We combine new and published satellite observations and the results of a coupled ice-ocean model to provide the first estimate of changes in the quantity of ice floating in the global oceans and the consequent sea level contribution. Rapid losses of Arctic sea ice and small Antarctic ice shelves are partially offset by thickening of Antarctic sea ice and large Antarctic ice shelves. Altogether, 746 +/- 127 km(3) yr(-1) of floating ice was lost between 1994 and 2004, a value that exceeds considerably the reduction in grounded ice over the same period. Although the losses are equivalent to a small (49 +/- 8 μm yr(-1)) rise in mean sea level, there may be large regional variations in the degree of ocean freshening and mixing. Ice shelves at the Antarctic Peninsula and in the Amundsen Sea, for example, have lost 481 +/- 38 km(3) yr(-1)

Improving the spatial distribution of modeled Arctic sea ice thickness

The spatial distribution of ice thickness/draft in the Arctic Ocean is examined using a sea ice model. A comparison of model predictions with submarine observations of sea ice draft made during cruises between 1987 and 1997 reveals that the model has the same deficiencies found in previous studies, namely ice that is too thick in the Beaufort Sea and too thin near the North Pole. We find that increasing the large scale shear strength of the sea ice leads to substantial improvements in the model's spatial distribution of sea ice thickness, and simultaneously improves the agreement between modeled and ERS-derived 1993–2001 mean winter ice thickness

Consistent and contrasting decadal Arctic sea ice thickness predictions from a highly optimized sea ice model

[1] Decadal hindcast simulations of Arctic Ocean sea ice thickness made by a modern dynamic-thermodynamic sea ice model and forced independently by both the ERA-40 and NCEP/NCAR reanalysis data sets are compared for the first time. Using comprehensive data sets of observations made between 1979 and 2001 of sea ice thickness, draft, extent, and speeds, we find that it is possible to tune model parameters to give satisfactory agreement with observed data, thereby highlighting the skill of modern sea ice models, though the parameter values chosen differ according to the model forcing used. We find a consistent decreasing trend in Arctic Ocean sea ice thickness since 1979, and a steady decline in the Eastern Arctic Ocean over the full 40-year period of comparison that accelerated after 1980, but the predictions of Western Arctic Ocean sea ice thickness between 1962 and 1980 differ substantially. The origins of differing thickness trends and variability were isolated not to parameter differences but to differences in the forcing fields applied, and in how they are applied. It is argued that uncertainty, differences and errors in sea ice model forcing sets complicate the use of models to determine the exact causes of the recently reported decline in Arctic sea ice thickness, but help in the determination of robust features if the models are tuned appropriately against observations

A multithickness sea ice model accounting for sliding friction

A multithickness sea ice model explicitly accounting for the ridging and sliding friction contributions to sea ice stress is developed. Both ridging and sliding contributions depend on the deformation type through functions adopted from the Ukita and Moritz kinematic model of floe interaction. In contrast to most previous work, the ice strength of a uniform ice sheet of constant ice thickness is taken to be proportional to the ice thickness raised to the 3/2 power, as is revealed in discrete element simulations by Hopkins. The new multithickness sea ice model for sea ice stress has been implemented into the Los Alamos “CICE” sea ice model code and is shown to improve agreement between model predictions and observed spatial distribution of sea ice thickness in the Arctic

The long-term and interannual variability of summer fresh water storage over the eastern Siberian shelf: Implication for climatic change

A time series of summer fresh water content anomalies (FWCA) over the Laptev and East Siberian sea shelves was constructed from historical hydrographic records for the period from 1920 to 2005. Results from a multiple regression between FCWA and various atmospheric and oceanic indices show that the fresh water content on the shelves is mainly controlled by atmospheric vorticity on quasi-decadal timescales. When the vorticity of the atmosphere on the shelves is antycyclonic, approximately 500 km3 of fresh water migrates from the eastern Siberian shelf to the Arctic Ocean through the northeastern Laptev Sea. When the vorticity of the atmosphere is cyclonic, this fresh water remains on the southern Laptev and East Siberian sea shelves. This FWCA represents approximately 35% of the total fresh water inflow provided by river discharge and local sea-ice melt, and is about ten times larger than the standard deviation of the Lena River summer long-term mean discharge. However, the large interannual and spatial variability in the fresh water content of the shelves, as well as the spatial coverage of the hydrographic data, makes it difficult to detect the long-term tendency of fresh water storage associated with climate change

Utilizing CryoSat-2 sea ice thickness to initialize a coupled ice-ocean modeling system

Two CryoSat-2 sea ice thickness products derived with independent algorithms are used to initialize a coupled ice-ocean modeling system in which a series of reanalysis studies are performed for the period of March 15, 2014–September 30, 2015. Comparisons against moored upward looking sonar, drifting ice mass balance buoy, and NASA Operation IceBridge ice thickness data show that the modeling system exhibits greatly reduced bias using the satellite-derived ice thickness data versus the operational model run without these data. The model initialized with CryoSat-2 ice thickness exhibits skill in simulating ice thickness from the initial period to up to 6 months. We find that the largest improvements in ice thickness occur over multi-year ice. Based on the data periods examined here, we find that for the 18-month study period, when compared with upward looking sonar measurements, the CryoSat-2 reanalyses show significant improvement in bias (0.47–0.75) and RMSE (0.89–1.04) versus the control run without these data (1.44 and 1.60, respectively). An ice drift comparison reveals little change in ice velocity statistics for the Pan Arctic region; however some improvement is seen during the summer/autumn months in 2014 for the Bering/Beaufort/Chukchi and Greenland/Norwegian Seas. These promising results suggest that such a technique should be used to reinitialize operational sea ice modeling systems

Sensitivity of simulated regional Arctic climate to the choice of coupled model domain

The climate over the Arctic has undergone changes in recent decades. In order to evaluate the coupled response of the Arctic system to external and internal forcing, our study focuses on the estimation of regional climate variability and its dependence on large-scale atmospheric and regional ocean circulations. A global ocean–sea ice model with regionally high horizontal resolution is coupled to an atmospheric regional model and global terrestrial hydrology model. This way of coupling divides the global ocean model setup into two different domains: one coupled, where the ocean and the atmosphere are interacting, and one uncoupled, where the ocean model is driven by prescribed atmospheric forcing and runs in a so-called stand-alone mode. Therefore, selecting a specific area for the regional atmosphere implies that the ocean–atmosphere system can develop ‘freely’ in that area, whereas for the rest of the global ocean, the circulation is driven by prescribed atmospheric forcing without any feedbacks. Five different coupled setups are chosen for ensemble simulations. The choice of the coupled domains was done to estimate the influences of the Subtropical Atlantic, Eurasian and North Pacific regions on northern North Atlantic and Arctic climate. Our simulations show that the regional coupled ocean–atmosphere model is sensitive to the choice of the modelled area. The different model configurations reproduce differently both the mean climate and its variability. Only two out of five model setups were able to reproduce the Arctic climate as observed under recent climate conditions (ERA-40 Reanalysis). Evidence is found that the main source of uncertainty for Arctic climate variability and its predictability is the North Pacific. The prescription of North Pacific conditions in the regional model leads to significant correlation with observations, even if the whole North Atlantic is within the coupled model domain. However, the inclusion of the North Pacific area into the coupled system drastically changes the Arctic climate variability to a point where the Arctic Oscillation becomes an ‘internal mode’ of variability and correlations of year-to-year variability with observational data vanish. In line with previous studies, our simulations provide evidence that Arctic sea ice export is mainly due to ‘internal variability’ within the Arctic region. We conclude that the choice of model domains should be based on physical knowledge of the atmospheric and oceanic processes and not on ‘geographic’ reasons. This is particularly the case for areas like the Arctic, which has very complex feedbacks between components of the regional climate system