Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Clay”

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 793 records · Page 44

Determination of fundamental mechanical properties of biomass using the cubical triaxial tester to model biomass flow

The flowability of biomass is an indicator of how amenable biomass is to handling. The flowability measurements are often determined with empirical or tertiary experiments. This status quo impedes employing engineering principles to advance the design and operation of biomass handling systems. It is imperative to establish an experimental protocol minimizing empirical aspects of flowability characterization to account for the variability of biomass. This study demonstrates the operational principles of the large chamber cubical triaxial tester and the procedure of triaxial tests to quantify bulk flow behaviors of two milled biomass feedstocks, namely corn stover 2 mm (CS) and Douglas fir 1 mm (DF). The Mohr-Coulomb (MC), Drucker-Prager (DP) and modified Cam Clay (mCC) models are calibrated for both biomasses. Analysis indicates that CS will exhibit a cohesive flow with a larger cohesion coefficient of MC (3.8 ± 3.5 kPa) than DF (0.42 ± 0.9 kPa) and a larger d value of DP (6.9 kPa) than DF (0.0 kPa), respectively. CS exhibits a higher spring-back index of mCC (0.39 ± 0.05) than DF (0.27 ± 0.05) also suggesting handling issues, which agrees with the experiences in the industry. This study demonstrates the capability of a CTT in a quantitative investigation of biomass handling characteristics.

09 BIOMASS FUELS↗

Continental-scale patterns of extracellular enzyme activity in the subsoil: an overlooked reservoir of microbial activity

Abstract Chemical stabilization of microbial-derived products such as extracellular enzymes (EE) onto mineral surfaces has gained attention as a possibly important mechanism leading to the persistence of soil organic carbon (SOC). While the controls on EE activities and their stabilization in the surface soil are reasonably well-understood, how these activities change with soil depth and possibly diverge from those at the soil surface due to distinct physical, chemical, and biotic conditions remains unclear. We assessed EE activity to a depth of 1 m (10 cm increments) in 19 soil profiles across the Critical Zone Observatory Network, which represents a wide range of climates, soil orders, and vegetation types. For all EEs, activities per mass of soil correlated positively with microbial biomass (MB) and SOC, and all three of these variables decreased logarithmically with depth ( p < 0.05). Across all sites, over half of the potential EE activities per mass soil consistently occurred below 20 cm for all measured EEs. Activities per unit MB or SOC were substantially higher at depth (soils below 20 cm accounted for 80% of whole-profile EE activity), suggesting an accumulation of stabilized (i.e. mineral sorbed) EEs in subsoil horizons. The pronounced enzyme stabilization in subsurface horizons was corroborated by mixed-effects models that showed a significant, positive relationship between clay concentration and MB-normalized EE activities in the subsoil. Furthermore, the negative relationships between soil C, N, and P and C-, N-, and P-acquiring EEs found in the surface soil decoupled below 20 cm, which could have also been caused by EE stabilization. This finding suggests that EEs may not reflect soil nutrient availabilities deeper in the soil profile. Taken together, our results suggest that deeper soil horizons hold a significant reservoir of EEs, and that the controls of subsoil EEs differ from their surface soil counterparts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimating perturbed stress from 3-D borehole displacements induced by fluid injection in fractured or faulted shales

SUMMARY Hydrofracturing stress measurements in fractured and anisotropic shales are notoriously difficult, because opening of existing geological features tends to prevent the creation of a pure hydraulic fracture perpendicular to the least compressive principal stress. Here we show how adding 3-D borehole-displacement measurements while conducting the hydraulic injection test helps to better constrain the principal stress orientations and magnitudes. We developed a 3-D fully coupled hydromechanical numerical model to analyse the displacement, fluid pressure and injection flow-rate data measured during an injection pressure-step-rate test conducted to activate a faulted borehole interval in the Mont Terri Opalinus Clay (Switzerland). We find that injected fluids can only penetrate the fault when it is at or above the Coulomb failure pressure. Borehole displacement orientations are sensitive to a ∼15° variation in the stress–tensor orientation and a 1 MPa stress magnitude variation. Although some dispersion occurs while rupture is propagating along the fault plane ∼4 m away from the borehole, the maximum density of displacement orientations consistently informs about the stress orientation. Thus, an extended injection step-rate approach coupled with an accurate in situ measurement of the borehole wall displacements can be used to better constrain the local stress field perturbations in fractured shales and in heterogeneous rock in general.

