Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗
NASA NTRS1988

Sulfide mineralization: Its role in chemical weathering of Mars

Pyrrhotite-pentlandite assemblages in mafic and ultramafic igneous rocks may have contributed significantly to the chemical weathering reactions that produced degradation products in the Martian regolith. By analogy with terrestrial processes, a model is proposed whereby supergene alteration of these primary Fe-Ni sulfides on Mars has generated secondary sulfides (e.g., pyrite) below the water table and produced acidic groundwater containing high concentrations of dissolved Fe, Ni and sulfate ions. The low pH solutions also initiated weathering reactions of igneous feldspars and ferromagnesian silicates to form clay silicate and ferric oxyhydroxide phases. Near-surface oxidation and hydrolysis of ferric sulfato- and hydroxo-complex ions and sols formed gossans above the water table consisting of poorly crystalline hydrated ferric sulfates (e.g., jarosite), oxides (ferrihydrite, goethite) and silica (opal). Underlying groundwater, now permafrost, contains hydroxo sulfato complexes of Fe, Al, Mg, Ni, etc., which may be stabilized in frozen acidic solutions beneath the surface of Mars. Sublimation of permafrost may replenish colloidal ferric oxides, sulfates and phyllosilicates during dust storms on Mars.

Burns, Roger G.↗
NASA NTRS1991

Interaction of the sea breeze with a river breeze in an area of complex coastal heating

The interaction of the sea-breeze circulation with a river-breeze circulation in an area of complex coastal heating (east coast of Florida) was studied using a 3D finite-element mesoscale model. The model simulations are compared with temperature and wind fields observed on a typical fall day during the Kennedy Space Center Atmospheric Boundary Layer Experiment. The results from numerical experiments designed to isolate the effect of the river breeze indicate that the convergence in the sea-breeze front is suppressed when it passes over the cooler surface of the rivers.

Zhong, Shiyuan↗
NASA NTRS1994

What would we miss if we characterized the Moon and Mars with just planetary meteorites, remote mapping, and robotic landers?

Exploration of the Moon and planets began with telescopic studies of their surfaces, continued with orbiting spacecraft and robotic landers, and will culminate with manned exploration and sample return. For the Moon and Mars we also have accidental samples provided by impacts on their surfaces, the lunar and martian meteorites. How much would we know about the lunar surface if we only had lunar meteorites, orbital spacecraft, and robotic exploration, and not the Apollo and Luna returned samples? What does this imply for Mars? With martian meteorites and data from Mariner, Viking, and the future Pathfinder missions, how much could we learn about Mars? The basis of most of our detailed knowledge about the Moon is the Apollo samples. They provide ground truth for the remote mapping, timescales for lunar processes, and samples from the lunar interior. The Moon is the foundation of planetary science and the basis for our interpretation of the other planets. Mars is similar to the Moon in that impact and volcanism are the dominant processes, but Mars' surface has also been affected by wind and water, and hence has much more complex surface geology. Future geochemical or mineralogical mapping of Mars' surface should be able to tell us whether the dominant rock types of the ancient southern highlands are basaltic, anorthositic, granitic, or something else, but will not be able to tell us the detailed mineralogy, geochemistry, or age. Without many more martian meteorites or returned samples we will not know the diversity of martian rocks, and therefore will be limited in our ability to model martian geological evolution.

Lindstrom, M. M.↗
NASA NTRS1995

Updated Panel-Method Computer Program

Panel code PMARC_12 (Panel Method Ames Research Center, version 12) computes potential-flow fields around complex three-dimensional bodies such as complete aircraft models. Contains several advanced features, including internal mathematical modeling of flow, time-stepping wake model for simulating either steady or unsteady motions, capability for Trefftz computation of drag induced by plane, and capability for computation of off-body and on-body streamlines, and capability of computation of boundary-layer parameters by use of two-dimensional integral boundary-layer method along surface streamlines. Investigators interested in visual representations of phenomena, may want to consider obtaining program GVS (ARC-13361), General visualization System. GVS is Silicon Graphics IRIS program created to support scientific-visualization needs of PMARC_12. GVS available separately from COSMIC. PMARC_12 written in standard FORTRAN 77, with exception of NAMELIST extension used for input.

