Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrologic Connectivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Data for: Spatial access and resource limitations control carbon mineralization in soils

This dataset contains data and code used for the paper "Spatial access and resource limitations control carbon mineralization in soils", https://doi.org/10.1016/j.soilbio.2021.108427. Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale water-soluble soil organic matter (SOM) mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble carbon (acetate) and nitrogen (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Soil respiration was measured during the incubation, after which pore water was extracted from the saturated soils and analyzed for water soluble organic carbon concentrations and characterization. Our results showed that carbon (C) amendments increased the cumulative carbon dioxide (CO2) evolved from the soil cores, suggesting that the soils were C-limited. Drought and rewetting increased soil respiration, and there was a greater abundance of complex aromatic molecules in pore waters sampled from these soils. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We had hypothesized that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (pore water samples), drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules, suggesting that fresh inputs stimulated the C-and N-limited microbial populations present deeper in the soil profile. Our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/kaizadp/TES_spatial_access_2021 (DOI: 10.5281/zenodo.5522938).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Editorial: Linking Hydrological and Biogeochemical Processes in Riparian Corridors

The riparian corridor is a key component of the critical zone and an essential component of watershed systems. According to Merriam-Webster Dictionary, the word riparian is derived from the Latin word riparius, meaning “existing alongside a river.” Riparian corridors typically extend from a few meters to hundreds of meters adjacent to a river and are marked by rich biodiversity, vegetation, and intense biogeochemical activity. They act as integrators of watershed processes and constitute the primary pathways for the subsurface geochemical exports from the watershed. Although riparian corridors comprise only 2–10% of a watershed's area, as much as 90–98% of biogeochemical processing in watersheds occurs in this region, thereby affecting the subsurface geochemical exports and downstream river water quality. Indeed, the riparian corridor is a good example of the Pareto principle. This outsized contribution occurs at the interface between aquatic (river) and terrestrial (land) environments, where interactions between hydrologic and biogeochemical processes are intensified. Variations in the river corridor over time can thus also have outsize impacts. Therefore, it is important to understand the hydrological and biogeochemical linkages in riparian corridors to determine water availability and quality for sustainable management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Exploring Hydrologic Model Process Connectivity at the Continental Scale Through an Information Theory Approach

Exploring water fluxes between hydrological model (HM) components is essential to assess and improve model realism. Many classical metrics for HM diagnosis rely solely on streamflow and hence provide limited insights into model performance across processes. This study applies an information theory measure known as “transfer entropy” (TE) to systematically quantify the transfer of information among major HM components. To test and demonstrate the benefits of TE, we use the Framework for Understanding Structural Errors (FUSE) model to mimic and compare four commonly used HM structures, VIC, PRMS, SACRAMENTO, and TOPMODEL, across 671 catchments spanning a variety of hydrologic regimes in the conterminous United States. We explore connections between HM components and catchment landscape characteristics (e.g., climate, topography, soil, and vegetation) and characterize their nonlinear associations using distance correlation and Spearman correlation coefficients. Our results indicate that while the information transferred from precipitation to runoff is similar across model structures (likely as a result of calibration), the information transferred among other components can vary significantly from a FUSE structure to another. We find that aridity, precipitation duration and frequency, snow fraction, mean elevation, forest area, and leaf area index are often significantly associated with TE between the main HM components. Here, we propose that the presence of meaningful nonlinear associations can be used to diagnose process representation in HMs. Our results highlight the necessity to enhance the conventional streamflow-only calibration approach for a more realistic representation of water dynamics in the models.

large-sample hydrology↗
NASA NTRS2020

Lambayeque Water Resources: Assessing Hydrologic Patterns Using NASA Earth Observations to Address Tree Mortality in Peru’s Coastal Mesquite Forests

The mesquite (Prosopis sp.) forests in Northwestern Peru have had a significant increase in tree mortality in the past fifty years. Within this time frame, 17% of the forest extent was lost and the forest saw an average annual declination rate of 0.33% (Ektvedt, Vetaas, & Lundberg, 2012). These habitats support the region's rich biodiversity and play an important role in local community economies. Several hydrologic causes for mesquite mortality have been hypothesized, but local researchers lack spatially comprehensive techniques to address the problem. While in situ research is currently being utilized in an attempt to explain this recent anomaly, landscape-level visualizations through remote sensing have not been produced to find connections between hydrologic trends and the health of Northwestern Peru’s mesquite forests. Climate Hazards Group InfraRed Precipitation with Station data and Global Land Data Assimilation System data were analyzed to assess hydrologic patterns pertaining to precipitation and soil moisture in the Lambayeque region of Peru. These data were then paired with vegetation indices derived from Terra Moderate Resolution Imaging Spectroradiometer (MODIS) and Suomi National Polar-orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) to display how changing hydrologic patterns relate to the health of mesquite trees. These data were then compiled on a monthly basis over the 30-year study period to create a time series product, which was later referenced with background research to find the likely causes of recent forest decline. The results will assist in the Lambayeque Regional Government’s understanding of recent biological declination and will help them design strategies to mitigate forest decline in Lambayeque, Peru and surrounding regions.

