USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Advances in Modeling Earth Systems”

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 19 recordsLinked to original sources

Mean squared error, deconstructed

As science becomes increasingly cross-disciplinary and scientific models become increasingly cross-coupled, standardized practices of model evaluation are more important than ever. For normally distributed data, mean squared error (MSE) is ideal as an objective measure of model performance, but it gives little insight into what aspects of model performance are “good” or “bad.” This apparent weakness has led to a myriad of specialized error metrics, which are sometimes aggregated to form a composite score. Such scores are inherently subjective, however, and while their components may be interpretable, the composite itself is not. We contend that, a better approach to model benchmarking and interpretation is to decompose MSE into interpretable components. To demonstrate the versatility of this approach, we outline some fundamental types of decomposition and apply them to predictions at 1,021 streamgages across the conterminous United States from three streamflow models. Through this demonstration, we hope to show that each component in a decomposition represents a distinct concept, like “season” or “variability,” and that simple decompositions can be combined to represent more complex concepts, like “seasonal variability,” creating an expressive language through which to interrogate models and data.

Journal of Advances in Earth Systems Modeling

The GFDL global atmosphere and land model AM4.0/LM4.0: 2. Model description, sensitivity studies, and tuning strategies

In Part 2 of this two‐part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part 1. Part 2 provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.

Journal of Advances in Modeling Earth Systems

The GFDL global atmosphere and land model AM4.0/LM4.0: 1. Simulation characteristics with prescribed SSTs

In this two‐part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100 km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part 1, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode—with prescribed sea surface temperatures (SSTs) and sea‐ice distribution—is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top‐of‐atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part 2, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.

Journal of Advances in Modeling Earth Systems

SPEAR: The next generation GFDL modeling system for seasonal to multidecadal prediction and projection

