Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explainable”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2025

Environmental Controls on Water Vapor Deuterium Excess in the Coastal Boundary Layer: An Information Theory Perspective

We use information theory to quantify the environmental controls on water vapor deuterium excess (D-excess) in coastal Southern California from June 2023 through February 2024. Using Shannon entropy, mutual information (MI), and joint mutual information, metrics that capture both linear and nonlinear relationships, we identify the most informative variables and variable combinations governing D-excess across contrasting marine and continental regimes. Relative humidity with respect to sea surface temperature (RHS) is consistently the strongest individual predictor, explaining up to 27% of D-excess variability during marine conditions but only 10% in continental air masses. The Relative humidity(RHS) + sea surface temperature (SST) combination demonstrates synergistic effects, where their joint influence (explaining up to 36% of D-excess variability) exceeds what either variable achieves individually, confirming their coupled influence on deuterium excess. Wind direction complements RHS most effectively during continental conditions. The best three-variable combination (RHS + SST + Planetary Boundary Layer height) explains 38% of D-excess variability in marine air, while no combination exceeds 20% explanatory power during continental periods. Information theory shows that heteroscedasticity in D-excess relationships indicates regime shifts in controlling processes and quantifies fundamental constraints on predictor variables: some environmental factors like surface pressure or water vapor flux contain insufficient information content to explain D-excess variability regardless of their physical relevance. These results highlight the different predictability limits between marine and continental regimes, challenging the adequacy of linear models and providing a rigorous framework for quantifying the information content of isotope-climate relationships with implications for both modern and paleoclimate applications.

information theory↗
DOE OSTI2026

The impact of optical measurement techniques on measured aerosol particle size distributions

Ambient aerosol particle size distributions measured by the Ultra-High Sensitivity Aerosol Spectrometer (UHSAS) at various sites around the world exhibit modes at optical diameters near 600 nm and 850 nm. These modes are not present in concurrent measurements with the Grimm 11-D Optical Particle Counter (OPC). Here, in this study, we argue that these modes result from the optical measurement technique itself, and we explain why they appear in measurements by the UHSAS but not in those by the Grimm OPC. We construct computer models of the UHSAS and Grimm (“digital UHSAS” and “digital Grimm”) and use these to investigate the size distribution that would result from measurements of artificial aerosol particle size distributions that do not contain modes. The appearance of modes for the structureless incoming size distributions sampled by the digital UHSAS is explained by the nonlinear behavior of partial scattering cross sections of uniform spherical particles as a function of their diameter. The absence of modes in the digital Grimm is explained by the coarser size resolution of that instrument. Detailed analysis of the relationship between optical and geometric diameters for uniform spherical particles reveals two important results. First, these diameters generally have different numerical values for the same particle, and second, the relationship is nonlinear; thus, widths of size bins in terms of optical diameter differ from those in geometric diameter. These results explain the modes observed in the ambient size distributions and highlight concerns with attempts to create a merged size distribution by combining measurements from different instruments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Economic losses from extreme weather in the U.S. Gulf Coast region: spatially differential contributions of climate hazard and socioeconomic exposure and vulnerability

Worldwide economic losses from extreme weather events (EWE) have increased over recent decades, with significant geographic heterogeneity in damages. The IPCC defines the risk from EWE as a function of the climate hazard, socioeconomic exposure, and vulnerability. Although these three drivers vary at fine spatial scales, spatial variability largely has been overlooked in assessments of the drivers of economic loss from EWE. Using cluster analysis, we developed a novel socio-climate hazard typology (SCT) that integrates locally defined climate hazard and socioeconomic exposure and social vulnerability typologies. The results identified 838 unique SCT types impacted by EWE across the Gulf Coastal United States during 1981–2010. We regressed the SCT types and their constituent hazard and socioeconomic components against the cumulative economic loss (1981–2010) from EWE for each SCT type. Across the landscape, economic damages of SCT types were determined by unique, spatially explicit combinations of different risk factors, even in explaining the same level of economic loss. For example, multi-billion-dollar damages in the central Gulf Coast and peninsular Florida were explained by different drivers of risk, with damages in the former explained by additive interaction between climate hazard and multiplicative interaction between climate hazard and socioeconomic exposure and vulnerability, and in the latter explained by socioeconomic exposure and vulnerability. These results highlight the need to diagnose additive and multiplicative interactions among drivers of EWE risk in a spatially explicit context.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Above-ground net primary productivity in regenerating seasonally dry tropical forest: Contributions of rainfall, forest age and soil

