The Bohai Sea, located in the innermost part of the East Asian marginal seas, is of particular interest in the studies of relative sea level (RSL) and glacial isostatic adjustment (GIA) due to its apparent far-field position. This study analyzed six cores and one pre-existing archaeological site from the southern Bohai Sea coast, generating 11 sea-level index points (SLIPs) from supratidal, upper tidal, and mid-lower tidal flat sediments using foraminifera tests and sedimentary analysis. All SLIPs were corrected for possible self-compaction, long-term tectonic effect, and the lowering effect due to water extraction. This enabled a high-quality reconstruction of RSL changes over the age range of 9000–3000 cal a BP. The RSL rose rapidly from about −17.19 ± 1.32 m to 1.76 ± 1.32 m MSL between around 9000 and 7000 cal a BP. However, the RSL rise rates quickly decreased from around 8–10 mm/a before 8000 cal a BP to ∼4mm/a by 7000 cal a BP. After 7000 cal a BP, the RSL continued to rise, reaching its peak of 2.44 ± 1.34 m MSL around 6000 cal a BP followed by a gradual decline to the present height at 3000–4000 cal a BP, with a maximum falling rate of 1.2 mm/a around 5000 cal a BP. Comparisons of the reconstructed RSL in this study with those from other non-deltaic regions, including the inner Hangzhou Bay and the surrounding coasts of western Bohai Sea, as well as with GIA predictions by ICE6G, ANU, and region-specific ‘final model’ reveal: (1) a significant and long-term early-to-mid Holocene continental levering effect on the southern Bohai Sea coast; (2) a close fit with the region-specific ‘final model’, albeit with SLIPs that fall slightly below predictions from 9000 to 8000 cal a BP and form a distinct mid-Holocene highstand above model projections after 7000 cal a BP; (3) these data-model misfits may result, in the former case, from the GIA model's Holocene ice melting history not embedding accelerated sea-level rise events, and in the latter, from the presence of a relatively strong upper mantle at the southern Bohai Sea coast, rather than the extremely weak upper mantle assumed for the Bohai Sea region.

The Ganymede Laser Altimeter (GALA) on the Jupiter Icy Moons Explorer (JUICE) mission, is in charge of a comprehensive geodetic mapping of Europa, Ganymede, and Callisto on the basis of Laser range measurements. While multiple topographic profiles will be obtained for Europa and Callisto during flybys, GALA will provide a high-resolution global shape model of Ganymede while in orbit around this moon based on at least 600 million range measurements from altitudes of 500 km and 200 km above the surface. By measuring the diurnal tidal deformation of Ganymede, which crucially depends on the decoupling of the outer ice shell from the deeper interior by a liquid water ocean, GALA will obtain evidence for (or against) a subsurface ocean on Ganymede and will provide constraints on the ice shell thickness above the ocean. In combination with other instruments, it will characterize the morphology of surface units on Ganymede, Europa, and Callisto providing not only topography but also measurements of surface roughness on the scale of the laser footprint, i.e. at a scale of about 50 m from 500 km altitude, and albedo values at the laser wavelength of 1064 nm. GALA is a single-beam laser altimeter, operating at a nominal frequency of 30 Hz, with a capability of reaching up to 48 Hz. It uses a Nd:YAG laser to generate pulses with pulse lengths of 5.5 ± 2.5 ns. The return pulse is detected by an Avalanche Photo Diode (APD) with 100 MHz bandwidth and the signal is digitized at a sampling rate of 200 MHz providing range measurements with a sub-sample resolution of 0.1 m. Research institutes and industrial partners from Germany, Japan, Switzerland and Spain collaborated to build the instrument. JUICE, conducted under responsibility of the European Space Agency (ESA), was successfully launched in April 2023 and is scheduled for arrival at the Jupiter system in July 2031. The nominal science mission including multiple close flybys at Europa, Ganymede, and Callisto, as well as the final Ganymede orbit phase will last from 2031 to 2035. In May 2023 GALA has completed its Near-Earth Commissioning, showing full functionality of all units. Here we summarize the scientific objectives, instrument design and implementation, performance, and operational aspects of GALA.