Water Resources↗
NASA NTRS2020

Lambayeque Water Resources: Assessing Hydrologic Patterns Using NASA Earth Observations to Address Tree Mortality in Peru’s Coastal Mesquite Forests

The mesquite (Prosopi ssp.) forests in Northwestern Peru have had a significant increase in tree mortality in the past fifty years. Within this time frame, 17% of the forest extent was lost and the forest saw an average annual declination rate of 0.33% (Ektvedt, Vetaas, & Lundberg, 2012). These habitats support the region's rich biodiversity and play an important role in local community economies. Several hydrologic causes for mesquite mortality have been hypothesized, but local researchers lack spatially comprehensive techniques to address the problem. While in situ research is currently being utilized in an attempt to explain this recent anomaly, landscape-level visualizations through remote sensing have not been produced to find connections between hydrologic trends and the health of Northwestern Peru’s mesquite forests. Climate Hazards Group InfraRed Precipitation with Station dataand Global Land Data Assimilation Systemdata were analyzed to assess hydrologic patterns pertaining to precipitation and soil moisture in the Lambayeque region of Peru. These data were then paired with vegetation indices derived from Terra Moderate Resolution Imaging Spectroradiometer (MODIS) and Suomi National Polar-orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) to display how changing hydrologic patterns relate to the health of mesquite trees. These data were then compiled on a monthly basis over the 30-year study period to create a time series product, which was later referenced with background research to find the likely causes of recent forest decline. The results will assist in the Lambayeque Regional Government’s understanding of recent biological declination and will help them design strategies to mitigate forest decline in Lambayeque, Peru and surrounding regions.

Water Resources↗
DOE OSTI2024

Old-Aged groundwater contributes to mountain hillslope hydrologic dynamics

Understanding connectivity between the soil and deeper bedrock groundwater is needed to accurately predict a watershed’s response to perturbation, such as drought. Yet, the bedrock groundwater dynamics in mountainous environments are typically under-constrained and excluded from watershed hydrologic models. Here, we investigate the role of groundwater characterized with decadal and longer water ages on the hydrologic and mass-transport processes within a steep snow-dominated mountain hillslope in the Central Rocky Mountains (USA). We quantify subsurface and surface water mass-balance, groundwater flowpaths, and age distributions using the ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models, which are compared to hydrometric and environmental tracer observations. An ensemble of models with varied soil and hydrogeologic parameters reproduces observed groundwater levels and century-scale mean ages inferred from environmental tracers. The numerical models suggest soil water near the toe of the hillslope contains considerable (>60 % of the mass-flux) contributions from bedrock flowpaths characterized with water ages >10 years. Flowpath connectivity between the deeper bedrock and soil systems is present throughout the year, highlighting the potentially critical role of groundwater with old ages on processes such as evapotranspiration and streamflow generation. The coupled numerical model and groundwater age observations show the bedrock groundwater system influences the hillslope hydrodynamics and should be considered in mountain watershed conceptual and numerical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The hidden influence of terrestrial groundwater on salt marsh function and resilience

Salt marshes are hotspots of nutrient processing and carbon sequestration. So far, studies addressing spatiotemporal variability in and drivers of salt marsh biogeochemical function, carbon storage and resilience have focused on ocean-driven surface hydrologic influences, neglecting effects of terrestrial hydrology through subsurface connections. Here we evaluate drivers of salt marsh redox potential, a proxy for biogeochemical state, through wavelet analyses and information theory using data from seven marshes. The results point to terrestrial groundwater level as a dominant control on redox variability across all sites. Because redox is a key driver of biogeochemical processes, and specifically oxidation of organic matter that sequesters carbon and maintains marsh elevation, these terrestrial influences are critical to understanding marsh function and evolution. Further, the newly identified links between onshore groundwater levels and marsh redox conditions shift the traditional paradigm and suggest that terrestrial hydrology is a primary control on salt marsh carbon sequestration potential and resilience.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Data-driven modeling to enhance municipal water demand estimates in response to dynamic climate conditions

