Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Atmospheres”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53Linked to original sources
USGS2020

A lacustrine paleoenvironment recorded at Vera Rubin ridge, Gale crater: Overview of the sedimentology and stratigraphy observed by the Mars Science Laboratory Curiosity rover

For ~500 Martian solar days (sols), the Mars Science Laboratory team explored Vera Rubin ridge (VRR), a topographic feature on the northwest slope of Aeolis Mons. Here we review the sedimentary facies and stratigraphy observed during sols 1,800–2,300, covering more than 100 m of stratigraphic thickness. Curiosity's traverse includes two transects across the ridge, which enables investigation of lateral variability over a distance of ~300 m. Three informally named stratigraphic members of the Murray formation are described: Blunts Point, Pettegrove Point, and Jura, with the latter two exposed on VRR. The Blunts Point member, exposed just below the ridge, is characterized by a recessive, fine-grained facies that exhibits extensive planar lamination and is crosscut by abundant curvi-planar veins. The Pettegrove Point member is more resistant, fine-grained, thinly planar laminated, and contains a higher abundance of diagenetic concretions. Conformable above the Pettegrove Point member is the Jura member, which is also fine-grained and parallel stratified, but is marked by a distinct step in topography, which coincides with localized meter-scale inclined strata, a thinly and thickly laminated facies, and occasional crystal molds. All members record low-energy lacustrine deposition, consistent with prior observations of the Murray formation. Uncommon outcrops of low-angle stratification suggest possible subaqueous currents, and steeply inclined beds may be the result of slumping. Collectively, the rocks exposed at VRR provide additional evidence for a long-lived lacustrine environment (in excess of 10 6 years via comparison to terrestrial records of sedimentation), which extends our understanding of the duration of habitable conditions in Gale crater.

Journal of Geophysical Research: Planets↗
USGS2020

The chemostratigraphy of the Murray formation and role of diagenesis at Vera Rubin ridge in Gale crater, Mars, as observed by the ChemCam instrument

Geochemical results are presented from Curiosity’s exploration of the Vera Rubin ridge (VRR), in addition to the full chemostratigraphy of the predominantly lacustrine mudstone Murray formation up to and including VRR. VRR is a prominent ridge flanking Aeolis Mons (informally Mt. Sharp), the central mound in Gale crater, Mars, and was a key area of interest for the Mars Science Laboratory mission. ChemCam data show that VRR is overall geochemically similar to lower-lying members of the Murray formation, even though the top of VRR is characterized by strong hematite spectral signatures as observed from orbit. While overall geochemically similar, VRR is characterized by a prominent decrease in Li abundance and Chemical Index of Alteration across the ridge. This decrease follows the morphology of the ridge rather than elevation and is inferred to reflect a non-depositionally controlled decrease in clay mineral abundance in VRR rocks. Additionally, a notable enrichment in Mn above baseline levels is observed on VRR. While not supporting a single model, the results suggest that VRR rocks were likely affected by multiple episodes of post-depositional groundwater interactions that made them more erosionally resistant than surrounding Murray rocks – thus resulting in the modern-day ridge after subsequent erosion.

Journal of Geophysical Research: Planets↗
USGS2019

Combining numerical and statistical models to predict storm-induced dune erosion

Dune erosion is an important aspect to consider when assessing coastal flood risk, as dune elevation loss makes the protected areas more susceptible to flooding. However, most advanced dune erosion numerical models are computationally expensive, which hinders their application in early-warning systems. Based on a combination of probabilistic and process-based numerical modeling, we develop an efficient statistical tool to predict dune erosion during storms. The analysis focuses on Dauphin Island, AL in the northern Gulf of Mexico, where we combine synthetic sea storms with a calibrated and validated XBeach model to develop and test a range of different surrogate models for their ability to predict barrier-island geometric parameters under storm conditions. Surrogate models are developed by combining the oceanographic forcing from 100 optimally sampled sea-storm events covering the entire multivariate parameter space (used as XBeach input) and associated changes in the dune system (XBeach output). We test four surrogate models using a k-fold approach for validation. All models perform well in predicting changes in dune elevation, barrier-island area, and width, but are less accurate in predicting alterations in the cross-shore locations of dune morphological features. Multivariate adaptive regression splines are identified as the best surrogate model based on their fast development and good performance, attaining a modified Mielke index of 0.81 for dune crest height. As demonstrated at Dauphin Island, our approach shows potential to be used in an operational framework to predict dune response (in particular crest elevation change) when water level and wave forecasts are available.