Identifying factors controlling forest productivity is critical to understanding forest-climate change feedbacks, modelling vegetation dynamics and carbon finance schemes. However, little research has focused on productivity in regenerating tropical forests which are expanding in their fraction of global area have an order of magnitude larger carbon uptake rates relative to older forest. Here We examined above-ground net primary productivity (ANPP) and its components (wood production and litterfall) over 10 years in forest plots that vary in successional age, soil characteristics and species composition using band dendrometers and litterfall traps in regenerating seasonally dry tropical forests in northwestern Costa Rica. We show that the components of ANPP are differentially driven by age and annual rainfall and that local soil variation is important. Total ANPP was explained by a combination of age, annual rainfall and soil variation. Wood production comprised 35% of ANPP on average across sites and years, and was explained by annual rainfall but not forest age. Conversely, litterfall increased with forest age and soil fertility yet was not affected by annual rainfall. In this region, edaphic variability is highly correlated with plant community composition. Thus, variation in ecosystem processes explained by soil may also be partially explained by species composition. These results suggest that future changes in annual rainfall can alter the secondary forest carbon sink, but this effect will be buffered by the litterfall flux which varies little among years. In determining the long-term strength of the secondary forest carbon sink, both rainfall and forest age will be critical variables to track. We also conclude that detailed understanding of local site variation in soils and plant community may be required to accurately predict the impact of changing rainfall on forest carbon uptake. Synthesis. We show that in seasonally dry tropical forest, annual rainfall has a positive relationship with the growth of above-ground woody tissues of trees and that droughts lead to significant reductions in above-ground productivity. These results provide evidence for climate change—carbon cycle feedbacks in the seasonal tropics and highlight the value of longitudinal data on forest regeneration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Microbial community response to a decade of simulated global changes depends on the plant community

Global changes such as increased drought and atmospheric nitrogen deposition perturb both the microbial and plant communities that mediate terrestrial ecosystem functioning. However, few studies consider how microbial responses to global changes may be influenced by interactions with plant communities. To begin to address the role of microbial–plant interactions, we tested the hypothesis that the response of microbial communities to global change depends on the plant community. We characterized bacterial and fungal communities from 395 plant litter samples taken from the Loma Ridge Global Change Experiment, a decade-long global change experiment in Southern California that manipulates rainfall and nitrogen levels across two adjacent ecosystems, a grassland and a coastal sage scrubland. The differences in bacterial and fungal composition between ecosystems paralleled distinctions in plant community composition. In addition to the direct main effects, the global change treatments altered microbial composition in an ecosystem-dependent manner, in support of our hypothesis. The interaction between the drought treatment and ecosystem explained nearly 5% of the variation in bacterial community composition, similar to the variation explained by the ecosystem-independent effects of drought. Unexpectedly, we found that the main effect of drought was approximately four times as strong on bacterial composition as that of nitrogen addition, which did not alter fungal or plant composition. Overall, the findings underscore the importance of considering plant–microbe interactions when considering the transferability of the results of global change experiments across ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗
DOE OSTI2026

Developing a GUI for the Robotic Test Stand

The introduction of this poster explains the technology behind DUNE’s far and near detectors and how passing neutrinos generate electrons that drift into a wire grid. I then explain how 3 ASICs manage signals received from electron interception. Next, the poster states how COLDATA chips are undergoing quality control by a Robotic Test Stand using a state machine. I further explained how earlier tests were done via a command line script and the necessity to implement a user-friendly Graphical User Interface with new features a command line can’t implement. For the implementation section, tools and methods for implementation are listed such as Python, tkinter, and GitHub as well as how multithreading and queue implementation was necessary for GUI functionality. Then, I elaborated on the GUIs new features. Finally, I explain how the GUI will be distributed across multiple institutions and future changes planned for the GUI. Photos of the RTS, far detector cave, diagram of anode assembly plane, COLDATA chips, set up tab, result tab, and legacy command line interface are shown.

Gutierrez Villanueva, Jaziel [DuPage Coll.]↗
NASA NTRS1983

Solar system history as recorded in the Saturnian ring structure

Holberg's analysis of the Voyager Saturn photographs in reflected and transparent light, and occultation data of stars seen through the rings are discussed. A hyperfine structure, with 10,000 ringlets can be explained by the Baxter-Thompson negative diffusion. This gives the ringlets a stability which makes it possible to interpret them as fossils, which originated at cosmogonic times. It is shown that the bulk structure can be explained by the combined cosmogonic shadows of the satellites Mimas, Janus and the Shepherd satellites. This structure originated at the transition from the plasma phase to the planetesimal phase. The shadows are not simple void regions but exhibit a characteristic signature. Parts of the fine structure, explained by Holberg as resonances with satellites, are interpreted as cosmogonic shadow effects. However, there are a number of ringlets which can neither be explained by cosmogonic nor by resonance effects. Analysis of ring data can reconstruct the plasma-planetesimal transition with an accuracy of a few percent.

