Pan-Arctic Ice Ocean Modeling And Assimilation System Research Articles

Landfast ice thickness in the Canadian Arctic Archipelago from observations and models

Abstract. Observed and modelled landfast ice thickness variability and trends spanning more than 5 decades within the Canadian Arctic Archipelago (CAA) are summarized. The observed sites (Cambridge Bay, Resolute, Eureka and Alert) represent some of the Arctic's longest records of landfast ice thickness. Observed end-of-winter (maximum) trends of landfast ice thickness (1957–2014) were statistically significant at Cambridge Bay (−4.31 ± 1.4 cm decade−1), Eureka (−4.65 ± 1.7 cm decade−1) and Alert (−4.44 ± 1.6 cm −1) but not at Resolute. Over the 50+-year record, the ice thinned by ∼ 0.24–0.26 m at Cambridge Bay, Eureka and Alert with essentially negligible change at Resolute. Although statistically significant warming in spring and fall was present at all sites, only low correlations between temperature and maximum ice thickness were present; snow depth was found to be more strongly associated with the negative ice thickness trends. Comparison with multi-model simulations from Coupled Model Intercomparison project phase 5 (CMIP5), Ocean Reanalysis Intercomparison (ORA-IP) and Pan-Arctic Ice-Ocean Modeling and Assimilation System (PIOMAS) show that although a subset of current generation models have a "reasonable" climatological representation of landfast ice thickness and distribution within the CAA, trends are unrealistic and far exceed observations by up to 2 orders of magnitude. ORA-IP models were found to have positive correlations between temperature and ice thickness over the CAA, a feature that is inconsistent with both observations and coupled models from CMIP5.

Modeling the seasonal evolution of the Arctic sea ice floe size distribution

Abstract To better simulate the seasonal evolution of sea ice in the Arctic, with particular attention to the marginal ice zone, a sea ice model of the distribution of ice thickness, floe size, and enthalpy was implemented into the Pan-arctic Ice–Ocean Modeling and Assimilation System (PIOMAS). Theories on floe size distribution (FSD) and ice thickness distribution (ITD) were coupled in order to explicitly simulate multicategory FSD and ITD distributions simultaneously. The expanded PIOMAS was then used to estimate the seasonal evolution of the Arctic FSD in 2014 when FSD observations are available for model calibration and validation. Results indicate that the simulated FSD, commonly described equivalently as cumulative floe number distribution (CFND), generally follows a power law across space and time and agrees with the CFND observations derived from TerraSAR-X satellite images. The simulated power-law exponents also correlate with those derived using MODIS images, with a low mean bias of –2%. In the marginal ice zone, the modeled CFND shows a large number of small floes in winter because of stronger winds acting on thin, weak first-year ice in the ice edge region. In mid-spring and summer, the CFND resembles an upper truncated power law, with the largest floes mostly broken into smaller ones; however, the number of small floes is lower than in winter because floes of small sizes or first-year ice are easily melted away. In the ice pack interior there are fewer floes in late fall and winter than in summer because many of the floes are “welded” together into larger floes in freezing conditions, leading to a relatively flat CFND with low power-law exponents. The simulated mean floe size averaged over all ice-covered areas shows a clear annual cycle, large in winter and smaller in summer. However, there is no obvious annual cycle of mean floe size averaged over the marginal ice zone. The incorporation of FSD into PIOMAS results in reduced ice thickness, mainly in the marginal ice zone, which improves the simulation of ice extent and yields an earlier ice retreat.

Improving Arctic Sea Ice Prediction Using PIOMAS Initial Sea Ice Thickness in a Coupled Ocean–Atmosphere Model