Alabama↗
USGS2019

Widespread initiation, reactivation, and acceleration of landslides in the northern California Coast Ranges due to extreme rainfall

Episodically to continuously active slow-moving landslides are driven by precipitation. Climate change, which is altering both the frequency and magnitude of precipitation world21 wide, is therefore predicted to have a major impact on landslides. Here we examine the behavior of hundreds of slow-moving landslides in northern California in response to large changes in annual precipitation that occurred between 2016 and 2018. We quantify the landslide displacement using repeat-pass radar interferometry and pixel offset tracking techniques on a novel dataset from the airborne NASA/JPL Uninhabited Aerial Vehicle Synthetic Aperture Radar. We found that 312 landslides were moving due to extreme rainfall during 2017, compared to 119 during 2016, which was the final year of a historic multi-year drought. However, with a return to below average rainfall in 2018, only 146 landslides remained in motion. The increased number of landslides during 2017 was primarily accommodated by landslides that were smaller than the landslides that remained active between 2016 and 2018. Furthermore, by examining a subset of 51 landslides, we found that 49 had increased velocities during 2017 when compared to 2016. Our results show that slow-moving landslides are sensitive to large changes in annual precipitation, particularly the smaller and thinner landslides that likely experience larger basal pore water pressure changes. Based on climate model predictions for the next century in California, which include increases in annual precipitation and increases in the frequency of dry-to-wet extremes, we hypothesize that there will be an overall increase in landslide activity.

California↗
USGS2020

Causes of variability in suspended‐sand concentration evaluated using measurements in the Colorado River in Grand Canyon

Rivers commonly exhibit substantial variability in suspended‐sand concentration, even at constant water discharge. Here we derive an approach for evaluating how much of this variability arises from mean bed‐sand grain size. We apply this approach to the Colorado River in Grand Canyon, where discharge‐independent concentration of suspended sand varies by more than a factor of 23 ( N = 1.4 × 10 6 ). Theory predicts that where concentration is controlled by bed‐sand grain size, concentration and grain size in suspension will be inversely correlated (i.e., coarsening of the bed causes suspended sand to become coarser in grain size and lower in concentration). Although the observed correlation is negative, riverbed grain size accounts for only 40% of the variability in concentration. The residuals vary by an order of magnitude; they arise from other processes, such as changes in topography or distribution of sand that cause shear stress to change at constant discharge, changes in the fine tail of bed‐sand grain sizes or changing bedforms. Both bed sand and the other factors influence concentration for durations from less than 1 day to several years. Predictions of concentration based on bed‐sand grain size ( N = 4 × 10 4 ) are less accurate than predictions based on suspended‐sand grain size, probably because suspended sand is a natural integrator of sand‐transporting processes, giving more weight to those areas of the bed that exchange more sand with the flow. Although the causes of variability vary from one river to another, the approach illustrated here is applicable to any river in which concentration varies at constant water discharge.

Arizona↗
USGS2020

Water tracks enhance water flow above permafrost in upland Arctic Alaska hillslopes

Upland permafrost regions occupy approximately one third of the Arctic landscape. In upland regions, hydrologic fluxes are influenced by water tracks, curvilinear features on hillslopes that preferentially fill with and route water in response to snowmelt and rainfall when the soil above continuous permafrost thaws in the summer. As continued warming of the Arctic may alter hydrologic cycling leading to increased frequency of extreme hydrologic events like drought and flooding as well as modification of biogeochemical cycling, it is imperative to untangle the interplay between precipitation, runoff, and subsurface flow as water is routed from upland Arctic regions to the Arctic Ocean. This study quantifies how ground surface temperatures affect groundwater discharge from hillslopes with water tracks in the upland Arctic by employing a three-dimensional, physically based subsurface flow model with variable saturation and freeze and thaw capabilities that is calibrated to field measurements from the Upper Kuparuk River watershed on the North Slope of Alaska, USA. Model analysis indicates that higher ground surface temperatures along water track hillslopes promote increases in groundwater discharge where water tracks act as conduits for large recharge events and continue to discharge groundwater into the autumn after the adjacent hillslope has frozen. Simulating the conditions that distinguish water tracks from their hillslope watersheds changes subsurface water storage and ground thermal responses but does not alter the total magnitude of groundwater discharge outside of parameter uncertainty. These findings suggest that water tracks play a complex and critical role in hydrologic cycles of the upland Arctic.

