Search

Detailed imaging of accretionary wedges reveals splay fault networks that could pose a significant tsunami hazard. However, the dynamics of multiple splay fault activation during megathrust earthquakes and the consequent effects on tsunami generation are not well understood. We use a 2‐D dynamic rupture model with complex topo‐bathymetry and six curved splay fault geometries constrained from realistic tectonic loading modeled by a geodynamic seismic cycle model with consistent initial stress and strength conditions. We find that all splay faults rupture coseismically. While the largest splay fault slips due to a complex rupture branching process from the megathrust, all other splay faults are activated either top down or bottom up by dynamic stress transfer induced by trapped seismic waves. We ascribe these differences to local non‐optimal fault orientations and variable along‐dip strength excess. Generally, rupture on splay faults is facilitated by their favorable stress orientations and low strength excess as a result of high pore‐fluid pressures. The ensuing tsunami modeled with non‐linear 1‐D shallow water equations consists of one high‐amplitude crest related to rupture on the longest splay fault and a second broader wave packet resulting from slip on the other faults. This results in two episodes of flooding and a larger run‐up distance than the single long‐wavelength (300 km) tsunami sourced by the megathrust‐only rupture. Since splay fault activation is determined by both variable stress and strength conditions and dynamic activation, considering both tectonic and earthquake processes is relevant for understanding tsunamigenesis., Plain Language Summary: In subduction zones, where one tectonic plate moves beneath another, earthquakes can occur on many different faults. Splay faults are relatively steep faults that branch off the largest fault (the megathrust) in a subduction zone. As they are steeper than the megathrust, the same amount of movement on them could result in more vertical displacement of the seafloor. Therefore, splay faults are thought to play an important role in the generation of tsunamis. Here, we use computer simulations to study if an earthquake can break multiple splay faults at once and how this affects the resulting tsunami. We find that multiple splay faults can indeed fail during a single earthquake due to the stress changes from trapped seismic waves, which promote rupture on splay faults. Rupture on splay faults results in larger seafloor displacements with smaller wavelengths, so the ensuing tsunami is bigger and results in two main flooding episodes at the coast. Our results show that it is important to consider rupture on splay faults when assessing tsunami hazard., Key Points: Multiple splay faults can be activated during a single earthquake by megathrust slip and dynamic stress transfer due to trapped waves. Splay fault activation is facilitated by their favorable orientation with respect to the local stress field and their closeness to failure. Long‐term geodynamic stresses and fault geometries affect dynamic splay fault rupture and the subsequent tsunami., Volkswagen Foundation (VolkswagenStiftung) http://dx.doi.org/10.13039/501100001663, Royal Society (The Royal Society) http://dx.doi.org/10.13039/501100000288, EC | H2020 | H2020 Priority Excellent Science | H2020 European Research Council (ERC) http://dx.doi.org/10.13039/100010663, Deutsche Forschungsgemeinschaft (DFG) http://dx.doi.org/10.13039/501100001659, National Science Foundation (NSF) http://dx.doi.org/10.13039/100000001, https://github.com/TUM-I5/SWE, https://doi.org/10.5281/zenodo.6969455