We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR ( S eamless System for P rediction and EA rth System R esearch) is built from component models recently developed at GFDL—the AM4 atmosphere model, MOM6 ocean code, LM4 land model, and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0° (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1° to 0.25°. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM4 models but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM4 models.

Journal of Advances in Modeling Earth Systems

Quantitative attribution of major driving forces on soil organic carbon dynamics

Soil organic carbon (SOC) storage plays a major role in the global carbon cycle and is affected by many factors including land use/management changes (e.g., biofuel production-oriented changes). However, the contributions of various factors to SOC changes are not well understood and quantified. This study was designed to investigate the impacts of changing farming practices, initial SOC levels, and biological enhancement of grain production on SOC dynamics and to attribute the relative contributions of major driving forces (CO 2 enrichment and farming practices) using a fractional factorial modeling design. The case study at a crop site in Iowa in the United States demonstrated that the traditional corn-soybean (CS) rotation could still accumulate SOC over this century (from 4.2 to 6.8 kg C/m 2 ) under the current condition; whereas the continuous-corn (CC) system might have a higher SOC sequestration potential than CS. In either case, however, residue removal could reduce the sink potential substantially. Long-term simulation results also suggested that the equilibrium SOC level may vary greatly (∼5.7 to ∼11 kg C/m 2 ) depending on cropping systems and management practices, and projected growth enhancement could make the magnitudes higher (∼7.8 to ∼13 kg C/m 2 ). Importantly, the factorial design analysis indicated that residue management had the most significant impact (contributing 49.4%) on SOC changes, followed by CO 2 Enrichment (37%), Tillage (6.2%), the combination of CO 2 Enrichment-Residue removal (5.8%), and Fertilization (1.6%). In brief, this study is valuable for understanding the major forces driving SOC dynamics of agroecosystems and informative for decision-makers when seeking the enhancement of SOC sequestration potential and sustainability of biofuel production, especially in the Corn Belt region of the United States.

Journal of Advances in Modeling Earth Systems

The GFDL Earth System Model Version 4.1 (GFDL-ESM 4.1): Overall coupled model description and simulation characteristics

We describe the baseline coupled model configuration and simulation characteristics of GFDL's Earth System Model Version 4.1 (ESM4.1), which builds on component and coupled model developments at GFDL over 2013–2018 for coupled carbon-chemistry-climate simulation contributing to the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's CM4.0 development effort that focuses on ocean resolution for physical climate, ESM4.1 focuses on comprehensiveness of Earth system interactions. ESM4.1 features doubled horizontal resolution of both atmosphere (2° to 1°) and ocean (1° to 0.5°) relative to GFDL's previous-generation coupled ESM2-carbon and CM3-chemistry models. ESM4.1 brings together key representational advances in CM4.0 dynamics and physics along with those in aerosols and their precursor emissions, land ecosystem vegetation and canopy competition, and multiday fire; ocean ecological and biogeochemical interactions, comprehensive land-atmosphere-ocean cycling of CO 2 , dust and iron, and interactive ocean-atmosphere nitrogen cycling are described in detail across this volume of JAMES and presented here in terms of the overall coupling and resulting fidelity. ESM4.1 provides much improved fidelity in CO 2 and chemistry over ESM2 and CM3, captures most of CM4.0's baseline simulations characteristics, and notably improves on CM4.0 in (1) Southern Ocean mode and intermediate water ventilation, (2) Southern Ocean aerosols, and (3) reduced spurious ocean heat uptake. ESM4.1 has reduced transient and equilibrium climate sensitivity compared to CM4.0. Fidelity concerns include (1) moderate degradation in sea surface temperature biases, (2) degradation in aerosols in some regions, and (3) strong centennial scale climate modulation by Southern Ocean convection.

Journal of Advances in Modeling Earth Systems (JAM

Development and application of an Infragravity Wave (InWave) driver to simulate nearshore processes

Infragravity waves are key components of the hydro-sedimentary processes in coastal areas, especially during extreme storms. Accurate modeling of coastal erosion and breaching requires consideration of the effects of infragravity waves. Here, we present InWave, a new infragravity wave driver of the Coupled Ocean-Atmopshere-Waves-Sediment Transport (COAWST) modeling system. InWave computes the spatial and temporal variation of wave energy at the wave group scale and the associated incoming bound infragravity wave. Wave group-varying forces drive free infragravity wave growth and propagation within the hydrodynamic model of the coupled modeling system, which is the Regional Ocean Modeling System (ROMS) in this work. Since ROMS is a three-dimensional model, this coupling allows for the combined formation of undertow currents and infragravity waves. We verified the coupled InWave-ROMS with one idealized test case, one laboratory experiment, and one field experiment. The coupled modeling system correctly reproduced the propagation of gravity wave energy with acceptable numerical dissipation. It also captured the transfer of energy from the gravity band to the infragravity band, and within the different infragravity bands in the surf zone, the measured three-dimensional flow structure, and dune morphological evolution satisfactorily. The idealized case demonstrated that the infragravity wave variance depends on the directional resolution and horizontal grid resolution, which are known challenges with the approach taken here. The addition of InWave to COAWST enables novel investigation of nearshore hydro-sedimentary dynamics driven by infragravity waves using the strengths of the other modeling components, namely the three-dimensional nature of ROMS and the sediment transport routines.

Journal of Advances in Modeling Earth Systems

Parameter optimization, sensitivity, and uncertainty analysis of an ecosystem model at a forest flux tower site in the United States

Ecosystem models are useful tools for understanding ecological processes and for sustainable management of resources. In biogeochemical field, numerical models have been widely used for investigating carbon dynamics under global changes from site to regional and global scales. However, it is still challenging to optimize parameters and estimate parameterization uncertainty for complex process-based models such as the Erosion Deposition Carbon Model (EDCM), a modified version of CENTURY, that consider carbon, water, and nutrient cycles of ecosystems. This study was designed to conduct the parameter identifiability, optimization, sensitivity, and uncertainty analysis of EDCM using our developed EDCM-Auto, which incorporated a comprehensive R package—Flexible Modeling Framework (FME) and the Shuffled Complex Evolution (SCE) algorithm. Using a forest flux tower site as a case study, we implemented a comprehensive modeling analysis involving nine parameters and four target variables (carbon and water fluxes) with their corresponding measurements based on the eddy covariance technique. The local sensitivity analysis shows that the plant production-related parameters (e.g., PPDF1 and PRDX) are most sensitive to the model cost function. Both SCE and FME are comparable and performed well in deriving the optimal parameter set with satisfactory simulations of target variables. Global sensitivity and uncertainty analysis indicate that the parameter uncertainty and the resulting output uncertainty can be quantified, and that the magnitude of parameter-uncertainty effects depends on variables and seasons. This study also demonstrates that using the cutting-edge R functions such as FME can be feasible and attractive for conducting comprehensive parameter analysis for ecosystem modeling.

Massachusetts

Subsurface redox interactions regulate ebullitive methane flux in heterogeneous Mississippi River deltaic wetland

As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive-transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient-driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

Journal of Advances in Modeling Earth Systems

Modular compositional learning improves 1D hydrodynamic lake model performance by merging process-based modeling with deep learning

Hybrid Knowledge-Guided Machine Learning (KGML) models, which are deep learning models that utilize scientific theory and process-based model simulations, have shown improved performance over their process-based counterparts for the simulation of water temperature and hydrodynamics. We highlight the modular compositional learning (MCL) methodology as a novel design choice for the development of hybrid KGML models in which the model is decomposed into modular sub-components that can be process-based models and/or deep learning models. We develop a hybrid MCL model that integrates a deep learning model into a modularized, process-based model. To achieve this, we first train individual deep learning models with the output of the process-based models. In a second step, we fine-tune one deep learning model with observed field data. In this study, we replaced process-based calculations of vertical diffusive transport with deep learning. Finally, this fine-tuned deep learning model is integrated into the process-based model, creating the hybrid MCL model with improved overall projections for water temperature dynamics compared to the original process-based model. We further compare the performance of the hybrid MCL model with the process-based model and two alternative deep learning models and highlight how the hybrid MCL model has the best performance for projecting water temperature, Schmidt stability, buoyancy frequency, and depths of different isotherms. Modular compositional learning can be applied to existing modularized, process-based model structures to make the projections more robust and improve model performance by letting deep learning estimate uncertain process calculations.

Journal of Advances in Modeling Earth Systems

Using open-science workflow tools to produce SCEC CyberShake physics-based probabilistic seismic hazard models

The Statewide (formerly Southern) California Earthquake Center (SCEC) conducts multidisciplinary earthquake system science research that aims to develop predictive models of earthquake processes, and to produce accurate seismic hazard information that can improve societal preparedness and resiliency to earthquake hazards. As part of this program, SCEC has developed the CyberShake platform, which calculates physics-based probabilistic seismic hazard analysis (PSHA) models for regions with high-quality seismic velocity and fault models. The CyberShake platform implements a sophisticated computational workflow that includes over 15 individual codes written by 6 developers. These codes are heterogeneous, ranging from short-running high-throughput serial CPU codes to large, long-running, parallel GPU codes. Additionally, CyberShake simulation campaigns are computationally extensive, typically producing tens of terabytes of meaningful scientific data and metadata over several months of around-the-clock execution on leadership-class supercomputers. To meet the needs of the CyberShake platform, we have developed an extreme-scale workflow stack, including the Pegasus Workflow Management System, HTCondor, Globus, and custom tools. We present this workflow software stack and identify how the CyberShake platform and supporting tools enable us to meet a variety of challenges that come with large-scale simulations, such as automated remote job submission, data management, and verification and validation. This platform enabled us to perform our most recent simulation campaign, CyberShake Study 22.12, from December 2022 to April 2023. During this time, our workflow tools executed approximately 32,000 jobs, and used up to 73% of the Summit system at Oak Ridge Leadership Computing Facility. Our workflow tools managed about 2.5 PB of total temporary and output data, and automatically staged 19 million output files totaling 74 TB back to archival storage on the University of Southern California's Center for Advanced Research Computing systems, including file-based relational data and large binary files to efficiently store millions of simulated seismograms. CyberShake extreme-scale workflows have generated simulation-based probabilistic seismic hazard models that are being used by seismological, engineering, and governmental communities.

Frontiers Earth Science Journal

Ultraviolet and visible remote sensing of volcanic gases

As magma rises in volcanic systems, volatile species exsolve from the silicate melt and are emitted as gases into the atmosphere. Measuring the magnitude and composition of gas emissions from volcanoes provides insights into processes occurring deep within the Earth and helps constrain the impact of volcanic degassing on atmospheric chemistry. Optical remote sensing techniques allow volcanic gas emissions to be characterized without the need to access hazardous areas near active volcanic vents. This paper reviews the state of the art in ultraviolet and visible volcanic gas remote sensing from the ground, air, and space. Special attention is given to discussing the physics of atmospheric radiative transfer on which these techniques are based. The functionality and limitations of different remote sensing instruments are examined, making clear that the ideal choice of instrumentation will depend on the volcanic system to which it is applied and the sought measurement parameters. Common algorithms for determining trace gas column densities, gas burdens, and volcanic emission rates from measurements of spectral radiance are outlined and compared, showing how some algorithms attempt to model the physics of the measurement while others maximize sensitivity. Several examples demonstrate how remote sensing measurements continue to advance our understanding of volcanic systems and their impact on the atmosphere. Finally, a few promising directions of inquiry are suggested that could lead to improvements in remote sensing instrumentation and analysis techniques. By combining spectroscopic and imaging techniques, improving our understanding of atmospheric radiative transfer, expanding the suite of target gases, and increasing the coverage and frequency of observations, we stand to significantly improve our ability to detect and quantify volcanic gas emissions and gain new insights into important Earth-system processes.

Journal of Volcanology and Geothermal Research

Lessons from (triggered) tremor

I test a “clock-advance” model that implies triggered tremor is ambient tremor that occurs at a sped-up rate as a result of loading from passing seismic waves. This proposed model predicts that triggering probability is proportional to the product of the ambient tremor rate and a function describing the efficacy of the triggering wave to initiate a tremor event. Using data mostly from Cascadia, I have compared qualitatively a suite of teleseismic waves that did and did not trigger tremor with ambient tremor rates. Many of the observations are consistent with the model if the efficacy of the triggering wave depends on wave amplitude. One triggered tremor observation clearly violates the clock-advance model. The model prediction that larger triggering waves result in larger triggered tremor signals also appears inconsistent with the measurements. I conclude that the tremor source process is a more complex system than that described by the clock-advance model predictions tested. Results of this and previous studies also demonstrate that (1) conditions suitable for tremor generation exist in many tectonic environments, but, within each, only occur at particular spots whose locations change with time; (2) any fluid flow must be restricted to less than a meter; (3) the degree to which delayed failure and secondary triggering occurs is likely insignificant; and 4) both shear and dilatational deformations may trigger tremor. Triggered and ambient tremor rates correlate more strongly with stress than stressing rate, suggesting tremor sources result from time-dependent weakening processes rather than simple Coulomb failure.

Cascadia

Characterization of errors in a coupled snow hydrology-microwave emission model

Traditional approaches to the direct estimation of snow properties from passive microwave remote sensing have been plagued by limitations such as the tendency of estimates to saturate for moderately deep snowpacks and the effects of mixed land cover within remotely sensed pixels. An alternative approach is to assimilate satellite microwave emission observations directly, which requires embedding an accurate microwave emissions model into a hydrologic prediction scheme, as well as quantitative information of model and observation errors. In this study a coupled snow hydrology [Variable Infiltration Capacity (VIC)] and microwave emission [Dense Media Radiative Transfer (DMRT)] model are evaluated using multiscale brightness temperature (TB) measurements from the Cold Land Processes Experiment (CLPX). The ability of VIC to reproduce snowpack properties is shown with the use of snow pit measurements, while TB model predictions are evaluated through comparison with Ground-Based Microwave Radiometer (GBMR), air-craft [Polarimetric Scanning Radiometer (PSR)], and satellite [Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E)] TB measurements. Limitations of the model at the point scale were not as evident when comparing areal estimates. The coupled model was able to reproduce the TB spatial patterns observed by PSR in two of three sites. However, this was mostly due to the presence of relatively dense forest cover. An interesting result occurs when examining the spatial scaling behavior of the higher-resolution errors; the satellite-scale error is well approximated by the mode of the (spatial) histogram of errors at the smaller scale. In addition, TB prediction errors were almost invariant when aggregated to the satellite scale, while forest-cover fractions greater than 30% had a significant effect on TB predictions. ?? 2008 American Meteorological Society.