Ashby, Dale L.↗
NASA NTRS2005

Chimera Grid Tools

Chimera Grid Tools (CGT) is a software package for performing computational fluid dynamics (CFD) analysis utilizing the Chimera-overset-grid method. For modeling flows with viscosity about geometrically complex bodies in relative motion, the Chimera-overset-grid method is among the most computationally cost-effective methods for obtaining accurate aerodynamic results. CGT contains a large collection of tools for generating overset grids, preparing inputs for computer programs that solve equations of flow on the grids, and post-processing of flow-solution data. The tools in CGT include grid editing tools, surface-grid-generation tools, volume-grid-generation tools, utility scripts, configuration scripts, and tools for post-processing (including generation of animated images of flows and calculating forces and moments exerted on affected bodies). One of the tools, denoted OVERGRID, is a graphical user interface (GUI) that serves to visualize the grids and flow solutions and provides central access to many other tools. The GUI facilitates the generation of grids for a new flow-field configuration. Scripts that follow the grid generation process can then be constructed to mostly automate grid generation for similar configurations. CGT is designed for use in conjunction with a computer-aided-design program that provides the geometry description of the bodies, and a flow-solver program.

Chan, William M.↗
NASA NTRS2019

Reduced Order Modeling for Transonic Aeroservoelastic Control Law Development

As aircraft become more flexible, aeroelastic considerations become increasingly important and complex, particularly for transonic flight where nonlinearities in the flow render linear analysis tools less effective. In order to analyze these aeroelastic interactions between the fluid and the structure efficiently, reduced order models (ROMs) are sometimes generated from and used in place of computational fluid dynamics solutions. In this paper, several aerodynamic ROMs are generated and coupled with structural models to form aeroelastic ROMs. The aerodynamic ROMs generated here include the effects of control surface motion. Hence, the aeroelastic ROMs presented here are appropriate for use in aeroservoelastic applications and are intended to be used for aeroservoelastic control law development. These ROMs are used to simulate a number of test cases with and without control surface involvement. Results show that several of the ROMs generated in the paper are able to predict results similar to solutions of higher-order computational methods.

Waite, Josiah M.↗
DOE OSTI2022

Thermo-Mechanical Distortion of Tungsten-Coated Steel During High Heat Flux Testing Using Plasma Arc Lamps

An experimental setup and a test section were designed and fabricated for high heat flux testing (HHFT) of neutron-irradiated specimens using water-wall plasma arc lamps. Because of the radiological considerations and limitations of reactor irradiation, the size of the test articles was limited to disks less than 10 mm in diameter. The specimen was clamped onto an actively cooled block, and clamping allowed the insertion of several thermocouples on the back surface of the specimen through a copper (Cu) block. Five vacuum plasma sprayed tungsten (W)–coated F82H steel specimens were subjected to HHFT. Surface profilometry measurements, which were conducted after HHFT, revealed central bowing of the top W surface. This type of residual distortion occurred for all of the specimens, and the larger the specimens were, the larger was the distortion.In an attempt to understand specimen distortion and address the science questions related to the testing of subsize specimens during HHFT, a simplified thermo-mechanical model was developed. By using a measured temperature in the Cu as an isothermal boundary condition, the model eliminated the need for coupling cooling fluid flow models with stress models, greatly simplifying the analysis. The main variable in the proposed model is hC, i.e., the thermal contact conductance between the F82H and the Cu washer. Inelastic properties, including hardening properties, were considered for F82H steel and Cu. Numerical simulation results demonstrated a buildup of residual deformation during HHFT and a very complex state of stress and deformation during typical heat flux (HF) cycling. Additionally, hoop stress evolution during a high heat flux cycle reveals that F82H at an interface with W would be mainly in compression during HF application and experienced a transition to a tension state during cooldown. Also, specimen distortion evolves during each HF cycle, as the specimen bows downward during HF application and upward during the cooldown period between HF cycles. The final specimen distortion, i.e., upward bowing of the specimen center, was qualitatively predicted for hC values of 4000 to 5000 W/(m 2 ·K). This hC range of values, for which bulging is obtained, is at the lower spectrum of the range of values for hC, consistent with the low thermal contact conductance expected from the unpolished F82H surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

