The magmatic system under Hunga volcano before and after the 15 January 2022 eruption

One of the largest explosive eruptions instrumentally recorded occurred at Hunga volcano on 15 January 2022. The magma plumbing system under this volcano is unexplored because of inherent difficulties caused by its submarine setting. We use marine gravity data derived from satellite altimetry combined with multibeam bathymetry to model the architecture and dynamics of the magmatic system before and after the January 2022 eruption. We provide geophysical evidence for substantial high-melt content magma accumulation in three reservoirs at shallow depths (2 to 10 kilometers) under the volcano. We estimate that less than ~30% of the existing magma was evacuated by the main eruptive phases, enough to trigger caldera collapse. The eruption and caldera collapse reorganized magma storage, resulting in an increased connectivity between the two spatially distinct reservoirs. Modeling global satellite altimetry-derived gravity data at undersea volcanoes offer a promising reconnaissance tool to probe the subsurface for eruptible magma.

Satellite Altimetry: Achievements and Future Trends by a Scientometrics Analysis

Scientometric reviews, facilitated by computational and visual analytical approaches, allow researchers to gain a thorough understanding of research trends and areas of concentration from a large number of publications. With the fast development of satellite altimetry, which has been effectively applied to a wide range of research topics, it is timely to summarize the scientific achievements of the previous 50 years and identify future trends in this field. A comprehensive overview of satellite altimetry was presented using a total of 8541 publications from the Web of Science Core Collection covering the years from 1970 to 2021. We begin by presenting the fundamental statistical results of the publications, such as the annual number of papers, study categories, countries/regions, afflictions, journals, authors, and keywords, in order to provide a comprehensive picture of satellite altimetry research. We discuss the co-occurrence of the authors in order to reveal the global collaboration network of satellite altimetry research. Finally, we utilised co-citation networks to detect the development trend and associated crucial publications for various specific topics. The findings show that satellite altimetry research has been changed immensely during the last half-century. The United States, France, China, England, and Germany made the most significant contributions in the field of satellite altimetry. The analysis reveals a clear link between technology advancements and the trend in satellite altimetry research. As a result, wide swath altimetry, GNSS-reflectometry, laser altimetry, terrestrial hydrology, and deep learning are among the most frontier study subjects. The findings of this work could guide a thorough understanding of satellite altimetry’s overall development and research front.

Seafloor Topography Estimation from Gravity Anomaly and Vertical Gravity Gradient Using Nonlinear Iterative Least Square Method

Currently, seafloor topography inversion based on satellite altimetry gravity data provides the principal means to predict the global seafloor topography. Researchers often use sea surface geoid height or gravity anomaly to predict sea depth in the space domain. In this paper, a comprehensive discussion on seafloor topography inversion formulas in the space domain is presented using sea surface geoid height, gravity anomaly and introduces an approach that uses vertical gravity gradient. This would be the first study to estimate seafloor topography by vertical gravity gradient in the space domain. Further, a nonlinear iterative least-square inversion process is discussed. Using the search area for the Malaysia Airlines Flight MH370 as study site, we used the DTU17 gravity anomaly model and SIO V29.1 vertical gravity gradient to generate the seafloor topography. The results of the proposed bathymetric models were analyzed and compared with the DTU18 and SIO V20.1 bathymetric models. The experimental results show that the gravity anomaly and vertical gravity gradient in the study area are strongly correlated with the seafloor topography in the 20–200 km wavelength range. The optimal initial iteration values for seafloor topography variance and correlation length are 0.6365 km2 and 10.5′, respectively. Shipborne measurements from SONAR data were used as external checkpoints to evaluate the bathymetric models. The results show that the RMS for BAT_VGG_ILS (inversion model constructed by vertical gravity gradient) is smaller than for BAT_GA_ILS (inversion model constructed by gravity anomaly) and BAT_GA_VGG_ILS (inversion model constructed by gravity anomaly and vertical gravity gradient). The relative accuracy of the DTU18 bathymetry model was 9.27%, while the relative accuracy of the proposed seafloor models was higher than 4%. Within the 200 m difference range, the proportion of checkpoints for BAT_VGG_ILS was close to 95%, about 80% for BAT_GA_ILS and BAT_GA_VGG_ILS, and less than 50% for the DTU18. The results show that the nonlinear iterative least square method in the space domain is feasible.