Journal of Hydrometeorology

Can beaches survive climate change?

Anthropogenic climate change is driving sea level rise, leading to numerous impacts on the coastal zone, such as increased coastal flooding, beach erosion, cliff failure, saltwater intrusion in aquifers, and groundwater inundation. Many beaches around the world are currently experiencing chronic erosion as a result of gradual, present-day rates of sea level rise (about 3 mm/year) and human-driven restrictions in sand supply (e.g., harbor dredging and river damming). Accelerated sea level rise threatens to worsen coastal erosion and challenge the very existence of natural beaches throughout the world. Understanding and predicting the rates of sea level rise and coastal erosion depends on integrating data on natural systems with computer simulations. Although many computer modeling approaches are available to simulate shoreline change, few are capable of making reliable long-term predictions needed for full adaption or to enhance resilience. Recent advancements have allowed convincing decadal to centennial-scale predictions of shoreline evolution. For example, along 500 km of the Southern California coast, a new model featuring data assimilation predicts that up to 67% of beaches may completely erode by 2100 without large-scale human interventions. In spite of recent advancements, coastal evolution models must continue to improve in their theoretical framework, quantification of accuracy and uncertainty, computational efficiency, predictive capability, and integration with observed data, in order to meet the scientific and engineering challenges produced by a changing climate.

