TY - JOUR A1 - Kreuzer, Moritz A1 - Reese, Ronja A1 - Huiskamp, Willem Nicholas A1 - Petri, Stefan A1 - Albrecht, Torsten A1 - Feulner, Georg A1 - Winkelmann, Ricarda T1 - Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOMS (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain JF - Geoscientific model development : an interactive open access journal of the European Geosciences Union N2 - The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Pots-dam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3 degrees respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice-ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups. Y1 - 2021 U6 - https://doi.org/10.5194/gmd-14-3697-2021 SN - 1991-959X SN - 1991-9603 VL - 14 IS - 6 SP - 3697 EP - 3714 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Zeitz, Maria A1 - Haacker, Jan M. A1 - Donges, Jonathan A1 - Albrecht, Torsten A1 - Winkelmann, Ricarda T1 - Dynamic regimes of the Greenland Ice Sheet emerging from interacting melt-elevation and glacial isostatic adjustment feedbacks JF - Earth system dynamics N2 - The stability of the Greenland Ice Sheet under global warming is governed by a number of dynamic processes and interacting feedback mechanisms in the ice sheet, atmosphere and solid Earth. Here we study the long-term effects due to the interplay of the competing melt-elevation and glacial isostatic adjustment (GIA) feedbacks for different temperature step forcing experiments with a coupled ice-sheet and solid-Earth model. Our model results show that for warming levels above 2 degrees C, Greenland could become essentially ice-free within several millennia, mainly as a result of surface melting and acceleration of ice flow. These ice losses are mitigated, however, in some cases with strong GIA feedback even promoting an incomplete recovery of the Greenland ice volume. We further explore the full-factorial parameter space determining the relative strengths of the two feedbacks: our findings suggest distinct dynamic regimes of the Greenland Ice Sheets on the route to destabilization under global warming - from incomplete recovery, via quasi-periodic oscillations in ice volume to ice-sheet collapse. In the incomplete recovery regime, the initial ice loss due to warming is essentially reversed within 50 000 years, and the ice volume stabilizes at 61 %-93 % of the present-day volume. For certain combinations of temperature increase, atmospheric lapse rate and mantle viscosity, the interaction of the GIA feedback and the melt-elevation feedback leads to self-sustained, long-term oscillations in ice-sheet volume with oscillation periods between 74 000 and over 300 000 years and oscillation amplitudes between 15 %-70 % of present-day ice volume. This oscillatory regime reveals a possible mode of internal climatic variability in the Earth system on timescales on the order of 100 000 years that may be excited by or synchronized with orbital forcing or interact with glacial cycles and other slow modes of variability. Our findings are not meant as scenario-based near-term projections of ice losses but rather providing insight into of the feedback loops governing the "deep future" and, thus, long-term resilience of the Greenland Ice Sheet. Y1 - 2022 U6 - https://doi.org/10.5194/esd-13-1077-2022 SN - 2190-4979 SN - 2190-4987 VL - 13 IS - 3 SP - 1077 EP - 1096 PB - Copernicus Publ. CY - Göttingen ER - TY - JOUR A1 - Albrecht, Torsten A1 - Winkelmann, Ricarda A1 - Levermann, Anders T1 - Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM) BT - part 2: parameter ensemble analysis JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2020 U6 - https://doi.org/10.5194/tc-14-633-2020 SN - 1994-0416 SN - 1994-0424 VL - 14 IS - 2 SP - 633 EP - 656 PB - Copernicus Publ. CY - Göttingen ER - TY - JOUR A1 - Schlemm, Tanja A1 - Feldmann, Johannes A1 - Winkelmann, Ricarda A1 - Levermann, Anders T1 - Stabilizing effect of melange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2022 U6 - https://doi.org/10.5194/tc-16-1979-2022 SN - 1994-0416 SN - 1994-0424 VL - 16 IS - 5 SP - 1979 EP - 1996 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Feldmann, Johannes A1 - Reese, Ronja A1 - Winkelmann, Ricarda A1 - Levermann, Anders T1 - Shear-margin melting causes stronger transient ice discharge than ice-stream melting in idealized simulations JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2022 U6 - https://doi.org/10.5194/tc-16-1927-2022 SN - 1994-0416 SN - 1994-0424 VL - 16 IS - 5 SP - 1927 EP - 1940 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Wunderling, Nico A1 - Willeit, Matteo A1 - Donges, Jonathan A1 - Winkelmann, Ricarda T1 - Global warming due to loss of large ice masses and Arctic summer sea ice JF - Nature Communications N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1038/s41467-020-18934-3 SN - 2041-1723 VL - 11 IS - 1 PB - Nature Publishing Group CY - Berlin ER - TY - JOUR A1 - Zeitz, Maria A1 - Reese, Ronja A1 - Beckmann, Johanna A1 - Krebs-Kanzow, Uta A1 - Winkelmann, Ricarda T1 - Impact of the melt-albedo feedback on the future evolution of the Greenland Ice Sheet with PISM-dEBM-simple JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - Surface melting of the Greenland Ice Sheet contributes a large amount to current and future sea level rise. Increased surface melt may lower the reflectivity of the ice sheet surface and thereby increase melt rates: the so-called melt-albedo feedback describes this self-sustaining increase in surface melting. In order to test the effect of the melt-albedo feedback in a prognostic ice sheet model, we implement dEBM-simple, a simplified version of the diurnal Energy Balance Model dEBM, in the Parallel Ice Sheet Model (PISM). The implementation includes a simple representation of the melt-albedo feedback and can thereby replace the positive-degree-day melt scheme. Using PISM-dEBM-simple, we find that this feedback increases ice loss through surface warming by 60 % until 2300 for the high-emission scenario RCP8.5 when compared to a scenario in which the albedo remains constant at its present-day values. With an increase of 90 % compared to a fixed-albedo scenario, the effect is more pronounced for lower surface warming under RCP2.6. Furthermore, assuming an immediate darkening of the ice surface over all summer months, we estimate an upper bound for this effect to be 70 % in the RCP8.5 scenario and a more than 4-fold increase under RCP2.6. With dEBM-simple implemented in PISM, we find that the melt-albedo feedback is an essential contributor to mass loss in dynamic simulations of the Greenland Ice Sheet under future warming. Y1 - 2021 U6 - https://doi.org/10.5194/tc-15-5739-2021 SN - 1994-0416 SN - 1994-0424 VL - 15 IS - 12 SP - 5739 EP - 5764 PB - Copernicus CY - Katlenburg-Lindau ER - TY - JOUR A1 - Klose, Ann Kristin A1 - Wunderling, Nico A1 - Winkelmann, Ricarda A1 - Donges, Jonathan T1 - What do we mean, 'tipping cascade'? JF - Environmental research letters : ERL N2 - 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. KW - tipping cascade KW - domino effect KW - tipping interactions KW - cascading regime KW - shifts KW - early warning indicators Y1 - 2021 U6 - https://doi.org/10.1088/1748-9326/ac3955 SN - 1748-9326 VL - 16 IS - 12 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Durand, Gael A1 - van den Broeke, Michiel R. A1 - Le Cozannet, Goneri A1 - Edwards, Tamsin L. A1 - Holland, Paul R. A1 - Jourdain, Nicolas C. A1 - Marzeion, Ben A1 - Mottram, Ruth A1 - Nicholls, Robert J. A1 - Pattyn, Frank A1 - Paul, Frank A1 - Slangen, Aimee B. A. A1 - Winkelmann, Ricarda A1 - Burgard, Clara A1 - van Calcar, Caroline J. A1 - Barre, Jean-Baptiste A1 - Bataille, Amelie A1 - Chapuis, Anne T1 - Sea-Level rise: from global perspectives to local services JF - Frontiers in Marine Science N2 - Coastal areas are highly diverse, ecologically rich, regions of key socio-economic activity, and are particularly sensitive to sea-level change. Over most of the 20th century, global mean sea level has risen mainly due to warming and subsequent expansion of the upper ocean layers as well as the melting of glaciers and ice caps. Over the last three decades, increased mass loss of the Greenland and Antarctic ice sheets has also started to contribute significantly to contemporary sea-level rise. The future mass loss of the two ice sheets, which combined represent a sea-level rise potential of similar to 65 m, constitutes the main source of uncertainty in long-term (centennial to millennial) sea-level rise projections. Improved knowledge of the magnitude and rate of future sea-level change is therefore of utmost importance. Moreover, sea level does not change uniformly across the globe and can differ greatly at both regional and local scales. The most appropriate and feasible sea level mitigation and adaptation measures in coastal regions strongly depend on local land use and associated risk aversion. Here, we advocate that addressing the problem of future sea-level rise and its impacts requires (i) bringing together a transdisciplinary scientific community, from climate and cryospheric scientists to coastal impact specialists, and (ii) interacting closely and iteratively with users and local stakeholders to co-design and co-build coastal climate services, including addressing the high-end risks. KW - sea-level rise KW - Antarctic KW - Greenland KW - glaciers KW - local impact Y1 - 2022 U6 - https://doi.org/10.3389/fmars.2021.709595 SN - 2296-7745 VL - 8 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Pang, Peter Tsun Ho A1 - Dietrich, Tim A1 - Tews, Ingo A1 - Van Den Broeck, Chris T1 - Parameter estimation for strong phase transitions in supranuclear matter using gravitational-wave astronomy JF - Physical review research N2 - At supranuclear densities, explored in the core of neutron stars, a strong phase transition from hadronic matter to more exotic forms of matter might be present. To test this hypothesis, binary neutron-star mergers offer a unique possibility to probe matter at densities that we cannot create in any existing terrestrial experiment. In this work, we show that, if present, strong phase transitions can have a measurable imprint on the binary neutron-star coalescence and the emitted gravitational-wave signal. We construct a new parametrization of the supranuclear equation of state that allows us to test for the existence of a strong phase transition and extract its characteristic properties purely from the gravitational-wave signal of the inspiraling neutron stars. We test our approach using a Bayesian inference study simulating 600 signals with three different equations of state and find that for current gravitational-wave detector networks already 12 events might be sufficient to verify the presence of a strong phase transition. Finally, we use our methodology to analyze GW170817 and GW190425 but do not find any indication that a strong phase transition is present at densities probed during the inspiral. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevResearch.2.033514 SN - 2643-1564 VL - 2 IS - 3 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Bulla, Mattia A1 - Coughlin, Michael W. A1 - Dhawan, Suhail A1 - Dietrich, Tim T1 - Multi-messenger constraints on the Hubble constant through combination of gravitational waves, gamma-ray bursts and kilonovae from neutron star mergers JF - Universe : open access journal N2 - 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. KW - gravitational waves KW - stars: neutron KW - stars: binaries KW - cosmology: cosmological parameters KW - cosmology: distance scale KW - cosmology: cosmic background radiation Y1 - 2022 U6 - https://doi.org/10.3390/universe8050289 SN - 2218-1997 VL - 8 IS - 5 PB - MDPI CY - Basel ER - TY - JOUR A1 - Smirnov, Artem G. A1 - Kronberg, Elena A. A1 - Daly, Patrick W. A1 - Aseev, Nikita A1 - Shprits, Yuri A1 - Kellerman, Adam C. T1 - Adiabatic Invariants Calculations for Cluster Mission: A Long-Term Product for Radiation Belts Studies JF - Journal of Geophysical Research: Space Physics N2 - The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes. KW - L-Asterisk KW - magnetosphere KW - electrons KW - model Y1 - 2019 VL - 125 IS - 2 PB - John Wiley & Sons, Inc. CY - New Jersey ER - TY - JOUR A1 - Smirnov, Artem A1 - Shprits, Yuri A1 - Allison, Hayley A1 - Aseev, Nikita A1 - Drozdov, Alexander A1 - Kollmann, Peter A1 - Wang, Dedong A1 - Saikin, Anthony T1 - Storm-Time evolution of the Equatorial Electron Pitch Angle Distributions in Earth's Outer Radiation Belt JF - Frontiers in astronomy and space sciences N2 - In this study we analyze the storm-time evolution of equatorial electron pitch angle distributions (PADs) in the outer radiation belt region using observations from the Magnetic Electron Ion Spectrometer (MagEIS) instrument aboard the Van Allen Probes in 2012-2019. The PADs are approximated using a sum of the first, third and fifth sine harmonics. Different combinations of the respective coefficients refer to the main PAD shapes within the outer radiation belt, namely the pancake, flat-top, butterfly and cap PADs. We conduct a superposed epoch analysis of 129 geomagnetic storms and analyze the PAD evolution for day and night MLT sectors. PAD shapes exhibit a strong energy-dependent response. At energies of tens of keV, the PADs exhibit little variation throughout geomagnetic storms. Cap PADs are mainly observed at energies < 300 keV, and their extent in L shrinks with increasing energy. The cap distributions transform into the pancake PADs around the main phase of the storm on the nightside, and then come back to their original shapes during the recovery phase. At higher energies on the dayside, the PADs are mainly pancake during pre-storm conditions and become more anisotropic during the main phase. The quiet-time butterfly PADs can be observed on the nightside at L> 5.6. During the main phase, butterfly PADs have stronger 90 degrees-minima and can be observed at lower L-shells (down to L = 5), then transitioning into flat-top PADs at L similar to 4.5 - 5 and pancake PADs at L < 4.5. The resulting PAD coefficients for different energies, locations and storm epochs can be used to test the wave models and physics-based radiation belt codes in terms of pitch angle distributions. KW - pitch angle KW - pitch angle distributions KW - electrons KW - radiation belts KW - magnetosphere KW - van alien probes Y1 - 2022 U6 - https://doi.org/10.3389/fspas.2022.836811 SN - 2296-987X VL - 9 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Shprits, Yuri A1 - Menietti, J. D. A1 - Drozdov, Alexander A1 - Horne, Richard B. A1 - Woodfield, Emma E. A1 - Groene, J. B. A1 - de Soria-Santacruz, M. A1 - Averkamp, T. F. A1 - Garrett, H. A1 - Paranicas, C. A1 - Gurnett, Don A. T1 - Strong whistler mode waves observed in the vicinity of Jupiter’s moons JF - Nature Communications N2 - Understanding of wave environments is critical for the understanding of how particles are accelerated and lost in space. This study shows that in the vicinity of Europa and Ganymede, that respectively have induced and internal magnetic fields, chorus wave power is significantly increased. The observed enhancements are persistent and exceed median values of wave activity by up to 6 orders of magnitude for Ganymede. Produced waves may have a pronounced effect on the acceleration and loss of particles in the Jovian magnetosphere and other astrophysical objects. The generated waves are capable of significantly modifying the energetic particle environment, accelerating particles to very high energies, or producing depletions in phase space density. Observations of Jupiter’s magnetosphere provide a unique opportunity to observe how objects with an internal magnetic field can interact with particles trapped in magnetic fields of larger scale objects. Y1 - 2018 U6 - https://doi.org/10.1038/s41467-018-05431-x SN - 2041-1723 VL - 9 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Adolfs, Marjolijn A1 - Hoque, Mohammed Mainul A1 - Shprits, Yuri T1 - Storm-time relative total electron content modelling using machine learning techniques JF - Remote sensing N2 - Accurately predicting total electron content (TEC) during geomagnetic storms is still a challenging task for ionospheric models. In this work, a neural-network (NN)-based model is proposed which predicts relative TEC with respect to the preceding 27-day median TEC, during storm time for the European region (with longitudes 30 degrees W-50 degrees E and latitudes 32.5 degrees N-70 degrees N). The 27-day median TEC (referred to as median TEC), latitude, longitude, universal time, storm time, solar radio flux index F10.7, global storm index SYM-H and geomagnetic activity index Hp30 are used as inputs and the output of the network is the relative TEC. The relative TEC can be converted to the actual TEC knowing the median TEC. The median TEC is calculated at each grid point over the European region considering data from the last 27 days before the storm using global ionosphere maps (GIMs) from international GNSS service (IGS) sources. A storm event is defined when the storm time disturbance index Dst drops below 50 nanotesla. The model was trained with storm-time relative TEC data from the time period of 1998 until 2019 (2015 is excluded) and contains 365 storms. Unseen storm data from 33 storm events during 2015 and 2020 were used to test the model. The UQRG GIMs were used because of their high temporal resolution (15 min) compared to other products from different analysis centers. The NN-based model predictions show the seasonal behavior of the storms including positive and negative storm phases during winter and summer, respectively, and show a mixture of both phases during equinoxes. The model's performance was also compared with the Neustrelitz TEC model (NTCM) and the NN-based quiet-time TEC model, both developed at the German Aerospace Agency (DLR). The storm model has a root mean squared error (RMSE) of 3.38 TEC units (TECU), which is an improvement by 1.87 TECU compared to the NTCM, where an RMSE of 5.25 TECU was found. This improvement corresponds to a performance increase by 35.6%. The storm-time model outperforms the quiet-time model by 1.34 TECU, which corresponds to a performance increase by 28.4% from 4.72 to 3.38 TECU. The quiet-time model was trained with Carrington averaged TEC and, therefore, is ideal to be used as an input instead of the GIM derived 27-day median. We found an improvement by 0.8 TECU which corresponds to a performance increase by 17% from 4.72 to 3.92 TECU for the storm-time model using the quiet-time-model predicted TEC as an input compared to solely using the quiet-time model. KW - ionosphere KW - relative total electron content KW - geomagnetic storms KW - neural KW - networks KW - NTCM KW - European storm-time model Y1 - 2022 U6 - https://doi.org/10.3390/rs14236155 SN - 2072-4292 VL - 14 IS - 23 PB - MDPI CY - Basel ER - TY - JOUR A1 - Prol, Fabricio S. A1 - Smirnov, Artem G. A1 - Hoque, M. Mainul A1 - Shprits, Yuri T1 - Combined model of topside ionosphere and plasmasphere derived from radio-occultation and Van Allen Probes data JF - Scientific reports N2 - In the last years, electron density profile functions characterized by a linear dependence on the scale height showed good results when approximating the topside ionosphere. The performance above 800 km, however, is not yet well investigated. This study investigates the capability of the semi-Epstein functions to represent electron density profiles from the peak height up to 20,000 km. Electron density observations recorded by the Van Allen Probes were used to resolve the scale height dependence in the plasmasphere. It was found that the linear dependence of the scale height in the topside ionosphere cannot be directly used to extrapolate profiles above 800 km. We find that the dependence of scale heights on altitude is quadratic in the plasmasphere. A statistical model of the scale heights is therefore proposed. After combining the topside ionosphere and plasmasphere by a unified model, we have obtained good estimations not only in the profile shapes, but also in the Total Electron Content magnitude and distributions when compared to actual measurements from 2013, 2014, 2016 and 2017. Our investigation shows that Van Allen Probes can be merged to radio-occultation data to properly represent the upper ionosphere and plasmasphere by means of a semi-Epstein function. Y1 - 2022 U6 - https://doi.org/10.1038/s41598-022-13302-1 SN - 2045-2322 VL - 12 IS - 1 PB - Macmillan Publishers Limited, part of Springer Nature CY - London ER - TY - JOUR A1 - Haas, Bernhard A1 - Shprits, Yuri A1 - Allison, Hayley A1 - Wutzig, Michael A1 - Wang, Dedong T1 - Which parameter controls ring current electron dynamics JF - Frontiers in astronomy and space sciences N2 - Predicting the electron population of Earth's ring current during geomagnetic storms still remains a challenging task. In this work, we investigate the sensitivity of 10 keV ring current electrons to different driving processes, parameterised by the Kp index, during several moderate and intense storms. Results are validated against measurements from the Van Allen Probes satellites. Perturbing the Kp index allows us to identify the most dominant processes for moderate and intense storms respectively. We find that during moderate storms (Kp < 6) the drift velocities mostly control the behaviour of low energy electrons, while loss from wave-particle interactions is the most critical parameter for quantifying the evolution of intense storms (Kp > 6). Perturbations of the Kp index used to drive the boundary conditions at GEO and set the plasmapause location only show a minimal effect on simulation results over a limited L range. It is further shown that the flux at L & SIM; 3 is more sensitive to changes in the Kp index compared to higher L shells, making it a good proxy for validating the source-loss balance of a ring current model. KW - ring current KW - magnetosphere KW - electron lifetimes KW - electrons KW - van allen probes (RBSP) KW - ring current model KW - verb Y1 - 2022 U6 - https://doi.org/10.3389/fspas.2022.911002 SN - 2296-987X VL - 9 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Woodfield, Emma E. A1 - Horne, Richard B. A1 - Glauert, S. A. A1 - Menietti, J. D. A1 - Shprits, Yuri A1 - Kurth, William S. T1 - Formation of electron radiation belts at Saturn by Z-mode wave acceleration JF - Nature Communications N2 - At Saturn electrons are trapped in the planet’s magnetic field and accelerated to relativistic energies to form the radiation belts, but how this dramatic increase in electron energy occurs is still unknown. Until now the mechanism of radial diffusion has been assumed but we show here that in-situ acceleration through wave particle interactions, which initial studies dismissed as ineffectual at Saturn, is in fact a vital part of the energetic particle dynamics there. We present evidence from numerical simulations based on Cassini spacecraft data that a particular plasma wave, known as Z-mode, accelerates electrons to MeV energies inside 4 RS (1 RS = 60,330 km) through a Doppler shifted cyclotron resonant interaction. Our results show that the Z-mode waves observed are not oblique as previously assumed and are much better accelerators than O-mode waves, resulting in an electron energy spectrum that closely approaches observed values without any transport effects included. Y1 - 2018 U6 - https://doi.org/10.1038/s41467-018-07549-4 SN - 2041-1723 VL - 9 PB - Nature Publ. Group CY - London ER -