Meltwater Pulses (MWPs) from the last deglaciation are highly relevant today. Studying their timing and characteristics offers valuable insights into past periods of accelerated sea-level response under intense climate forcing and provides an analog for potential future scenarios. The western Bohai Sea coast (WBSC), located deep within the Asian continent and far from major glaciation centers and subduction zones, exhibits sensitive relative sea level (RSL) responses to any high latitude ice-sheet meltwater influx, offering valuable data to refine the chronology and magnitude of MWPs. This paper presents 25 high-confidence early Holocene sea-level index points (SLIPs) from high marsh vegetation samples, generated through analysis of sediments, microfossil foraminifera, and radiocarbon dating of 14 cores from the WBSC. The altitudes of SLIPs were adjusted for factors such as self-compaction, long-term tectonic subsidence, and ground lowering due to water extraction. The early Holocene RSL history in the WBSC reflects a rapid rise in sea level during ∼9800–6900 cal a BP, occurring within the global context of global sea-level rise. During ∼9800–8100 cal a BP, the RSL rise was faster, averaging 7–8 mm/a, characterized by two distinct stepwise increases. A rapid sea-level rise event before 9500 cal a BP was shown by a sharp RSL rise of at least 5.4 ± 0.63 m, within the narrow timeframe of 9603 ± 288 cal a BP (1σ). The second MWP took place during 8500–8100 cal a BP, adding a total rise of 2.92 ± 0.93 m on the top of background RSL rise, with a two-phase structure where the main phase produced 2.55 ± 0.69 m within a narrow window around 8505 ± 135 cal BP (1σ). Although the rate of RSL rise decreased rapidly to 5.0 mm/a by 7000 cal a BP after ∼8100 cal a BP, sea-level data from the WBSC still show an additional rise of 1.45 ± 0.64 m during 7593 ± 78 cal a BP (1σ). These three distinct sea-level jump events, happened before 9500, ∼8500, and ∼7600 cal a BP, respectively, provide indications of three episodes of enhanced meltwater discharge in the early Holocene, with the latter two closely linked to the decaying history of the Laurentide Ice Sheet.

We present an overview of the operations and engineering interface for Planetary Radio Interferometry and Doppler Experiment (PRIDE) radio astronomy observations as a scientific component of the ESA’s Jupiter Icy Moons Explorer (JUICE) mission, as well as other prospective planetary and space science missions. The article discusses advanced scheduling and planning methods that make it possible to create observing schedules for observations of specific spacecraft in concurrence with observations of natural radio sources. In order to put this into practice and find suitable natural background calibrator sources for PRIDE of JUICE mission, we developed planning and scheduling software. The conventional scheduling software for natural celestial radio sources is not set up to include spacecraft as observation targets in the necessary control files. Therefore, difficulties already arise during observation planning. We report on the development of new and the adaptation of existing routines used in astrophysical and geodetic VLBI for satellite scheduling and planning. The analysis of the PRIDE science observations led to improved observational planning, and the mission’s scheduling methodologies were studied using a systems engineering approach. In addition, we highlighted the new procedures, like finding charts for selecting calibrator radio sources over a range of frequency bands and the outcomes of those strategies for science operation planning. A simulation of the flyby of Venus during the cruise phase of the JUICE spacecraft, based on the Tudat software, is also presented, resulting in a promising opportunity to test PRIDE techniques and evaluate the improvements that PRIDE observables can make to natural bodies’ ephemerides. The first K _a-band (32 GHz) observations of the ESA’s BepiColombo by a radio telescope in the VLBI network, which employs a similar radio communications system as JUICE, were also demonstrated as a test case. The primary objective of these activities is to serve as a practice run for the upcoming operational PRIDE JUICE operations. We showcase the capabilities of the planning and scheduling software for other space missions.