Recurrent and large forest fires negatively impact ecosystem, air quality, and human health. Moderate Resolution Imaging Spectroradiometer fire product is used to identify forest fires over central India domain, an extremely fire prone region. The study finds that from 2001 to 2020, ∼70% of yearly forest fires over the region occurred during March (1,857.5 counts/month) and April (922.8 counts/month). Some years such as 2009, 2012, and 2017 show anomalously high forest fires. The role of persistent warmer temperatures and multiple climate extremes in increasing forest fire activity over central India is comprehensively investigated. Warmer period from 2006 to 2020 showed doubling and tripling of forest fire activity during forest fire (February–June; FMAMJ) and non‐fire (July–January; JASONDJ) seasons, respectively. From 2015 JASONDJ to 2018 FMAMJ, central India experienced a severe heatwave, a rare drought and an extremely strong El Niño, the combined effect of which is linked to increased forest fires. Further, the study assesses quinquennial spatiotemporal changes in forest fire characteristics such as fire count density and average fire intensity. Deciduous forests of Jagdalpur‐Gadchiroli Range and Indravati National Park in Chhattisgarh state are particularly fire prone (>61 fire counts/grid) during FMAMJ and many forest fires are of high intensity (>45 MW). Statistical associations link high near surface air temperature and low precipitation during FMAMJ to significantly high soil temperature, low soil moisture content, low evapotranspiration and low normalized difference vegetation index. This creates a significantly drier environment, conducive for high forest fire activity in the region., Key Points Compared to 2001–2005, central India domain forest fire activity during 2006–2020 doubled in forest fire season and tripled in non‐fire seasonRole of persistent warmer temperatures and multiple climate extremes in increasing forest fire activity over central India is highlightedDeciduous forests of Jagdalpur‐Gadchiroli Range and Indravati National Park are extremely fire prone and fires are of high intensity

Along ultraslow spreading ridges melt is distributed unequally, but melt focusing guides melt away from amagmatic segments toward volcanic centers. An interplay of tectonism and magmatism is thought to control melt ascent, but the detailed process of melt extraction is not yet understood. We present a detailed image of the seismic velocity structure of the Logachev volcanic center and adjacent region along the Knipovich Ridge. With travel times of P‐ and S‐waves of 3,959 earthquakes we performed a local earthquake tomography. We simultaneously inverted for source locations, velocity structure and the Vp/Vs‐ratio. An extensive low velocity anomaly coincident with high Vp/Vs‐ratios >1.9 lies underneath the volcanic center at depths of 10 km below sea level in an aseismic area. More shallow, tightly clustered earthquake swarms connect the anomaly to a shallow anomaly with high Vp/Vs‐ratio beneath the basaltic seafloor. We consider the deep low‐velocity anomaly to represent an area of partial melt from which melts ascent vertically to the surface and northwards into the adjacent segment. By comparing tomographic studies of the Logachev and Southwest Indian Ridge Segment‐8 volcano we conclude that volcanic centers of ultraslow spreading ridges host spatially confined, circular partial melt areas below 10 km depth, in contrast to the shallow extended melt lenses along fast spreading ridges. Lateral feeding over distances of 35 km is possible at orthogonal spreading segments, but limited at the obliquely spreading Knipovich Ridge., Plain Language Summary: Mid‐ocean ridges mark the tectonic plate boundaries, where the plates drift apart. Fresh magma rises into the gap and builds new seafloor. The slower the plates drift apart, the less magma is present underneath the ridge. At very slow spreading ridges there is not enough magma to build new seafloor along the entire length of the ridge. Rather, melt is guided toward individual volcanic centers spaced at about 100 km, where melt accumulates and ascents. In our study we try to find melt storage areas and ascent paths of such a volcanic center. With velocities of different seismic wave types from earthquakes we map the velocity structure of the area underneath the major Logachev volcanic center. Lower velocities indicate an area partly including melt at depths of more than 10 km, far deeper than at mid‐ocean ridges with sufficient melt supply. From the deep magma reservoir, many earthquake swarms map the long ascent path of melt to the surface. The interplay of magmatic and tectonic activity is important here. In a comparison with results from another volcanic center, we find that lateral magma feeding is possible in orthogonal spreading, but limited in oblique spreading, as at the Knipovich Ridge., Key Points: Active volcanic centers at ultraslow spreading ridges host deeper and more confined partial melt areas than faster spreading ridges. Earthquake swarms delineate melt ascent paths from the partial melt area to the surface. Lateral feeding at shallow depths into subordinate segments is prevented by ridge obliquity., Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Following extreme drought during the 2019–2020 bushfire summer, the eastern part of Australia suffered from a week‐long intense rainfall and extensive flooding in March 2021. Understanding how much water storage changes in response to these climate extremes is critical for developing timely water management strategies. To quantify prompt water storage changes associated with the 2021 March flooding, we processed the low‐latency (1–3 days), high‐precision intersatellite laser ranging measurements from GRACE Follow‐On spacecraft and determined instantaneous gravity changes along spacecraft orbital passes. Such new data processing detected an abrupt surge of water storage approaching 60–70 trillion liters (km3 of water) over a week in the region, which concurrently caused land subsidence of ∼5 mm measured by a network of ground GPS stations. This was the highest speed of ground water recharge ever recorded in the region over the last two decades. Compared to the condition in February 2020, the amount of recharged water was similar but the recharge speed was much faster in March 2021. While these two events together replenished the region up to ∼80% of the maximum storage over the last two decades, the wet antecedent condition of soils in 2021 was distinctly different from the dry conditions in 2020 and led to generating extensive runoff and flooding in 2021., Key Points New GRACE Follow‐On gravity data processing quantified prompt changes in water storage by the heavy rainfall and flooding in March 2021The eastern Australia experienced the highest speed of ground water recharge and wet antecedent soils were responsible for extensive flooding in 2021The study demonstrated the feasibility of rapid processing (1–3 days) for immediate assessment of mass changes from extreme events