Alaska↗
USGS2020

Seasonal subsurface thaw dynamics of an aufeis feature inferred from geophysical methods

Aufeis are sheets of ice unique to cold regions that originate from repeated flooding and freezing events during the winter. They have hydrological importance associated with summer flows and possibly winter insulation, but little is known about the seasonal dynamics of the unfrozen sediment layer beneath them. This layer may support perennial groundwater flow in regions with otherwise continuous permafrost. For this study, ground penetrating radar (GPR) were collected in September 2016 (maximum thaw) and April 2017 (maximum frozen) at the Kuparuk aufeis field on the North Slope of Alaska. Supporting surface nuclear magnetic resonance (NMR) data were collected during the maximum frozen campaign. These point-in-time geophysical data sets were augmented by continuous subsurface temperature data and periodic Structure-from-Motion digital elevation models (DEM) collected seasonally. GPR and difference DEM data showed maximum ice thicknesses of up to 6 m over the sediment surface. Below the ice, GPR and NMR identified regions of permafrost and regions of seasonally frozen sediment (i.e., the active layer) underlain by a substantial perennially unfrozen zone or “lateral talik” that ranged from 0 m to over 13 m thick. The seasonally frozen cobble layer above the talik was typically 3 to 5 m thick, with freezing apparently enabled by relatively high thermal diffusivity of the overlying ice and rock cobbles. The large talik beneath the aufeis and active layer suggests that year-round groundwater flow and coupled heat transport occurs beneath much of the feature. Highly permeable alluvial material and discrete zones of apparent groundwater upwelling indicated by geophysical and ground temperature data allows direct connection between the aufeis and the talik below.

Alaska↗
USGS2020

Reconstructing the velocity and deformation of a rapid landslide using multiview video

Noncontact measurements of spatially varied ground surface deformation during landslide motion can provide important constraints on landslide mechanics. Here, we present and test a new method for extracting measurements of rapid landslide surface displacement and velocity (accelerations of approximately 1 m/s 2 ) using sequences of stereo images obtained from a pair of inexpensive, stationary 4K video cameras with nominal frame rates of 29.97 Hz. The method combines elements of Structure from Motion with those of optical flow to extract data on 3‐D evolution of the ground surface during slope failure. We apply the method to an experiment at the U.S. Geological Survey debris‐flow flume in which a high‐speed, liquefying landslide was triggered by gradually adding water to a 6‐m 3 prism of loosely packed sediment on a 31° slope. Strip‐scanning lidar measurements made during the experiment corroborate our video‐based measurements, but the latter encompassed the entire landslide surface and were much lower in cost. Our video‐based measurements enabled computation of depth‐integrated landslide dilation/contraction rates. The range of computed rates was within the ranges inferred from independent measurements of evolving pore water pressures and reasonable estimates of the hydraulic permeability of the sediment. Dilation and contraction rates play a crucial role in landslide mechanics. The dilation and contraction we observe contradict the incompressible flow assumption used in many studies that have employed noncontact methods to infer landslide properties.

Oregon↗
USGS2020

The influence of frost weathering on the debris flow sediment supply in an alpine basin

Rocky, alpine mountains are prone to mass wasting from debris flows. The Chalk Cliffs study area (central Colorado, USA) produces debris flows annually. These debris flows are triggered when overland flow driven by intense summer convective storms mobilizes large volumes of sediment within the channel network. Understanding the debris flow hazard in this, and similar alpine settings, requires determining the magnitude of sed- iment accumulation between debris flow seasons, and identifying the control on sediment production. To address these knowledge gaps, we measured changes in sediment produc- tion using a sediment retention fence to quantify how sedimentation was influenced by temperature at the plot scale. These measurements were extrapolated to a larger area, where we extended the sediment fence results to explore how rockfall sedimentation con- tributed to channel refilling between active debris flow periods. This work shows debris flow channel refilling is correlated with low temperatures and time in the frost-cracking window, implicating frost weathering mechanisms as a key driver of sedimentation. This sediment production process resulted in a large amount of sediment accumulation dur- ing a single winter season in our study reach (up to 0.4 m in some locations). Using these observations, we develop a channel refilling model that generally describes the mass bal- ance of debris flow watersheds in alpine areas.