Planetary Radio Interferometry and Doppler Experiment (PRIDE) is a multi-purpose experimental technique aimed at enhancing the science return of planetary missions. The technique exploits the science payload and spacecraft service systems without requiring a dedicated onboard instrumentation or imposing on the existing instrumentation any special for PRIDE requirements. PRIDE is based on the near-field phase-referencing Very Long Baseline Interferometry (VLBI) and evaluation of the Doppler shift of the radio signal transmitted by spacecraft by observing it with multiple Earth-based radio telescopes. The methodology of PRIDE has been developed initially at the Joint Institute for VLBI ERIC (JIVE) for tracking the ESA’s Huygens Probe during its descent in the atmosphere of Titan in 2005. From that point on, the technique has been demonstrated for various planetary and other space science missions. The estimates of lateral position of the target spacecraft are done using the phase-referencing VLBI technique. Together with radial Doppler estimates, these observables can be used for a variety of applications, including improving the knowledge of the spacecraft state vector. The PRIDE measurements can be applied to a broad scope of research fields including studies of atmospheres through the use of radio occultations, the improvement of planetary and satellite ephemerides, as well as gravity field parameters and other geodetic properties of interest, and estimations of interplanetary plasma properties. This paper presents the implementation of PRIDE as a component of the ESA’s Jupiter Icy Moons Explorer (JUICE) mission.

Long-term monitoring of global mass transport within the Earth system improves our ability to mitigate natural hazards and better understand their relations to climate change. Satellite gravity is widely used to monitor surface mass variations for its unprecedented spatial and temporal coverage. However, the gravity data contain signals from visco-elastic deformation in response to past ice sheet melting, preventing us from extracting signals of present-day surface mass trend (PDMT) directly. Here we present a global inversion scheme that separates PDMT and visco-elastic glacial isostatic adjustment (GIA) signatures by combining satellite gravimetry with satellite altimetry and ground observations. Our inversion provides global dual data coverage that enables a robust separation of PDMT and GIA spherical harmonic coefficients. It has the advantage of providing estimates of Earth's long wavelength deformation signatures and their uncertainties. Our GIA result, along with its uncertainty estimates, can be used in future GRACE processing to better assess the impact of GIA on surface mass change. Our GIA estimates include a rapid GIA uplift in the Southeast Alaska and the Amundsen Sea Embayment, due to the visco-elastic response to recent glacial unloading. We estimate the average surface mass change rate from 2002–2010 to be −203 ± 3 GT·a⁻¹ in Greenland, −126 ± 18 GT·a⁻¹ in Antarctica and, −62 ± 5 GT·a⁻¹ in Alaska. The GIA low degree spherical harmonic coefficients are sensitive to rheological properties in Earth's deep interior. Our low-degree GIA estimates include geocenter motion and (Formula presented.) which provide unique constraints to understand Earth's lower mantle and ice history.

In this study, we focus on improved constraint of the glacial isostatic adjustment (GIA) signal at present-day, and its role as a contributor to present-day sea level budgets. The main study area extends from the coastal regions of northwestern Europe to northern Europe. Both Holocene relative sea level (RSL) data as well as vertical land motion (VLM) data are incorporated as constraints in a semi-empirical GIA model. 71 geological rates of GIA-driven RSL change are inferred from Holocene proxy data and 108 rates of vertical land motion from GNSS provide an additional measure of regional GIA deformation. Within the study area, the geological RSL data complement the spatial gaps of the VLM data and vice versa. Both data sets are inverted in a semi-empirical GIA model to yield updated estimates of regional present-day GIA deformations. A regional validation using tide gauges is presented for the North Sea, where the GIA signal may be complicated by lateral variations in Earth structure and existing predictions of regional and global GIA models show discrepancies. The model validation in the North Sea region suggests that geological data are needed to fit independent estimates of GIA-related RSL change inferred from tide gauge rates, indicating that geological rates from Holocene data do provide an important additional constraint for data-driven approaches to GIA estimation.

