Although iceberg models have been used for decades, they have received far more widespread attention in recent years, due in part to efforts to explicitly represent icebergs in climate models. This calls for increased scrutiny of all aspects of typical iceberg models. An important component of iceberg models is the representation of iceberg capsizing, or rolling. Rolling occurs spontaneously when the ratio of iceberg width to height falls below a critical threshold. Here we examine previously proposed representations of this threshold, and we find that there have been crucial flaws in the representation of rolling in many modeling studies to date. We correct these errors and identify an accurate model representation of iceberg rolling. Next, we assess how iceberg rolling influences simulation results in a hierarchy of models. Rolling is found to substantially prolong the lifespan of individual icebergs and allow them to drift farther offshore. However, rolling occurs only after large icebergs have lost most of their initial volume, and it thus has a relatively small impact on the large-scale freshwater distribution in comprehensive model simulations. The results suggest that accurate representations of iceberg rolling may be of particular importance for operational forecast models of iceberg drift, as well as for regional changes in high-resolution climate model simulations. (C) 2017 Elsevier Ltd. All rights reserved.

Observations indicate that the Arctic sea ice cover is rapidly retreating while the Antarctic sea ice cover is steadily expanding. State-of-the-art climate models, by contrast, typically simulate a moderate decrease in both the Arctic and Antarctic sea ice covers. However, in each hemisphere there is a small subset of model simulations that have sea ice trends similar to the observations. Based on this, a number of recent studies have suggested that the models are consistent with the observations in each hemisphere when simulated internal climate variability is taken into account. Here sea ice changes during 1979-2013 are examined in simulations from the most recent Coupled Model Intercomparison Project (CMIP5) as well as the Community Earth System Model Large Ensemble (CESM-LE), drawing on previous work that found a close relationship in climate models between global-mean surface temperature and sea ice extent. All of the simulations with 1979-2013 Arctic sea ice retreat as fast as observations are found to have considerably more global warming than observations during this time period. Using two separate methods to estimate the sea ice retreat that would occur under the observed level of global warming in each simulation in both ensembles, it is found that simulated Arctic sea ice retreat as fast as observations would occur less than 1% of the time. This implies that the models are not consistent with the observations. In the Antarctic, simulated sea ice expansion as fast as observations is found to typically correspond with too little global warming, although these results are more equivocal. As a result, the simulations do not capture the observed asymmetry between Arctic and Antarctic sea ice trends. This suggests that the models may be getting the right sea ice trends for the wrong reasons in both polar regions.

Sea ice features a dense inner pack ice zone surrounded by a marginal ice zone (MIZ) in which the sea ice properties are modified by interaction with the ice-free open ocean. The width of the MIZ is a fundamental length scale for polar physical and biological dynamics. Several different criteria for establishing MIZ boundaries have emerged in the literature-wave penetration, floe size, sea ice concentration, etc.-and a variety of definitions for the width between the MIZ boundaries have been published. Here, three desirable mathematical properties for defining MIZ width are proposed: invariance with respect to translation and rotation on the sphere; uniqueness at every point in the MIZ; and generality, including nonconvex shapes. The previously published streamline definition is shown to satisfy all three properties, where width is defined as the arc length of a streamline through the solution to Laplaces's equation within the MIZ boundaries, while other published definitions each satisfy only one of the desired properties. When defining MIZ spatial average width from streamline results, the rationale for averaging with respect to distance along both MIZ boundaries was left implicit in prior studies. Here it is made rigorous by developing and applying the mathematics of an analytically tractable idealization of MIZ geometry-the eccentric annulus. Finally, satellite-retrieved Arctic sea ice concentrations are used to investigate how well streamline-based MIZ spatial average width is approximated by alternative definitions that lack desirable mathematical properties or local width values but offer computational efficiency.