Journal of Geophysical Research F: Earth Surface

Establishing the Antarctic Dome C community reference standard site towards consistent measurements from Earth observation satellites

Establishing satellite measurement consistency by using common desert sites has become increasingly more important not only for climate change detection but also for quantitative retrievals of geophysical variables in satellite applications. Using the Antarctic Dome C site (75°06′S, 123°21′E, elevation 3.2 km) for satellite radiometric calibration and validation (Cal/Val) is of great interest owing to its unique location and characteristics. The site surface is covered with uniformly distributed permanent snow, and the atmospheric effect is small and relatively constant. In this study, the long-term stability and spectral characteristics of this site are evaluated using well-calibrated satellite instruments such as the Moderate Resolution Imaging Spectroradiometer (MODIS) and Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Preliminary results show that despite a few limitations, the site in general is stable in the long term, the bidirectional reflectance distribution function (BRDF) model works well, and the site is most suitable for the Cal/Val of reflective solar bands in the 0.4–1.0 µm range. It was found that for the past decade, the reflectivity change of the site is within 1.35% at 0.64 µm, and interannual variability is within 2%. The site is able to resolve calibration biases between instruments at a level of ~1%. The usefulness of the site is demonstrated by comparing observations from seven satellite instruments involving four space agencies, including OrbView-2–SeaWiFS, Terra–Aqua MODIS, Earth Observing 1 (EO-1) – Hyperion, Meteorological Operational satellite programme (MetOp) – Advanced Very High Resolution Radiometer (AVHRR), Envisat Medium Resolution Imaging Spectrometer (MERIS) – dvanced Along-Track Scanning Radiometer (AATSR), and Landsat 7 Enhanced Thematic Mapper Plus (ETM+). Dome C is a promising candidate site for climate quality calibration of satellite radiometers towards more consistent satellite measurements, as part of the framework for climate change detection and data quality assurance for the Global Earth Observation System of Systems (GEOSS).

Canadian Journal of Remote Sensing

Advances in a distributed approach for ocean model data interoperability

An infrastructure for earth science data is emerging across the globe based on common data models and web services. As we evolve from custom file formats and web sites to standards-based web services and tools, data is becoming easier to distribute, find and retrieve, leaving more time for science. We describe recent advances that make it easier for ocean model providers to share their data, and for users to search, access, analyze and visualize ocean data using MATLAB® and Python®. These include a technique for modelers to create aggregated, Climate and Forecast (CF) metadata convention datasets from collections of non-standard Network Common Data Form (NetCDF) output files, the capability to remotely access data from CF-1.6-compliant NetCDF files using the Open Geospatial Consortium (OGC) Sensor Observation Service (SOS), a metadata standard for unstructured grid model output (UGRID), and tools that utilize both CF and UGRID standards to allow interoperable data search, browse and access. We use examples from the U.S. Integrated Ocean Observing System (IOOS®) Coastal and Ocean Modeling Testbed, a project in which modelers using both structured and unstructured grid model output needed to share their results, to compare their results with other models, and to compare models with observed data. The same techniques used here for ocean modeling output can be applied to atmospheric and climate model output, remote sensing data, digital terrain and bathymetric data.

Journal of Marine Science and Engineering

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