Tidal heating can play an important role in the formation and evolution of subsurface oceans of outer-planet moons. Up until now tidal heating has only been studied in subsurface oceans of spatially uniform thickness. We develop a numerical model to consider oceans of spatially variable thickness. We use the Laplace Tidal Equations for the ocean and model the ice shell using membrane theory. The problem is solved using the commercial Finite Element software Comsol Multiphysics®. We use this new model to study the tidal response of Enceladus' ocean with a twofold objective: to understand how ocean thickness variations modify the tidal response of a subsurface ocean and to assess if tidal dissipation in an Enceladan ocean with varying ocean thickness can explain the high heat flux emanating from Enceladus' South Polar Terrain and the perdurance of a subsurface ocean. We consider the effect of meridional ocean thickness changes of spherical harmonic degree two and three as suggested by topography and gravimetry data. We observe that an ocean with degree two topography responds with the same eigenmodes as an ocean of constant thickness but resonances occur for thicker oceans. However, resonant ocean thicknesses are still thin compared to current estimates for Enceladus ocean thickness. Rossby-Haurwitz waves, excited by the obliquity tide for thick oceans of constant thickness, are not excited at the tidal frequency when oceans of variable thickness are considered. This result implies that the role of the obliquity tide in ocean tidal-dissipation might have been overestimated for Enceladus and other icy worlds. An antisymmetric, degree-three ocean thickness variation mixes the ocean modes excited in a constant thickness ocean by the eccentricity and obliquity tide.

Tidal dissipation makes Jupiter's moon Io the most volcanically active body in the solar system. Most of the heat generated in the interior is lost through volcanic activity. In this study, we aim to answer the questions: Can convection and melt migration in the mantle explain the spatial characteristics of Io's observed volcanic pattern? And, if so, what constraints does this place on the viscosity and thickness of the convective layer? We examine three different spatial characteristics of Io's volcanic activity: (i) The presence of global volcanism, (ii) the presence of large-scale variations in Io's volcanic activity, and (iii) the number of Io's volcanic systems. Our study relies on the assumptions that melt in the mantle controls Io's global volcanism, that the large-scale variations of Io's volcanic activity are caused by nonuniform tidal heating, and that the spatial density of volcanoes correlates with the spatial density of convective anomalies in the mantle. The results show that the observed small and large-scale characteristics of Io's volcanic pattern can be explained by sublithospheric anomalies influenced and caused by convective flow. Solutions that allow for active volcanism and Io's specific large-scale variations in volcanic activity range from a thick mantle of a high viscosity ((Formula presented.) Pa s) to a thin asthenosphere of a low viscosity ((Formula presented.) Pa s). Provided that Io's volcanoes are induced by convective anomalies in the mantle, we find that more than 80% of Io's internal heat is transported by magmatic processes and that Io's upper mantle needs to be thicker than 50 km.

Changes in ocean properties and circulation lead to a spatially non-uniform pattern of ocean dynamic sea-level change (DSLC). The projections of ocean dynamic sea level presented in the IPCC AR5 were constructed with global climate models (GCMs) from the Coupled Model Intercomparison Project 5 (CMIP5). Since CMIP5 GCMs have a relatively coarse resolution and exclude tides and surges it is unclear whether they are suitable for providing DSLC projections in shallow coastal regions such as the Northwestern European Shelf (NWES). One approach to addressing these shortcomings is dynamical downscaling – i.e. using a high-resolution regional model forced with output from GCMs. Here we use the regional shelf seas model AMM7 to show that, depending on the driving CMIP5 GCM, dynamical downscaling can have a large impact on DSLC simulations in the NWES region. For a business-as-usual greenhouse gas concentration scenario, we find that downscaled simulations of twenty-first century DSLC can be up to 15.5 cm smaller than DSLC in the GCM simulations along the North Sea coastline owing to unresolved processes in the GCM. Furthermore, dynamical downscaling affects the simulated time of emergence of sea-level change (SLC) above sea-level variability, and can result in differences in the projected change of the amplitude of the seasonal cycle of sea level of over 0.3 mm/yr. We find that the difference between GCM and downscaled results is of similar magnitude to the uncertainty of CMIP5 ensembles used for previous DSLC projections. Our results support a role for dynamical downscaling in future regional sea-level projections to aid coastal decision makers.