The discharge of icebergs into the polar oceans is expected to increase over the coming century, which raises the importance of accurate representations of icebergs in global climate models (GCMs) used for future projections. Here we analyze the prospects for interactive icebergs in GCMs by forcing an iceberg drift and decay model with circulation and temperature fields from (i) state-of-the-art GCM output and (ii) an observational state estimate. The spread of meltwater is found to be smaller for the GCM than for the observational state estimate, despite a substantial high wind bias in the GCM—a bias that is similar to most current GCMs. We argue that this large-scale reduction in the spread of meltwater occurs primarily due to localized differences in ocean currents, which may be related to the coarseness of the horizontal resolution in the GCM. The high wind bias in the GCM is shown to have relatively little impact on the meltwater distribution, despite Arctic iceberg drift typically being dominated by the wind forcing. We find that this is due to compensating effects between faster drift under stronger winds and larger wind-driven wave erosion. These results may have implications for future changes in the Atlantic meridional overturning circulation simulated with iceberg-enabled GCMs.

The downward trend in Arctic sea ice extent is one of the most dramatic signals of climate change during recent decades. Comprehensive climate models have struggled to reproduce this trend, typically simulating a slower rate of sea ice retreat than has been observed. However, this bias has been widely noted to have decreased in models participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5) compared with the previous generation of models (CMIP3). Here simulations are examined from both CMIP3 and CMIP5. It is found that simulated historical sea ice trends are influenced by volcanic forcing, which was included in all of the CMIP5 models but in only about half of the CMIP3 models. The volcanic forcing causes temporary simulated cooling in the 1980s and 1990s, which contributes to raising the simulated 1979-2013 global-mean surface temperature trends to values substantially larger than observed. It is shown that this warming bias is accompanied by an enhanced rate of Arctic sea ice retreat and hence a simulated sea ice trend that is closer to the observed value, which is consistent with previous findings of an approximately linear relationship between sea ice extent and global-mean surface temperature. Both generations of climate models are found to simulate Arctic sea ice that is substantially less sensitive to global warming than has been observed. The results imply that much of the difference in Arctic sea ice trends between CMIP3 and CMIP5 occurred because of the inclusion of volcanic forcing, rather than improved sea ice physics or model resolution.

Previous studies have suggested that the global ocean density stratification below approximate to 3000m is approximately set by its direct connection to the Southern Ocean surface density, which in turn is constrained by the atmosphere. Here the role of Southern Ocean surface forcing in glacial-interglacial stratification changes is investigated using a comprehensive climate model and an idealized conceptual model. Southern Ocean surface forcing is found to control the global deep ocean stratification up to approximate to 2000m, which is much shallower than previously thought and contrary to the expectation that the North Atlantic surface forcing should strongly influence the ocean at intermediate depths. We show that this is due to the approximately fixed surface freshwater fluxes, rather than a fixed surface density distribution in the Southern Ocean as was previously assumed. These results suggest that Southern Ocean surface freshwater forcing controls glacial-interglacial stratification changes in much of the deep ocean.

Understanding the causes of recent climatic trends and variability in the high-latitude Southern Hemisphere is hampered by a short instrumental record. Here, we analyse recent atmosphere, surface ocean and sea-ice observations in this region and assess their trends in the context of palaeoclimate records and climate model simulations. Over the 36-year satellite era, significant linear trends in annual mean sea-ice extent, surface temperature and sea-level pressure are superimposed on large interannual to decadal variability. Most observed trends, however, are not unusual when compared with Antarctic palaeoclimate records of the past two centuries. With the exception of the positive trend in the Southern Annular Mode, climate model simulations that include anthropogenic forcing are not compatible with the observed trends. This suggests that natural variability overwhelms the forced response in the observations, but the models may not fully represent this natural variability or may overestimate the magnitude of the forced response.

Uncovering universal early warning signals for critical transitions has become a coveted goal in diverse scientific disciplines, ranging from climate science to financial mathematics. There has been a flurry of recent research proposing such signals, with increasing autocorrelation and increasing variance being among the most widely discussed candidates. A number of studies have suggested that increasing autocorrelation alone may suffice to signal an impending transition, although some others have questioned this. Here we consider variance and autocorrelation in the context of sea ice loss in an idealized model of the global climate system. The model features no bifurcation, nor increased rate of retreat, as the ice disappears. Nonetheless, the autocorrelation of summer sea ice area is found to increase in a global warming scenario. The variance, by contrast, decreases. A simple physical mechanism is proposed to explain the occurrence of increasing autocorrelation but not variance when there is no approaching bifurcation. Additionally, a similar mechanism is shown to allow an increase in both indicators with no physically attainable bifurcation. This implies that relying on autocorrelation and variance as early warning signals can raise false alarms in the climate system, warning of "tipping points" that are not actually there.

