Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “The Earth Observer”

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 1,351 records · Page 75Linked to original sources
USGS2018

Unraveling the dynamics of magmatic CO2 degassing at Mammoth Mountain, California

The accumulation of magmatic CO 2 beneath low-permeability barriers may lead to the formation of CO 2 -rich gas reservoirs within volcanic systems. Such accumulation is often evidenced by high surface CO 2 emissions that fluctuate over time. The temporal variability in surface degassing is believed in part to reflect a complex interplay between deep magmatic degassing and the permeability of degassing pathways. A better understanding of the dynamics of CO 2 degassing is required to improve monitoring and hazards mitigation in these systems. Owing to the availability of long-term records of CO 2 emissions rates and seismicity, Mammoth Mountain in California constitutes an ideal site towards such predictive understanding. Mammoth Mountain is characterized by intense soil CO 2 degassing (up to ∼1000 t d −1 ) and tree kill areas that resulted from leakage of CO 2 from a CO 2 -rich gas reservoir located in the upper ∼4 km. The release of CO 2 -rich fluids from deeper basaltic intrusions towards the reservoir induces seismicity and potentially reactivates faults connecting the reservoir to the surface. While this conceptual model is well-accepted, there is still a debate whether temporally variable surface CO 2 fluxes directly reflect degassing of intrusions or variations in fault permeability. Here, we report the first large-scale numerical model of fluid and heat transport for Mammoth Mountain. We discuss processes (i) leading to the initial formation of the CO 2 -rich gas reservoir prior to the occurrence of high surface CO 2 degassing rates and (ii) controlling current CO 2 degassing at the surface. Although the modeling settings are site-specific, the key mechanisms discussed in this study are likely at play at other volcanic systems hosting CO 2 -rich gas reservoirs. In particular, our model results illustrate the role of convection in stripping a CO 2 -rich gas phase from a rising hydrothermal fluid and leading to an accumulation of a large mass of CO 2 (∼10 7 –10 8 t) in a shallow gas reservoir. Moreover, we show that both, short-lived (months to years) and long-lived (hundreds of years) events of magmatic fluid injection can lead to critical pressures within the reservoir and potentially trigger fault reactivation. Our sensitivity analysis suggests that observed temporal fluctuations in surface degassing are only indirectly controlled by variations in magmatic degassing and are mainly the result of temporally variable fault permeability. Finally, we suggest that long-term CO 2 emission monitoring, seismic tomography and coupled thermal–hydraulic–mechanical modeling are important for CO 2 -related hazard mitigation.

California↗
USGS2005

Forecasting the evolution of seismicity in southern California: Animations built on earthquake stress transfer

We develop a forecast model to reproduce the distibution of main shocks, aftershocks and surrounding seismicity observed during 1986-200 in a 300 ?? 310 km area centered on the 1992 M = 7.3 Landers earthquake. To parse the catalog into frames with equal numbers of aftershocks, we animate seismicity in log time increments that lengthen after each main shock; this reveals aftershock zone migration, expansion, and densification. We implement a rate/state algorithm that incorporates the static stress transferred by each M ??? 6 shock and then evolves. Coulomb stress changes amplify the background seismicity, so small stress changes produce large changes in seismicity rate in areas of high background seismicity. Similarly, seismicity rate declines in the stress shadows are evident only in areas with previously high seismicity rates. Thus a key constituent of the model is the background seismicity rate, which we smooth from 1981 to 1986 seismicity. The mean correlation coefficient between observed and predicted M ??? 1.4 shocks (the minimum magnitude of completeness) is 0.52 for 1986-2003 and 0.63 for 1992-2003; a control standard aftershock model yields 0.54 and 0.52 for the same periods. Four M ??? 6.0 shocks struck during the test period; three are located at sites where the expected seismicity rate falls above the 92 percentile, and one is located above the 75 percentile. The model thus reproduces much, but certainly not all, of the observed spatial and temporal seismicity, from which we infer that the decaying effect of stress transferred by successive main shocks influences seismicity for decades. Finally, we offer a M ??? 5 earthquake forecast for 2005-2015, assigning probabilities to 324 10 ?? 10 km cells.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Planetary sensor models interoperability using the community sensor model specification