Satellite and recent Earth-based observations of Io's surface reveal a specific spatial pattern of persisting hotspots and sudden high-intensity events. Io's major heat producing mechanism is tidal dissipation, which is thought to be non-uniformly distributed within Io's mantle and asthenosphere. The question arises to what extent Io's non-homogeneous heat production can cause long-wavelength variations in the interior and volcanic activity at the surface. We investigate dissipation patterns resulting from two different initially spherical symmetric visco-elastic rheological structures, which are consistent with geodetic observations. The spatial distributions of the time-averaged tidal heat production are computed by a finite element model. Whereas for the first rheological structure heat is produced only in the upper viscous layer (asthenosphere-heating model), the second rheological structure results in a more evenly distributed dissipation pattern (mixed-heating model) with tidal heating occurring in the deep mantle and the asthenosphere. To relate the heat production to the interior temperature and melt distribution, we use steady-state scaling laws of mantle convection and a simple melt migration model. The resulting long-wavelength thermal heterogeneities strongly depend on the initial tidal dissipation pattern, the thickness of the convective layer, the mantle viscosity, and the ratio between magmatic and convective heat transport. While for the asthenosphere-heating model a strong lateral temperature signal with up to 190 K peak-to-peak difference can remain, convection within a thick convective layer, as for the mixed-heating model, can reduce the lateral temperature variation to <1 K, if the mantle viscosity is sufficiently low. Models with a dominating magma heat transport preserve the long-wavelength pattern of tidal dissipation much better and are favoured, because they are better to explain Io's thick crust. The approach presented here can also be applied to investigate the effect of an arbitrary interior heating pattern on Io's volcanic activity pattern.

A growing number of satellites in the outer solar system likely have global oceans beneath their outer icy shells. While the presence of liquid water makes these ocean worlds compelling astrobiological targets, the exchange of heat and materials between the deep interior and the surface also plays a critical role in promoting habitable environments. In this article, we combine geophysical, geochemical, and geological observations of the Jovian satellites Europa, Ganymede, and Callisto as well as the Saturnian satellites Enceladus and Titan to summarize our current state of understanding of their interiors and surface exchange processes. Potential mechanisms for driving exchange processes upward from the ocean floor and downward from the satellite surface are then reviewed, which are primarily based on numerical models of ice shell and ocean dynamics and complemented by terrestrial analog studies. Future missions to explore these exo-oceans will further revolutionize our understanding of ice-ocean exchange processes and their implications for the habitability of these worlds.

Some of the moons of the outer solar system harbour subsurface liquid oceans. Tidal dissipation plays an important role in preventing these oceans from freezing. In the past, most studies considered only tidal dissipation in the solid layers of these bodies (rock and ice). Recently, new studies considering tidal dissipation in the oceans of these moons have appeared. All of them make use of the shallow water approximation. However, the use of this approximation might not be adequate. Here we consider the linear non-hydrostatic three dimensional response of these oceans to tidal forcing with the full Coriolis force. To do so we consider an ocean of homogeneous density contained within a perfectly spherical shell and neglect the effect of the ice shell. We force the ocean with a time changing tidal potential and observe patterns of periodic inertial waves that take energy from the global tidal forcing and focus it along thin shear layers that propagate in the fluid. We focus on Europa and Enceladus, showing that inertial waves result in fluid flows of significant amplitude (a few cm/s). Nevertheless, we find that under the previously mentioned assumptions tidal dissipation due to inertial waves is several orders of magnitude smaller than Europa's radiogenic heating and Enceladus’ observed heat flux. Finally, we propose additional dissipation mechanisms that might play a relevant role in Europa and Enceladus and could be further investigated.

Existing lineaments on the surface of the Jovian moon Europa are thought to be the result of ongoing brittle crack formation in the elastic regime. Arcuate features are called cycloids and can be modeled using linear elastic fracture mechanics. Here, we build on existing terrestrial models of rift propagation and extend them to cycloids on the moon. We propose that these cracks tend to grow as a series of nearly instantaneous events, spaced by periods of inactivity. The behavior is similar to what is observed on Antarctic ice shelves, where rifts can remain dormant for years. We argue that dormant periods between growth events could explain the presence of cycloids on Europa even without invoking secular motion of the crust. Furthermore, being able to model propagation events and their timing should help future missions exploring the moon.