Guglielmi, Yves↗

K in an Urban World: New Contexts for Hydraulic Conductivity

Abstract Hydraulic conductivity ( K ) is a key hydrologic parameter widely recognized to be difficult to estimate and constrain, with little consistent assessment in disturbed, urbanized soils. To estimate K , it is either measured, or simulated by pedotransfer functions, which relate K to easily measured soil properties. We measured K in urbanized soils by double‐ring infiltrometer ( K dring ), near‐saturated tension infiltrometry ( K minidisk ), and constant head borehole permeametry ( K borehole ), along with other soil properties across the major soil orders in 12 United States cities. We compared measured K with that predicted from the pedotransfer function, ROSETTA. We found that regardless of soil texture, K dring was consistently larger than K minidisk ; with the latter having slightly less sample variance. K borehole was dependent upon specific subsurface conditions, and contrary to common expectations, did not always decrease with depth. Based on either soil textural class, or percent textural separates (sand, silt clay), ROSETTA did not accurately predict measured K for surface nor subsurface soils. We go on to discuss how K varies in urban landscapes, the role of measurement methods and artifacts in the perception of this metric, and implications for hydrologic modeling. Overall, we aim to inspire consistency and coherence when addressing K ‐related challenges in sustainable urban water management.

Shuster, W. D.↗

Comparing Experimental Results for Large Particle Separation from Non-Newtonian Slurries Using Full and Tapered Bump Arrays

To separate particles from carrier fluids arrays of staggered posts inserted into the flow stream can be employed. The smaller particles follow the flow stream while the larger particles move laterally to one side. During flow, the concentration of particulate increases in one direction and decreases in the other crosswise direction. This technique separates the particles from the carrier fluid. Experiments using a slurry with non-Newtonian rheology were conducted to evaluate large particle separation from the non-Newtonian carrier fluid using bump arrays. The bentonite kaolin clay slurry with non-Newtonian rheology was spiked with large diameter inert glass particles. Experiments with flow through lateral displacement arrays evaluated the performance of particle separation using two array configurations: a full array with all posts and a tapered array with posts removed in a triangular portion of the array past the rows of anticipated particle of separation. The posts are staggered posts to promote particle segregation to one side of the flow channel. Using this type of array for removing particles from non-Newtonian yield stress slurries is novel and we present unique results. In this paper, the performance of full and tapered arrays are compared. These results have many industrial applications including removing particles from slurries of nuclear waste. Large particle removal is an important step in waste processing.

slurry, solids removal, deterministic lateral disp↗

Biogenic-to-lithogenic handoff of particulate Zn affects the Zn cycle in the Southern Ocean

Zinc (Zn) is vital to marine organisms. Its active uptake by phytoplankton results in a substantial depletion of dissolved Zn, and Zn bound to particulate organic matter replenishes dissolved Zn in the ocean through remineralization. However, we found that particulate Zn changes from Zn bound to phosphoryls in cells to recalcitrant inorganic pools that include biogenic silica, clays, and iron, manganese, and aluminum oxides in the Southern Ocean water column. The abundances of inorganic pools increase with depth and are the only phases preserved in sediments. Changes in the particulate-Zn speciation influence Zn bioavailability and explain the decoupling of Zn and phosphorus and the correlation of Zn and silicon in the water column. Here, these findings reveal a new dimension to the ocean Zn cycle, implicating an underappreciated role of inorganic Zn particles and their impact on biological productivity.

Duan, J.↗

Genomic fingerprints of the world’s soil ecosystems