Record lows in Arctic sea ice extent have been making frequent headlines in recent years. The change in albedo when sea ice is replaced by open water introduces a nonlinearity that has sparked an ongoing debate about the stability of the Arctic sea ice cover and the possibility of Arctic "tipping points.'' Previous studies identified instabilities for a shrinking ice cover in two types of idealized climate models: (i) annual-mean latitudinally varying diffusive energy balance models (EBMs) and (ii) seasonally varying single-column models (SCMs). The instabilities in these low-order models stand in contrast with results from comprehensive global climate models (GCMs), which typically do not simulate any such instability. To help bridge the gap between low-order models and GCMs, an idealized model is developed that includes both latitudinal and seasonal variations. The model reduces to a standard EBM or SCM as limiting cases in the parameter space, thus reconciling the two previous lines of research. It is found that the stability of the ice cover vastly increases with the inclusion of spatial communication via meridional heat transport or a seasonal cycle in solar forcing, being most stable when both are included. If the associated parameters are set to values that correspond to the current climate, the ice retreat is reversible and there is no instability when the climate is warmed. The two parameters have to be reduced by at least a factor of 3 for instability to occur. This implies that the sea ice cover may be substantially more stable than has been suggested in previous idealized modeling studies.

A fine-resolution (1/10°) ocean/sea ice model configured in the Community Earth System Model framework is compared with observations and studied to determine the basin-scale and local balances controlling the variability of sea ice anomalies from the mean seasonal cycle in the Bering Sea for the time period 1980–1989. The model produces variations in total Bering Sea ice area anomalies that are highly correlated with observations. Surface air temperature, which is specified from reanalysis atmospheric forcing, strongly controls the ice volume variability in this simulation. The thermodynamic ice volume change is dominated by surface energy flux via atmosphere-ice sensible heat flux, except near the southern ice edge where it is largely controlled by ocean-ice heat flux. While thermodynamic processes dominate the variations in ice volume in the Bering Sea on the large scale, dynamic processes are important on the local scale near ice margins (both oceanic and land), where dynamic and thermodynamic ice volume changes have opposite signs and nearly cancel each other. Ice motion is generally consistent with winds driving the flow, except near certain straits in the north where ice motion largely follows ocean currents. Two key climate events, strong ice growth with cold air temperature and northerly wind anomalies in February 1984 and weak ice growth with warm air temperature and southerly wind anomalies in February 1989, are studied here in detail. While the processes controlling the ice changes are generally similar to those in other years, these large events help reveal the characteristic spatial patterns of ice growth/melt and transport anomalies.

The seasonal cycle of sea ice variability in the Bering Sea, together with the thermodynamic and dynamic processes that control it, are examined in a fine resolution (1/10°) global coupled ocean/sea-ice model configured in the Community Earth System Model (CESM) framework. The ocean/sea-ice model consists of the Los Alamos National Laboratory Parallel Ocean Program (POP) and the Los Alamos Sea Ice Model (CICE). The model was forced with time-varying reanalysis atmospheric forcing for the time period 1970–1989. This study focuses on the time period 1980–1989. The simulated seasonal-mean fields of sea ice concentration strongly resemble satellite-derived observations, as quantified by root-mean-square errors and pattern correlation coefficients. The sea ice energy budget reveals that the seasonal thermodynamic ice volume changes are dominated by the surface energy flux between the atmosphere and the ice in the northern region and by heat flux from the ocean to the ice along the southern ice edge, especially on the western side. The sea ice force balance analysis shows that sea ice motion is largely associated with wind stress. The force due to divergence of the internal ice stress tensor is large near the land boundaries in the north, and it is small in the central and southern ice-covered region. During winter, which dominates the annual mean, it is found that the simulated sea ice was mainly formed in the northern Bering Sea, with the maximum ice growth rate occurring along the coast due to cold air from northerly winds and ice motion away from the coast. South of St Lawrence Island, winds drive the model sea ice southwestward from the north to the southwestern part of the ice-covered region. Along the ice edge in the western Bering Sea, model sea ice is melted by warm ocean water, which is carried by the simulated Bering Slope Current flowing to the northwest, resulting in the S-shaped asymmetric ice edge. In spring and fall, similar thermodynamic and dynamic patterns occur in the model, but with typically smaller magnitudes and with season-specific geographical and directional differences.