Rising sea levels due to climate change can have severe consequences for coastal populations and ecosystems all around the world. Understanding and projecting sea-level rise is especially important for low-lying countries such as the Netherlands. It is of specific interest for vulnerable ecological and morphodynamic regions, such as the Wadden Sea UNESCO World Heritage region. Here we provide an overview of sea-level projections for the 21st century for the Wadden Sea region and a condensed review of the scientific data, understanding and uncertainties underpinning the projections. The sea-level projections are formulated in the framework of the geological history of the Wadden Sea region and are based on the regional sea-level projections published in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5). These IPCC AR5 projections are compared against updates derived from more recent literature and evaluated for the Wadden Sea region. The projections are further put into perspective by including interannual variability based on long-Term tide-gauge records from observing stations at Den Helder and Delfzijl. We consider three climate scenarios, following the Representative Concentration Pathways (RCPs), as defined in IPCC AR5: The RCP2.6 scenario assumes that greenhouse gas (GHG) emissions decline after 2020; the RCP4.5 scenario assumes that GHG emissions peak at 2040 and decline thereafter; and the RCP8.5 scenario represents a continued rise of GHG emissions throughout the 21st century. For RCP8.5, we also evaluate several scenarios from recent literature where the mass loss in Antarctica accelerates at rates exceeding those presented in IPCC AR5. For the Dutch Wadden Sea, the IPCC AR5-based projected sea-level rise is 0.07±0.06m for the RCP4.5 scenario for the period 2018-30 (uncertainties representing 5-95%), with the RCP2.6 and RCP8.5 scenarios projecting 0.01m less and more, respectively. The projected rates of sea-level change in 2030 range between 2.6mma^-1 for the 5th percentile of the RCP2.6 scenario to 9.1mma^-1 for the 95th percentile of the RCP8.5 scenario. For the period 2018-50, the differences between the scenarios increase, with projected changes of 0.16±0.12m for RCP2.6, 0.19±0.11m for RCP4.5 and 0.23±0.12m for RCP8.5. The accompanying rates of change range between 2.3 and 12.4mma^-1 in 2050. The differences between the scenarios amplify for the 2018-2100 period, with projected total changes of 0.41±0.25m for RCP2.6, 0.52±0.27m for RCP4.5 and 0.76±0.36m for RCP8.5. The projections for the RCP8.5 scenario are larger than the high-end projections presented in the 2008 Delta Commission Report (0.74m for 1990-2100) when the differences in time period are considered. The sea-level change rates range from 2.2 to 18.3mma^-1 for the year 2100. We also assess the effect of accelerated ice mass loss on the sea-level projections under the RCP8.5 scenario, as recent literature suggests that there may be a larger contribution from Antarctica than presented in IPCC AR5 (potentially exceeding 1m in 2100). Changes in episodic extreme events, such as storm surges, and periodic (tidal) contributions on (sub-)daily timescales, have not been included in these sea-level projections. However, the potential impacts of these processes on sea-level change rates have been assessed in the report.

Sea-level indicators dated to the Last Interglacial, or Marine Isotope Stage (MIS) 5e, have a twofold value. First, they can be used to constrain the melting of Greenland and Antarctic Ice Sheets in response to global warming scenarios. Second, they can be used to calculate the vertical crustal rates at active margins. For both applications, the contribution of glacio- and hydro-isostatic adjustment (GIA) to vertical displacement of sea-level indicators must be calculated. In this paper, we re-assess MIS 5e sea-level indicators at 11 Mediterranean sites that have been generally considered tectonically stable or affected by mild tectonics. These are found within a range of elevations of 2–10 m above modern mean sea level. Four sites are characterized by two separate sea-level stands, which suggest a two-step sea-level highstand during MIS 5e. Comparing field data with numerical modeling we show that (i) GIA is an important contributor to the spatial and temporal variability of the sea-level highstand during MIS 5e, (ii) the isostatic imbalance from the melting of the MIS 6 ice sheet can produce a >2.0 m sea-level highstand, and (iii) a two-step melting phase for the Greenland and Antarctic Ice Sheets reduces the differences between observations and predictions. Our results show that assumptions of tectonic stability on the basis of the MIS 5e records carry intrinsically large uncertainties, stemming either from uncertainties in field data and GIA models. The latter are propagated to either Holocene or Pleistocene sea-level reconstructions if tectonic rates are considered linear through time.