Alfven, H.↗
NASA NTRS1983

Genesis stories for the millisecond pulsar

Theoretical models proposed to explain the origin of the millisecond pulsar (MP) PSR 1937+214 are reviewed, examining their ability to explain its low surface dipole magnetic field (B), its low birth temperature (less than 10 to the 8th K), the absence of a companion or remnant, and its low velocity perpendicular to the Galactic plane. The models discussed are a single isolated explosion forming a rapidly spinning neutron star, spin-up of a dead pulsar by accretion from a companion, collapse of an accreting spinning white dwarf, and fusion of a tight binary composed of two old neutron stars. Although all of the models have difficulties in explaining one or more of the MP characteristics, the second model is found to be most probable in the light of present knowledge. The lack of a companion is explained by its tidal disruption after it had fed the accreting pre-pulsar for 1 Gyr or more and its mass had decreased to about 0.01 solar mass. Neutron stars accreting in this way have been observed in Galactic-bulge X-ray sources.

Ruderman, M. A.↗
NASA NTRS1983

Solar system history as recorded in the Saturnian ring structure

Holberg's analysis of the Voyager Saturn photographs in reflected and transparent light, and occultation data of stars seen through the rings are discussed. A hyperfine structure with 10,000 ringlets can be explained by the Baxter-Thompson negative diffusion. This gives the ringlets a stability which makes it possible to interpret them as fossils which originated at cosmogonic times. It is shown that the bulk structure can be explained by the combined cosmogonic shadows of the satellites Mimas and Janus and the Shepherd satellites. This structure originated at the transition from the plasma phase to the planetesimal phase. The shadows are not simple void regions but exhibit a characteristic signature. Parts of the fine structure, explained by Holberg as resonances with satellites, are interpreted as cosmogonic shadow effects. However, there are a number of ringlets which can neither be explained by cosmogonic nor by resonance effects. Analysis of ring data can reconstruct the plasma-planetesimal transition with an accuracy of a few percent. Previously announced in STAR as N84-12013

Alfven, H.↗
NASA NTRS1991

Ancient fluvial processes in the equatorial highlands of Mars

Martian highland craters typically lack ejecta deposits, have no noticeable rim, and are flat floored. In addition, crater size frequency distribution curves show that highland craters have depleted populations less than 20 km in diameter. A variety of processes have been suggested to explain these observations including deposition of aeolian or volcanic materials up to the crater rim crests, thermal creep, terrain softening, and mass wasting. However, none of these processes adequately explains both the crater morphology and population distribution. In order to explain both the Martian highland crater morphology and population distribution, a fluvial process is proposed which is capable of removing the loose crater rim material. The resulting effect is to decrease the crater diameter, thereby causing the population curves to bendover. The eroded material is redistributed, burying or partially burying smaller diameter craters before complete erosion. This material may also be deposited into local topographic lows, creating the depositional basins observed. A fluvial process explains both sets of observations: crater morphology and crater population distribution curves.

Craddock, Robert A.↗
NASA NTRS1992

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗
NASA NTRS1993

Helical jets and the misalignment distribution for core-dominated radio sources