Altered precipitation and temperature patterns from a changing climate will affect supply, demand, and overall municipal water system operations throughout the arid western U.S. While supply forecasts leverage hydrological models to connect climate influences with surface water availability, demand forecasts typically estimate water use independent of climate and other externalities. Stemming from an increased focus on seasonal water demand management, we use the Salt Lake City, Utah municipal water system as a test bed to assess model accuracy versus complexity trade-offs between simple climate-independent econometric-based models and complex climate-sensitive data-driven models to average to extreme wet and dry climate conditions—representative of a new climate normal. Here, the climate-independent model displayed low performance during extreme dry conditions with predictions exceeding 90% and 40% of the observed monthly and seasonal volumetric demands, respectively, which we attribute to insufficient model complexity. The climate-sensitive models displayed greater accuracy in all conditions, with an ordinary least squares model demonstrating a measurable reduction in prediction bias (3.4% vs. -27.3%) and RMSE (74.0 lpcd vs. 294 lpcd) compared to the climate-independent model. The climate-sensitive workflow increased model accuracy and characterized climate-demand interactions, demonstrating a novel tool to enhance water system management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2020

Surface Atmosphere Integrated Field Laboratory

To help improve understanding and modeling of processes connected to mountain hydrology, the Atmospheric Radiation Measurement (ARM) user facility will operate the Surface Atmosphere Integrated Field Laboratory (SAIL) field campaign in central Colorado. From September 2021 to June 2023, an ARM Mobile Facility will collect measurements with a focus on the 300-square-kilometer (116-square-mile) East River Watershed within the Upper Colorado River Basin. SAIL will combine atmospheric measurements from ARM, a U.S. Department of Energy (DOE) scientific user facility, with existing surface and subsurface data from partner research organizations. Researchers will develop detailed measurements of mountainous water-cycle processes as they pertain to the Colorado River, which supplies water for 40 million people in the American West. Through SAIL, researchers from DOE national laboratories, universities, and research centers and agencies will enable an atmosphere-through-bedrock understanding of mountainous water cycles. By observing what is happening above and below ground to influence hydrology at various scales across the East River Watershed, SAIL will provide insights into how Upper Colorado River watersheds interact with the atmosphere to produce water.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

The scale and persistence of soil moisture anomalies as simulated in a global model

Short term variability of climate is intimately connected with soil moisture variability. Soil moisture provides the storage and subsequent return to the atmosphere, through evaporation and transpiration, of precipitation anomalies over land. Global Circulation Model (GCM) simulations enable consistent identification of correlations and dynamical connections between the hydrologic variables, many of which are incompletely observed. One way to facilitate understanding with these increasingly intricate models is to perform sensitivity studies in which a boundary condition or process is prescribed. In this study we will report on a sensitivity study in which a GCM with a sophisticated land surface representation is used to investigate soil moisture variability in the model climate. The simulations to be used in this study were made at R15 resolution (approximately 4.5 deg latitude x 7.5 deg longitude) with prescribed sea surface temperatures (SST) in the GENESIS model (Thompson and Pollard, 1994), which is coupled to a Land Surface Transfer model (LSX) at 2 deg x 2 deg resolution (Pollard and Thompson, 1994). All the results represented here were taken from the monthly averages of the model results. The LSX model accounts for the physical effects of vegetation with two layers specified at each grid point. Vegetation attributes such as leaf area indices, fractional cover, leaf albedos, etc., were taken from the global dataset in Dorman and Sellers (1989). A six-layer soil model extends from the surface to 4.25 m depth. SST's were prescribed in two ten year experiments using monthly SST values with the daily value being interpolated from the nearest two months. In the first experiment monthly climatological values were used, and in the second, the Atmospheric Model Intercomparison Project (AMIP) observed SST's for the years 1979 through 1988 were used (Gates, 1992). Thus, the former experiment gives a measure of the intrinsic model variability, to be compared with that of the latter experiment, which includes month-to-month variability due to ocean forcing.

Fitzjarrald, D.↗
DOE OSTI2021

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

WHONDRS-GUI: a web application for global survey of surface water metabolites