A Geometry Based Infra-Structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis suites (computational fluid dynamics (CFD), structural analysis, heat transfer and etc.) are: (1) Surface Generation -- usually by employing a Computer Assisted Design (CAD) system; (2) Grid Generation -- preparing the volume for the simulation; (3) Flow Solver -- producing the results at the specified operational point; (4) Post-processing Visualization -- interactively attempting to understand the results. For structural analysis, integrated systems can be obtained from a number of commercial vendors. These vendors couple directly to a number of CAD systems and are executed from within the CAD Graphical User Interface (GUI). It should be noted that the structural analysis problem is more tractable than CFD; there are fewer mesh topologies used and the grids are not as fine (this problem space does not have the length scaling issues of fluids). For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. In most cases, the output from a CAD system could go to Initial Graphics Exchange Specification (IGES) or Standard Exchange Program (STEP) files. The output from Grid Generators and Solvers do not really have standards though there are a couple of file formats that can be used for a subset of the gridding (i.e. PLOT3D data formats). The user would have to patch up the data or translate from one format to another to move to the next step. Sometimes this could take days. Specifically the problems with this procedure are:(1) File based -- Information flows from one step to the next via data files with formats specified for that procedure. File standards, when they exist, are wholly inadequate. For example, geometry from CAD systems (transmitted via IGES files) is defined as disjoint surfaces and curves (as well as masses of other information of no interest for the Grid Generator). This is particularly onerous for modern CAD systems based on solid modeling. The part was a proper solid and in the translation to IGES has lost this important characteristic. STEP is another standard for CAD data that exists and supports the concept of a solid. The problem with STEP is that a solid modeling geometry kernel is required to query and manipulate the data within this type of file. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. Adroit multi-block methods are not far behind. This means that a million node steady-state solution can be computed on the order of hours (using current high performance computers) starting from this 'good' geometry. Unfortunately, the geometry usually transmitted from the CAD system is not 'good' in the grid generator sense. The grid generator needs smooth closed solid geometry. It can take a week (or more) of interaction with the CAD output (sometimes by hand) before the process can begin. One way Communication. (3) One-way Communication -- All information travels on from one phase to the next. This makes procedures like node adaptation difficult when attempting to add or move nodes that sit on bounding surfaces (when the actual surface data has been lost after the grid generation phase). Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive. There is also no way to easily deal with this system in a modular manner. One can only replace the grid generator, for example, if the software reads and writes the same files. Instead of the serial approach to analysis as described above, CAPRI takes a geometry centric approach. This makes the actual geometry (not a discretized version) accessible to all phases of the analysis. The connection to the geometry is made through an Application Programming Interface (API) and NOT a file system. This API isolates the top-level applications (grid generators, solvers and visualization components) from the geometry engine. Also this allows the replacement of one geometry kernel with another, without effecting these top-level applications. For example, if UniGraphics is used as the CAD package then Parasolid (UG's own geometry engine) can be used for all geometric queries so that no solid geometry information is lost in a translation. This is much better than STEP because when the data is queried, the same software is executed as used in the CAD system. Therefore, one analyzes the exact part that is in the CAD system. CAPRI uses the same idea as the commercial structural analysis codes but does not specify control. Software components of the CAD system are used, but the analysis suite, not the CAD operator, specifies the control of the software session. This also means that the license issues (may be) minimized and individuals need not have to know how to operate a CAD system in order to run the suite.

Haimes, Robert↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2020

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Imaging interfacial topography with coherent x-ray reflectivity

The use of coherent x-ray reflectivity to image interfaces is illustrated using model calculations for three cases: a thin film, a semi-infinite substrate, and a coherently strained epitaxial thin film. The intensities and phases of the coherent scattering signals are calculated based on a fully atomistic model (where the phases are assumed to have been recovered by a phasing algorithm). Here, we describe the characteristics of the effective density that are obtained by direct inverse Fourier transformation of the complex structure factors, using a window function to define the data range. These effective densities include significant artifacts like oscillations and negative values derived from the finite vertical momentum transfer of the window function. Two approaches are described to extract the surface topography from the effective densities for three-dimensional visualization: (i) multiplying the effective density by a phase factor and (ii) analyzing the data as a hybrid structure factor F(x, Q z ). The explicit surface sensitivity and specificity of measurements is illustrated when the film or substrate Bragg peaks are excluded. These considerations show that images of interfacial topography can be obtained using currently available experimental capabilities with parameters corresponding to epitaxial thin films.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A Bimodal Diagnostic Cloud Fraction Parameterization. Part II: Evaluation and Resolution Sensitivity

A wide range of approaches exists to account for subgrid cloud variability in regional simulations of the atmosphere. This paper addresses the following questions: 1) Is there still benefit in representing subgrid variability of cloud in convection-permitting simulations? 2) What is the sensitivity to the cloud fraction parameterization complexity? 3) Are current cloud fraction parameterizations scale-aware across convection-permitting resolutions? These questions are addressed for regional simulations of a 6-week observation campaign in the U.S. southern Great Plains. Particular attention is given to a new diagnostic cloud fraction scheme with a bimodal subgrid saturation-departure PDF, described in Part I. The model evaluation is performed using ground-based remote sensing synergies, satellite-based retrievals, and surface observations. It is shown that not using a cloud fraction parameterization results in underestimated cloud frequency and water content, even for stratocumulus. The use of a cloud fraction parameterization does not guarantee improved cloud property simulations, however. Diagnostic and prognostic cloud schemes with a symmetric subgrid saturation-departure PDF underestimate cloud fraction and cloud optical thickness, and hence overestimate surface shortwave radiation. These schemes require empirical bias-correction techniques to improve the cloud cover. The new cloud fraction parameterization, introduced in Part I, improves cloud cover, liquid water content, cloud-base height, optical thickness, and surface radiation compared to schemes reliant on a symmetric PDF. Furthermore, cloud parameterizations using turbulence-based, rather than prescribed constant subgrid variances, are shown to be more scale-aware across convection-permitting resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Learning Forecasts of Rare Stratospheric Transitions from Short Simulations