Despite the explosion of soil metagenomic data, we lack a synthesized understanding of patterns in the distribution and functions of soil microorganisms. These patterns are critical to predictions of soil microbiome responses to climate change and resulting feedbacks that regulate greenhouse gas release from soils. To address this gap, we assay 1,512 manually curated soil metagenomes using complementary annotation databases, read-based taxonomy, and machine learning to extract multidimensional genomic fingerprints of global soil microbiomes. Our objective is to uncover novel biogeographical patterns of soil microbiomes across environmental factors and ecological biomes with high molecular resolution. We reveal shifts in the potential for (i) microbial nutrient acquisition across pH gradients; (ii) stress-, transport-, and redox-based processes across changes in soil bulk density; and (iii) greenhouse gas emissions across biomes. We also use an unsupervised approach to reveal a collection of soils with distinct genomic signatures, characterized by coordinated changes in soil organic carbon, nitrogen, and cation exchange capacity and in bulk density and clay content that may ultimately reflect soil environments with high microbial activity. Genomic fingerprints for these soils highlight the importance of resource scavenging, plant-microbe interactions, fungi, and heterotrophic metabolisms. Across all analyses, we observed phylogenetic coherence in soil microbiomes—more closely related microorganisms tended to move congruently in response to soil factors. Collectively, the genomic fingerprints uncovered here present a basis for global patterns in the microbial mechanisms underlying soil biogeochemistry and help beget tractable microbial reaction networks for incorporation into process-based models of soil carbon and nutrient cycling.

59 BASIC BIOLOGICAL SCIENCES↗

Wide-blocky veins explained by dependency of crystal growth rate on fracture surface type: Insights from phase-field modeling

Vein microstructures contain a wealth of information on coupled chemical and mechanical processes of fracturing, fluid transport, and crystal growth. Numerical simulations have been used for exploring the factors controlling the development of vein microstructures; however, they have not been quantitatively validated against natural veins. Here we combined phase-field modeling with microtextural analysis of previously unexplained wide-blocky calcite veins in natural limestone and of the fresh fracture surface in this limestone. Results show that the wide-blocky vein textures can only be reproduced if ~10%–20% of crystals grow faster than the rest. This fraction corresponds to the amount of transgranularly broken grains that were observed on the experimental fracture surfaces, which are dominantly intergranular. We hypothesize that transgranular fractures allow faster growth of vein minerals due to the lack of clay coatings and other nucleation discontinuities that are common along intergranular cracks. Our simulation results show remarkable similarity to the natural veins and reproduce the nonlinear relationship between vein crystal width and vein aperture. This allows accurate simulations of crystal growth processes and related permeability evolution in fractured rocks.

Geology↗

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics↗

Validation of Play Fairway Analysis of the geothermal potential of Camas Prairie, south-central Idaho, by an exploration well

Abstract Play Fairway Analysis (PFA) methodology was adapted for geothermal exploration at Camas Prairie, Idaho. Geophysical data, structural and geologic mapping, volcanic rock ages and vent locations, and the distribution of thermal springs and wells all indicated a relatively high geothermal potential along the southern margin of the Prairie. An exploration well (USU Camas-1) was drilled to a depth of 618.3 m to validate the PFA. A permeable zone was encountered at ~ 357.5 m with a maximum measured temperature of ~ 80 °C, which was suppressed following the injection of cold water. A moderate transmissivity of ~ 0.25–1 cm 2 /s estimated from an injection test as well a seasonal artesian flow at ~ 0.7 L/s corroborate the presence of a permeable zone. The existence of a lacustrine clay seal was confirmed near the bottom of the basin-fill sediment occupying the upper 314 m of the well. Geothermometers suggest the USU Camas-1 well water equilibrated at a reservoir temperature of ~ 120 °C. Based on the locations of both thermal and cold wells, geothermal fluids appear to be flowing upward along one or both of two fault systems. The presence of young basalts and elevated helium isotope ratios suggest that the heat source of Camas Prairie is magmatic. However, the faults may be acting as a conduit for geothermal fluids to rise from great depth without a shallow magmatic source being present. Camas Prairie is a promising area for geothermal development, but the relatively low reservoir temperatures indicate this resource may not be suitable for electric generation. Perhaps the best use would be for heating.

Lachmar, Thomas E. (ORCID:0000000249640119)↗

Acoustic velocity and permeability of acidized and propped fractures in shale