Vertical Deflections and Gravity Disturbances Derived from HY-2A Data

The first Chinese altimetry satellite, Haiyang-2A (HY-2A), which was launched in 2011, has provided a large amount of sea surface heights which can be used to derive marine gravity field. This paper derived the vertical deflections and gravity disturbances using HY-2A observations for the major area of the whole Earth’s ocean from 60°S and 60°N. The results showed that the standard deviations (STD) of vertical deflections differences were 1.1 s and 3.5 s for the north component and the east component between HY-2A’s observations and those from EGM2008 and EIGEN-6C4, respectively. This indicates the accuracy of the east component was poorer than that of the north component. In order to clearly demonstrate contribution of HY-2A’s observations to gravity disturbances, reference models and the commonly used remove-restore method were not adopted in this study. Therefore, the results can be seen as ‘pure’ signals from HY-2A. Assuming the values from EGM2008 were the true values, the accuracy of the gravity disturbances was about −1.1 mGal in terms of mean value of the errors and 8.0 mGal in terms of the STD. This shows systematic errors if only HY-2A observations were used. An index of STD showed that the accuracy of HY-2A was close to the theoretical accuracy according to the vertical deflection products. To verify whether the systematic errors of gravity field were from the long wavelengths, the long-wavelength parts of HY-2A’s gravity disturbance with wavelengths larger than 500 km were replaced by those from EGM2008. By comparing with ‘pure’ HY-2A version of gravity disturbance, the accuracy of the new version products was improved largely. The systematic errors no longer existed and the error STD was reduced to 6.1 mGal.

Integrated geophysical investigation of recent earth tremors in Nigeria using aeromagnetic and gravity data

Understanding the tectonics at active sites of recent seismic activities in Nigeria is fundamental towards disaster mitigation and emergency planning. We apply geophysical techniques of gravity and magnetic methods to investigate crustal depths, subsurface geologic faults and fractures, and the overall subsurface interaction at Mpape region and environs. Estimated depths to the bottom of magnetic crust (basal depths) range between 11.0 and 11.4 km at the Mpape region and decrease further southward towards Guabe town. This signifies the depth range of the active crust within the region. Comparative deeper basal depths (15.0-16.2 km) were obtained at locations farther from Mpape-Guabe towns at Nasarawa, Rubochi, and Fuka regions, showing a more stable region away from Mpape region. Computed Moho depths from gravity data show deeper depths at the Mpape region (~ 34.1 km) suggesting that the active crust exists in the upper crust. Two-dimensional modeling analysis along a profile taken across the Mpape region shows a conspicuous subsurface basement intrusion at the Mpape region with deep faults and fractures reaching depths of 7-14 km. Shallow basal depths at the Mpape region resulting from significant subsurface intrusion and concentrated subsurface faults at the intruded region may be responsible for the instability of the Mpape region. The most affected area is located within the Mpape-Guabe towns. We recommend the establishment of seismic monitoring facilities in this area for effective monitoring and evaluation.

Effects of Interferometric Radar Altimeter Errors on Marine Gravity Field Inversion

The traditional altimetry satellite, which is based on pulse-limited radar altimeter, only measures ocean surface heights along tracks; hence, leads to poorer accuracy in the east component of the vertical deflections compared to the north component, which in turn limits the final accuracy of the marine gravity field inversion. Wide-swath altimetry using radar interferometry can measure ocean surface heights in two dimensions and, thus, can be used to compute vertical deflections in an arbitrary direction with the same accuracy. This paper aims to investigate the impact of Interferometric Radar Altimeter (InRA) errors on gravity field inversion. The error propagation between gravity anomalies and InRA measurements is analyzed, and formulas of their relationship are given. By giving a group of possible InRA parameters, numerical simulations are conducted to analyze the accuracy of gravity anomaly inversion. The results show that the accuracy of the gravity anomalies is mainly influenced by the phase errors of InRA; and the errors of gravity anomalies have a linear approximation relationship with the phase errors. The results also show that the east component of the vertical deflections has almost the same accuracy as the north component.

Inversion and Validation of Improved Marine Gravity Field Recovery in South China Sea by Incorporating HY-2A Altimeter Waveform Data