Intense glaciation during the end of Cryogenian time (similar to 635 million years ago) marks the coldest climate state in Earth history - a time when glacial deposits accumulated at low, tropical paleolatitudes. The leading idea to explain these deposits, the snowball Earth hypothesis, predicts globally frozen surface conditions and subfreezing temperatures, with global climate models placing surface temperatures in the tropics between -20 degrees C and -60 degrees C. However, precise paleosurface temperatures based upon geologic constraints have remained elusive and the global severity of the glaciation undetermined. Here we make new geologic observations of tropical periglacial, aeolian and fluvial sedimentary structures formed during the end-Cryogenian, Marinoan glaciation in South Australia; these observations allow us to constrain ancient surface temperatures. We find periglacial sand wedges and associated deformation suggest that ground temperatures were sufficiently warm to allow for ductile deformation of a sandy regolith. The wide range of deformation structures likely indicate the presence of a paleoactive layer that penetrated 2-4 m below the ground surface. These observations, paired with a model of ground temperature forced by solar insolation, constrain the local mean annual surface temperature to within a few degrees of freezing. This temperature constraint matches well with our observations of fluvial deposits, which require temperatures sufficiently warm for surface runoff. Although this estimate coincides with one of the coldest near sea-level tropical temperatures in Earth history, if these structures represent peak Marinaon glacial conditions, they do not support the persistent deep freeze of the snowball Earth hypothesis. Rather, surface temperatures near 0 degrees C allow for regions of seasonal surface melting, atmosphere-ocean coupling and possible tropical refugia for early metazoans. If instead these structures formed during glacial onset or deglaciation, then they have implications for the timescale and character for the transition into or out of a snowball state. (C) 2014 Elsevier B.V. All rights reserved.

The transient response of the Atlantic Meridional Overturning Circulation (AMOC) to a deglacial ice sheet retreat is studied using the Community Climate System Model version 3 (CCSM3), with a focus on orographic effects rather than meltwater discharge. It is found that the AMOC weakens significantly (41%) in response to the deglacial ice sheet retreat. The AMOC weakening follows the decrease of the Northern Hemisphere ice sheet volume linearly, with no evidence of abrupt thresholds. A wind-driven mechanism is proposed to explain the weakening of the AMOC: lowering the Northern Hemisphere ice sheets induces a northward shift of the westerlies, which causes a rapid eastward sea ice transport and expanded sea ice cover over the subpolar North Atlantic; this expanded sea ice insulates the ocean from heat loss and leads to suppressed deep convection and a weakened AMOC. A sea ice-ocean positive feedback could be further established between the AMOC decrease and sea ice expansion.

Recent estimates indicate that the Antarctic sea ice cover is expanding at a statistically significant rate with a magnitude one-third as large as the rapid rate of sea ice retreat in the Arctic. However, during the mid-2000s, with several fewer years in the observational record, the trend in Antarctic sea ice extent was reported to be considerably smaller and statistically indistinguishable from zero. Here, we show that much of the increase in the reported trend occurred due to the previously undocumented effect of a change in the way the satellite sea ice observations are processed for the widely used Bootstrap algorithm data set, rather than a physical increase in the rate of ice advance. Specifically, we find that a change in the intercalibration across a 1991 sensor transition when the data set was reprocessed in 2007 caused a substantial change in the long-term trend. Although our analysis does not definitively identify whether this change introduced an error or removed one, the resulting difference in the trends suggests that a substantial error exists in either the current data set or the version that was used prior to the mid-2000s, and numerous studies that have relied on these observations should be reexamined to determine the sensitivity of their results to this change in the data set. Furthermore, a number of recent studies have investigated physical mechanisms for the observed expansion of the Antarctic sea ice cover. The results of this analysis raise the possibility that much of this expansion may be a spurious artifact of an error in the processing of the satellite observations.