Colorado↗
USGS2020

Does channel narrowing by floodplain growth necessarily indicate sediment surplus? Lessons from sediment‐transport analyses in the Green and Colorado rivers, Canyonlands, Utah

Analyses of suspended sediment transport provide valuable insight into the role that sediment supply plays in causing geomorphic change. The sediment supply within a river system evolves depending on the discharge, flood frequency and duration, changes in sediment input, and ecohydraulic conditions that modify sediment transport processes. Changes in supply can be evaluated through analyses of coupled changes in suspended sediment concentration and grain size. The concentration of sand in transport in the Green and Colorado Rivers is most strongly controlled by discharge and the bed sand grain size distribution. Since the 1950s, sand loads have decreased in response to declines in peak discharge in the Green River and coarsening of the bed sand in the Colorado River. However, changes in the bed sand grain size distribution are associated with large changes in suspended sand concentration in both rivers; concentration varies by a factor of ~3 in the Green River and a factor of ~8 in the Colorado River, depending on the bed sand grain size distribution. Analyses of hysteresis in suspended sediment measurements show that sediment depletion during annual floods is most strongly controlled by flood duration, with peak discharge being nearly equally important in the Green River. Despite channel narrowing in both rivers, periods of bed sand coarsening and sediment depletion during annual floods indicate that these rivers are not necessarily in sediment surplus. Channel narrowing appears to be strongly controlled by short‐term declines in flood magnitude and the ecohydraulic effects of vegetation and may not be indicative of the long‐term sediment budget.

Utah↗
USGS2020

Sea‐level rise will drive divergent sediment transport patterns on fore reefs and reef flats, potentially causing erosion on atoll islands

Atoll reef islands primarily consist of unconsolidated sediment, and their ocean‐facing shorelines are maintained by sediment produced and transported across their reefs. Changes in incident waves can alter cross‐shore sediment exchange and, thus, affect the sediment budget and morphology of atoll reef islands. Here we investigate the influence of sea level rise and projected wave climate change on wave characteristics and cross‐shore sediment transport across an atoll reef at Kwajalein Island, Republic of the Marshall Islands. Using a phase‐resolving model, we quantify the influence on sediment transport of quantities not well captured by wave‐averaged models, namely, wave asymmetry and skewness and flow acceleration. Model results suggest that for current reef geometry, sea level, and wave climate, potential bedload transport is directed onshore, decreases from the fore reef to the beach, and is sensitive to the influence of flow acceleration. We find that a projected 12% decrease in annual wave energy by 2100 CE has negligible influence on reef flat hydrodynamics. However, 0.5–2.0 m of sea level rise increases wave heights, skewness, and shear stress on the reef flat and decreases wave skewness and shear stress on the fore reef. These hydrodynamic changes decrease potential sediment inputs onshore from the fore reef where coral production is greatest but increase potential cross‐reef sediment transport from the outer reef flat to the beach. Assuming sediment production on the fore reef remains constant or decreases due to increasing ocean temperatures and acidification, these processes have the potential to decrease net sediment delivery to atoll islands, causing erosion.

Journal of Geophysical Research – Earth Surface↗
USGS2019

Carbon chemistry of intact versus chronically drained peatlands in the southeastern USA

The Great Dismal Swamp (GDS) is a large temperate swamp in Virginia/North Carolina with peat soils historically resistant to microbial decomposition. However, this peatland has been subject to ~200 years of disturbance during which extensive drainage, fire suppression, and wide-spread logging have increased decomposition and dramatically decreased the distribution of Atlantic white cedar (AWC). The purpose of this study was to determine the impact of long-term drainage and AWC loss on the carbon chemistry of GDS peats. Peat cores were collected from three drained GDS vegetation communities (pocosin, AWC, and red maple-black gum) and compared to cores collected from an intact, undrained AWC peatland at the Alligator River National Wildlife Refuge (AR) in North Carolina, USA. The AR peats had higher lignin content in the deeper peat intervals, and lignin content and % organic carbon were largely invariant with depth compared to the GDS peats. The concentrations of syringyl group phenols were greater in the surface layers of GDS peats, likely reflecting the selective removal of AWC and transition from gymnosperms to angiosperms. Acid to aldehyde ratios for vanillyl and syringyl group phenols indicated that the GDS peats were more decomposed, particularly at depth, and that this occurred under aerobic conditions. Moreover, solid state 13C NMR confirmed a coincident loss of carbohydrates and increase in recalcitrant byproducts of carbohydrate degradation with depth. These data indicate that long-term drainage has accelerated the decomposition of peat at the GDS, reducing the capacity and stability of the carbon sink.

