Refine
Year of publication
Language
- English (75)
Is part of the Bibliography
- yes (75)
Keywords
- climate change (7)
- climate impacts (3)
- model (3)
- sea-level rise (3)
- tropical cyclones (3)
- chains (2)
- costs (2)
- damage (2)
- impact (2)
- meteorological extremes (2)
Institute
- Institut für Physik und Astronomie (75) (remove)
We present a minimal conceptual model for the Atlantic meridional overturning circulation which incorporates the advection of salinity and the basic dynamics of the oceanic pycnocline. Four tracer transport processes following Gnanadesikan in Science 283(5410):2077-2079, (1999) allow for a dynamical adjustment of the oceanic pycnocline which defines the vertical extent of a mid-latitudinal box. At the same time the model captures the salt-advection feedback (Stommel in Tellus 13(2):224-230, (1961)). Due to its simplicity the model can be solved analytically in the purely wind- and purely mixing-driven cases. We find the possibility of abrupt transition in response to surface freshwater forcing in both cases even though the circulations are very different in physics and geometry. This analytical approach also provides expressions for the critical freshwater input marking the change in the dynamics of the system. Our analysis shows that including the pycnocline dynamics in a salt-advection model causes a decrease in the freshwater sensitivity of its northern sinking up to a threshold at which the circulation breaks down. Compared to previous studies the model is restricted to the essential ingredients. Still, it exhibits a rich behavior which reaches beyond the scope of this study and might be used as a paradigm for the qualitative behaviour of the Atlantic overturning in the discussion of driving mechanisms.
In recent decades, the Greenland Ice Sheet has been losing mass and has thereby contributed to global sea-level rise. The rate of ice loss is highly relevant for coastal protection worldwide. The ice loss is likely to increase under future warming. Beyond a critical temperature threshold, a meltdown of the Greenland Ice Sheet is induced by the self-enforcing feedback between its lowering surface elevation and its increasing surface mass loss: the more ice that is lost, the lower the ice surface and the warmer the surface air temperature, which fosters further melting and ice loss. The computation of this rate so far relies on complex numerical models which are the appropriate tools for capturing the complexity of the problem. By contrast we aim here at gaining a conceptual understanding by deriving a purposefully simple equation for the self-enforcing feedback which is then used to estimate the melt time for different levels of warming using three observable characteristics of the ice sheet itself and its surroundings. The analysis is purely conceptual in nature. It is missing important processes like ice dynamics for it to be useful for applications to sea-level rise on centennial timescales, but if the volume loss is dominated by the feedback, the resulting logarithmic equation unifies existing numerical simulations and shows that the melt time depends strongly on the level of warming with a critical slow-down near the threshold: the median time to lose 10% of the present-day ice volume varies between about 3500 years for a temperature level of 0.5 degrees C above the threshold and 500 years for 5 degrees C. Unless future observations show a significantly higher melting sensitivity than currently observed, a complete meltdown is unlikely within the next 2000 years without significant ice-dynamical contributions.
Both ice sheets in Greenland and Antarctica are discharging ice into the ocean. In many regions along the coast of the ice sheets, the icebergs calve into a bay. If the addition of icebergs through calving is faster than their transport out of the embayment, the icebergs will be frozen into a melange with surrounding sea ice in winter. In this case, the buttressing effect of the ice melange can be considerably stronger than any buttressing by mere sea ice would be. This in turn stabilizes the glacier terminus and leads to a reduction in calving rates. Here we propose a simple parametrization of ice melange buttressing which leads to an upper bound on calving rates and can be used in numerical and analytical modelling.
Over large coastal regions in Greenland and Antarctica the ice sheet calves directly into the ocean. In contrast to ice-shelf calving, an increase in calving from grounded glaciers contributes directly to sea-level rise. Ice cliffs with a glacier freeboard larger than approximate to 100 m are currently not observed, but it has been shown that such ice cliffs are increasingly unstable with increasing ice thickness. This cliff calving can constitute a self-amplifying ice loss mechanism that may significantly alter sea-level projections both of Greenland and Antarctica. Here we seek to derive a minimalist stress-based parametrization for cliff calving from grounded glaciers whose freeboards exceed the 100m stability limit derived in previous studies. This will be an extension of existing calving laws for tidewater glaciers to higher ice cliffs. To this end we compute the stress field for a glacier with a simplified two-dimensional geometry from the two-dimensional Stokes equation. First we assume a constant yield stress to derive the failure region at the glacier front from the stress field within the glacier. Secondly, we assume a constant response time of ice failure due to exceedance of the yield stress. With this strongly constraining but very simple set of assumptions we propose a cliff-calving law where the calving rate follows a power-law dependence on the freeboard of the ice with exponents between 2 and 3, depending on the relative water depth at the calving front. The critical freeboard below which the ice front is stable decreases with increasing relative water depth of the calving front. For a dry water front it is, for example, 75 m. The purpose of this study is not to provide a comprehensive calving law but to derive a particularly simple equation with a transparent and minimalist set of assumptions.
Indian monsoon rainfall is vital for a large share of the world's population. Both reliably projecting India's future precipitation and unraveling abrupt cessations of monsoon rainfall found in paleorecords require improved understanding of its stability properties. While details of monsoon circulations and the associated rainfall are complex, full-season failure is dominated by large-scale positive feedbacks within the region. Here we find that in a comprehensive climate model, monsoon failure is possible but very rare under pre-industrial conditions, while under future warming it becomes much more frequent. We identify the fundamental intraseasonal feedbacks that are responsible for monsoon failure in the climate model, relate these to observational data, and build a statistically predictive model for such failure. This model provides a simple dynamical explanation for future changes in the frequency distribution of seasonal mean all-Indian rainfall. Forced only by global mean temperature and the strength of the Pacific Walker circulation in spring, it reproduces the trend as well as the multidecadal variability in the mean and skewness of the distribution, as found in the climate model. The approach offers an alternative perspective on large-scale monsoon variability as the result of internal instabilities modulated by pre-seasonal ambient climate conditions.
Abrupt monsoon transitions as seen in paleorecords can be explained by moisture-advection feedback
(2016)
Earth’s surface temperature will continue to rise for another 20 to 30 years even with the strongest carbon emission reduction currently considered. The associated changes in rainfall patterns can result in an increased flood risk worldwide. We compute the required increase in flood protection to keep high-end fluvial flood risk at present levels. The analysis is carried out worldwide for subnational administrative units. Most of the United States, Central Europe, and Northeast and West Africa, as well as large parts of India and Indonesia, require the strongest adaptation effort. More than half of the United States needs to at least double their protection within the next two decades. Thus, the need for adaptation to increased river flood is a global problem affecting industrialized regions as much as developing countries.
Anthropogenic carbon emissions lock in long-term sea-level rise that greatly exceeds projections for this century, posing profound challenges for coastal development and cultural legacies. Analysis based on previously published relationships linking emissions to warming and warming to rise indicates that unabated carbon emissions up to the year 2100 would commit an eventual global sea-level rise of 4.3-9.9 m. Based on detailed topographic and population data, local high tide lines, and regional long-term sea-level commitment for different carbon emissions and ice sheet stability scenarios, we compute the current population living on endangered land at municipal, state, and national levels within the United States. For unabated climate change, we find that land that is home to more than 20 million people is implicated and is widely distributed among different states and coasts. The total area includes 1,185-1,825 municipalities where land that is home to more than half of the current population would be affected, among them at least 21 cities exceeding 100,000 residents. Under aggressive carbon cuts, more than half of these municipalities would avoid this commitment if the West Antarctic Ice Sheet remains stable. Similarly, more than half of the US population-weighted area under threat could be spared. We provide lists of implicated cities and state populations for different emissions scenarios and with and without a certain collapse of the West Antarctic Ice Sheet. Although past anthropogenic emissions already have caused sea-level commitment that will force coastal cities to adapt, future emissions will determine which areas we can continue to occupy or may have to abandon.
We present climatic consequences of the Representative Concentration Pathways (RCPs) using the coupled climate model CLIMBER-3 alpha, which contains a statistical-dynamical atmosphere and a three-dimensional ocean model. We compare those with emulations of 19 state-of-the-art atmosphere-ocean general circulation models (AOGCM) using MAGICC6. The RCPs are designed as standard scenarios for the forthcoming IPCC Fifth Assessment Report to span the full range of future greenhouse gas (GHG) concentrations pathways currently discussed. The lowest of the RCP scenarios, RCP3-PD, is projected in CLIMBER-3 alpha to imply a maximal warming by the middle of the 21st century slightly above 1.5 degrees C and a slow decline of temperatures thereafter, approaching today's level by 2500. We identify two mechanisms that slow down global cooling after GHG concentrations peak: The known inertia induced by mixing-related oceanic heat uptake; and a change in oceanic convection that enhances ocean heat loss in high latitudes, reducing the surface cooling rate by almost 50%. Steric sea level rise under the RCP3-PD scenario continues for 200 years after the peak in surface air temperatures, stabilizing around 2250 at 30 cm. This contrasts with around 1.3 m of steric sea level rise by 2250, and 2 m by 2500, under the highest scenario, RCP8.5. Maximum oceanic warming at intermediate depth (300-800 m) is found to exceed that of the sea surface by the second half of the 21st century under RCP3-PD. This intermediate-depth warming persists for centuries even after surface temperatures have returned to present-day values, with potential consequences for marine ecosystems, oceanic methane hydrates, and ice-shelf stability. Due to an enhanced land-ocean temperature contrast, all scenarios yield an intensification of monsoon rainfall under global warming.
Large-scale transitions in societies are associated with both individual behavioural change and restructuring of the social network. These two factors have often been considered independently, yet recent advances in social network research challenge this view. Here we show that common features of societal marginalization and clustering emerge naturally during transitions in a co-evolutionary adaptive network model. This is achieved by explicitly considering the interplay between individual interaction and a dynamic network structure in behavioural selection. We exemplify this mechanism by simulating how smoking behaviour and the network structure get reconfigured by changing social norms. Our results are consistent with empirical findings: The prevalence of smoking was reduced, remaining smokers were preferentially connected among each other and formed increasingly marginalized clusters. We propose that self-amplifying feedbacks between individual behaviour and dynamic restructuring of the network are main drivers of the transition. This generative mechanism for co-evolution of individual behaviour and social network structure may apply to a wide range of examples beyond smoking.
The future evolution of the Antarctic Ice Sheet represents the largest uncertainty in sea-level projections of this and upcoming centuries. Recently, satellite observations and high-resolution simulations have suggested the initiation of an ice-sheet instability in the Amundsen Sea sector of West Antarctica, caused by the last decades' enhanced basal ice-shelf melting. Whether this localized destabilization will yield a full discharge of marine ice from West Antarctica, associated with a global sea-level rise of more than 3 m, or whether the ice loss is limited by ice dynamics and topographic features, is unclear. Here we show that in the Parallel Ice Sheet Model, a local destabilization causes a complete disintegration of the marine ice in West Antarctica. In our simulations, at 5-km horizontal resolution, the region disequilibrates after 60 y of currently observed melt rates. Thereafter, the marine ice-sheet instability fully unfolds and is not halted by topographic features. In fact, the ice loss in Amundsen Sea sector shifts the catchment's ice divide toward the Filchner-Ronne and Ross ice shelves, which initiates grounding-line retreat there. Our simulations suggest that if a destabilization of Amundsen Sea sector has indeed been initiated, Antarctica will irrevocably contribute at least 3 m to global sea-level rise during the coming centuries to millennia.
Consequences of twenty-first-century policy for multi-millennial climate
and sea-level change
(2016)
Most of the policy debate surrounding the actions needed to mitigate and adapt to anthropogenic climate change has been framed by observations of the past 150 years as well as climate and sea-level projections for the twenty-first century. The focus on this 250-year window, however, obscures some of the most profound problems associated with climate change. Here, we argue that the twentieth and twenty-first centuries, a period during which the overwhelming majority of human-caused carbon emissions are likely to occur, need to be placed into a long-term context that includes the past 20 millennia, when the last Ice Age ended and human civilization developed, and the next ten millennia, over which time the projected impacts of anthropogenic climate change will grow and persist. This long-term perspective illustrates that policy decisions made in the next few years to decades will have profound impacts on global climate, ecosystems and human societies - not just for this century, but for the next ten millennia and beyond.
Projections of changes in Antarctic Ice Sheet (AIS) surface mass balance indicate a negative contribution to sea level because of the expected increase in precipitation due to the higher moisture holding capacity of warmer air(1). Observations over the past decades, however, are unable to constrain the relation between temperature and accumulation changes because both are dominated by strong natural variability(2-5). Here we derive a consistent continental-scale increase in accumulation of approximately 5 +/- 1% K-1, through the assessment of ice-core data (spanning the large temperature change during the last deglaciation, 21,000 to 10,000 years ago), in combination with palaeo-simulations, future projections by 35 general circulation models (GCMs), and one high-resolution future simulation. The ice-core data and modelling results for the last deglaciation agree, showing uniform local sensitivities of similar to 6% K-1. The palaeo-simulation allows for a continental-scale aggregation of accumulation changes reaching 4.3% K-1. Despite the different timescales, these sensitivities agree with the multi-model mean of 6.1 +/- 2.6% K-1 (GCMprojections) and the continental-scale sensitivity of 4.9% K-1 (high-resolution future simulation). Because some of the mass gain of the AIS is offset by dynamical losses induced by accumulation(6,7), we provide a response function allowing projections of sea-level fall in terms of continental-scale accumulation changes that compete with surface melting and dynamical losses induced by other mechanisms(6,8,9).
The possibility of an impact of global warming on the Indian monsoon is of critical importance for the large population of this region. Future projections within the Coupled Model Intercomparison Project Phase 3 (CMIP-3) showed a wide range of trends with varying magnitude and sign across models. Here the Indian summer monsoon rainfall is evaluated in 20 CMIP-5 models for the period 1850 to 2100. In the new generation of climate models, a consistent increase in seasonal mean rainfall during the summer monsoon periods arises. All models simulate stronger seasonal mean rainfall in the future compared to the historic period under the strongest warming scenario RCP-8.5. Increase in seasonal mean rainfall is the largest for the RCP-8.5 scenario compared to other RCPs. Most of the models show a northward shift in monsoon circulation by the end of the 21st century compared to the historic period under the RCP-8.5 scenario. The interannual variability of the Indian monsoon rainfall also shows a consistent positive trend under unabated global warming. Since both the long-term increase in monsoon rainfall as well as the increase in interannual variability in the future is robust across a wide range of models, some confidence can be attributed to these projected trends.
Even if greenhouse gas emissions were stopped today, sea level would continue to rise for centuries, with the long-term sea-level commitment of a 2 degrees C warmer world significantly exceeding 2 m. In view of the potential implications for coastal populations and ecosystems worldwide, we investigate, from an ice-dynamic perspective, the possibility of delaying sea-level rise by pumping ocean water onto the surface of the Antarctic ice sheet. We find that due to wave propagation ice is discharged much faster back into the ocean than would be expected from a pure advection with surface velocities. The delay time depends strongly on the distance from the coastline at which the additional mass is placed and less strongly on the rate of sea-level rise that is mitigated. A millennium-scale storage of at least 80% of the additional ice requires placing it at a distance of at least 700 km from the coastline. The pumping energy required to elevate the potential energy of ocean water to mitigate the currently observed 3 mmyr(-1) will exceed 7% of the current global primary energy supply. At the same time, the approach offers a comprehensive protection for entire coastlines particularly including regions that cannot be protected by dikes.
Tropical cyclones range among the costliest of all meteorological events worldwide and planetary scale warming provides more energy and moisture to these storms. Modelling the national and global economic repercussions of 2017's Hurricane Harvey, we find a qualitative change in the global economic response in an increasingly warmer world.
While the United States were able to balance regional production failures by the original 2017 hurricane, this option becomes less viable under future warming.
In our simulations of over 7000 regional economic sectors with more than 1.8 million supply chain connections, the US are not able to offset the losses by use of national efforts with intensifying hurricanes under unabated warming.
At a certain warming level other countries have to step in to supply the necessary goods for production, which gives US economic sectors a competitive disadvantage. In the highly localized mining and quarrying sector-which here also comprises the oil and gas production industry-this disadvantage emerges already with the original Hurricane Harvey and intensifies under warming.
Eventually, also other regions reach their limit of what they can offset.
While we chose the example of a specific hurricane impacting a specific region, the mechanism is likely applicable to other climate-related events in other regions and other sectors.
It is thus likely that the regional economic sectors that are best adapted to climate change gain significant advantage over their competitors under future warming.
Direct observations, satellite measurements and paleo records reveal strong variability in the Atlantic subpolar gyre on various time scales. Here we show that variations of comparable amplitude can only be simulated in a coupled climate model in the proximity of a dynamical threshold. The threshold and the associated dynamic response is due to a positive feedback involving increased salt transport in the subpolar gyre and enhanced deep convection in its centre. A series of sensitivity experiments is performed with a coarse resolution ocean general circulation model coupled to a statistical-dynamical atmosphere model which in itself does not produce atmospheric variability. To simulate the impact of atmospheric variability, the model system is perturbed with freshwater forcing of varying, but small amplitude and multi-decadal to centennial periodicities and observational variations in wind stress. While both freshwater and wind-stress-forcing have a small direct effect on the strength of the subpolar gyre, the magnitude of the gyre's response is strongly increased in the vicinity of the threshold. Our results indicate that baroclinic self-amplification in the North Atlantic ocean can play an important role in presently observed SPG variability and thereby North Atlantic climate variability on multi-decadal scales.
Assessing global impacts of unexpected meteorological events in an increasingly connected world economy is important for estimating the costs of climate change. We show that since the beginning of the 21st century, the structural evolution of the global supply network has been such as to foster an increase of climate-related production losses. We compute first- and higher-order losses from heat stress-induced reductions in productivity under changing economic and climatic conditions between 1991 and 2011. Since 2001, the economic connectivity has augmented in such a way as to facilitate the cascading of production loss. The influence of this structural change has dominated over the effect of the comparably weak climate warming during this decade. Thus, particularly under future warming, the intensification of international trade has the potential to amplify climate losses if no adaptation measures are taken.
The Indian summer monsoon shapes the livelihood of a large share of the world's population. About 80% of annual precipitation over India occurs during the monsoon season from June through September. Next to its seasonal mean rainfall, the day-to-day variability is crucial for the risk of flooding, national water supply, and agricultural productivity. Here we show that the latest ensemble of climate model simulations, prepared for the AR-5 of the Intergovernmental Panel on Climate Change, consistently projects significant increases in day-to-day rainfall variability under unmitigated climate change. The relative increase by the period 2071-2100 with respect to the control period 1871-1900 ranges from 13% to 50% under the strongest scenario (Representative Concentration Pathways, RCP-8.5), in the 10 models with the most realistic monsoon climatology; and 13% to 85% when all the 20 models are considered. The spread across models reduces when variability increase per degree of global warming is considered, which is independent of the scenario in most models, and is 8% +/- 4%/K on average. This consistent projection across 20 comprehensive climate models provides confidence in the results and suggests the necessity of profound adaptation measures in the case of unmitigated climate change.
Thawing of permafrost and the associated release of carbon constitutes a positive feedback in the climate system, elevating the effect of anthropogenic GHG emissions on global-mean temperatures. Multiple factors have hindered the quantification of this feedback, which was not included in climate carbon-cycle models which participated in recent model intercomparisons (such as the Coupled Carbon Cycle Climate Model Intercomparison Project - (CMIP)-M-4). There are considerable uncertainties in the rate and extent of permafrost thaw, the hydrological and vegetation response to permafrost thaw, the decomposition timescales of freshly thawed organic material, the proportion of soil carbon that might be emitted as carbon dioxide via aerobic decomposition or as methane via anaerobic decomposition, and in the magnitude of the high latitude amplification of global warming that will drive permafrost degradation. Additionally, there are extensive and poorly characterized regional heterogeneities in soil properties, carbon content, and hydrology. Here, we couple a new permafrost module to a reduced complexity carbon-cycle climate model, which allows us to perform a large ensemble of simulations. The ensemble is designed to span the uncertainties listed above and thereby the results provide an estimate of the potential strength of the feedback from newly thawed permafrost carbon. For the high CO2 concentration scenario (RCP8.5), 33-114 GtC (giga tons of Carbon) are released by 2100 (68% uncertainty range). This leads to an additional warming of 0.04-0.23 degrees C. Though projected 21st century permafrost carbon emissions are relatively modest, ongoing permafrost thaw and slow but steady soil carbon decomposition means that, by 2300, about half of the potentially vulnerable permafrost carbon stock in the upper 3 m of soil layer (600-1000 GtC) could be released as CO2, with an extra 1-4% being released as methane. Our results also suggest that mitigation action in line with the lower scenario RCP3-PD could contain Arctic temperature increase sufficiently that thawing of the permafrost area is limited to 9-23% and the permafrost-carbon induced temperature increase does not exceed 0.04-0.16 degrees C by 2300.
Recent observations and modeling studies emphasize the crucial role of fracture mechanics for the stability of ice shelves and thereby the evolution of ice sheets. Here we introduce a macroscopic fracture-density field into a prognostic continuum ice-flow model and compute its evolution incorporating the initiation and growth of fractures as well as their advection with two-dimensional ice flow. To a first approximation, fracture growth is assumed to depend on the spreading rate only, while fracture initiation is defined in terms of principal stresses. The inferred fracture-density fields compare well with observed elongate surface structures. Since crevasses and other deep-reaching fracture structures have been shown to influence the overall ice-shelf dynamics, we propose the fracture-density field introduced here be used as a measure for ice softening and decoupling of the ice flow in fracture-weakened zones. This may yield more accurate and realistic velocity patterns in prognostic simulations. Additionally, the memory of past fracture events links the calving front to the upstream dynamics. Thus the fracture-density field proposed here may be employed in fracture-based calving parameterizations. The aim of this study is to introduce the field and investigate which of the observed surface structures can be reproduced by the simplest physically motivated fracture source terms.
Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by the dynamic effects of fracture processes within the protective ice shelves, leading to accelerated ice flow and hence to a sea-level contribution. In order to account for the macroscopic effect of fracture processes on large-scale viscous ice dynamics (i.e., ice-shelf scale) we apply a continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM) and compare the results to observations. To this end we introduce a higher order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in the fracture-initiation threshold. As a result of prognostic flow simulations, sharp across-flow velocity gradients appear in fracture-weakened regions. These modeled gradients compare well in magnitude and location with those in observed flow patterns. This model framework is in principle expandable to grounded ice streams and provides simple means of investigating climate-induced effects on fracturing (e. g., hydro fracturing) and hence on the ice flow. It further constitutes a physically sound basis for an enhanced fracture-based calving parameterization.
Here we report on a cyclic, physical ice-discharge instability in the Parallel Ice Sheet Model, simulating the flow of a three-dimensional, inherently buttressed ice-sheet-shelf system which periodically surges on a millennial timescale. The thermomechanically coupled model on 1 km horizontal resolution includes an enthalpy-based formulation of the thermodynamics, a nonlinear stress-balance-based sliding law and a very simple subglacial hydrology. The simulated unforced surging is characterized by rapid ice streaming through a bed trough, resulting in abrupt discharge of ice across the grounding line which is eventually calved into the ocean. We visualize the central feedbacks that dominate the subsequent phases of ice buildup, surge and stabilization which emerge from the interaction between ice dynamics, thermodynamics and the subglacial till layer. Results from the variation of surface mass balance and basal roughness suggest that ice sheets of medium thickness may be more susceptible to surging than relatively thin or thick ones for which the surge feedback loop is damped. We also investigate the influence of different basal sliding laws (ranging from purely plastic to nonlinear to linear) on possible surging. The presented mechanisms underlying our simulations of self-maintained, periodic ice growth and destabilization may play a role in large-scale ice-sheet surging, such as the surging of the Laurentide Ice Sheet, which is associated with Heinrich events, and ice-stream shutdown and reactivation, such as observed in the Siple Coast region of West Antarctica.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (approximate to 210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation-age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model-data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016).
This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4 +/- 4.1m (or 6.5 +/- 2.0 x 10(6) km(3)), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.
Increasing Earth’s surface air temperature yields an intensification of its hydrological cycle. As a consequence, the risk of river floods will increase regionally within the next two decades due to the atmospheric warming caused by past anthropogenic greenhouse gas emissions. The direct economic losses caused by these floods can yield regionally heterogeneous losses and gains by propagation within the global trade and supply network. Here we show that, in the absence of large-scale structural adaptation, the total economic losses due to fluvial floods will increase in the next 20 years globally by 17% despite partial compensation through market adjustment within the global trade network. China will suffer the strongest direct losses, with an increase of 82%. The United States is mostly affected indirectly through its trade relations. By contrast to the United States, recent intensification of the trade relations with China leaves the European Union better prepared for the import of production losses in the future.
Damage due to tropical cyclones accounts for more than 50% of all meteorologically-induced economic losses worldwide. Their nominal impact is projected to increase substantially as the exposed population grows, per capita income increases, and anthropogenic climate change manifests. So far, historical losses due to tropical cyclones have been found to increase less than linearly with a nation's affected gross domestic product (GDP). Here we show that for the United States this scaling is caused by a sub-linear increase with affected population while relative losses scale super-linearly with per capita income. The finding is robust across a multitude of empirically derived damage models that link the storm's wind speed, exposed population, and per capita GDP to reported losses. The separation of both socio-economic predictors strongly affects the projection of potential future hurricane losses. Separating the effects of growth in population and per-capita income, per hurricane losses with respect to national GDP are projected to triple by the end of the century under unmitigated climate change, while they are estimated to decrease slightly without the separation.
Changes in ice discharge from Antarctica constitute the largest uncertainty in future sea-level projections, mainly because of the unknown response of its marine basins(1). Most of West Antarctica's marine ice sheet lies on an inland-sloping bed(2) and is thereby prone to a marine ice sheet instability(3-5). A similar topographic configuration is found in large parts of East Antarctica, which holds marine ice equivalent to 19 m of global sea-level rise(6), that is, more than five times that of West Antarctica. Within East Antarctica, the Wilkes Basin holds the largest volume of marine ice that is fully connected by subglacial troughs. This ice body was significantly reduced during the Pliocene epoch(7). Strong melting underneath adjacent ice shelves with similar bathymetry(8) indicates the ice sheet's sensitivity to climatic perturbations. The stability of the Wilkes marine ice sheet has not been the subject of any comprehensive assessment of future sea level. Using recently improved topographic data(6) in combination with ice-dynamic simulations, we show here that the removal of a specific coastal ice volume equivalent to less than 80 mm of global sea-level rise at the margin of the Wilkes Basin destabilizes the regional ice flow and leads to a self-sustained discharge of the entire basin and a global sea-level rise of 3-4 m. Our results are robust with respect to variation in ice parameters, forcing details and model resolution as well as increased surface mass balance, indicating that East Antarctica may become a large contributor to future sea-level rise on timescales beyond a century.
Ten ice-sheet models are used to study sensitivity of the Greenland and Antarctic ice sheets to prescribed changes of surface mass balance, sub-ice-shelf melting and basal sliding. Results exhibit a large range in projected contributions to sea-level change. In most cases, the ice volume above flotation lost is linearly dependent on the strength of the forcing. Combinations of forcings can be closely approximated by linearly summing the contributions from single forcing experiments, suggesting that nonlinear feedbacks are modest. Our models indicate that Greenland is more sensitive than Antarctica to likely atmospheric changes in temperature and precipitation, while Antarctica is more sensitive to increased ice-shelf basal melting. An experiment approximating the Intergovernmental Panel on Climate Change's RCP8.5 scenario produces additional first-century contributions to sea level of 22.3 and 8.1 cm from Greenland and Antarctica, respectively, with a range among models of 62 and 14 cm, respectively. By 200 years, projections increase to 53.2 and 26.7 cm, respectively, with ranges of 79 and 43 cm. Linear interpolation of the sensitivity results closely approximates these projections, revealing the relative contributions of the individual forcings on the combined volume change and suggesting that total ice-sheet response to complicated forcings over 200 years can be linearized.
Anthropogenic climate change is likely to cause continuing global sea level rise(1), but some processes within the Earth system may mitigate the magnitude of the projected effect. Regional and global climate models simulate enhanced snowfall over Antarctica, which would provide a direct offset of the future contribution to global sea level rise from cryospheric mass loss(2,3) and ocean expansion(4). Uncertainties exist in modelled snowfall(5), but even larger uncertainties exist in the potential changes of dynamic ice discharge from Antarctica(1,6) and thus in the ultimate fate of the precipitation-deposited ice mass. Here we show that snowfall and discharge are not independent, but that future ice discharge will increase by up to three times as a result of additional snowfall under global warming. Our results, based on an ice-sheet model(7) forced by climate simulations through to the end of 2500 (ref. 8), show that the enhanced discharge effect exceeds the effect of surface warming as well as that of basal ice-shelf melting, and is due to the difference in surface elevation change caused by snowfall on grounded versus floating ice. Although different underlying forcings drive ice loss from basal melting versus increased snowfall, similar ice dynamical processes are nonetheless at work in both; therefore results are relatively independent of the specific representation of the transition zone. In an ensemble of simulations designed to capture ice-physics uncertainty, the additional dynamic ice loss along the coastline compensates between 30 and 65 per cent of the ice gain due to enhanced snowfall over the entire continent. This results in a dynamic ice loss of up to 1.25 metres in the year 2500 for the strongest warming scenario. The reported effect thus strongly counters a potential negative contribution to global sea level by the Antarctic Ice Sheet.
The Sea-level Response to Ice Sheet Evolution (SeaRISE) effort explores the sensitivity of the current generation of ice sheet models to external forcing to gain insight into the potential future contribution to sea level from the Greenland and Antarctic ice sheets. All participating models simulated the ice sheet response to three types of external forcings: a change in oceanic condition, a warmer atmospheric environment, and enhanced basal lubrication. Here an analysis of the spatial response of the Greenland ice sheet is presented, and the impact of model physics and spin-up on the projections is explored. Although the modeled responses are not always homogeneous, consistent spatial trends emerge from the ensemble analysis, indicating distinct vulnerabilities of the Greenland ice sheet. There are clear response patterns associated with each forcing, and a similar mass loss at the full ice sheet scale will result in different mass losses at the regional scale, as well as distinct thickness changes over the ice sheet. All forcings lead to an increased mass loss for the coming centuries, with increased basal lubrication and warmer ocean conditions affecting mainly outlet glaciers, while the impacts of atmospheric forcings affect the whole ice sheet.
The initiation of a marine ice-sheet instability (MISI) is generally discussed from the ocean side of the ice sheet. It has been shown that the reduction in ice-shelf buttressing and softening of the coastal ice can destabilize a marine ice sheet if the bedrock is sloping upward towards the ocean. Using a conceptional flow-line geometry, we investigate the possibility of whether a MISI can be triggered from the direction of the ice divide as opposed to coastal forcing and explore the interaction between connected basins. We find that the initiation of a MISI in one basin can induce a destabilization in the other. The underlying mechanism of basin interaction is based on dynamic thinning and a consecutive motion of the ice divide which induces a thinning in the adjacent basin and a successive initiation of the instability. Our simplified and symmetric topographic setup allows scaling both the geometry and the transition time between both instabilities. We find that the ice profile follows a universal shape that is scaled with the horizontal extent of the ice sheet and that the same exponent of 1/2 applies for the scaling relation between central surface elevation and horizontal extent as in the pure shallow ice approximation (Vialov profile). Altering the central bed elevation, we find that the extent of grounding-line retreat in one basin determines the degree of interaction with the other. Different scenarios of basin interaction are discussed based on our modeling results as well as on a conceptual flux-balance analysis. We conclude that for the three-dimensional case, the possibility of drainage basin interaction on timescales on the order of 1 kyr or larger cannot be excluded and hence needs further investigation.
Increasing greenhouse gas emissions are likely to impact not only natural systems but economies worldwide. If these impacts alter future economic development, the financial losses will be significantly higher than the mere direct damages. So far, potentially aggravating investment responses were considered negligible. Here we consistently incorporate an empirically derived temperature-growth relation into the simple integrated assessment model DICE. In this framework we show that, if in the next eight decades varying temperatures impact economic growth as has been observed in the past three decades, income is reduced by similar to 20% compared to an economy unaffected by climate change. Hereof similar to 40% are losses due to growth effects of which similar to 50% result from reduced incentive to invest. This additional income loss arises from a reduced incentive for future investment in anticipation of a reduced return and not from an explicit climate protection policy. Under economically optimal climate-change mitigation, however, optimal investment would only be reduced marginally as mitigation efforts keep returns high.
The Atlantic subpolar gyre (SPG) is one of the main drivers of decadal climate variability in the North Atlantic. Here we analyze its dynamics in pre-industrial control simulations of 19 different comprehensive coupled climate models. The analysis is based on a recently proposed description of the SPG dynamics that found the circulation to be potentially bistable due to a positive feedback mechanism including salt transport and enhanced deep convection in the SPG center. We employ a statistical method to identify multiple equilibria in time series that are subject to strong noise and analyze composite fields to assess whether the bistability results from the hypothesized feedback mechanism. Because noise dominates the time series in most models, multiple circulation modes can unambiguously be detected in only six models. Four of these six models confirm that the intensification is caused by the positive feedback mechanism.
Antarctica’s contribution to global sea-level rise has recently been increasing1. Whether its ice discharge will become unstable and decouple from anthropogenic forcing2,3,4 or increase linearly with the warming of the surrounding ocean is of fundamental importance5. Under unabated greenhouse-gas emissions, ocean models indicate an abrupt intrusion of warm circumpolar deep water into the cavity below West Antarctica’s Filchner–Ronne ice shelf within the next two centuries6,7. The ice basin’s retrograde bed slope would allow for an unstable ice-sheet retreat8, but the buttressing of the large ice shelf and the narrow glacier troughs tend to inhibit such instability9,10,11. It is unclear whether future ice loss will be dominated by ice instability or anthropogenic forcing. Here we show in regional and continental-scale ice-sheet simulations, which are capable of resolving unstable grounding-line retreat, that the sea-level response of the Filchner–Ronne ice basin is not dominated by ice instability and follows the strength of the forcing quasi-linearly. We find that the ice loss reduces after each pulse of projected warm water intrusion. The long-term sea-level contribution is approximately proportional to the total shelf-ice melt. Although the local instabilities might dominate the ice loss for weak oceanic warming12, we find that the upper limit of ice discharge from the region is determined by the forcing and not by the marine ice-sheet instability.
The Atlantic meridional overturning circulation (AMOC) carries large amounts of heat into the North Atlantic influencing climate regionally as well as globally. Palaeo-records and simulations with comprehensive climate models suggest that the positive salt-advection feedback may yield a threshold behaviour of the system. That is to say that beyond a certain amount of freshwater flux into the North Atlantic, no meridional overturning circulation can be sustained. Concepts of monitoring the AMOC and identifying its vicinity to the threshold rely on the fact that the volume flux defining the AMOC will be reduced when approaching the threshold. Here we advance conceptual models that have been used in a paradigmatic way to understand the AMOC, by introducing a density-dependent parameterization for the Southern Ocean eddies. This additional degree of freedom uncovers a mechanism by which the AMOC can increase with additional freshwater flux into the North Atlantic, before it reaches the threshold and collapses: an AMOC that is mainly wind-driven will have a constant upwelling as long as the Southern Ocean winds do not change significantly. The downward transport of tracers occurs either in the northern sinking regions or through Southern Ocean eddies. If freshwater is transported, either atmospherically or via horizontal gyres, from the low to high latitudes, this would reduce the eddy transport and by continuity increase the northern sinking which defines the AMOC until a threshold is reached at which the AMOC cannot be sustained. If dominant in the real ocean this mechanism would have significant consequences for monitoring the AMOC.
Non-linear intensification of Sahel rainfall as a possible dynamic response to future warming
(2017)
There is growing empirical evidence that anthropogenic climate change will substantially affect the electric sector. Impacts will stem both from the supply sidethrough the mitigation of greenhouse gasesand from the demand sidethrough adaptive responses to a changing environment. Here we provide evidence of a polarization of both peak load and overall electricity consumption under future warming for the worlds third-largest electricity marketthe 35 countries of Europe. We statistically estimate country-level doseresponse functions between daily peak/total electricity load and ambient temperature for the period 2006-2012. After removing the impact of nontemperature confounders and normalizing the residual load data for each country, we estimate a common doseresponse function, which we use to compute national electricity loads for temperatures that lie outside each countrys currently observed temperature range. To this end, we impose end-of-century climate on todays European economies following three different greenhouse-gas concentration trajectories, ranging from ambitious climate-change mitigationin line with the Paris agreementto unabated climate change. We find significant increases in average daily peak load and overall electricity consumption in southern and western Europe (similar to 3 to similar to 7% for Portugal and Spain) and significant decreases in northern Europe (similar to-6 to similar to-2% for Sweden and Norway). While the projected effect on European total consumption is nearly zero, the significant polarization and seasonal shifts in peak demand and consumption have important ramifications for the location of costly peak-generating capacity, transmission infrastructure, and the design of energy-efficiency policy and storage capacity.
The Paris Climate Agreement aims to keep temperature rise well below 2 degrees C. This implies mitigation costs as well as avoided climate damages. Here we show that independent of the normative assumptions of inequality aversion and time preferences, the agreement constitutes the economically optimal policy pathway for the century. To this end we consistently incorporate a damage-cost curve reproducing the observed relation between temperature and economic growth into the integrated assessment model DICE. We thus provide an intertemporally optimizing cost-benefit analysis of this century's climate problem. We account for uncertainties regarding the damage curve, climate sensitivity, socioeconomic future, and mitigation costs. The resulting optimal temperature is robust as can be understood from the generic temperature-dependence of the mitigation costs and the level of damages inferred from the observed temperature-growth relationship. Our results show that the politically motivated Paris Climate Agreement also represents the economically favourable pathway, if carried out properly.
Forests play a key role in a bio-based economy by providing renewable materials, mitigating climate change, and accommodating biodiversity. However, forests experience massive increases in stresses in their ecological and socioeconomic environments, threatening forest ecosystem services supply. Alleviating those stresses is hampered by conflicting and disconnected governance arrangements, competing interests and claims, and rapid changes in technology and social demands. Identifying which stresses threaten forest ecosystem services supply and which factors hamper their alleviation requires stakeholders' perceptions. Stakeholder-oriented stress tests for the supply of forest ecosystem services are therefore necessary but are not yet available. This perspective presents a roadmap to develop a stress test tailored to multiple stakeholders' needs and demands across spatial scales. We provide the Cascade and Resilience Rosetta, with accompanying performance- and resilience indicators, as tools to facilitate development of the stress test. The application of the stress test will facilitate the transition toward a bio-based economy in which healthy and diverse forests provide sustainable and resilient ecosystem services.
After the United Kingdom has left the European Union it remains unclear whether the two parties can successfully negotiate and sign a trade agreement within the transition period. Ongoing negotiations, practical obstacles and resulting uncertainties make it highly unlikely that economic actors would be fully prepared to a “no-trade-deal” situation. Here we provide an economic shock simulation of the immediate aftermath of such a post-Brexit no-trade-deal scenario by computing the time evolution of more than 1.8 million interactions between more than 6,600 economic actors in the global trade network. We find an abrupt decline in the number of goods produced in the UK and the EU. This sudden output reduction is caused by drops in demand as customers on the respective other side of the Channel incorporate the new trade restriction into their decision-making. As a response, producers reduce prices in order to stimulate demand elsewhere. In the short term consumers benefit from lower prices but production value decreases with potentially severe socio-economic consequences in the longer term.
After the United Kingdom has left the European Union it remains unclear whether the two parties can successfully negotiate and sign a trade agreement within the transition period. Ongoing negotiations, practical obstacles and resulting uncertainties make it highly unlikely that economic actors would be fully prepared to a “no-trade-deal” situation. Here we provide an economic shock simulation of the immediate aftermath of such a post-Brexit no-trade-deal scenario by computing the time evolution of more than 1.8 million interactions between more than 6,600 economic actors in the global trade network. We find an abrupt decline in the number of goods produced in the UK and the EU. This sudden output reduction is caused by drops in demand as customers on the respective other side of the Channel incorporate the new trade restriction into their decision-making. As a response, producers reduce prices in order to stimulate demand elsewhere. In the short term consumers benefit from lower prices but production value decreases with potentially severe socio-economic consequences in the longer term.
The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66% range: 0.02-0.14 m; 90% range: 0.0-0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66% range: 0.04-0.21 m; 90% range: 0.01-0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these values increase to 0.09 m (66% range: 0.04-0.17 m; 90% range: 0.02-0.25 m) for RCP-2.6 and 0.15 m (66% range: 0.07-0.28 m; 90% range: 0.04-0.43 m) for RCP-8.5. All probability distributions are highly skewed towards high values. The applied ice-sheet models are coarse resolution with limitations in the representation of grounding-line motion. Within the constraints of the applied methods, the uncertainty induced from different ice-sheet models is smaller than that induced by the external forcing to the ice sheets.
The radiation model is a parameter-free model of human mobility that has been applied primarily for short-distance moves, such as commuting. When applied to migration, it underestimates the number of long-range moves, such as between different US states. Here we show that it additionally suffers from a conceptual inconsistency that can have substantial numerical effects on long-distance moves.
We propose a modification of the radiation model that introduces a dependence on the angle between any two alternative potential destinations, accounting for the possibility that migrants may have preferences about the approximate direction of their move.
We demonstrate that this modification mitigates the conceptual inconsistency and improves the model fit to observational migration data, without introducing any fitting parameters.
Regional and sectoral disaggregation of multi-regional input-output tables - a flexible algorithm
(2015)
A common shortcoming of available multi-regional input-output (MRIO) data sets is their lack of regional and sectoral detail required for many research questions (e.g. in the field of disaster impact analysis). We present a simple algorithm to refine MRIO tables regionally and/or sectorally. By the use of proxy data, each MRIO flow in question is disaggregated into the corresponding sub-flows. This downscaling procedure is complemented by an adjustment rule ensuring that the sub-flows match the superordinate flow in sum. The approximation improves along several iteration steps. The algorithm unfolds its strength through the flexible combination of multiple, possibly incomplete proxy data sources. It is also flexible in a sense that any target sector and region resolution can be chosen. As an exemplary case we apply the algorithm to a regional and sectoral refinement of the Eora MRIO database.
Due to climate change the frequency and character of precipitation are changing as the hydrological cycle intensifies. With regards to snowfall, global warming has two opposing influences; increasing humidity enables intense snowfall, whereas higher temperatures decrease the likelihood of snowfall. Here we show an intensification of extreme snowfall across large areas of the Northern Hemisphere under future warming. This is robust across an ensemble of global climate models when they are bias-corrected with observational data. While mean daily snowfall decreases, both the 99th and the 99.9th percentiles of daily snowfall increase in many regions in the next decades, especially for Northern America and Asia. Additionally, the average intensity of snowfall events exceeding these percentiles as experienced historically increases in many regions. This is likely to pose a challenge to municipalities in mid to high latitudes. Overall, extreme snowfall events are likely to become an increasingly important impact of climate change in the next decades, even if they will become rarer, but not necessarily less intense, in the second half of the century.
Recently a multitude of empirically derived damage models have been applied to project future tropical cyclone (TC) losses for the United States. In their study (Geiger et al 2016 Environ. Res. Lett. 11 084012) compared two approaches that differ in the scaling of losses with socio-economic drivers: the commonly-used approach resulting in a sub-linear scaling of historical TC losses with a nation's affected gross domestic product (GDP), and the disentangled approach that shows a sub-linear increase with affected population and a super-linear scaling of relative losses with per capita income. Statistics cannot determine which approach is preferable but since process understanding demands that there is a dependence of the loss on both GDP per capita and population, an approach that accounts for both separately is preferable to one which assumes a specific relation between the two dependencies. In the accompanying comment, Rybski et al argued that there is no rigorous evidence to reach the conclusion that high-income does not protect against hurricane losses. Here we affirm that our conclusion is drawn correctly and reply to further remarks raised in the comment, highlighting the adequateness of our approach but also the potential for future extension of our research.
The dynamics of the North Atlantic subpolar gyre (SPG) are assessed under present and glacial boundary conditions by investigating the SPG sensitivity to surface wind-stress changes in a coupled climate model. To this end, the gyre transport is decomposed in Ekman, thermohaline, and bottom transports. Surface wind-stress variations are found to play an important indirect role in SPG dynamics through their effect on water-mass densities. Our results suggest the existence of two dynamically distinct regimes of the SPG, depending on the absence or presence of deep water formation (DWF) in the Nordic Seas and a vigorous Greenland-Scotland ridge (GSR) overflow. In the first regime, the GSR overflow is weak and the SPG strength increases with wind-stress as a result of enhanced outcropping of isopycnals in the centre of the SPG. As soon as a vigorous GSR overflow is established, its associated positive density anomalies on the southern GSR slope reduce the SPG strength. This has implications for past glacial abrupt climate changes, insofar as these can be explained through latitudinal shifts in North Atlantic DWF sites and strengthening of the North Atlantic current. Regardless of the ultimate trigger, an abrupt shift of DWF into the Nordic Seas could result both in a drastic reduction of the SPG strength and a sudden reversal in its sensitivity to wind-stress variations. Our results could provide insight into changes in the horizontal ocean circulation during abrupt glacial climate changes, which have been largely neglected up to now in model studies.