We report that from geochemical reactions to proppant emplacement, hydraulic fracturing induces various chemomechanical fracture alterations in shale reservoirs. Hydraulic fracturing through the injection of a vast amount and variety of fluids and proppants has substantial impacts on fluid flow and hydrocarbon production. There is a strong need to improve our understanding on how fracture alterations affect flow pathways within the stimulated rock volume and develop monitoring tools. We have conducted time-lapse rock-physics experiments on clay-rich (carbonate-poor) Marcellus shales to characterize the acoustic velocity and permeability responses to fracture acidizing and propping. Acoustic P- and S-wave velocities and fracture permeability were measured before and after laboratory-induced fracture alterations along with microstructural imaging through X-ray computed tomography and scanning electron microscopy. Our experiments indicate that the S-wave velocity is an important geophysical observable, particularly the S-wave polarized perpendicular to fractures because it is sensitive to fracture stiffness. The acidizing and propping of a fracture decrease its elastic stiffness. This effect is stronger for acidizing, so it is possible that proppant monitoring will be masked by chemical alteration except when propping is highly efficient (i.e., most fractures are propped). However, fracture permeability is undermined by the softening of fracture surfaces due to acidizing, while being greatly enhanced by propping. These contrasting effects on fluid flow in combination with similar seismic attributes indicate the importance of experiments to improve existing rock-physics models, which must include changes to the rock frame. Such improvements are necessary for a correct interpretation of seismic velocity monitoring of flow pathways in stimulated shales.

58 GEOSCIENCES↗

Predicting variations of the least principal stress with depth: Application to unconventional oil and gas reservoirs using a log-based viscoelastic stress relaxation model

Knowledge of layer-to-layer variations of the least principal stress, S hmin , with depth is essential for optimization of multi-stage hydraulic fracturing in unconventional reservoirs. Utilizing a geomechanical model based on viscoelastic stress relaxation in relatively clay rich rocks, we present a new method for predicting continuous S hmin variations with depth. The method utilizes geophysical log data and S hmin measurements from routine diagnostic fracture injection tests (DFITs) at several depths for calibration. We consider a case study in the Wolfcamp formation in the Midland Basin, where both geophysical logs and values of S hmin from DFITs are available. We compute a continuous stress profile as a function of the well logs that fits all of the DFITs well. We utilized several machine learning technologies, such as bootstrap aggregation (or bagging), to improve the generalization of the model and demonstrate that the excellent fit between predicted and observed stress values is not the result of over-fitting the calibration points. The model is then validated by accurately predicting hold-out stress measurements from four wells within the study area and, without recalibration, accurately predicting stress as a function of depth in an offset pad about 6 miles away.

58 GEOSCIENCES↗

<i>Finite Element Modeling of Biomass Hopper Flow</i>

The hopper is one of the widely used biomass handling devices that channels the bulk biomass from storage to the subsequent handling systems. Although Jenike's approach has been successfully used for hoppers handling grains and other agricultural produces for decades, designing a hopper, ensuring reliable biomass handling is found to be challenging. This study aims to address this engineering problem with alternative constitutive material models concerning the flow behavior of bulk solids. Finite element modeling is an approach that allows for implementing different material models. Underlying constitutive theories of different material models assist in investigating the mechanical behavior of a particulate system, including the flow characteristics of ground biomass. This study demonstrates hopper flow models of two types of biomass, i.e., ground corn stover and Douglas fir, with the Mohr-Coulomb model, which is based on Jenike's approach, modified Cam-Clay model, and Drucker-Prager/Cap model. The modeling results are compared with hopper flow experiments to highlights the advantages and shortfalls of each constitutive biomass flow model.

09 BIOMASS FUELS↗

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

