SUMMARY: The joint ESA/NASA Mass-change And Geosciences International Constellation (MAGIC) has the objective to extend time-series from previous gravity missions, including an improvement of accuracy and spatio-temporal resolution. The long-term monitoring of Earth’s gravity field carries information on mass change induced by water cycle, climate change and mass transport processes between atmosphere, cryosphere, oceans and solid Earth. MAGIC will be composed of two satellite pairs flying in different orbit planes. The NASA/DLR-led first pair (P1) is expected to be in a near-polar orbit around 500 km of altitude; while the second ESA-led pair (P2) is expected to be in an inclined orbit of 65°–70° at approximately 400 km altitude. The ESA-led pair P2 Next Generation Gravity Mission shall be launched after P1 in a staggered manner to form the MAGIC constellation. The addition of an inclined pair shall lead to reduction of temporal aliasing effects and consequently of reliance on de-aliasing models and post-processing. The main novelty of the MAGIC constellation is the delivery of mass-change products at higher spatial resolution, temporal (i.e. subweekly) resolution, shorter latency and higher accuracy than the Gravity Recovery and Climate Experiment (GRACE) and Gravity Recovery and Climate Experiment Follow-On (GRACE-FO). This will pave the way to new science applications and operational services. In this paper, an overview of various fields of science and service applications for hydrology, cryosphere, oceanography, solid Earth, climate change and geodesy is provided. These thematic fields and newly enabled applications and services were analysed in the frame of the initial ESA Science Support activities for MAGIC. The analyses of MAGIC scenarios for different application areas in the field of geosciences confirmed that the double-pair configuration will significantly enlarge the number of observable mass-change phenomena by resolving smaller spatial scales with an uncertainty that satisfies evolved user requirements expressed by international bodies such as IUGG. The required uncertainty levels of dedicated thematic fields met by MAGIC unfiltered Level-2 products will benefit hydrological applications by recovering more than 90 per cent of the major river basins worldwide at 260 km spatial resolution, cryosphere applications by enabling mass change signal separation in the interior of Greenland from those in the coastal zones and by resolving small-scale mass variability in challenging regions such as the Antarctic Peninsula, oceanography applications by monitoring meridional overturning circulation changes on timescales of years and decades, climate applications by detecting amplitude and phase changes of Terrestrial Water Storage after 30 yr in 64 and 56 per cent of the global land areas and solid Earth applications by lowering the Earthquake detection threshold from magnitude 8.8 to magnitude 7.4 with spatial resolution increased to 333 km.

The AlpArray Gravity Research Group (AAGRG), as part of the European AlpArray program, focuses on the compilation of a homogeneous surface-based gravity data set across the Alpine area. In 2017 10 European countries in the Alpine realm agreed to contribute with gravity data for a new compilation of the Alpine gravity field in an area spanning from 2 to 23^∘ E and from 41 to 51^∘ N. This compilation relies on existing national gravity databases and, for the Ligurian and the Adriatic seas, on shipborne data of the Service Hydrographique et Océanographique de la Marine and of the Bureau Gravimétrique International. Furthermore, for the Ivrea zone in the Western Alps, recently acquired data were added to the database. This first pan-Alpine gravity data map is homogeneous regarding input data sets, applied methods and all corrections, as well as reference frames. Here, the AAGRG presents the data set of the recalculated gravity fields on a 4 km × 4 km grid for public release and a 2 km × 2 km grid for special request. The final products also include calculated values for mass and bathymetry corrections of the measured gravity at each grid point, as well as height. This allows users to use later customized densities for their own calculations of mass corrections. Correction densities used are 2670 kg m⁻³ for landmasses, 1030 kg m⁻³ for water masses above the ellipsoid and −1640 kg m⁻³ for those below the ellipsoid and 1000 kg m⁻³ for lake water masses. The correction radius was set to the Hayford zone O₂ (167 km). The new Bouguer anomaly is station completed (CBA) and compiled according to the most modern criteria and reference frames (both positioning and gravity), including atmospheric corrections. Special emphasis was put on the gravity effect of the numerous lakes in the study area, which can have an effect of up to 5 mGal for gravity stations located at shorelines with steep slopes, e.g., for the rather deep reservoirs in the Alps. The results of an error statistic based on cross validations and/or “interpolation residuals” are provided for the entire database. As an example, the interpolation residuals of the Austrian data set range between about −8 and +8 mGal and the cross-validation residuals between −14 and +10 mGal; standard deviations are well below 1 mGal. The accuracy of the newly compiled gravity database is close to ±5 mGal for most areas. A first interpretation of the new map shows that the resolution of the gravity anomalies is suited for applications ranging from intra-crustal- to crustal-scale modeling to interdisciplinary studies on the regional and continental scales, as well as applications as joint inversion with other data sets. The data are published with the DOI https://doi.org/10.5880/fidgeo.2020.045 (Zahorec et al., 2021) via GFZ Data Services.

