@article{SteffenRoeckstromRichardsonetal.2018, author = {Steffen, Will and R{\"o}ckstrom, Johan and Richardson, Katherine and Lenton, Timothy M. and Folke, Carl and Liverman, Diana and Summerhayes, Colin P. and Barnosky, Anthony D. and Cornell, Sarah E. and Crucifix, Michel and Donges, Jonathan and Fetzer, Ingo and Lade, Steven J. and Scheffer, Marten and Winkelmann, Ricarda and Schellnhuber, Hans Joachim}, title = {Trajectories of the Earth System in the Anthropocene}, series = {Proceedings of the National Academy of Sciences of the United States of America}, volume = {115}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {33}, publisher = {National Acad. of Sciences}, address = {Washington}, issn = {0027-8424}, doi = {10.1073/pnas.1810141115}, pages = {8252 -- 8259}, year = {2018}, abstract = {We explore the risk that self-reinforcing feedbacks could push the Earth System toward a planetary threshold that, if crossed, could prevent stabilization of the climate at intermediate temperature rises and cause continued warming on a "Hothouse Earth" pathway even as human emissions are reduced. Crossing the threshold would lead to a much higher global average temperature than any interglacial in the past 1.2 million years and to sea levels significantly higher than at any time in the Holocene. We examine the evidence that such a threshold might exist and where it might be. If the threshold is crossed, the resulting trajectory would likely cause serious disruptions to ecosystems, society, and economies. Collective human action is required to steer the Earth System away from a potential threshold and stabilize it in a habitable interglacial-like state. Such action entails stewardship of the entire Earth System-biosphere, climate, and societies-and could include decarbonization of the global economy, enhancement of biosphere carbon sinks, behavioral changes, technological innovations, new governance arrangements, and transformed social values.}, language = {en} } @article{ReeseWinkelmannGudmundsson2018, author = {Reese, Ronja and Winkelmann, Ricarda and Gudmundsson, Gudmundur Hilmar}, title = {Grounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {12}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {10}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-12-3229-2018}, pages = {3229 -- 3242}, year = {2018}, abstract = {Currently, several large-scale ice-flow models impose a condition on ice flux across grounding lines using an analytically motivated parameterisation of grounding-line flux. It has been suggested that employing this analytical expression alleviates the need for highly resolved computational domains around grounding lines of marine ice sheets. While the analytical flux formula is expected to be accurate in an unbuttressed flow-line setting, its validity has hitherto not been assessed for complex and realistic geometries such as those of the Antarctic Ice Sheet. Here the accuracy of this analytical flux formula is tested against an optimised ice flow model that uses a highly resolved computational mesh around the Antarctic grounding lines. We find that when applied to the Antarctic Ice Sheet the analytical expression provides inaccurate estimates of ice fluxes for almost all grounding lines. Furthermore, in many instances direct application of the analytical formula gives rise to unphysical complex-valued ice fluxes. We conclude that grounding lines of the Antarctic Ice Sheet are, in general, too highly buttressed for the analytical parameterisation to be of practical value for the calculation of grounding-line fluxes.}, language = {en} } @article{CiemerRehmKurthsetal.2020, author = {Ciemer, Catrin and Rehm, Lars and Kurths, J{\"u}rgen and Donner, Reik Volker and Winkelmann, Ricarda and Boers, Niklas}, title = {An early-warning indicator for Amazon droughts exclusively based on tropical Atlantic sea surface temperatures}, series = {Environmental Research Letters}, volume = {15}, journal = {Environmental Research Letters}, number = {9}, publisher = {IOP - Institute of Physics Publishing}, address = {Bristol}, pages = {10}, year = {2020}, abstract = {Droughts in tropical South America have an imminent and severe impact on the Amazon rainforest and affect the livelihoods of millions of people. Extremely dry conditions in Amazonia have been previously linked to sea surface temperature (SST) anomalies in the adjacent tropical oceans. Although the sources and impacts of such droughts have been widely studied, establishing reliable multi-year lead statistical forecasts of their occurrence is still an ongoing challenge. Here, we further investigate the relationship between SST and rainfall anomalies using a complex network approach. We identify four ocean regions which exhibit the strongest overall SST correlations with central Amazon rainfall, including two particularly prominent regions in the northern and southern tropical Atlantic. Based on the time-dependent correlation between SST anomalies in these two regions alone, we establish a new early-warning method for droughts in the central Amazon basin and demonstrate its robustness in hindcasting past major drought events with lead-times up to 18 months.}, language = {en} } @article{SeroussiNowickiPayneetal.2020, author = {Seroussi, Helene and Nowicki, Sophie and Payne, Antony J. and Goelzer, Heiko and Lipscomb, William H. and Abe-Ouchi, Ayako and Agosta, Cecile and Albrecht, Torsten and Asay-Davis, Xylar and Barthel, Alice and Calov, Reinhard and Cullather, Richard and Dumas, Christophe and Galton-Fenzi, Benjamin K. and Gladstone, Rupert and Golledge, Nicholas R. and Gregory, Jonathan M. and Greve, Ralf and Hattermann, Tore and Hoffman, Matthew J. and Humbert, Angelika and Huybrechts, Philippe and Jourdain, Nicolas C. and Kleiner, Thomas and Larour, Eric and Leguy, Gunter R. and Lowry, Daniel P. and Little, Chistopher M. and Morlighem, Mathieu and Pattyn, Frank and Pelle, Tyler and Price, Stephen F. and Quiquet, Aurelien and Reese, Ronja and Schlegel, Nicole-Jeanne and Shepherd, Andrew and Simon, Erika and Smith, Robin S. and Straneo, Fiammetta and Sun, Sainan and Trusel, Luke D. and Van Breedam, Jonas and van de Wal, Roderik S. W. and Winkelmann, Ricarda and Zhao, Chen and Zhang, Tong and Zwinger, Thomas}, title = {ISMIP6 Antarctica}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {9}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-3033-2020}, pages = {3033 -- 3070}, year = {2020}, abstract = {Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015-2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between 7:8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between 6 :1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.}, language = {en} } @article{ReeseLevermannAlbrechtetal.2020, author = {Reese, Ronja and Levermann, Anders and Albrecht, Torsten and Seroussi, Helene and Winkelmann, Ricarda}, title = {The role of history and strength of the oceanic forcing in sea level projections from Antarctica with the Parallel Ice Sheet Model}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {9}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-3097-2020}, pages = {3097 -- 3110}, year = {2020}, abstract = {Mass loss from the Antarctic Ice Sheet constitutes the largest uncertainty in projections of future sea level rise. Ocean-driven melting underneath the floating ice shelves and subsequent acceleration of the inland ice streams are the major reasons for currently observed mass loss from Antarctica and are expected to become more important in the future. Here we show that for projections of future mass loss from the Antarctic Ice Sheet, it is essential (1) to better constrain the sensitivity of sub-shelf melt rates to ocean warming and (2) to include the historic trajectory of the ice sheet. In particular, we find that while the ice sheet response in simulations using the Parallel Ice Sheet Model is comparable to the median response of models in three Antarctic Ice Sheet Intercomparison projects - initMIP, LARMIP-2 and ISMIP6 - conducted with a range of ice sheet models, the projected 21st century sea level contribution differs significantly depending on these two factors. For the highest emission scenario RCP8.5, this leads to projected ice loss ranging from 1:4 to 4:0 cm of sea level equivalent in simulations in which ISMIP6 ocean forcing drives the PICO ocean box model where parameter tuning leads to a comparably low sub-shelf melt sensitivity and in which no surface forcing is applied. This is opposed to a likely range of 9:1 to 35:8 cm using the exact same initial setup, but emulated from the LARMIP-2 experiments with a higher melt sensitivity, even though both projects use forcing from climate models and melt rates are calibrated with previous oceanographic studies. Furthermore, using two initial states, one with a previous historic simulation from 1850 to 2014 and one starting from a steady state, we show that while differences between the ice sheet configurations in 2015 seem marginal at first sight, the historic simulation increases the susceptibility of the ice sheet to ocean warming, thereby increasing mass loss from 2015 to 2100 by 5\% to 50 \%. Hindcasting past ice sheet changes with numerical models would thus provide valuable tools to better constrain projections. Our results emphasize that the uncertainty that arises from the forcing is of the same order of magnitude as the ice dynamic response for future sea level projections.}, language = {en} } @article{ZeitzLevermannWinkelmann2020, author = {Zeitz, Maria and Levermann, Anders and Winkelmann, Ricarda}, title = {Sensitivity of ice loss to uncertainty in flow law parameters in an idealized one-dimensional geometry}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {10}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-3537-2020}, pages = {3537 -- 3550}, year = {2020}, abstract = {Acceleration of the flow of ice drives mass losses in both the Antarctic and the Greenland Ice Sheet. The projections of possible future sea-level rise rely on numerical ice-sheet models, which solve the physics of ice flow, melt, and calving. While major advancements have been made by the ice-sheet modeling community in addressing several of the related uncertainties, the flow law, which is at the center of most process-based ice-sheet models, is not in the focus of the current scientific debate. However, recent studies show that the flow law parameters are highly uncertain and might be different from the widely accepted standard values. Here, we use an idealized flow-line setup to investigate how these uncertainties in the flow law translate into uncertainties in flow-driven mass loss. In order to disentangle the effect of future warming on the ice flow from other effects, we perform a suite of experiments with the Parallel Ice Sheet Model (PISM), deliberately excluding changes in the surface mass balance. We find that changes in the flow parameters within the observed range can lead up to a doubling of the flow-driven mass loss within the first centuries of warming, compared to standard parameters. The spread of ice loss due to the uncertainty in flow parameters is on the same order of magnitude as the increase in mass loss due to surface warming. While this study focuses on an idealized flow-line geometry, it is likely that this uncertainty carries over to realistic three-dimensional simulations of Greenland and Antarctica.}, language = {en} } @article{ReeseGudmundssonLevermannetal.2017, author = {Reese, Ronja and Gudmundsson, Gudmundur Hilmar and Levermann, Anders and Winkelmann, Ricarda}, title = {The far reach of ice-shelf thinning in Antarctica}, series = {Nature climate change}, volume = {8}, journal = {Nature climate change}, number = {1}, publisher = {Nature Publ. Group}, address = {London}, issn = {1758-678X}, doi = {10.1038/s41558-017-0020-x}, pages = {53 -- 57}, year = {2017}, abstract = {Floating ice shelves, which fringe most of Antarctica's coastline, regulate ice flow into the Southern Ocean1,2,3. Their thinning4,5,6,7 or disintegration8,9 can cause upstream acceleration of grounded ice and raise global sea levels. So far the effect has not been quantified in a comprehensive and spatially explicit manner. Here, using a finite-element model, we diagnose the immediate, continent-wide flux response to different spatial patterns of ice-shelf mass loss. We show that highly localized ice-shelf thinning can reach across the entire shelf and accelerate ice flow in regions far from the initial perturbation. As an example, this 'tele-buttressing' enhances outflow from Bindschadler Ice Stream in response to thinning near Ross Island more than 900 km away. We further find that the integrated flux response across all grounding lines is highly dependent on the location of imposed changes: the strongest response is caused not only near ice streams and ice rises, but also by thinning, for instance, well-within the Filchner-Ronne and Ross Ice Shelves. The most critical regions in all major ice shelves are often located in regions easily accessible to the intrusion of warm ocean waters10,11,12, stressing Antarctica's vulnerability to changes in its surrounding ocean.}, language = {en} } @article{AlbrechtWinkelmannLevermann2020, author = {Albrecht, Torsten and Winkelmann, Ricarda and Levermann, Anders}, title = {Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {2}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-633-2020}, pages = {633 -- 656}, year = {2020}, abstract = {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.}, language = {en} } @article{GarbeAlbrechtLevermannetal.2020, author = {Garbe, Julius and Albrecht, Torsten and Levermann, Anders and Donges, Jonathan and Winkelmann, Ricarda}, title = {The hysteresis of the Antarctic Ice Sheet}, series = {Nature : the international weekly journal of science}, volume = {585}, journal = {Nature : the international weekly journal of science}, number = {7826}, publisher = {Macmillan Publishers Limited}, address = {Berlin}, issn = {0028-0836}, doi = {10.1038/s41586-020-2727-5}, pages = {538 -- 544}, year = {2020}, abstract = {More than half of Earth's freshwater resources are held by the Antarctic Ice Sheet, which thus represents by far the largest potential source for global sea-level rise under future warming conditions(1). Its long-term stability determines the fate of our coastal cities and cultural heritage. Feedbacks between ice, atmosphere, ocean, and the solid Earth give rise to potential nonlinearities in its response to temperature changes. So far, we are lacking a comprehensive stability analysis of the Antarctic Ice Sheet for different amounts of global warming. Here we show that the Antarctic Ice Sheet exhibits a multitude of temperature thresholds beyond which ice loss is irreversible. Consistent with palaeodata(2)we find, using the Parallel Ice Sheet Model(3-5), that at global warming levels around 2 degrees Celsius above pre-industrial levels, West Antarctica is committed to long-term partial collapse owing to the marine ice-sheet instability. Between 6 and 9 degrees of warming above pre-industrial levels, the loss of more than 70 per cent of the present-day ice volume is triggered, mainly caused by the surface elevation feedback. At more than 10 degrees of warming above pre-industrial levels, Antarctica is committed to become virtually ice-free. The ice sheet's temperature sensitivity is 1.3 metres of sea-level equivalent per degree of warming up to 2 degrees above pre-industrial levels, almost doubling to 2.4 metres per degree of warming between 2 and 6 degrees and increasing to about 10 metres per degree of warming between 6 and 9 degrees. Each of these thresholds gives rise to hysteresis behaviour: that is, the currently observed ice-sheet configuration is not regained even if temperatures are reversed to present-day levels. In particular, the West Antarctic Ice Sheet does not regrow to its modern extent until temperatures are at least one degree Celsius lower than pre-industrial levels. Our results show that if the Paris Agreement is not met, Antarctica's long-term sea-level contribution will dramatically increase and exceed that of all other sources.
Modelling shows that the Antarctic Ice Sheet exhibits multiple temperature thresholds beyond which ice loss would become irreversible, and once melted, the ice sheet can regain its previous mass only if the climate cools well below pre-industrial temperatures.}, language = {en} } @article{SchlemmFeldmannWinkelmannetal.2022, author = {Schlemm, Tanja and Feldmann, Johannes and Winkelmann, Ricarda and Levermann, Anders}, title = {Stabilizing effect of melange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {16}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-16-1979-2022}, pages = {1979 -- 1996}, year = {2022}, abstract = {Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice melange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6 m within 100 a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice melange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-melange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet.}, language = {en} } @phdthesis{Bastian2023, author = {Bastian, Martin}, title = {An emergent machine learning approach for seasonal cyclone activity forecasts}, school = {Universit{\"a}t Potsdam}, pages = {135}, year = {2023}, abstract = {Seasonal forecasts are of great interest in many areas. Knowing the amount of precipitation for the upcoming season in regions of water scarcity would facilitate a better water management. If farmers knew the weather conditions of the upcoming summer at sowing time, they could select those cereal species that are best adapted to these conditions. This would allow farmers to improve the harvest and potentially even reduce the amount of pesticides used. However, the undoubted advantages of seasonal forecasts are often opposed by their high degree of uncertainty. The great challenge of generating seasonal forecasts with lead times of several months mainly originates from the chaotic nature of the earth system. In a chaotic system, even tiny differences in the initial conditions can lead to strong deviations in the system's state in the long run. In this dissertation we propose an emergent machine learning approach for seasonal forecasting, called the AnlgModel. The AnlgModel combines the analogue method with myopic feature selection and bootstrapping. To benchmark the abilities of the AnlgModel we apply it to seasonal cyclone activity forecasts in the North Atlantic and Northwest Pacific. The AnlgModel demonstrates competitive hindcast skills with two operational forecasts and even outperforms these for long lead times. In the second chapter we comprehend the forecasting strategy of the Anlg-Model. We thereby analyse the analogue selection process for the 2017 North Atlantic and the 2018 Northwest Pacific seasonal cyclone activity. The analysis shows that those climate indices which are known to influence the seasonal cyclone activity, such as the Ni{\~n}o 3.4 SST, are correctly represented among the selected analogues. Furthermore the selected analogues reflect large-scale climate patterns that were identified by expert reports as being determinative for these particular seasons. In the third chapter we analyse the features that are used by the AnlgModel for its predictions. We therefore inspect the feature relevance (FR). The FR patterns learned by the AnlgModel show a high congruence with the predictor regions used by the operational forecasts. However, the AnlgModel also discovered new features, such as the SST anomaly in the Gulf of Guinea during November. This SST pattern exhibits a remarkably high predictive potential for the upcoming Atlantic hurricane activity. In the final chapter we investigate potential mechanisms, that link two of these regions with high feature relevance to the Atlantic hurricane activity. We mainly focus on ocean surface transport. The ocean surface flow paths are calculated using Lagrangian particle analysis. We demonstrate that the FR patterns in the region of the Canary islands do not correspond with ocean surface transport. It is instead likely that these FR patterns fingerprint a wind transport of latent heat. The second region to be studied is situated in the Gulf of Guinea. Our analysis shows that the FR patterns seen there do fingerprint ocean surface transport. However, our simulations also show that at least one other mechanism is involved in linking the Gulf of Guinea SST anomaly in November to the hurricane activity of the upcoming season. In this work the AnlgModel does not only demonstrate its outstanding forecast skills but also shows its capabilities as research tool for detecting oceanic and atmospheric mechanisms.}, language = {en} } @article{FeldmannReeseWinkelmannetal.2022, author = {Feldmann, Johannes and Reese, Ronja and Winkelmann, Ricarda and Levermann, Anders}, title = {Shear-margin melting causes stronger transient ice discharge than ice-stream melting in idealized simulations}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {16}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-16-1927-2022}, pages = {1927 -- 1940}, year = {2022}, abstract = {Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet, the melting increases the ice discharge into the ocean. Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an isothermal, inherently buttressed ice-sheet-shelf system. In the idealized numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's (1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and (2) shear margins, where the ice flow is slower. The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction. Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming.}, language = {en} } @article{WunderlingWilleitDongesetal.2020, author = {Wunderling, Nico and Willeit, Matteo and Donges, Jonathan and Winkelmann, Ricarda}, title = {Global warming due to loss of large ice masses and Arctic summer sea ice}, series = {Nature Communications}, volume = {11}, journal = {Nature Communications}, number = {1}, publisher = {Nature Publishing Group}, address = {Berlin}, issn = {2041-1723}, doi = {10.1038/s41467-020-18934-3}, pages = {14}, year = {2020}, abstract = {Several large-scale cryosphere elements such as the Arctic summer sea ice, the mountain glaciers, the Greenland and West Antarctic Ice Sheet have changed substantially during the last century due to anthropogenic global warming. However, the impacts of their possible future disintegration on global mean temperature (GMT) and climate feedbacks have not yet been comprehensively evaluated. Here, we quantify this response using an Earth system model of intermediate complexity. Overall, we find a median additional global warming of 0.43 degrees C (interquartile range: 0.39-0.46 degrees C) at a CO2 concentration of 400 ppm. Most of this response (55\%) is caused by albedo changes, but lapse rate together with water vapour (30\%) and cloud feedbacks (15\%) also contribute significantly. While a decay of the ice sheets would occur on centennial to millennial time scales, the Arctic might become ice-free during summer within the 21st century. Our findings imply an additional increase of the GMT on intermediate to long time scales. The disintegration of cryosphere elements such as the Arctic summer sea ice, mountain glaciers, Greenland and West Antarctica is associated with temperature and radiative feedbacks. In this work, the authors quantify these feedbacks and find an additional global warming of 0.43 degrees C.}, language = {en} } @article{KloseWunderlingWinkelmannetal.2021, author = {Klose, Ann Kristin and Wunderling, Nico and Winkelmann, Ricarda and Donges, Jonathan}, title = {What do we mean, 'tipping cascade'?}, series = {Environmental research letters : ERL}, volume = {16}, journal = {Environmental research letters : ERL}, number = {12}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1748-9326}, doi = {10.1088/1748-9326/ac3955}, pages = {11}, year = {2021}, abstract = {Based on suggested interactions of potential tipping elements in the Earth's climate and in ecological systems, tipping cascades as possible dynamics are increasingly discussed and studied. The activation of such tipping cascades would impose a considerable risk for human societies and biosphere integrity. However, there are ambiguities in the description of tipping cascades within the literature so far. Here we illustrate how different patterns of multiple tipping dynamics emerge from a very simple coupling of two previously studied idealized tipping elements. In particular, we distinguish between a two phase cascade, a domino cascade and a joint cascade. A mitigation of an unfolding two phase cascade may be possible and common early warning indicators are sensitive to upcoming critical transitions to a certain degree. In contrast, a domino cascade may hardly be stopped once initiated and critical slowing down-based indicators fail to indicate tipping of the following element. These different potentials for intervention and anticipation across the distinct patterns of multiple tipping dynamics should be seen as a call to be more precise in future analyses of cascading dynamics arising from tipping element interactions in the Earth system.}, language = {en} } @article{KunertPangTewsetal.2022, author = {Kunert, Nina and Pang, Peter T. H. and Tews, Ingo and Coughlin, Michael W. and Dietrich, Tim}, title = {Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources}, series = {Physical review D}, volume = {105}, journal = {Physical review D}, number = {6}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0010}, doi = {10.1103/PhysRevD.105.L061301}, pages = {7}, year = {2022}, abstract = {With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections is combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyze them with a set of four different waveform models. For our analysis with 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to 1250 m (2\% at 90\% credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.}, language = {en} } @misc{BullaCoughlinDhawanetal.2022, author = {Bulla, Mattia and Coughlin, Michael W. and Dhawan, Suhail and Dietrich, Tim}, title = {Multi-messenger constraints on the Hubble constant through combination of gravitational waves, gamma-ray bursts and kilonovae from neutron star mergers}, series = {Universe : open access journal}, volume = {8}, journal = {Universe : open access journal}, number = {5}, publisher = {MDPI}, address = {Basel}, issn = {2218-1997}, doi = {10.3390/universe8050289}, pages = {21}, year = {2022}, abstract = {The simultaneous detection of gravitational waves and light from the binary neutron star merger GW170817 led to independent measurements of distance and redshift, providing a direct estimate of the Hubble constant H-0 that does not rely on a cosmic distance ladder, nor assumes a specific cosmological model. By using gravitational waves as "standard sirens", this approach holds promise to arbitrate the existing tension between the H-0 value inferred from the cosmic microwave background and those obtained from local measurements. However, the known degeneracy in the gravitational-wave analysis between distance and inclination of the source led to a H-0 value from GW170817 that was not precise enough to resolve the existing tension. In this review, we summarize recent works exploiting the viewing-angle dependence of the electromagnetic signal, namely the associated short gamma-ray burst and kilonova, to constrain the system inclination and improve on H-0. We outline the key ingredients of the different methods, summarize the results obtained in the aftermath of GW170817 and discuss the possible systematics introduced by each of these methods.}, language = {en} } @article{NiCaoShpritsetal.2018, author = {Ni, Binbin and Cao, Xing and Shprits, Yuri and Summers, Danny and Gu, Xudong and Fu, Song and Lou, Yuequn}, title = {Hot Plasma Effects on the Cyclotron-Resonant Pitch-Angle Scattering Rates of Radiation Belt Electrons Due to EMIC Waves}, series = {Geophysical research letters}, volume = {45}, journal = {Geophysical research letters}, number = {1}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2017GL076028}, pages = {21 -- 30}, year = {2018}, abstract = {To investigate the hot plasma effects on the cyclotron-resonant interactions between electromagnetic ion cyclotron (EMIC) waves and radiation belt electrons in a realistic magnetospheric environment, calculations of the wave-induced bounce-averaged pitch angle diffusion coefficients are performed using both the cold and hot plasma dispersion relations. The results demonstrate that the hot plasma effects have a pronounced influence on the electron pitch angle scattering rates due to all three EMIC emission bands (H+, He+, and O+) when the hot plasma dispersion relation deviates significantly from the cold plasma approximation. For a given wave spectrum, the modification of the dispersion relation by hot anisotropic protons can strongly increase the minimum resonant energy for electrons interacting with O+ band EMIC waves, while the minimum resonant energies for H+ and He+ bands are not greatly affected. For H+ band EMIC waves, inclusion of hot protons tends to weaken the pitch angle scattering efficiency of >5MeV electrons. The most crucial differences introduced by the hot plasma effects occur for >3MeV electron scattering rates by He+ band EMIC waves. Mainly due to the changes of resonant frequency and wave group velocity when the hot protons are included, the difference in scattering rates can be up to an order of magnitude, showing a strong dependence on both electron energy and equatorial pitch angle. Our study confirms the importance of including hot plasma effects in modeling the scattering of ultra-relativistic radiation belt electrons by EMIC waves.}, language = {en} } @article{BernardiBerdjaDaniGuzmanetal.2021, author = {Bernardi, Rafael L. and Berdja, Amokrane and Dani Guzman, Christian and Torres-Torriti, Miguel and Roth, Martin M.}, title = {Restoration of images with a spatially varying PSF of the T80-S telescope optical model using neural networks}, series = {Monthly notices of the Royal Astronomical Society}, volume = {510}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stab3400}, pages = {4284 -- 4294}, year = {2021}, abstract = {Most image restoration methods in astronomy rely upon probabilistic tools that infer the best solution for a deconvolution problem. They achieve good performances when the point spread function (PSF) is spatially invariant in the image plane. However, this condition is not always satisfied in real optical systems. We propose a new method for the restoration of images affected by static and anisotropic aberrations using Deep Neural Networks that can be directly applied to sky images. The network is trained using simulated sky images corresponding to the T80-S Telescope optical model, a 80-cm survey imager at Cerro Tololo (Chile), which are synthesized using a Zernike polynomial representation of the optical system. Once trained, the network can be used directly on sky images, outputting a corrected version of the image that has a constant and known PSF across its field of view. The method is to be tested on the T80-S Telescope. We present the method and results on synthetic data.}, language = {en} } @article{GuoNiFuetal.2022, author = {Guo, Yingjie and Ni, Binbin and Fu, Song and Wang, Dedong and Shprits, Yuri and Zhelavskaya, Irina and Feng, Minghang and Guo, Deyu}, title = {Identification of controlling geomagnetic and solar wind factors for magnetospheric chorus intensity using feature selection techniques}, series = {Journal of geophysical research : A, Space physics}, volume = {127}, journal = {Journal of geophysical research : A, Space physics}, number = {1}, publisher = {Wiley}, address = {Hoboken, NJ}, issn = {2169-9380}, doi = {10.1029/2021JA029926}, pages = {14}, year = {2022}, abstract = {Using over-5-year EMFISIS wave measurements from Van Allen Probes, we present a detailed survey to identify the controlling factors among the geomagnetic indices and solar wind parameters for the 1-min root mean square amplitudes of lower band chorus (LBC) and upper band chorus (UBC). A set of important features are automatically determined by feature selection techniques, namely, Random Forest and Maximum Relevancy Minimum Redundancy. Our analysis results indicate the AE index with zero-time-delay dominates the intensity evolution of LBC and UBC, consistent with the evidence that chorus waves prefer to occur and amplify during enhanced substorm periods. Regarding solar wind parameters, solar wind speed and IMF B-z are identified as the controlling factors for chorus wave intensity. Using the combination of all these important features, a predictive neural network model of chorus wave intensity is established to reconstruct the temporal variations of chorus wave intensity, for which application of Random Forest produces the overall best performance. Plain Language Summary Whistler mode chorus waves are electromagnetic waves observed in the low-density region near the geomagnetic equator outside the plasmapause. The dynamics of Earth's radiation belts are largely influenced by chorus waves owing to their dual contributions to both radiation belt electron acceleration and loss. In this study, we use feature selection techniques to identify the controlling geomagnetic and solar wind factors for magnetospheric chorus waves. Feature selection techniques implement the processes which can select the features most influential to the output. In this study, the inputs are geomagnetic indices and solar wind parameters and the output is the chorus wave intensity. The results indicate that AE index with zerotime delay dominates the chorus wave intensity. Furthermore, solar wind speed and IMF B-z are identified as the most important solar wind drivers for chorus wave intensity. On basis of the combination of all these important geomagnetic and solar wind controlling factors, we develop a neural network model of chorus wave intensity, and find that the model with the inputs identified using the Random Forest method produces the overall best performance.}, language = {en} } @phdthesis{Smirnov2023, author = {Smirnov, Artem}, title = {Understanding the dynamics of the near-earth space environment utilizing long-term satellite observations}, doi = {10.25932/publishup-61371}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-613711}, school = {Universit{\"a}t Potsdam}, pages = {xxxvi, 286}, year = {2023}, abstract = {The near-Earth space environment is a highly complex system comprised of several regions and particle populations hazardous to satellite operations. The trapped particles in the radiation belts and ring current can cause significant damage to satellites during space weather events, due to deep dielectric and surface charging. Closer to Earth is another important region, the ionosphere, which delays the propagation of radio signals and can adversely affect navigation and positioning. In response to fluctuations in solar and geomagnetic activity, both the inner-magnetospheric and ionospheric populations can undergo drastic and sudden changes within minutes to hours, which creates a challenge for predicting their behavior. Given the increasing reliance of our society on satellite technology, improving our understanding and modeling of these populations is a matter of paramount importance. In recent years, numerous spacecraft have been launched to study the dynamics of particle populations in the near-Earth space, transforming it into a data-rich environment. To extract valuable insights from the abundance of available observations, it is crucial to employ advanced modeling techniques, and machine learning methods are among the most powerful approaches available. This dissertation employs long-term satellite observations to analyze the processes that drive particle dynamics, and builds interdisciplinary links between space physics and machine learning by developing new state-of-the-art models of the inner-magnetospheric and ionospheric particle dynamics. The first aim of this thesis is to investigate the behavior of electrons in Earth's radiation belts and ring current. Using ~18 years of electron flux observations from the Global Positioning System (GPS), we developed the first machine learning model of hundreds-of-keV electron flux at Medium Earth Orbit (MEO) that is driven solely by solar wind and geomagnetic indices and does not require auxiliary flux measurements as inputs. We then proceeded to analyze the directional distributions of electrons, and for the first time, used Fourier sine series to fit electron pitch angle distributions (PADs) in Earth's inner magnetosphere. We performed a superposed epoch analysis of 129 geomagnetic storms during the Van Allen Probes era and demonstrated that electron PADs have a strong energy-dependent response to geomagnetic activity. Additionally, we showed that the solar wind dynamic pressure could be used as a good predictor of the PAD dynamics. Using the observed dependencies, we created the first PAD model with a continuous dependence on L, magnetic local time (MLT) and activity, and developed two techniques to reconstruct near-equatorial electron flux observations from low-PA data using this model. The second objective of this thesis is to develop a novel model of the topside ionosphere. To achieve this goal, we collected observations from five of the most widely used ionospheric missions and intercalibrated these data sets. This allowed us to use these data jointly for model development, validation, and comparison with other existing empirical models. We demonstrated, for the first time, that ion density observations by Swarm Langmuir Probes exhibit overestimation (up to ~40-50\%) at low and mid-latitudes on the night side, and suggested that the influence of light ions could be a potential cause of this overestimation. To develop the topside model, we used 19 years of radio occultation (RO) electron density profiles, which were fitted with a Chapman function with a linear dependence of scale height on altitude. This approximation yields 4 parameters, namely the peak density and height of the F2-layer and the slope and intercept of the linear scale height trend, which were modeled using feedforward neural networks (NNs). The model was extensively validated against both RO and in-situ observations and was found to outperform the International Reference Ionosphere (IRI) model by up to an order of magnitude. Our analysis showed that the most substantial deviations of the IRI model from the data occur at altitudes of 100-200 km above the F2-layer peak. The developed NN-based ionospheric model reproduces the effects of various physical mechanisms observed in the topside ionosphere and provides highly accurate electron density predictions. This dissertation provides an extensive study of geospace dynamics, and the main results of this work contribute to the improvement of models of plasma populations in the near-Earth space environment.}, language = {en} }