HaiYang-2A (HY-2A, where ‘Haiyang’ means ‘Ocean’ in Chinese) has provided reliable sea surface height observations for gravity with uniform ocean data coverage on a global scale for more than 8 years, particularly with denser across track sampling during the geodetic mission since March 2016. This paper aims at modeling and evaluating the regional marine gravity field at 1′×1′ resolution by incorporating HY-2A altimeter waveform data from 7 complete 168-day cycles in the geodetic mission phase. Initial evaluation indicates that, firstly, the measurements in the geodetic mission stay at a consistent accuracy level with observations at the start-of-life stage according to statistics of discrepancies at crossover points cycle by cycle. Secondly, range precision improvement can be achieved using a two-pass weighted least-squares retracker. Thirdly, a downsampling procedure combined with a low-pass filter is designed for HY-2A 20 Hz data to obtain 5 Hz measurements with enhanced precision. We calculate the 1′×1′ marine gravity field model over the South China Sea area by using the EGM2008 model as a reference field with the remove/restore method. The verifications with published models and shipborne gravimetric data show that HY-2A GM data is capable of improving marine gravity field modeling. Results show slightly higher accuracy than other models with similar input datasets but not including HY-2A. The accuracy is also compared with the latest DTU17 and SIO V27.1 model.

Geotectonic Controls on CO2 Formation and Distribution Processes in the Brazilian Pre-Salt Basins

Exploratory work for hydrocarbons along the southeastern Brazilian Margin discovered high concentrations of CO2 in several fields, setting scientific challenges to understand these accumulations. Despite significant progress in understanding the consequences of high CO2 in these reservoirs, the role of several variables that may control such accumulations of CO2 is still unclear. For example, significant differences in the percentages of CO2 have been found in reservoirs of otherwise similar prospects lying close to each other. In this paper, we present a hypothesis on how the rifting geodynamics are related to these CO2-rich accumulations. CO2-rich mantle material may be intruded into the upper crustal levels through hyper-stretched continental crust during rifting. Gravimetric and magnetic potential methods were used to identify major intrusive bodies, crustal thinning and other geotectonic elements of the southeastern Brazilian Margin. Modeling based on magnetic, gravity, and seismic data suggests a major intrusive magmatic body just below the reservoir where a high CO2 accumulation was found. Small faults connecting this magmatic body with the sedimentary section could be the fairway for the magmatic sourced gas rise to reservoirs. Mapping and understanding the crustal structure of sedimentary basins are shown to be important steps for “de-risking” the exploration process.

Evaluation of the Accuracy of Bathymetry on the Nearshore Coastlines of Western Korea from Satellite Altimetry, Multi-Beam, and Airborne Bathymetric LiDAR

Bathymetric mapping is traditionally implemented using shipborne single-beam, multi-beam, and side-scan sonar sensors. Procuring bathymetric data near coastlines using shipborne sensors is difficult, however, this type of data is important for maritime safety, marine territory management, climate change monitoring, and disaster preparedness. In recent years, the bathymetric light detection and ranging (LiDAR) technique has been tried to get seamless geospatial data from land to submarine topography. This paper evaluated the accuracy of bathymetry generated near coastlines from satellite altimetry-derived gravity anomalies and multi-beam bathymetry using a tuning density contrast of 5000 kg/m³ determined by the gravity-geologic method. Comparing with the predicted bathymetry of using only multi-beam depth data, 78% root mean square error from both multi-beam and airborne bathymetric LiDAR was improved in shallow waters of nearshore coastlines of the western Korea. As a result, the satellite-derived bathymetry estimated from the multi-beam and the airborne bathymetric LiDAR was enhanced to the accuracy of about 0.2 m.

Marine geophysics. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure

Gravity models are powerful tools for mapping tectonic structures, especially in the deep ocean basins where the topography remains unmapped by ships or is buried by thick sediment. We combined new radar altimeter measurements from satellites CryoSat-2 and Jason-1 with existing data to construct a global marine gravity model that is two times more accurate than previous models. We found an extinct spreading ridge in the Gulf of Mexico, a major propagating rift in the South Atlantic Ocean, abyssal hill fabric on slow-spreading ridges, and thousands of previously uncharted seamounts. These discoveries allow us to understand regional tectonic processes and highlight the importance of satellite-derived gravity models as one of the primary tools for the investigation of remote ocean basins.

Gravimetric structure for the abyssal mantle massif of Saint Peter and Saint Paul peridotite ridge, Equatorial Atlantic Ocean, and its relation to active uplift