Abstract Because sea ice thickness is known to influence future patterns of sea ice concentration, it is likely that an improved initialization of sea ice thickness in a coupled ocean–atmosphere model would improve Arctic sea ice cover forecasts. Here, two sea ice thickness datasets as possible candidates for forecast initialization were investigated: the Climate Forecast System Reanalysis (CFSR) and the Pan-Arctic Ice Ocean Modeling and Assimilation System (PIOMAS). Using Ice, Cloud, and Land Elevation Satellite (ICESat) data, it was shown that the PIOMAS dataset had a more realistic representation of sea ice thickness than CFSR. Subsequently, both March CFSR and PIOMAS sea ice thicknesses were used to initialize hindcasts using the Climate Forecast System, version 2 (CFSv2), model. A second set of model runs was also done in which the original model physics were modified to more physically reasonable settings—namely, increasing the number of marine stratus clouds in the Arctic region and enabling realistic representation of the ice–ocean heat flux. Hindcasts were evaluated using sea ice concentration observations from the National Aeronautics and Space Administration (NASA) Team and Bootstrap algorithms. Results show that using PIOMAS initial sea ice thickness in addition to the physics modifications yielded significant improvement in the prediction of September Arctic sea ice extent along with increased interannual predictive skill. Significant local improvements in sea ice concentration were also seen in distinct regions for the change to PIOMAS initial thickness or the physics adjustments, with the most improvement occurring when these changes were applied concurrently.

Assessment of sea ice simulations in the CMIP5 models

Abstract. The historical simulations of sea ice during 1979 to 2005 by the Coupled Model Intercomparison Project Phase 5 (CMIP5) are compared with satellite observations, Global Ice-Ocean Modeling and Assimilation System (GIOMAS) output data and Pan-Arctic Ice Ocean Modeling and Assimilation System (PIOMAS) output data in this study. Forty-nine models, almost all of the CMIP5 climate models and earth system models with historical simulation, are used. For the Antarctic, multi-model ensemble mean (MME) results can give good climatology of sea ice extent (SIE), but the linear trend is incorrect. The linear trend of satellite-observed Antarctic SIE is 1.29 (±0.57) × 105 km2 decade−1; only about 1/7 CMIP5 models show increasing trends, and the linear trend of CMIP5 MME is negative with the value of −3.36 (±0.15) × 105 km2 decade−1. For the Arctic, both climatology and linear trend are better reproduced. Sea ice volume (SIV) is also evaluated in this study, and this is a first attempt to evaluate the SIV in all CMIP5 models. Compared with the GIOMAS and PIOMAS data, the SIV values in both the Antarctic and the Arctic are too small, especially for the Antarctic in spring and winter. The GIOMAS Antarctic SIV in September is 19.1 × 103 km3, while the corresponding Antarctic SIV of CMIP5 MME is 13.0 × 103 km3 (almost 32% less). The Arctic SIV of CMIP5 in April is 27.1 × 103 km3, which is also less than that from PIOMAS SIV (29.5 × 103 km3). This means that the sea ice thickness simulated in CMIP5 is too thin, although the SIE is fairly well simulated.

Uncertainty in modeled Arctic sea ice volume

[1] Uncertainty in the Pan-Arctic Ice-Ocean Modeling and Assimilation System (PIOMAS) Arctic sea ice volume record is characterized. A range of observations and approaches, including in situ ice thickness measurements, ICESat retrieved ice thickness, and model sensitivity studies, yields a conservative estimate for October Arctic ice volume uncertainty of 1.35 × 103 km3 and an uncertainty of the ice volume trend over the 1979–2010 period of 1.0 × 103 km3 decade–1. A conservative estimate of the trend over this period is −2.8 × 103 km3 decade–1. PIOMAS ice thickness estimates agree well with ICESat ice thickness retrievals (<0.1 m mean difference) for the area for which submarine data are available, while difference outside this area are larger. PIOMAS spatial thickness patterns agree well with ICESat thickness estimates with pattern correlations of above 0.8. PIOMAS appears to overestimate thin ice thickness and underestimate thick ice, yielding a smaller downward trend than apparent in reconstructions from observations. PIOMAS ice volume uncertainties and trends are examined in the context of climate change attribution and the declaration of record minima. The distribution of 32 year trends in a preindustrial coupled model simulation shows no trends comparable to those seen in the PIOMAS retrospective, even when the trend uncertainty is accounted for. Attempts to label September minima as new record lows are sensitive to modeling error. However, the September 2010 ice volume anomaly did in fact exceed the previous 2007 minimum by a large enough margin to establish a statistically significant new record.