Abstract Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event’s progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology.

Finkel, Justin↗
DOE OSTI2020

BESS-STAIR: a framework to estimate daily, 30m, and all-weather crop evapotranspiration using multi-source satellite data for the US Corn Belt

Abstract. With increasing crop water demands and drought threats, mapping andmonitoring of cropland evapotranspiration (ET) at high spatial and temporalresolutions become increasingly critical for water management andsustainability. However, estimating ET from satellites for precise waterresource management is still challenging due to the limitations in bothexisting ET models and satellite input data. Specifically, the process of ETis complex and difficult to model, and existing satellite remote-sensing datacould not fulfill high resolutions in both space and time. To address theabove two issues, this study presents a new high spatiotemporal resolution ETmapping framework, i.e., BESS-STAIR, which integrates a satellite-drivenwater–carbon–energy coupled biophysical model, BESS (Breathing Earth SystemSimulator), with a generic and fully automated fusion algorithm, STAIR(SaTallite dAta IntegRation). In this framework, STAIR provides daily 30'mmultispectral surface reflectance by fusing Landsat and MODIS satellite datato derive a fine-resolution leaf area index and visible/near-infrared albedo,all of which, along with coarse-resolution meteorological and CO 2 data, are used to drive BESS to estimate gap-free 30 m resolution daily ET.We applied BESS-STAIR from 2000 through 2017 in six areas across the US CornBelt and validated BESS-STAIR ET estimations using flux-tower measurementsover 12 sites (85 site years). Results showed that BESS-STAIR daily ETachieved an overall R2=0.75, with root mean square error RMSE=0.93 mm d -1 and relative error RE =27.9 % when benchmarkedwith the flux measurements. In addition, BESS-STAIR ET estimations capturedthe spatial patterns, seasonal cycles, and interannual dynamics well indifferent sub-regions. The high performance of the BESS-STAIR frameworkprimarily resulted from (1) the implementation of coupled constraints onwater, carbon, and energy in BESS, (2) high-quality daily 30 m data from theSTAIR fusion algorithm, and (3) BESS's applicability under all-skyconditions. BESS-STAIR is calibration-free and has great potentials to be areliable tool for water resource management and precision agricultureapplications for the US Corn Belt and even worldwide given the globalcoverage of its input data.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Composite propellant combustion modeling with a porous plate burner