We present a data‐driven approach to clustering or grouping Global Navigation Satellite System (GNSS) stations according to observed velocities, displacements or other selected characteristics. Clustering GNSS stations provides useful scientific information, and is a necessary initial step in other analysis, such as detecting aseismic transient signals (Granat et al., 2013, https://doi.org/10.1785/0220130039). Desired features of the data can be selected for clustering, including some subset of displacement or velocity components, uncertainty estimates, station location, and other relevant information. Based on those selections, the clustering procedure autonomously groups the GNSS stations according to a selected clustering method. We have implemented this approach as a Python application, allowing us to draw upon the full range of open source clustering methods available in Python's scikit‐learn package (Pedregosa et al., 2011, https://doi.org/10.5555/1953048.2078195). The application returns the stations labeled by group as a table and color coded KML file and is designed to work with the GNSS information available from GeoGateway (Donnellan et al., 2021, https://doi.org/10.1007/s12145-020-00561-7; Heflin et al., 2020, https://doi.org/10.1029/2019ea000644) but is easily extensible. We demonstrate the methodology on California and western Nevada. The results show partitions that follow faults or geologic boundaries, including for recent large earthquakes and post‐seismic motion. The San Andreas fault system is most prominent, reflecting Pacific‐North American plate boundary motion. Deformation reflected as class boundaries is distributed north and south of the central California creeping section. For most models a cluster boundary connects the southernmost San Andreas fault with the Eastern California Shear Zone (ECSZ) rather than continuing through the San Gorgonio Pass., Key Points Unsupervised clustering methods provide a data‐driven way of analyzing and partitioning Global Navigation Satellite System observations of crustal deformationDeformation is distributed across the San Andreas fault system but is localized at the creeping section in central CaliforniaThe Southern San Andreas fault connects with the Eastern California Shear Zone rather than the rest of the San Andreas fault system

The effects of peat burning on organic‐rich agricultural soils of the Mezzano Lowland (NE Italy) were evaluated on soil profiles variously affected by smoldering. Profiles were investigated for pH, electrical conductivity, bulk density, elemental and isotopic composition of distinct carbon (and nitrogen) fractions. The results suggest that the horizons affected by carbon loss lie at depths 10–70 cm, where the highest temperatures are developed. We suggest that the exothermal oxidation of methane (mediated by biological activity) plays a significant role in the triggering mechanism. In the interested soils we estimated a potential loss of Soil Organic Carbon of approximately 110 kg m −2 within the first meter, corresponding to 580 kg CO2 m −3. The released greenhouse gas is coupled with a loss of soil structure and nutrients. Moreover, the process plausibly triggers mobility of metals bound in organometallic complexes. All these consequences negatively affect the environment, the agricultural activities and possibly also health of the local people., Key Points Peat burning in the Mezzano Lowland is a process active since wetland reclamation and variously affects the soil structure and carbon contentThe main carbon loss occurs in the depth interval 10–70 cm, where temperature up to 750°C is recordedThe process induces negative effects to surrounding environment, agriculture activities and possibly also human health

Interferometric synthetic aperture radar (InSAR) has been used to quantify a range of surface and near surface physical properties in permafrost landscapes. Most previous InSAR studies have utilized spaceborne InSAR platforms, but InSAR datasets over permafrost landscapes collected from airborne platforms have been steadily growing in recent years. Most existing algorithms dedicated toward retrieval of permafrost physical properties were originally developed for spaceborne InSAR platforms. In this study, which is the first in a two part series, we introduce a series of calibration techniques developed to apply a novel joint retrieval algorithm for permafrost active layer thickness retrieval to an airborne InSAR dataset acquired in 2017 by NASA's Uninhabited Aerial Vehicle Synthetic Aperture Radar over Alaska and Western Canada. We demonstrate how InSAR measurement uncertainties are mitigated by these calibration methods and quantify remaining measurement uncertainties with a novel method of modeling interferometric phase uncertainty using a Gaussian mixture model. Finally, we discuss the impact of native SAR resolution on InSAR measurements, the limitation of using few interferograms per retrieval, and the implications of our findings for cross‐comparison of airborne and spaceborne InSAR datasets acquired over Arctic regions underlain by permafrost., Key Points We develop and present several calibration and postprocessing methods for seasonal subsidence estimation from interferometric synthetic aperture radar deformationNovel methods for phase referencing and uncertainty quantification due to nonergodicity within the multilook window are proposedResidual sources of uncertainty in active layer thickness estimation are discussed and quantified

Seismological observations highlighted that earthquakes are often followed by changes in elastic properties around the fault zone. Here, we studied the origin of these variations using stick‐slip experiments on saw‐cut granite samples presenting different degrees of bulk damage (i.e., microcracks). Stick‐slip events were induced under triaxial compression configuration with continuous active ultrasonic measurements at confining pressures representative of upper crustal conditions (15–120 MPa). Both the P‐wave velocity (VP) and amplitude (AP) showed drops, concurrently with stress drops, and had a non‐monotonic dependence toward the fault's stress state. Our experimental results suggest that co‐seismic changes in VP were mostly controlled by the elastic re‐opening of microcracks in the bulk, rather than by co‐seismic damage or the formation of fault gouge. Co‐seismic changes in AP were controlled by a combination of elastic re‐opening of microcracks in the bulk and inelastic processes (i.e., co‐seismic damage and gouge formation and dilation)., Key Points Laboratory earthquakes are associated with seismic velocity and amplitude drops in crustal rocksVelocity drops are controlled by the stress acting on the fault and degree of damage in the fault wallAmplitude drops are partly controlled by the state of stress and partly affected by dissipative (inelastic) phenomena occurring on‐fault

The understanding of earthquake physics is hindered by the poor knowledge of fault strength and temperature evolution during seismic slip. Experiments reproducing seismic velocity (∼1 m/s) allow us to measure both the evolution of fault strength and the associated temperature increase due to frictional heating. However, temperature measurements were performed with techniques having insufficient spatial and temporal resolution. Here we conduct high velocity friction experiments on Carrara marble rock samples sheared at 20 MPa normal stress, velocity of 0.3 and 6 m/s, and 20 m of total displacement. We measured the temperature evolution of the fault surface at the acquisition rate of 1 kHz and over a spatial resolution of ∼40 µm with an optical fiber conveying the infrared radiation to a two‐color pyrometer. Temperatures up to 1,250°C and low coseismic fault shear strength are compatible with the activation of grain size dependent viscous creep., Key Points Optical fibers were used to measure temperature evolution during simulated seismic slipTemperatures were up to 1,250°C and compatible with coseismic grain size dependent creep“Frictional work” during simulated seismic slip was mostly dissipated as heat

Improved understanding of how our coasts will evolve over a range of time scales (years‐decades) is critical for effective and sustainable management of coastal infrastructure. A robust knowledge of the spatial, directional and temporal variability of the inshore wave climate is required to predict future coastal evolution and hence vulnerability. However, the variability of the inshore directional wave climate has received little attention, and an improved understanding could drive development of skillful seasonal or decadal forecasts of coastal response. We examine inshore wave climate at 63 locations throughout the United Kingdom and Ireland (1980–2017) and show that 73% are directionally bimodal. We find that winter‐averaged expressions of six leading atmospheric indices are strongly correlated (r = 0.60–0.87) with both total and directional winter wave power (peak spectral wave direction) at all studied sites. Regional inshore wave climate classification through hierarchical cluster analysis and stepwise multi‐linear regression of directional wave correlations with atmospheric indices defined four spatially coherent regions. We show that combinations of indices have significant skill in predicting directional wave climates (R 2 = 0.45–0.8; p, Key Points Over 70% of inshore wave climates analyzed throughout the United Kingdom and Ireland were directionally bimodalCombinations of winter atmospheric indices NAO, WEPA, SCAND, and EA are significantly correlated with directional wave climate in all regionsRegression models using multiple winter atmospheric indices enable skillful reconstructions of directional wave climate in all regions

Earthquakes preceding large events are commonly referred to as foreshocks. They are often considered as precursory phenomena reflecting the nucleation process of the main rupture. Such foreshock sequences may also be explained by cascades of triggered events. Recent advances in earthquake detection motivates a reevaluation of seismicity variations prior to mainshocks. Based on a highly complete earthquake catalog, previous studies suggested that mainshocks in Southern California are often preceded by anomalously elevated seismicity. In this study, we test the same catalog against the Epidemic Type Aftershock Sequence model that accounts for temporal clustering due to earthquake interactions. We find that 10/53 mainshocks are preceded by a significantly elevated seismic activity compared with our model. This shows that anomalous foreshock activity is relatively uncommon when tested against a model of earthquake interactions. Accounting for the recurrence of anomalies over time, only 3/10 mainshocks present a mainshock‐specific anomaly with a high predictive power., Key Points We further investigate previous claims of significantly elevated seismic activity prior to large earthquakes in Southern California10 out of 53 mainshocks are preceded by anomalously high seismicity, but only 3 of these anomalies are exclusively related to the mainshockThese selected foreshock sequences are likely due to additional pre‐slip, aseismic processes

The decline in global emissions of carbon dioxide due to the COVID‐19 pandemic provides a unique opportunity to investigate the sensitivity of the global carbon cycle and climate system to emissions reductions. Recent efforts to study the response to these emissions declines has not addressed their impact on the ocean, yet ocean carbon absorption is particularly susceptible to changing atmospheric carbon concentrations. Here, we use ensembles of simulations conducted with an Earth system model to explore the potential detection of COVID‐related emissions reductions in the partial pressure difference in carbon dioxide between the surface ocean and overlying atmosphere (ΔpCO2), a quantity that is regularly measured. We find a unique fingerprint in global‐scale ΔpCO2 that is attributable to COVID, though the fingerprint is difficult to detect in individual model realizations unless we force the model with a scenario that has four times the observed emissions reduction., Key Points COVID‐related emissions reductions will be imperceptible in surface ocean pH observationsThe CanESM5 COVID ensemble predicts a unique fingerprint of COVID‐related emissions reductions in global mean ΔpCO2 (pCO2oc ‐ pCO2atm)The fingerprint is potentially detectable in global‐scale observations of ΔpCO2, but only with large emissions reductions

We use Ice, Cloud, and land Elevation Satellite 2 (ICESat‐2) laser altimetry crossovers and repeat tracks collected over the North Slope of Alaska to estimate ground surface‐height change due to the seasonal freezing and thawing of the active layer. We compare these measurements to a time series of surface deformation from Sentinel‐1 interferometric synthetic aperture radar (InSAR) and demonstrate agreement between these independent observations of surface deformation at broad spatial scales. We observe a relationship between ICESat‐2‐derived surface subsidence/uplift and changes in normalized accumulated degree days, which is consistent with the thermodynamically driven seasonal freezing and thawing of the active layer. Integrating ICESat‐2 crossover estimates of surface‐height change yields an annual time series of surface‐height change that is sensitive to changes in snow cover during spring and thawing of the active layer throughout spring and summer. Furthermore, this time series exhibits temporal correlation with independent reanalysis datasets of temperature and snow cover, as well as an InSAR‐derived time series. ICESat‐2‐derived surface‐height change estimates can be significantly affected by short length‐scale topographic gradients and changes in snow cover and snow depth. We discuss optimal strategies of post‐processing ICESat‐2 data for permafrost applications, as well as the future potential of joint ICESat‐2 and InSAR investigations of permafrost surface‐dynamics., Key Points ICESat‐2 altimetry can resolve surface subsidence that is related to changes in snow‐cover depth and seasonal thawing of the active layerICESat‐2 measurements of surface‐height change are affected by along‐track topographic gradients and complex surface roughnessComplementary ICESat‐2 and InSAR datasets can be jointly leveraged for future studies in periglacial environments

Laboratory studies suggest that seismogenic rupture on faults in carbonate terrains can be explained by a transition from high friction, at low sliding velocities (V), to low friction due to rapid dynamic weakening as seismic slip velocities are approached. However, consensus on the controlling physical processes is lacking. We previously proposed a microphysically based model (the “Chen–Niemeijer–Spiers” [CNS] model) that accounts for the (rate‐and‐state) frictional behavior of carbonate fault gouges seen at low velocities characteristic of rupture nucleation. In the present study, we extend the CNS model to high velocities (1 mm/s ≤ V ≤ 10 m/s) by introducing multiple grain‐scale deformation mechanisms activated by frictional heating. As velocity and hence temperature increase, the model predicts a continuous transition in dominant deformation mechanisms, from frictional granular flow with partial accommodation by plasticity at low velocities and temperatures, to grain boundary sliding with increasing accommodation by solid‐state diffusion at high velocities and temperatures. Assuming that slip occurs in a localized shear band, within which grain size decreases with increasing velocity, the model results capture the main mechanical trends seen in high‐velocity friction experiments on room‐dry calcite‐rich rocks, including steady‐state and transient aspects, with reasonable quantitative agreement and without the need to invoke thermal decomposition or fluid pressurization effects. The extended CNS model covers the full spectrum of slip velocities from earthquake nucleation to seismic slip rates. Since it is based on realistic fault structure, measurable microstructural state variables, and established deformation mechanisms, it may offer an improved basis for extrapolating lab‐derived friction data to natural fault conditions., Key Points We extend a microphysical friction model from earthquake nucleation to seismic slip velocitiesThe extended model predicts steady‐state and transient frictionModel results reproduce high‐velocity friction experiments on simulated calcite fault gouges

Transient deformation associated with foreshocks activity has been observed before large earthquakes, suggesting the occurrence of a detectable preseismic slow slip during the initiation phase. A critical issue consists in discriminating the relative contributions from seismic and aseismic fault slip during the preparation phase of large earthquakes. We focus on the April–May 2017 Valparaíso earthquake sequence, which involved a M W = 6.9 earthquake preceded by intense foreshock activity. To assess the relative contribution of seismic and aseismic slip, we compare surface displacement predicted from foreshocks source models with transient motion measured prior to the mainshock. The comparison between observed and predicted displacements shows that only half of the total displacement can be explained by the contribution of foreshocks. This result suggests the presence of aseismic pre‐slip during an initiation phase preceding the mainshock., Key Points The 2017 Valparaíso M W = 6.9 earthquake presents a preseismic transient displacementWe evaluate the contribution of foreshock‐induced displacement to the preseismic Global Positioning System observationsResults suggest that 50% ± 11% of the preseismic displacement is caused by aseismic slip

A reduced‐physics model is employed at 1/25° to 1/100° global resolution to determine (a) if linear dynamics can reproduce the observed low‐mode M2 internal tide, (b) internal‐tide sensitivity to bathymetry, stratification, surface tides, and dissipation parameterizations, and (c) the amount of power transferred to the nonstationary internal tide. The simulations predict 200 GW of mode‐1 internal‐tide generation, consistent with a general circulation model and semianalytical theory. Mode‐1 energy is sensitive to damping, but a simulation using parameterizations for wave drag and wave‐mean interaction predicts 84% of satellite observed sea‐surface height amplitude variance on a 1° × 1° grid. The simulation energy balance indicates that 16% of stationary mode‐1 energy is scattered to modes 2–4 and negligible energy propagates onto the shelves. The remaining 84% of energy is lost through parameterizations for high‐mode scattering over rough topography (54%) and wave‐mean interactions that transfer energy to the nonstationary internal tide (29%)., Key Points Global mode‐1 internal tides are predictable from a linearized, reduced physics modelThe choice of damping parameterizations predetermines mode‐1 internal‐tide energySimulations with damping by wave drag and wave‐mean interactions largely agree with satellite data

The conventional sea level budget (SLB) equates changes in sea surface height with the sum of ocean mass and steric change, where solid‐Earth movements are included as corrections but limited to the impact of glacial isostatic adjustment. However, changes in ocean mass load also deform the ocean bottom elastically. Until the early 2000s, ocean mass change was relatively small, translating into negligible elastic ocean bottom deformation (OBD), hence neglected in the SLB equation. However, recently ocean mass has increased rapidly; hence, OBD is no longer negligible and likely of similar magnitude to the deep steric sea level contribution. Here, we use a mass‐volume framework, which allows the ocean bottom to respond to mass load, to derive a SLB equation that includes OBD. We discuss the theoretical appearance of OBD in the SLB equation and its implications for the global SLB., Key Points The conventional sea level budget equation does not include elastic ocean bottom deformation, implicitly assuming it is negligibleRecent increases in ocean mass yield global‐mean ocean bottom deformation of similar magnitude to the deep steric sea level contributionWe use a mass‐volume approach to derive and update the sea level budget equation to include ocean bottom deformation as a process

Gravity survey is one of the passive geophysical techniques commonly used to delineate geological formations, especially in determining basement rock and the overlying deposit. Geologically, the study area is made up of thick quaternary alluvium deposited on top of the older basement rock. The Muda River basin constitutes, approximately, of more than 300 m of thick quaternary alluvium overlying the unknown basement rock type. Previous studies, including drilling and geo-electrical resistivity surveys, were conducted in the area but none of them managed to conclusively determine the basement rock type and depth precisely. Hence, a regional gravity survey was conducted to determine the thickness of the quaternary sediments prior to assessing the sustainability of the Muda River basin. Gravity readings were made at 347 gravity stations spaced at 3–5 km intervals using Scintrex CG-3 covering an area and a perimeter of 9000 km2 and 730 km, respectively. The gravity data were then conventionally reduced for drift, free air, latitude, Bouguer, and terrain corrections. These data were then consequently analyzed to generate Bouguer, regional and total horizontal derivative (THD) anomaly maps for qualitative and quantitative interpretations. The Bouguer gravity anomaly map shows low gravity values in the north-eastern part of the study area interpreted as representing the Main Range granitic body, while relatively higher gravity values observed in the south-western part are interpreted as representing sedimentary rocks of Semanggol and Mahang formations. Patterns observed in the THD anomaly and Euler deconvolution maps closely resembled the presence of structural features such as fault lineaments dominantly trending along NW-SE and NE-SW like the trends of topographic lineaments in the study area. Based on power spectral analysis of the gravity data, the average depth of shallow body, representing alluvium, and deep body, representing underlying rock formations, are 0.5 km and 1.2 km, respectively. The thickness of Quaternary sediment and the depth of sedimentary formation can be more precisely estimated by other geophysical techniques such as the seismic reflection survey.

Since 1982, high precision gravity measurements have been routinely carried out on Mt. Vesuvius. The gravity network consists of selected sites most of them coinciding with, or very close to, leveling benchmarks to remove the effect of the elevation changes from gravity variations. The reference station is located in Napoli, outside the volcanic area. Since 1986, absolute gravity measurements have been periodically made on a station on Mt. Vesuvius, close to a permanent gravity station established in 1987, and at the reference in Napoli. The results of the gravity measurements since 1982 are presented and discussed. Moderate gravity changes on short-time were generally observed. On long-term significant gravity changes occurred and the overall fields displayed well defined patterns. Several periods of evolution may be recognized. Gravity changes revealed by the relative surveys have been confirmed by repeated absolute measurements, which also confirmed the long-term stability of the reference site. The gravity changes over the recognized periods appear correlated with the seismic crises and with changes of the tidal parameters obtained by continuous measurements. The absence of significant ground deformation implies masses redistribution, essentially density changes without significant volume changes, such as fluids migration at the depth of the seismic foci, i.e. at a few kilometers. The fluid migration may occur through pre-existing geological structures, as also suggested by hydrological studies, and/or through new fractures generated by seismic activity. This interpretation is supported by the analyses of the spatial gravity changes overlapping the most significant and recent seismic crises.

Changes in gravity and/or surface deformation are often associated with volcanic activity. Usually, bodies with simple geometry (e.g., point sources, prolate or oblate spheroids) are used to model these signals considering anomalous mass and/or pressure variations. We present a new method for the simultaneous, nonlinear inversion of gravity changes and surface deformation using bodies with a free geometry. Assuming simple homogenous elastic conditions, the method determines a general geometrical configuration of pressure and density sources. These sources are described as an aggregate of pressure and density point sources, fitting the whole data set (given some regularity conditions). The approach works in a growth step-by-step process that allows us to build very general geometrical configurations. The methodology is validated against an ellipsoidal body with anomalous pressure and a parallelepiped body with anomalous density, buried in an elastic medium. The simultaneous inversion of deformation and gravity values permits a good reconstruction of the assumed bodies. Finally, the inversion method is applied to the interpretation of gravity, leveling, and interferometric synthetic aperture radar (InSAR) data from the volcanic area of Campi Flegrei (Italy) for the period 1992–2000. For this period, a model with no significant mass change and an extended decreasing pressure region satisfactorily fits the data. The pressure source is located at about ∼1500 m depth, and it is interpreted as corresponding to the dynamics of the shallow (depth 1–2 km) hydrothermal system confined to the caldera fill materials.

EGU2020: Sharing Geoscience Online, 4-8 May 2020, The modes in which the lithosphere deforms during continental collision and the mechanisms involved are not well understood. While continental subduction and mantle delamination are often invoked in tectonophysical studies, these processes are difficult to be confirmed in more complex tectonic regions such as the Gibraltar Arc. We study the present-day density and compositional structure of the lithosphere along a transect running from S Iberia to N Africa crossing the western Gibraltar Arc. This region is located in the westernmost continental segment of the African-Eurasian plates, characterized by a diffuse transpressive plate boundary. An integrated and self-consistent geophysical-petrological methodology is used to model the lithosphere structure variations and the thermophysical properties of the upper mantle. The crustal structure is mainly constrained by seismic experiments and geological data, whereas the composition of the lithospheric mantle is constrained by xenolith data. The results show large lateral variations in the topography of the lithosphere-asthenosphere boundary (LAB). We distinguish different chemical lithospheric mantle domains that reproduce the main trends of the geophysical observables and the modelled P- and S-wave seismic velocities. A sublithospheric body colder than the surrounding mantle is needed beneath the Betics-Rif to adjust the measured potential fields. We link this body to the Iberian slab localized just to the east of the profile and having some effect on the geoid and Bouguer anomalies. Local isostasy allows explaining most of the topography, but an elastic thickness higher than 10 km is needed to explain local misfits between the Atlas and the Rif Mountains., This work has been supported by Spanish Ministry by the projects MITE (CGL2014-59516) and GeoCAM (PGC2018-095154-B-100).