This is a final technical report for the project: Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems (EGS). The goal is to demonstrate the feasibility of foam fracturing in EGS applications. The project, led by Oak Ridge National Laboratory (ORNL), was conducted in collaboration with Temple University. The report describes the research activities with Task 1 at ORNL: foam fracturing testing system development and experimental study on foam fracturing, and Task 2 at Temple University: foam testing and foam characterization. Main findings are: 1. A foam fracturing test system has been developed at ORNL, which can be used to perform foam fracturing under pressure up to 6,000 psi. The system monitors foam density during fracturing online and is capable of testing materials in both monotonic and cyclic (up to 50 Hz) injections. 2. Foam fracturing tests were carried out on Charcoal black granite specimens with a blind borehole to the middle length. Two diameters of blind borehole were tested; G2 series: 9.53 mm and G3 series: 4.76 mm. N2-in-water foam was used with AOS as a surfactant. 3. There was a hole-size effect on fracture initiation pressure. The effect is smaller in the case of foam, which was influenced by the high penetrability of gas in foam. Breakdown pressure showed a behavior just as that of fracture pressure; namely an increased value for small hole samples, while the effect in water fracture was more impressive than in foam fracture. 4. Water mass was reduced in foam fracturing within similar range of breakdown pressures. In G2 series, it was decreased from 10.44 g for water fracturing to 5.17 g, representing more than 50% water reduction. Therefore, there is the potential to reduce water use in EGS stimulation through foam fracturing. 5. Use of cyclic injection has the potential to reduce the breakdown pressure and seismicity in EGS application. Experiments using 4-s cycle period found that specimens can be fractured with a low number of cycles. The fatigue pressure was approximately 64 - 77% of monotonic breakdown pressure for water fracturing and 58 - 94% of the breakdown pressure for foam fracturing. 6. A foam stability testing system has been developed that can test foam at 220 Deg C to 2,000 psi. Tested components of candidate foams included two gases: N2 and CO2; 4 surfactants: AOS, SDS, NP-40 and CTAC; 5 stabilizing agents: guar, bentonite clay, borate salt, silica NPs, and GO. 7. N2 and AOS provided the most stable performance over the tested ranges. Furthermore, the AOS foam with stabilizing agents of guar and borate salt (crosslinker) offered the highest half-life of 20 minutes at 200 Deg C and 1,000 psi. 8. Arrhenius equation and modified power law have been demonstrated to fit well the half-time vs. temperature and pressure data, respectively. These relations can be useful to provide the suggestion for future foam stability study. This submission contains the supporting data developed during the project: 1) A final technical report 2) Granite fracturing data in monotonic and cyclic injections with water and N2 foam Foam performance data in various temperatures and pressures, including half-time, is submitted separately.

15 GEOTHERMAL ENERGY↗

Experimental Images and Videos of Foam Stability (Half-life)

The experimental data obtained in this project is the thermal stability data of various foams measured using the setup established at Temple University during this study. The setup is installed with a portable digital camera which can take images and videos of foam evolution at a given pressure and temperature condition. Consequently, the half-life data was recorded from the images/videos, which are used as a measure of the thermal stability for foams. Over the 3 years of this project, four different surfactants and five different stabilizing agents were studied. The surfactants are, Alfa Olefin Sulfonate (AOS), Sodium Dodecyl Sulphate (SDS), Tergitol (NP-40), and Cetyltrimethylammonium chloride (CTAC). The stabilizing agents are, guar gum, bentonite clay, crosslinking agents, silicon dioxide nanoparticles (60 to 70nm), and graphene oxide dispersions. Foam stability was evaluated at different temperatures between 100C and 200cC, while the foam generation pressure varied between atmospheric pressure (14.7 psi) and 1000 psi. The images are saved as .jpg file and videos are saved as .avi files.

15 GEOTHERMAL ENERGY↗

Wetting behaviors of fluoroterpolymer fiber films

Various aspects of electrospun fibers prepared from terpolymer of tetrafluoroethylene (TFE), hexafluoropropylene (HFP), and vinylidene fluoride (VDF) (THV)/acetone solutions at two applied voltages, THV/acetone solutions having Texas montmorillonite with two ratios, and THV/ethyl acetate solutions using two needle sizes are described. Fibers from THV/acetone and THV/ethyl acetate solutions showed shallow indentations and pores, respectively. The clay, functioning as electrospinning agent, did not influence the fiber morphology, but yielded narrower fiber diameter distribution and the thinnest fibers. Heterogeneous fiber diameter distribution and increase in the fiber diameters were observed by lowering the voltage for fibers of THV/acetone solutions. Fibers from THV/ethyl acetate solutions had the largest diameter and the broadest diameter distribution. Electrospun THV fibers having both hydrophobic characteristics with nearly 140° water contact angles and oleophilic properties with oil contact angles less than 45° might have applications in areas such as water/oil separation.

02 PETROLEUM↗

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