A burner is designed to model on a large scale (millimeters) the complex vapor phase processes that take place on awkwardly small scales (about a hundred microns) in the combustion of practical AP/ composite propellants at conventional pressures. Binder vapor evolution is modeled with gaseous fuel (ethane in the experiments reported) flow through a porous plate and oxidizer vapor with the flow of a gaseous oxidizer (air and enriched air in the experiments reported) through discrete holes in the porous plate. Measured flame standoff distance and surface temperature variations are consistent with theoretical predictions at atmospheric pressure. Data obtained at several atmospheres are discussed in the light of the current theories of propellant burning that place varying emphasis on the roles of chemical kinetics and fluid dynamic diffusion/mixing in the vapor phase of a burning composite propellant. The potential and proposed future applications of the burner are indicated.

Kumar, R. N.↗
NASA NTRS1984

A fundamental model and efficient inference for SAR ocean imagery

Employing a synthetic aperture radar (SAR) imaging model based on fundamental models of nonlinear hydrodynamics, electromagnetic scattering from a two-scale surface, and SAR imaging of a time-variant scene, the optimal (minimum mean-square error) estimates of the parameters of a sinusoidal, long gravity wave, and the short gravity wave ensemble are found in an efficient recursive form and their performance evaluated, generally by numerical simulation, in a one-dimensional stationary version. An application is made to Seasat-SAR complex imagery.

Harger, R. O.↗
NASA NTRS1987

Weathering of sulfides on Mars

Pyrrhotite-pentlandite assemblages in mafic and ultramafic igneous rocks may have contributed significantly to the chemical weathering reactions that produce degradation products in the Martian regolith. By analogy and terrestrial processes, a model is proposed whereby supergene alteration of these primary Fe-Ni sulfides on Mars has generated secondary sulfides (e.g., pyrite) below the water table and produced acidic groundwater containing high concentrations of dissolved Fe, Ni, and sulfate ions. The low pH solutions also initiated weathering reactions of igneous feldspars and ferromagnesian silicates to form clay silicate and ferric oxyhydroxide phases. Near-surface oxidation and hydrolysis of ferric sulfato-and hydroxo-complex ions and sols formed gossan above the water table consisting of poorly crystalline hydrated ferric sulfates (e.g., jarosite), oxides (ferrihydrite, goethite), and silica (opal). Underlying groundwater, now permafrost contains hydroxo sulfato complexes of Fe, Al, Mg, Ni, which may be stabilized in frozen acidic solutions beneath the surface of Mars. Sublimation of permafrost may replenish colloidal ferric oxides, sulfates, and phyllosilicates during dust storms on Mars.

Burns, Roger G.↗