Virginia, North Carolina↗
USGS2019

Complex response of sediment phosphorus to land use and management within a river network

Rivers affected by anthropogenic nutrient inputs can retain some of the phosphorus (P) load through sediment retention and burial. Determining the influence of land use and management on sediment P concentrations and P retention in fluvial ecosystems is challenging because of different stressors operating at multiple spatial and temporal scales. In this study, we sought to determine how land use and management influenced sediment P concentrations and P retention within a river network draining a watershed modified by mixed land use activities; the Fox River, Wisconsin, USA. Results show current land use had no detectable effect on sediment P concentrations, and only a small potential effect on P retention capacity. However, sites draining predominantly forested areas were associated with riverbed sediments less saturated in P, whereas sites draining mainly agricultural areas were more likely to release sediment-bound P. Current management actions, including the implementation of best management practices, had a small positive effect on P retention capacity but had no observed effect on sediment P concentrations. Our results suggest that fine riverbed sediment throughout the Fox River network is saturated with P and has the capability to release P when water column P concentrations are low. These P-saturated riverbed sediments represent a potential legacy effect of past land use disturbances on P dynamics in the Fox River network.

Wisconsin↗
USGS2019

Heterogeneity in hyporheic flow, pore water chemistry, and microbial community composition in an alpine streambed

The hyporheic zone, where surface water and groundwater mix, is an important microbial habitat where biogeochemical reactions influence water quality. We show that spatial variability in hyporheic flow in the East River near Crested Butte, Colorado, drives heterogeneity in streambed geochemical conditions and microbial community assemblages, but the diversity of assemblages remains nearly constant throughout the reach. In July 2018, we collected approximately 100 pore water samples at 20 cm depth and analyzed them for anions, cations, dissolved organic carbon, dissolved organic matter (DOM) quality, and basic water quality parameters. Vertical hydraulic head gradients were also measured to assess the potential for upward or downward flow, and heat tracing was used to quantify vertical flux rates at a subset of locations. We found that regions of the streambed that are more groundwater-dominated contain less dissolved oxygen, higher concentrations of reduced metals, and more microbially-processed, recalcitrant DOM, while more surface water-dominated locations contain higher dissolved oxygen concentrations and terrestrially-derived, labile DOM. 16S rRNA gene sequencing of extracted DNA revealed that microbial community composition varies with geochemical gradients related to hyporheic flow. These findings provide a better understanding of hyporheic controls on streambed biogeochemistry during the baseflow season, which is expected to lengthen with climate change in alpine watersheds due to earlier snowmelt onset and reduced snowpack.

Colorado↗
USGS2020

Increasing rates of carbon burial in southwest Florida coastal wetlands

Rates of organic carbon (OC) burial in some coastal wetlands appear to be greater in recent years than they were in the past. Possible explanations include ongoing mineralization of older OC or the influence of an unaccounted‐for artefact of the methods used to measure burial rates. Alternatively, the trend may represent real acceleration in OC burial. We quantified OC burial rates of mangrove and coastal freshwater marshes in southwest Florida through a comparison of rates derived from 210 Pb, 137 Cs, and surface marker horizons (MHs). Age/depth profiles of lignin: OC were used to assess whether down‐core remineralization had depleted the OC pool relative to lignin, and lignin phenols were used to quantify the variability of lignin degradation. Over the past 120 years, OC burial rates at seven sites increased by factors ranging from 1.4 to 6.2. We propose that these increases represent net acceleration. Change in relative sea‐level rise is the most likely large‐scale driver of acceleration, and sediment deposition from large storms can contribute to periodic increases. Mangrove sites had higher OC and lignin burial rates than marsh sites, indicating inherent differences in OC burial factors between the two habitat types. The higher OC burial rates in mangrove soils mean that their encroachment into coastal freshwater marshes has the potential to increase burial rates in those locations even more than might be expected from the acceleration trends. Regionally, these findings suggest that burial represents a substantially growing proportion of the coastal wetland carbon budget.

Florida↗
USGS2020

Impacts of seagrass dynamics on the coupled long-term evolution of barrier-marsh-bay systems