The type of a model which would explain the large population of orthogonally misaligned sources is considered. Two of the best examples of misaligned sources are investigated, namely, the quasar 3C 309.1 (Wilkinson et al., 1986; Kus et al., 1990, 1991) and the BL Lac object 1823 + 568 (Foley, 1982; O'Dea et al., 1988). It is shown that simple bend geometries cannot explain the misaligned sources, and the geometries required to give orthogonal misalignments are explained. A generally applicable model which might explain misaligned sources is proposed.

Conway, J. E.↗
NASA NTRS1993

ACARA user's manual

ACARA (Availability, Cost, and Resource Allocation) is a computer program which analyzes system availability, lifecycle cost (LCC), and resupply scheduling using Monte Carlo analysis to simulate component failure and replacement. This manual was written to: (1) explain how to prepare and enter input data for use in ACARA; (2) explain the user interface, menus, input screens, and input tables; (3) explain the algorithms used in the program; and (4) explain each table and chart in the output.

Stalnaker, Dale K.↗
NASA NTRS1994

Coupling of ELF/ULF energy from lightning and MeV particles to the middle atmosphere, inosphere, and global circuit

In an attempt to explain numerous atmospheric electrical phenomena, the elements of the global electrical circuit are reexamined. In addition to being a 'quasi-static 'DC' generator' and source of radiated energy at VLF and higher, the thunderstorm is found to be a pulse generator, with most of the external energy contained in ELF and ULF pulse currents to the ionosphere (and Earth). The pulse energy is found to deposit largely in the middle atmosphere above the thunderstorm. The VLF and above components are well understood, as are the ULF components due to the conductivity gradient. However, a previously poorly understood ELF component on the millsecond timescale, or 'slow tail,' contains a large fraction of the electrical energy. This component couples strongly to the ionosphere and also launches a unipolar transverse electromagnetic (TEM) wavelet in the radial Earth-ionosphere transmission line. The increase in charge with distance associated with such wavelets, and their ensemble sum at a point, may explain some large mesospheric 'DC' fields but there are still difficulties explaining other than rare occurrences, except for antipodal reconvergence. These millisecond duration unipolar wavelets also coupled to the ionosphere and may trigger other lightning at a distance. A schema is elucidated by which the charge of MeV particles deposited in the middle atmosphere persists for much longer than the local relaxation time. This also gives rise to unipolar waves of global extent which may explain lower-latitude field perturbations associated with solar/geomagnetic events.

Hale, Leslie C.↗
NASA NTRS1994

Flat spectrum T Tauri stars: The case for infall

We show that the mid- to far-infrared fluxes of 'flat spectrum' T Tauri stars can be explained by radiative equilibrium emission from infalling dusty envelopes. Infall eliminates the need for accretion disks with non-standard temperature distributions. The simplicity and power of this explanantion indicates that models employing 'active' disks, in which the temperature distribution is a parameterized power law, should be invoked with caution. Infall also naturally explains the scattered light nebulae detected around many flat spectrum sources. To match the observed spectra, material must fall onto a disk rather than the central star, as expected for collapse of a rotating molecular cloud. It may be necessary to invoke cavities in the envelopes to explain the strength of optical and near-infrared emission; these cavities could be produced by the powerful bipolar outflows commonly observed from young stars. If viewed along the cavity, a source may be lightly extincted at visual wavelengths, while still accreting substantial amounts of material from the envelope. Infall may also be needed to explain the infrared-bright companions of many optical T Tauri stars. This picture suggests that many of the flat spectrum sources are 'protostars'-young stellar objects surrounded by dust infalling envelopes of substantial mass.

Calvet, Nuria↗
NASA NTRS1994

Carbonate-sulfate volcanism on Venus?

Venusian canali, outflow channels, and associated volcanic deposits resemble fluvial landforms more than they resmeble volcanic features on Earth and Mars. Some canali have meandering habits and features indicative of channel migration that are very similar to meandering river channels and flood plains on Earth, venusian outflow channels closely resemble water-carved outflow channels on Mars and the Channeled Scabland in Washington, collapsed terrains at the sources of some venusian channels resemble chaotic terrains at the sources of martian outflow channels, venusian lava deltas are similar to bird's-foot deltas such as the Mississippi delta, and venusian valley networks indicate sapping. We have developed an alternative possibility that the lava had a water-like rheology and a melting point slightly greater than Venus' surface temperature, thus accounting for the unusual behavior of the lava. Unlike silicate lavas, some carbonatites (including carbonate-sulfate-rich liquids) have these properties; thus they can flow great distances while retaining a high fluidity, significant mechanical erosiveness, and substantial capacity to transport and deposit sediment. Venusian geochemistry and petrology are consistent with extensive eruptions of carbonatite lavas, which could have crustal and/or mantle origins. Venus' atmosphere (especially CO2, HCl, and HF abundances) and rocks may be in local chemical equilibrium, which suggests that the upper crust contains large amounts of calcite, anhydrite, and other salts. Chemical analyses indicate, according to some models, that Venusian rocks may contain 4-19% calcite and anhydrite. Mixtures of crustal salts could melt at temperatures a few tens to a few hundred Kelvins higher than Venus' surface temperature; hence, melting may be induced by modest endogenetic or impact heating. Salts may have many of the same geologic roles on Venus as water and ice have on Mars. A molten salt (carbonatite) 'aquifer' may exist beneath a few hundred meters to several kilometers of solidified salt-rich 'permafrost.' Many geologic features can be explained by carbonatite magmatism: (1) impact melting of crustal salts can explain crater outflows, (2) small, sustained eruptions from molten salt aquifers can explain sapping valleys, (3) large, sustained eruptions may explain canali and their flood plans, and (4) catastrophic outbursts amy have formed outflow channels and chaotic terrain. Landforms created by carbonate-rich lavas would be thermally stable on Venus' surface, though some minerals may weather to other solid substances.

Kargel, Jeffrey S.↗