In the past twenty years, satellite gravimetry missions have successfully provided data for the determination of the Earth static gravity field (GOCE) and its temporal variations (GRACE and GRACE-FO). In particular, the possibility to study the evolution in time of Earth masses allows us to monitor global parameters underlying climate changes, water resources, flooding, melting of ice masses and the corresponding global sea level rise, all of which are of paramount importance, providing basic data on, e.g. geodynamics, earthquakes, hydrology or ice sheets changes. Recently, a large interest has developed in novel technologies and quantum sensing, which promise higher sensitivity, drift-free measurements, and higher absolute accuracy for both terrestrial surveys and space missions, giving direct access to more precise long-term measurements. Looking at a time frame beyond the present decade, in the MOCAST+ study (MOnitoring mass variations by Cold Atom Sensors and Time measures) a satellite mission based on an “enhanced” quantum payload is proposed, with cold atom interferometers acting as gravity gradiometers, and atomic clocks for optical frequency measurements, providing observations of differences of the gravitational potential. The main outcomes are the definition of the accuracy level to be expected from this payload and the accuracy level needed to detect and monitor phenomena identified in the Scientific Challenges of the ESA Living Planet Program, in particular Cryosphere, Ocean and Solid Earth. In this paper, the proposed payload, mission profile and preliminary platform design are presented, with end-to-end simulation results and assessment of the impact on geophysical applications.

SUMMARY: Since the completion of the Gravity field and steady-state Ocean Circulation Explorer mission (GOCE), global gravity models of uniform quality and coverage are available. We investigate their potential of being useful tools for estimating the thermal structure of the continental lithosphere, through simulation and real-data test in Central-Eastern Europe across the Trans-European Suture Zone. Heat flow, measured near the Earth surface, is the result of the superposition of a complex set of contributions, one of them being the heat production occurring in the crust. The crust is enriched in radioactive elements respect to the underlying mantle and crustal thickness is an essential parameter in isolating the thermal contribution of the crust. Obtaining reliable estimates of crustal thickness through inversion of GOCE-derived gravity models has already proven feasible, especially when weak constraints from other observables are introduced. We test a way to integrate this in a geothermal framework, building a 3-D, steady state, solid Earth conductive heat transport model, from the lithosphere–asthenosphere boundary to the surface. This thermal model is coupled with a crust-mantle boundary depth resulting from inverse modelling, after correcting the gravity model for the effects of topography, far-field isostatic roots and sediments. We employ a mixed space- and spectral-domain based forward modelling strategy to ensure full spectral coherency between the limited spectral content of the gravity model and the reductions. Deviations from a direct crustal thickness to crustal heat production relationship are accommodated using a subsequent substitution scheme, constrained by surface heat flow measurements, where available. The result is a 3-D model of the lithosphere characterised in temperature, radiogenic heat and thermal conductivity. It provides added information respect to the lithospheric structure and sparse heat flow measurements alone, revealing a satisfactory coherence with the geological features in the area and their controlling effect on the conductive heat transport.

The quantum technology absolute gravimeters, gradiometers, and clocks are at the forefront of the instrumentation to be exploited in a future gravity mission (the QSG mission concept). Apart from the quantum payload, the mission design defines the choice of the number of satellites and the satellite orbit constellation, with the goal of optimizing the observation of the earth’s gravity field and reducing aliasing phenomena. Our goal is to define the realistic gravity field changes generated by glaciers and lakes and define the sensitivity of the quantum gravity mission for the detection of hydrologic and cryospheric mass changes. The analysis focuses on mass changes in the high mountains of Asia and the South American continent. The mass changes are based on terrestrial and satellite observations and are of a climatic origin. We show that compared to the existing GRACE-FO mission, a quantum gravity mission significantly improves the detection of the climatic mass gain of lakes and mass loss of glaciers, allowing for smaller mass features to be distinguished, and smaller mass losses to be detected. The greater signal is the seasonal signal with a yearly period, which would be detected at the 10 Gt level for areas > 8000 km². The yearly mass loss of the Patagonian glaciers can be detected at the 5 Gt/yr level, an improvement from the 10 Gt/yr detectable by GRACE-FO. Spatial resolution would also be improved, with an increase of about 50% in spatial frequency for the detection of the mass change rate of lakes and glaciers in Tibet. The improved spatial resolution enables an improved localization of the lakes and glaciers affected by climatic mass change. The results will contribute to defining the user requirements of the future QSG missions.