Background The Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) is a consortium that aims to understand complex hydrologic, biogeochemical, and microbial connections within river corridors experiencing perturbations such as dam operations, floods, and droughts. For one ongoing WHONDRS sampling campaign, surface water metabolite and microbiome samples are collected through a global survey to generate knowledge across diverse river corridors. Metabolomics analysis and a suite of geochemical analyses have been performed for collected samples through the Environmental Molecular Sciences Laboratory (EMSL). The obtained knowledge and data package inform mechanistic and data-driven models to enhance predictions of outcomes of hydrologic perturbations and watershed function, one of the most critical components in model-data integration. To support efforts of the multi-domain integration and make the ever-growing data package more accessible for researchers across the world, a Shiny/R Graphical User Interface (GUI) called WHONDRS-GUI was created. Results The web application can be run on any modern web browser without any programming or operational system requirements, thus providing an open, well-structured, discoverable dataset for WHONDRS. Together with a context-aware dynamic user interface, the WHONDRS-GUI has functionality for searching, compiling, integrating, visualizing and exporting different data types that can easily be used by the community. The web application and data package are available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1484811 , which enables users to simultaneously obtain access to the data and code and to subsequently run the web app locally. The WHONDRS-GUI is also available for online use at Shiny Server ( https://xmlin.shinyapps.io/whondrs/ ).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Modeling the Effects of Artificial Drainage on Agriculture-Dominated Watersheds Using a Fully Distributed Integrated Hydrology Model

In agriculture-dominated watersheds where natural drainage is poor, agricultural ditches (narrow engineered channels) and tile drains (perforated pipes) are widely employed to enhance surface and subsurface drainage, respectively. Despite their relatively small scale, these features exert substantial control over the hydro-biogeochemical function of watersheds and their effects need to be represented in the models. We introduce a novel strategy to incorporate the effects of artificial agricultural drainage into a fully distributed basin-scale integrated surface-subsurface hydrology models. In our approach, narrow agriculture ditches for surface drainage are resolved efficiently using ditch-aligned computational meshes that are hydrologically conditioned to ensure connectivity in the stream/ditch network. For tile drainage in the subsurface, we use the physically based Hooghoudt's drainage equation as a subgrid model and route the water drained through tiles to the nearest ditch. Without site-specific calibration, this model reproduced observed streamflow in the Portage River Watershed (>1,000 km 2 ) as recorded by a USGS gauge with good accuracy (normalized KGE = 0.81) and outperformed a calibrated SWAT model (normalized KGE = 0.68). Numerical experiments confirm that artificial drainage reduces surface inundations and effectively controls the water table. At the watershed scale, artificial drainage increases baseflow but has little effect on watershed discharges above the 90th percentile. The strong physical underpinnings and reduced need for calibration allow us to study the impacts of artificial drainage on distributed hydrological response in terms of fluxes and states and provide a platform for investigating watershed-scale nutrient transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Thermo-Hydro-Chemical Simulation of Mid-Ocean Ridge Hydrothermal Systems: Static 2D Models and Effects of Paleo-Seawater Chemistry

Decades of research have resulted in characterization of the ocean floor manifestations of mid-ocean ridge (MOR) hydrothermal systems, yet numerical models accounting for the connections between heat transfer, hydrology and geochemistry have been slow to develop. The Thermo-hydro-chemical code ToughReact can be used to describe the coupled effects of fluid flow, heat transfer, and fluid-rock chemical interactions that occur in MOR systems. We describe the results of 2-dimensional simulations of steady state flow in fractured diabase with mineral-fluid chemical reactions. Basal heating and specified permeability yield maximum temperature of 400°C. Total fluid flux and high fracture flow velocities are in accord with observations. Fluid chemistry, mineralogical changes and 87 Sr/ 86 Sr ratios can be compared to observations to assess and calibrate models. Simulated high temperature fracture fluids have Mg and SO 4 near zero, elevated Ca and 87 Sr/ 86 Sr of about 0.7040. Total alteration is 10%–50% for simple models of spreading. Anhydrite forms mainly near the base of the upwelling zone and results in substantial local fracture porosity reduction. A calibrated model is used to predict how Sr isotopes and other features of altered oceanic crust would be different in the Cretaceous (95 Ma) early Proterozoic (1,800 Ma) and Archean (3,800 Ma), when seawater may have had high Ca and Sr concentrations, lower pH, higher temperature, and lower Na, Mg, and SO 4 . The simulations are offered as a start on what ultimately may require a longer-term community effort to better understand the role of MOR thermo-hydro-chemical systems in Earth evolution.

58 GEOSCIENCES↗