This paper presents the photogrammetric foundations upon which the Community Sensor Model specification depends, describes common coordinate system and reference frame transformations that support conversion between image sensor (charge‐coupled device) coordinates to some arbitrary body coordinate, and describes the U.S. Geological Survey Astrogeology Community Sensor Model implementation ( https://github.com/USGS-Astrogeology/usgscsm ). We present a new image support data specification that provides the position, pointing, timing, and metadata information necessary to properly locate a pixel or observations location on a body and describe a system architecture designed to explicitly identify the responsibilities of software components within a larger pipeline or analytical environment. This paper concludes with a set of experiments that illustrate positional and pointing error in the sensor location and the impact on the computed surface location.

Earth and Space Science↗
USGS2023

Vulnerability of estuarine systems in the contiguous United States to water quality change under future climate and land-use

Changes in climate and land-use and land-cover (LULC) are expected to influence surface water runoff and nutrient characteristics of estuarine watersheds, but the extent to which estuaries are vulnerable to altered nutrient loading under future conditions is poorly understood. The present work aims to address this gap through the development of a new vulnerability assessment framework that accounts for (a) estuarine exposure to projected changes in total nitrogen (TN) and total phosphorus (TP) loads as a function of LULC and climate change under several scenarios, (b) sensitivity, and (c) adaptive capacity. The framework was applied to 112 estuaries and their contributing watersheds across the contiguous U.S., specifically to look at regional variability in estuarine vulnerability to nutrient loading. Study findings revealed that the largest increases in estuarine nutrient loads are expected in the North and South Atlantic regions and eastern Gulf of Mexico, while the lowest increases are expected in the North and South Pacific regions and the western Gulf of Mexico. However, the North Atlantic and the South Pacific had the highest adaptive capacity, which could potentially counteract the effects of LULC and climate change on nutrient loads. Strong variation in predicted estuarine nutrient loads was observed as a function of climate model projections, while projected LULC changes were more consistently associated with elevated loads. Our findings illustrate the benefits of integrating natural and socio-ecological factors to identify opportunities to develop adaptation plans and policies to mitigate ecological degradation in vitally important estuaries.

Earth's Future↗
USGS2016

Short-period volcanic gas precursors to phreatic eruptions: Insights from Poás Volcano, Costa Rica

Volcanic eruptions involving interaction with water are amongst the most violent and unpredictable geologic phenomena on Earth. Phreatic eruptions are exceptionally difficult to forecast by traditional geophysical techniques. Here we report on short-term precursory variations in gas emissions related to phreatic blasts at Poás volcano, Costa Rica, as measured with an in situ multiple gas analyzer that was deployed at the edge of the erupting lake. Gas emitted from this hyper-acid crater lake approaches magmatic values of SO 2 /CO 2 1–6 days prior to eruption. The SO 2 flux derived from magmatic degassing through the lake is measureable by differential optical absorption spectrometry (sporadic campaign measurements), which allows us to constrain lake gas output and input for the major gas species during eruptive and non-eruptive periods. We can further calculate power supply to the hydrothermal system using volatile mass balance and thermodynamics, which indicates that the magmatic heat flux into the shallow hydrothermal system increases from ∼27 MW during quiescence to ∼59 MW during periods of phreatic events. These transient pulses of gas and heat from the deeper magmatic system generate both phreatic eruptions and the observed short-term changes in gas composition, because at high gas flux scrubbing of sulfur by the hydrothermal system is both kinetically and thermodynamically inhibited whereas CO 2 gas is always essentially inert in hyperacid conditions. Thus, the SO 2 /CO 2 of lake emissions approaches magmatic values as gas and power supply to the sub-limnic hydrothermal system increase, vaporizing fluids and priming the hydrothermal system for eruption. Our results suggest that high-frequency real-time gas monitoring could provide useful short-term eruptive precursors at volcanoes prone to phreatic explosions.

Poás Volcano↗
USGS1984

Crustal structure of the Appalachian Highlands in Tennessee

Crustal structure of the southern Appalachians and adjacent Interior Low Plateaus in Tennessee is derived from seismic-refraction measurements observed by the U.S. Geological Survey in 1965 along reversed lines, normal (NW-SE) and parallel (NE-SW) to the structure of the Appalachian Highlands' major geologic divisions. Its easternmost part is located approximately 80 km southwest of the westernmost part of the COCORP seismic-reflection traverse within the Blue Ridge province. The velocity-depth models derived for both observational directions consist of three crustal layers with surprisingly high velocities, being about 6.1-6.2 km/s in the upper crust down to 7-10 km depth, 6.7-6.8 km/s for the middle crust between about 17 and 34 km and varying from 7.1 to 7.4 km/s for the lower crust at about 40-47 km depth. The boundaries between the three crustal layers as well as the crust-mantle boundary are transition zones of up to 11 km thickness. Similar to old orogens in other parts of the earth, the main result is a thick crust, at places in excess of 50 km, with high average velocity and a broad crust-mantle transition zone.

Tennessee↗
USGS2014

Characterization of very-long-period seismicity accompanying summit activity at Kīlauea Volcano, Hawai'i: 2007-2013

Eruptive activity returned to the summit region of Kīlauea Volcano, Hawai'i with the formation of the “Overlook crater” within the Halema'uma'u Crater in March 2008. The new crater continued to grow through episodic collapse of the crater walls and as of late 2013 had grown into an approximately elliptical opening with dimensions of ~ 160 × 215 m extending to a depth of ~ 200 m. Occasional weak explosive events and a persistent gas plume continued to occur through 2013. Lava was first observed in the new crater in September 2008, and through 2009 the lava level remained deep in the crater and was only occasionally observed. Since early 2010 a lava lake with fluctuating level within the Overlook crater has been nearly continuously present, and has reached to within 22 m of the Overlook crater rim. Volcanic activity at Kīlauea Volcano is episodic at all time scales and the characterization of very-long-period seismicity in the band 2–100 s for the years 2007–2013 illuminates a portion of this broad spectrum of volcanic behavior. Three types of very-long-period events have been observed over this time and each is associated with distinct processes. Type 1 events are associated with vigorous degassing and occurred primarily between 2007 and 2009. Type 2 events are associated with rockfalls onto the lava lake and occurred primarily after early 2010. Both of these event types are induced by pressure and momentum changes at the top of the magma column that are transmitted downward to a source centroid ~ 1 km below the northeast corner of the Halema'uma'u Crater where the energy couples to the solid Earth at a geometrical discontinuity in the underlying dike system. Type 3 events are not related to surficial phenomena but are associated with transients in mass transfer that occur within the dike system. Very-long-period tremor has also accompanied the return of eruptive activity, with increasing amplitude associated with hours- to months-long changes in gas emission rates and summit deformation.

Hawaii↗
USGS1991

Fractal patterns of fractures in granites

Fractal measurements using the Cantor's dust method in a linear one-dimensional analysis mode were made on the fracture patterns revealed on two-dimensional, planar surfaces in four granites. This method allows one to conclude that: (1)|The fracture systems seen on two-dimensional surfaces in granites are consistent with the part of fractal theory that predicts a repetition of patterns on different scales of observation, self similarity. Fractal analysis gives essentially the same values of D on the scale of kilometres, metres and centimetres (five orders of magnitude) using mapped, surface fracture patterns in a Sierra Nevada granite batholith (Mt. Abbot quadrangle, Calif.). (2)|Fractures show the same fractal values at different depths in a given batholith. Mapped fractures (main stage ore veins) at three mining levels (over a 700 m depth interval) of the Boulder batholith, Butte, Mont. show the same fractal values although the fracture disposition appears to be different at different levels. (3)|Different sets of fracture planes in a granite batholith, Central France, and in experimental deformation can have different fractal values. In these examples shear and tension modes have the same fractal values while compressional fractures follow a different fractal mode of failure. The composite fracture patterns are also fractal but with a different, median, fractal value compared to the individual values for the fracture plane sets. These observations indicate that the fractal method can possibly be used to distinguish fractures of different origins in a complex system. It is concluded that granites fracture in a fractal manner which can be followed at many scales. It appears that fracture planes of different origins can be characterized using linear fractal analysis.

Earth and Planetary Science Letters↗
USGS2012

A robust method to forecast volcanic ash clouds

Ash clouds emanating from volcanic eruption columns often form trails of ash extending thousands of kilometers through the Earth's atmosphere, disrupting air traffic and posing a significant hazard to air travel. To mitigate such hazards, the community charged with reducing flight risk must accurately assess risk of ash ingestion for any flight path and provide robust forecasts of volcanic ash dispersal. In response to this need, a number of different transport models have been developed for this purpose and applied to recent eruptions, providing a means to assess uncertainty in forecasts. Here we provide a framework for optimal forecasts and their uncertainties given any model and any observational data. This involves random sampling of the probability distributions of input (source) parameters to a transport model and iteratively running the model with different inputs, each time assessing the predictions that the model makes about ash dispersal by direct comparison with satellite data. The results of these comparisons are embodied in a likelihood function whose maximum corresponds to the minimum misfit between model output and observations. Bayes theorem is then used to determine a normalized posterior probability distribution and from that a forecast of future uncertainty in ash dispersal. The nature of ash clouds in heterogeneous wind fields creates a strong maximum likelihood estimate in which most of the probability is localized to narrow ranges of model source parameters. This property is used here to accelerate probability assessment, producing a method to rapidly generate a prediction of future ash concentrations and their distribution based upon assimilation of satellite data as well as model and data uncertainties. Applying this method to the recent eruption of Eyjafjallajökull in Iceland, we show that the 3 and 6 h forecasts of ash cloud location probability encompassed the location of observed satellite-determined ash cloud loads, providing an efficient means to assess all of the hazards associated with these ash clouds.

Journal of Geophysical Research D: Atmospheres↗
USGS2021

Geometry of the décollement below eastern Bangladesh and implications for seismic hazard

Eastern Bangladesh sits on the seismically active Chittagong-Myanmar fold and thrust belt (CMFB), a north-trending accretionary wedge on the eastern side of the India-Eurasia collision. Earthquakes on the basal décollement and associated thrusts within the CMFB present a hazard to this densely populated region. In this study, we interpret 28 seismic reflection profiles from both published and unpublished sources to constrain the depth of the basal décollement. To convert profiles from the time domain to the depth domain, we integrate sonic log and seismic stacking velocity data to generate time-velocity relationships for different parts of the CMFB. Our analysis reveals that the décollement is ∼9 km deep in northeast and southeast Bangladesh, but shallows to ∼5 km in east-central Bangladesh. The décollement has an area of 7.25 × 10 4 km 2 (∼150 × 450 km), making it capable of an Mw 8.5 earthquake. However, the warped geometry of this fault might act as a rupture barrier were a large earthquake to occur on the décollement. Our combined velocity and fault model lay the groundwork for future studies to address seismic segmentation, ground shaking, and rupture modeling in the CMFB. Finally, we use our compiled data set to analyze the evolution of fold kinematics in the CMFB. We observe that folding style and failure mode varies, from mainly ductile deformation in the foreland to mainly brittle in the hinterland. The dual-failure modes within the CMFB support the hypothesis that a region with ductile deformation may still be capable of seismic behavior.

Journal of Geophysical Research: Solid Earth↗
USGS2012

Seismic constraints on the nature of lower crustal reflectors beneath the extending Southern Transition Zone of the Colorado Plateau, Arizona

We determine the reflection polarity and exploit variations in P and S wave reflectivity and P wave amplitude versus offset (AVO) to constrain the origin of lower crustal reflectivity observed on new three-component seismic data recorded across the structural transition of the Colorado Plateau. The near vertical incidence reflection data were collected by Stanford University in 1989 as part of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that traversed the Arizona Transition Zone of the Colorado Plateau. The results of independent waveform modeling methods are consistent with much of the lower crustal reflectivity resulting from thin, high-impedance layers. The reflection polarity of the cleanest lower crustal events is positive, which implies that these reflections result from high-velocity contrasts, and the waveform character indicates that the reflectors are probably layers less than or approximately equal to 200 m thick. The lower crustal events are generally less reflective to incident S waves than to P waves, which agrees with the predicted behavior of high-velocity mafic layering. Analysis of the P wave AVO character of lower crustal reflections demonstrates that the events maintain a constant amplitude with offset, which is most consistent with a mafic-layering model. One exception is a high-amplitude (10 dB above background) event near the base of lower crustal reflectivity which abruptly decreases in amplitude at increasing offsets. The event has a pronounced S wave response, which along with its negative AVO trend is a possible indication of the presence of fluids in the lower crust. The Arizona Transition Zone is an active but weakly extended province, which causes us to discard models of lower crustal layering resulting from shearing because of the high degree of strain required to create such layers. Instead, we favor horizontal basaltic intrusions as the primary origin of high-impedance reflectors based on (1) The fact that most xenoliths in eruptive basalts of the Transition Zone are of mafic igneous composition, (2) indications that a pulse of magmatic activity crossed the Transition Zone in the late Tertiary period, and (3) the high regional heat flow observed in the Transition Zone. The apparent presence of fluids near the base of the reflective zone may indicate a partially molten intrusion. We present a mechanism by which magma can be trapped and be induced to intrude horizontally at rheologic contrasts in extending crust.

Arizona↗
USGS2012

Morphometric variability within the axial zone of the southern Juan de Fuca Ridge: Interpretation from Sea MARC II, Sea MARC I, and deep-sea photography

The morphometric characteristics of the axial regions of oceanic spreading centers are determined by (1) the type of volcanic flows, (2) the relation between primary volcanic relief (on a scale of a few meters to tens of meters) and degree of sediment cover, and (3) the extent of surficial expression and timing of tectonic disruption of the young oceanic crust. Even within a single, continuous, linear spreading-ridge segment with relatively uniform axial valley dimensions over a distance of 50 or more kilometers, such as along the southern Juan de Fuca Ridge, the changes in morphometric characteristics along axis within the youngest crust indicate distinct variation in tectonic and volcanic activity over short distances within short time periods. An integrated analysis of Sea MARC I, Sea MARC II, and photographic data for the southernmost continuous segment of the Juan de Fuca Ridge shows that generalizations about tectonic and volcanic processes at spreading ridges must consider both the temporal scale of processes as well as the physical scales of observations if predictive models are to be successful. Comparison of the morphometric expression within the major hydrothermal vent area and the rest of the southernmost ridge segment suggests that the mapped distribution of hydrothermal vents may reflect the extent of survey effort rather than uniqueness of geologic setting.

Journal of Geophysical Research B: Solid Earth↗
USGS2014

Fluctuations of water‐surface in observation‐wells and at stream gaging‐stations in the Mokelumne Area, California, during the earthquake of December 20, 1932

On December 20, 1932, much of the western United States felt a pronounced earthquake‐shock which caused noticeable fluctuations of the water‐surface in several observation‐wells and at two stream gaging‐stations in the Mokelumne Area, central California. These effects of the earthquake are reported briefly in this paper. According to Dr. Byerly (personal communication, February 9, 1933), in charge of the seismograph‐station of the University of California at Berkeley, the earthquake‐ vibrations at that place began at 10 h 11 m 00 s p.m., Pacific Standard Time, December 20, 1932, and continued three hours on the records of his most sensitive instruments. He reports further that the maximum double amplitude of the Earth‐motion at Berkeley was of the order of two mm and that the north‐south and east‐west components of the amplitude were approximately equal. In the Mokelumne Area, which centers about the City of Lodi, 57 miles north 70° east from the seismograph‐station at Berkeley, the pronounced initial Earth‐shock caused chandeliers to sway, upset some bric‐a‐brac, and in at least one building opened a few cracks in the interior plastering. At Stockton, 15 miles south of Lodi, the glass of at least one display‐window in a shop was cracked.

California↗
USGS2008

Investigation and hazard assessment of the 2003 and 2007 Staircase Falls rock falls, Yosemite National Park, California, USA

Since 1857 more than 600 rock falls, rock slides, debris slides, and debris flows have been documented in Yosemite National Park, with rock falls in Yosemite Valley representing the majority of the events. On 26 December 2003, a rock fall originating from west of Glacier Point sent approximately 200 m 3 of rock debris down a series of joint-controlled ledges to the floor of Yosemite Valley. The debris impacted talus near the base of Staircase Falls, producing fragments of flying rock that struck occupied cabins in Curry Village. Several years later on 9 June 2007, and again on 26 July 2007, smaller rock falls originated from the same source area. The 26 December 2003 event coincided with a severe winter storm and was likely triggered by precipitation and/or frost wedging, but the 9 June and 26 July 2007 events lack recognizable triggering mechanisms. We investigated the geologic and hydrologic factors contributing to the Staircase Falls rock falls, including bedrock lithology, weathering, joint spacing and orientations, and hydrologic processes affecting slope stability. We improved upon previous geomorphic assessment of rock-fall hazards, based on a shadow angle approach, by using STONE, a three-dimensional rock-fall simulation computer program. STONE produced simulated rock-fall runout patterns similar to the mapped extent of the 2003 and 2007 events, allowing us to simulate potential future rock falls from the Staircase Falls detachment area. Observations of recent rock falls, mapping of rock debris, and simulations of rock fall runouts beneath the Staircase Falls detachment area suggest that rock-fall hazard zones extend farther downslope than the extent previously defined by mapped surface talus deposits.

Natural Hazards and Earth System Sciences↗
USGS2020

Ancient Martian aeolian sand dune deposits recorded in the stratigraphy of Valles Marineris and implications for past climates

Aeolian sediment transport, deposition, and erosion have been ongoing throughout Mars's history. This record of widespread aeolian processes is preserved in landforms and geologic units that retain important clues about past environmental conditions including wind patterns. In this study we describe landforms within Melas Chasma, Valles Marineris, that occur in distinct groups with linear to crescentic shapes, arranged with a characteristic wavelength; some possess slope profiles analogous to modern sand dunes yet show evidence for lithification. Based on the features' dimensions, asymmetry, and spatial patterns relative to modern equivalents, we interpret these landforms to be two classes of aeolian bedforms: decameter‐scale megaripples and sand dunes. The presence of superposed erosional features and depositional units indicates that these landforms were cemented and likely ancient. Melas paleodunes are found atop Hesperian‐aged layered deposits, but we estimate them to be younger, likely lithified in the Amazonian period. Although a range of degradation was observed, some paleodunes are >10 m tall and maintain steep lee sides (>25°), an uncommon scenario for terrestrial examples as other geologic processes lead to dune obliteration. The preserved paleobedform geometries are largely consistent with those of modern aeolian indicators, suggesting no major shifts in wind regime or contributing boundary conditions. Finally, we propose that their appearance and context require sequential periods of dune migration, stabilization following catastrophic burial, cementation, differential erosion, exposure, and burial. The presence of wholly preserved duneforms appears to be more common on Mars compared to the Earth and may signal something important about Martian landscape evolution.

Journal of Geophysical Research: Planets↗
USGS2022

Orbital and in-situ investigation of periodic bedrock ridges in Glen Torridon, Gale Crater, Mars

Wind has been the dominant agent of landscape modification on Mars for the past ~3 billion years. Among the diversity of features formed by aeolian abrasion on the surface of Mars are periodic bedrock ridges (PBRs), a relatively recently recognized class of erosional bedforms on Mars for which Earth analogues are rare. Gale crater, the field site for NASA’s Mars Science Laboratory Curiosity rover since it landed there in 2012, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of periodic, linear bedrock ridges that occur within the Fe/Mg-smectite clay-bearing Glen Torridon region of Aeolis Mons (informally Mount Sharp). During Curiosity’s exploration of the Glen Torridon region between sols ~2300-3080, the rover drove through this field of ridges, providing the first opportunity for the in situ observation of these enigmatic erosional features. This study characterizes the Glen Torridon ridges using orbiter and rover data to determine their morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, the Glen Torridon ridges are interpreted to be consistent with an origin as wind-eroded periodic bedrock ridges carved during the most recent exhumation of Mount Sharp into the present-day mound. Although there is evidence for multidirectional winds in the Glen Torridon region based on the orientation of modern ripples, megaripples, TARs and other bedrock indicators, the consistent orientation of the Glen Torridon ridges, coupled with morphologic asymmetries within the ridges, support formation and elongation of the Glen Torridon PBRs forms parallel to a net regional northerly wind direction in and around Gale crater.

Journal of Geophysical Research E: Planets↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2012

Three-dimensional modeling of pull-apart basins: implications for the tectonics of the Dead Sea Basin

We model the three-dimensional (3-D) crustal deformation in a deep pull-apart basin as a result of relative plate motion along a transform system and compare the results to the tectonics of the Dead Sea Basin. The brittle upper crust is modeled by a boundary element technique as an elastic block, broken by two en echelon semi-infinite vertical faults. The deformation is caused by a horizontal displacement that is imposed everywhere at the bottom of the block except in a stress-free “shear zone” in the vicinity of the fault zone. The bottom displacement represents the regional relative plate motion. Results show that the basin deformation depends critically on the width of the shear zone and on the amount of overlap between basin-bounding faults. As the width of the shear zone increases, the depth of the basin decreases, the rotation around a vertical axis near the fault tips decreases, and the basin shape (the distribution of subsidence normalized by the maximum subsidence) becomes broader. In contrast, two-dimensional plane stress modeling predicts a basin shape that is independent of the width of the shear zone. Our models also predict full-graben profiles within the overlapped region between bounding faults and half-graben shapes elsewhere. Increasing overlap also decreases uplift near the fault tips and rotation of blocks within the basin. We suggest that the observed structure of the Dead Sea Basin can be described by a 3-D model having a large overlap (more than 30 km) that probably increased as the basin evolved as a result of a stable shear motion that was distributed laterally over 20 to 40 km.

Dead Sea Basin↗