The decline of Arctic sea ice has been documented in over 30 y of satellite passive microwave observations. The resulting darkening of the Arctic and its amplification of global warming was hypothesized almost 50 y ago but has yet to be verified with direct observations. This study uses satellite radiation budget measurements along with satellite microwave sea ice data to document the Arctic-wide decrease in planetary albedo and its amplifying effect on the warming. The analysis reveals a striking relationship between planetary albedo and sea ice cover, quantities inferred from two independent satellite instruments. We find that the Arctic planetary albedo has decreased from 0.52 to 0.48 between 1979 and 2011, corresponding to an additional 6.4 +/- 0.9 W/m(2) of solar energy input into the Arctic Ocean region since 1979. Averaged over the globe, this albedo decrease corresponds to a forcing that is 25% as large as that due to the change in CO2 during this period, considerably larger than expectations from models and other less direct recent estimates. Changes in cloudiness appear to play a negligible role in observed Arctic darkening, thus reducing the possibility of Arctic cloud albedo feedbacks mitigating future Arctic warming.

Orbital precession changes the seasonal distribution of insolation at a given latitude but not the annual mean. Hence, the correlation of paleoclimate proxies of annual-mean precipitation with orbital precession implies a nonlinear rectification in the precipitation response to seasonal solar forcing. It has previously been suggested that the relevant nonlinearity is that of the Clausius-Clapeyron relationship. Here it is argued that a different nonlinearity related to moisture advection by the atmospheric circulation is more important. When perihelion changes from one hemisphere's summer solstice to the other's in an idealized aquaplanet atmospheric general circulation model, annual-mean precipitation increases in the hemisphere with the brighter, warmer summer and decreases in the other hemisphere, in qualitative agreement with paleoclimate proxies that indicate such hemispherically antisymmetric climate variations. The rectification mechanism that gives rise to the precipitation changes is identified by decomposing the perturbation water vapor budget into "thermodynamic" and "dynamic" components. Thermodynamic changes (caused by changes in humidity with unchanged winds) dominate the hemispherically antisymmetric annual-mean precipitation response to precession in the absence of land-sea contrasts. The nonlinearity that enables the thermodynamic changes to affect annual-mean precipitation is a nonlinearity of moisture advection that arises because precession-induced seasonal humidity changes correlate with the seasonal cycle in low-level convergence. This interpretation is confirmed using simulations in which the Clausius-Clapeyron relationship is explicitly linearized. The thermodynamic mechanism also operates in simulations with an idealized representation of land, although in these simulations the dynamic component of the precipitation changes is also important, adding to the thermodynamic precipitation changes in some latitudes and offsetting it in others.

The response of the monsoonal and annual-mean Hadley circulation to orbital precession is examined in an idealized atmospheric general circulation model with an aquaplanet slab-ocean lower boundary. Contrary to expectations, the simulated monsoonal Hadley circulation is weaker when perihelion occurs at the summer solstice than when aphelion occurs at the summer solstice. The angular momentum balance and energy balance are examined to understand the mechanisms that produce this result. That the summer with stronger insolation has a weaker circulation is the result of an increase in the atmosphere's energetic stratification, the gross moist stability, which increases more than the amount required to balance the change in atmospheric energy flux divergence necessitated by the change in top-of-atmosphere net radiation. The solstice-season changes result in annual-mean Hadley circulation changes (e.g., changes in circulation strength).

The response of the monsoonal and annual-mean Hadley circulation to orbital precession is examined in an idealized atmospheric general circulation model with a simplified representation of land surface processes in subtropical latitudes. When perihelion occurs in the summer of a hemisphere with a subtropical continent, changes in the top-of-atmosphere energy balance, together with a poleward shift of the monsoonal circulation boundary, lead to a strengthening of the monsoonal circulation. Spatial variations in surface heat capacity determine whether radiative perturbations are balanced by energy storage or by atmospheric energy fluxes. Although orbital precession does not affect annual-mean insolation, the annual-mean Hadley circulation does respond to orbital precession because its sensitivity to radiative changes varies over the course of the year: the monsoonal circulation in summer is near the angular momentum-conserving limit and responds directly to radiative changes; whereas in winter, the circulation is affected by the momentum fluxes of extratropical eddies and is less sensitive to radiative changes.