This paper presents gravimetric and morphologic analyses based on the satellite-derived data set of EGM2008 and TOPEX for the area of the oceanic mantle massif of the Saint Peter and Saint Paul peridotite ridge, Equatorial Atlantic Ocean. The free-air anomaly indicates that the present plate boundary is not situated along the longitudinal graben which cuts peridotite ridge, but about 20 km to the north of it. The high Bouguer anomaly of the peridotite ridge suggests that it is constituted mainly by unserpentinised ultramafic rocks. The absence of isostatic compensation and low-degree serpentinisation of the ultramafic rocks indicate that the peridotite ridge is sustained mainly by active tectonic uplift. The unparallel relation between the transform fault and the relative plate motion generates near north-south compression and the consequent tectonic uplift. In this sense, the peridotite massif is a pressure ridge due to the strike-slip displacement of the Saint Paul Transform Fault.

Toward 1-mGal accuracy in global marine gravity from CryoSat-2, Envisat, and Jason-1

More than 60% of the Earth's land and shallow marine areas are covered by > 2 km of sediments and sedimentary rocks, with the thickest accumulations on rifted continental margins (Figure 1). Free-air marine gravity anomalies derived from Geosat and ERS-1 satellite altimetry (Fairhead et al., 2001; Sandwell and Smith, 2009; Andersen et al., 2009) outline most of these major basins with remarkable precision. Moreover, gravity and bathymetry data derived from altimetry are used to identify current and paleo-submarine canyons, faults, and local recent uplifts. These geomorphic features provide clues to where to look for large deposits of sediments. While current altimeter data delineate large offshore basins and major structures, they do not resolve some of the smaller geomorphic features and basins (Yale et al., 1998; Fairhead et al., 2001). Improved accuracy and resolution is desirable: to facilitate comparisons between continental margins; as an exploration tool and to permit extrapolation of known structures from well-surveyed areas; to follow fracture zones out of the deep-ocean basin into antecedent continental structures, to define and compare segmentation of margins along strike and identify the position of the continent-ocean boundary; and to study mass anomalies (e.g., sediment type and distribution) and isostatic compensation at continental margins. In this article, we assess the accuracy of a new global marine gravity model based on a wealth of new radar altimetry data and demonstrate that these gravity data are superior in quality to the majority of publicly available academic and government ship gravity data.

Australia–SE Asia collision: plate tectonics and crustal flow

The Sundaland core of SE Asia is a heterogeneous assemblage of Tethyan sutures and Gondwana fragments. Its complex basement structure was one major influence on Cenozoic tectonics; the rifting history of the north Australian margin was another. Fragments that rifted from Australia in the Jurassic collided with Sundaland in the Cretaceous and terminated subduction. From 90 to 45 Ma Sundaland was largely surrounded by inactive margins with localized strike-slip deformation, extension and subduction. At 45 Ma Australia began to move north, and subduction resumed beneath Sundaland. At 23 Ma the Sula Spur promontory collided with the Sundaland margin. From 15 Ma there was subduction hinge rollback into the Banda oceanic embayment, major extension, and later collision of the Banda volcanic arc with the southern margin of the embayment. However, this plate tectonic framework cannot be reduced to a microplate scale to explain Cenozoic deformation. Sundaland has a weak thin lithosphere, highly responsive to plate boundary forces and a hot weak deep crust has flowed in response to tectonic and topographic forces, and sedimentary loading. Gravity-driven movements of the upper crust, unusually rapid vertical motions, exceptionally high rates of erosion, and massive movements of sediment have characterized this region.

Contrasting décollement and prism properties over the Sumatra 2004-2005 earthquake rupture boundary

Styles of subduction zone deformation and earthquake rupture dynamics are strongly linked, jointly influencing hazard potential. Seismic reflection profiles across the trench west of Sumatra, Indonesia, show differences across the boundary between the major 2004 and 2005 plate interface earthquakes, which exhibited contrasting earthquake rupture and tsunami generation. In the southern part of the 2004 rupture, we interpret a negative-polarity sedimentary reflector approximately 500 meters above the subducting oceanic basement as the seaward extension of the plate interface. This predécollement reflector corresponds to unusual prism structure, morphology, and seismogenic behavior that are absent along the 2005 rupture zone. Although margins like the 2004 rupture zone are globally rare, our results suggest that sediment properties influence earthquake rupture, tsunami hazard, and prism development at subducting plate boundaries.