Seagrass provides a wide range of economically and ecologically valuable ecosystem services, with shoreline erosion control often listed as a key service, but can also alter the sediment dynamics and waves within back-barrier bays. Here we incorporate seagrass dynamics into an existing barrier-marsh exploratory model, GEOMBEST++, to examine the coupled interactions of the back-barrier bay with both adjacent (marsh) and non-adjacent (barrier island) subsystems. While seagrass reduces marsh edge erosion rates and increases progradation rates in many of our 288 model simulations, seagrass surprisingly increases marsh edge erosion rates when sediment export from the back-barrier basin is negligible because the ability of seagrass to reduce the volume of marsh sediment eroded matters little for back-barrier basins in which all sediment is conserved. Our model simulations also suggest that adding seagrass to the bay subsystem leads to increased deposition in the bay, reduced sediment available to the marsh, and enhanced marsh edge erosion until the bay reaches a new, shallower equilibrium depth. In contrast, removing seagrass liberates previously-sequestered sediment that is then delivered to the marsh, leading to enhanced marsh progradation. Lastly, we find that seagrass reduces barrier island migration rates in the absence of back-barrier marsh by filling accommodation space in the bay. These model observations suggest that seagrass meadows operate as dynamic sources and sinks of sediment that can influence the evolution of coupled marsh and barrier island landforms in unanticipated ways.

Virginia↗
USGS2021

Estimating the impact of seep methane oxidation on ocean pH and dissolved inorganic radiocarbon along the U.S. mid‐Atlantic Bight

Ongoing ocean warming can release methane (CH 4 ) currently stored in ocean sediments as free gas and gas hydrates. Once dissolved in ocean waters, this CH 4 can be oxidized to carbon dioxide (CO 2 ). While it has been hypothesized that the CO 2 produced from aerobic CH 4 oxidation could enhance ocean acidification, a previous study conducted in Hudson Canyon shows that CH 4 oxidation has a small short‐term influence on ocean pH and dissolved inorganic radiocarbon. Here we expand upon that investigation to assess the impact of widespread CH 4 seepage on CO 2 chemistry and possible accumulation of this carbon injection along 234 km of the U.S. Mid‐Atlantic Bight. Consistent with the estimates from Hudson Canyon, we demonstrate that a small fraction of ancient CH 4 ‐derived carbon is being assimilated into the dissolved inorganic radiocarbon (mean fraction of 0.5 ± 0.4 %). The areas with the highest fractions of ancient carbon coincide with elevated CH 4 concentration and active gas seepage. This suggests that aerobic CH 4 oxidation has a greater influence on the dissolved inorganic pool in areas where CH 4 concentrations are locally elevated, instead of displaying a cumulative effect downcurrent from widespread groupings of CH 4 seeps. An upper limit approximation of the input rate of ancient‐derived DIC into the waters overlying the northern U. S Mid‐Atlantic Bight further suggests that oxidation of ancient CH 4 ‐derived carbon is not negligible on the global scale and could contribute to deep‐water acidification over longer time scales.

Journal of Geophysical Research- Biogeosciences↗
USGS2020

Combined seismic and geodetic analysis before, during and after the 2018 Mt. Etna eruption

In December 2018, Etna volcano experienced one of the largest episodes of unrest since the installation of geophysical monitoring networks in 1970. The unrest culminated in a short eruption with a small volume of lava erupted, a significant seismic crisis and deformation of the entire volcanic edifice of magnitude never recorded before at Mount Etna. Here we describe the evolution of the 2018 eruptive cycle from the analysis of seismic and geodetic data collected in the months preceding, during, and following the intrusion. We model the space-time evolution of high-rate deformation data starting from the active source previously identified from deformation data and the propagation of seismicity in a 3-D velocity model. The intrusion model suggests emplacement of two dikes: a smaller dike located beneath the eruptive fissure and a second, deeper dike between 1 and 5 km below sea level that opened ~2 m. The rise and eruption of magma from the shallower dike did not interrupt the pressurization of a long-lasting deeper reservoir (~6 km) that induced continuous inflation and intense deformation of the eastern flank. Shortly after the intrusion, on 26 December 2018, a M L 4.8 earthquake occurred near Pisano, destroying buildings and roads in two villages. We propose a time-dependent intrusion model that supports the hypothesis of the inflation inducing flank deformation and that this process has been active since September 2018.

Mt. Etna↗