Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface parameters”

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 235 records · Page 13

Automated Redox Titrations via Interdigitated Electrode Arrays: Application to the Mediated Electron Transfer Interrogation of Charge and Rate on Electrodeposited Polymers

Mediated electron transfer (MET) plays a crucial role in energy storage and conversion technologies such as redox targeting flow batteries (RTFBs), yet its experimental investigation often requires labor-intensive and low-throughput setups. To address this, we developed a microfabricated interdigitated electrode array (IDA) platform that enables automated, high-throughput electrochemical redox titration measurement to be performed to study the MET process. Our redox titration method enables simultaneous measurement of the charge capacity and rate of MET processes on a material or surface. Automated redox titration (ART) facilitates systematic investigation of the MET process across a broad parameter space, exemplified through the study of polypyrrole (PPy) and a pyrene-4,5,9,10-tetrone azo group-based polymer (PTAP), both redox-active polymers relevant to various energy storage applications. Using PPy as a model material, 500 redox titration measurements were conducted within 50 h, varying the electrode gap widths, polymer charging potentials, voltammetric scan rates, and electrolyte concentrations. Finite-element simulations confirmed the electrochemical responses and elucidated the kinetics of the MET reactions. Our automated methodology was further tested with PTAP, revealing a surprising charging potential dependence on the rate of MET. The automation, flexibility, and scalability of our redox titration platform pave the way not only for advanced studies of MET processes relevant to RTFBs, but also with implications in the understanding of next-generation energy storage materials, molecular electrocatalysis, and biosensing.

electrochemical analysis↗

Laser Powder Bed Fusion Additive Manufacture Nb1Zr Development

Next generation fission and fusion nuclear reactors require materials that can withstand operating temperatures greater than 500 °C, neutron irradiation doses of up to 200 displacements per atom (dpa), and potentially corrosive coolants such as the alkali liquid metals sodium, lithium, and NaK (Na33K eutectic alloy). Refractory alloys, such as Nb1Zr (Nb-1wt%Zr) and Molybdenum alloy TZM (Mo-0.5wt%Ti-0.08wt%Zr) have been traditionally considered viable candidates for advanced fission and fusion reactor concepts. However, it is relatively difficult to generate complex geometries of interest from these alloys using traditional manufacturing methods. In addition, there needs to be a concentrated effort to address refractory metal challenges at elevated temperature operation. In order to generate complex geometries of interest, modern manufacturing techniques are considered to increase the technological readiness level (TRL), cost-effectiveness, and schedule savings. This work focused on the continued development of laser powder bed fusion (L-PBF) additive manufacturing (AM) to improve both design flexibility, evaluate microstructure and properties, and ultimately accelerate the TRL and qualification of these processes and alloys for components to potentially be put into service. Niobium alloy Nb1Zr was identified through a down-selection process outlined in previous reports as a candidate to develop in L-PBF AM. Historically, Nb1Zr had been explored for high temperature fast spectrum fission reactors for both terrestrial and space applications. Molybdenum alloy TZM has also been considered for these reactor concepts due to exceptional high-temperature strength, creep resistance, and stability under irradiation. L-PBF AM of TZM has previously been investigated at LANL under the Microreactor program, NASA, ORNL, and in academia. However, due to the crack prone nature of TZM, L-PBF AM of TZM resulted in significant microcracking and additional development is required to pursue viable maturation. Other AM methods have been found to be more successful in printing TZM, and those alternatives approaches are discussed in this effort. The efforts detailed in this report focused on continued development of Nb1Zr through L-PBF and development of TZM via L-PBF and electron powder bed fusion (E-PBF). The objective of this work was to further the development of these AM techniques for the chosen refractory alloys, elucidating and addressing associated challenges through characterization of several demonstration builds. At LANL, Nb1Zr builds were completed using an EOS M290 and M400 machines, and a refractory alloy-dedicated L-PBF system, the Xact Metal XM200G, was installed. The XM200G primary purpose was to do the Nb1Zr parameter development process; however, due to difficulties associated with the machine installation and qualification process, it was decided to pivot development to the larger M400 and M290 machines. Although the supply of Nb1Zr powder was limited, it was sufficient to generate sub-scale metallographic specimens for the purpose of parameter development. This was first accomplished on the EOS M400 then the M290 due to machine schedule availability. Further development of TZM has been initiated at the University of Texas El Paso (UTEP) under contract with LANL to use both a heated build envelope L-PBF machine and E-PBF machine that have been found in the literature to mitigate microcracking. UTEP was provided with TZM powder and build plates to support parallel TZM parameter development across both machines. As part of the contract, UTEP will also be conducting microstructural characterization once optimized process parameters have been identified. The optimized process parameters for each machine will be used to generate a series of metallographic, mechanical, and surface finish specimens for subsequent characterization and testing. In the next section, we provide a detailed discussion of the methodology used for investigating the feasibility of leveraging these alloys for use in advanced reactor applications.

36 MATERIALS SCIENCE↗

Stellar Population Properties in the Stellar Streams around SPRC047

Abstract We have investigated the properties (e.g., age, metallicity) of the stellar populations of a ringlike tidal stellar stream (or streams) around the edge-on galaxy SPRC047 (z= 0.031) using spectral energy distribution (SED) fits to integrated broadband aperture flux densities. We used visual images in six different bands and Spitzer/IRAC 3.6μm data. We have attempted to derive best-fit stellar population parameters (metallicity, age) in three noncontiguous segments of the stream. Due to the very low surface brightness of the stream, we have performed a deconvolution with a Richardson–Lucy–type algorithm of the low spatial resolution 3.6μm IRAC image, thereby reducing the effect of the point-spread function aliasedemissionfrom the bright edge-on central galaxy at the locations of our three stream segments. Our SED fits that used several different star formation (SF) history priors, from an exponentially decaying SF burst to continuous SF, indicate that the age–metallicity–dust degeneracy is not resolved, most likely because of inadequate wavelength coverage and low signal-to-noise ratios of the low surface brightness features. We also discuss how future deep visual–near-infrared observations, combined with absolute flux calibration uncertainties at or below the 1% level, complemented by equally well absolute flux-calibrated observations in ultraviolet and mid-infrared bands, would improve the accuracy of broadband SED fitting results for low surface brightness targets, such as stellar streams around nearby galaxies that are not resolved into stars.

Astronomy & Astrophysics↗

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗

High and Ultra-High Temperature Reaction Kinetics by Single Nanoparticle Mass Spectrometry

Methodology is presented for non-destructive, optically-detected single nanoparticle (NP) mass spectrometry, with the goal of extracting surface reaction kinetics for single NPs at high temperatures. Methods for determining the NP charge, mass, and temperature as a function of time are discussed, and the data are used to extract both the absolute kinetics for mass change, as well as the efficiencies of the surface processes that cause them. Factors that contribute to the uncertainties in absolute and relative mass determination, and in the resulting kinetic parameters, are discussed. The method allows the NP-to-NP variations in initial reactivity to be measured directly, along with the time evolution of reactivity resulting from NP structural/compositional changes that occur under reaction conditions. The strengths and limitations of single nanoparticle mass spectrometry as a high temperature surface kinetics tool are discussed in the context of sublimation and O2 oxidation kinetics experiments for single hafnium (Hf) NPs at temperatures ranging above 2400 K. The Hf oxidation kinetics are compared to analogous oxidation experiments for silicon, graphite, and carbon black NPs. In all four cases, the oxidation chemistry was dominated by processes that result in net mass loss, and the distinct mechanisms responsible are discussed. All four NPs also eventually passivated, i.e., the efficiencies for oxidative etching decreased by at least two orders of magnitude, relative to the initial efficiencies. Furthermore, the passivation mechanisms, which are quite different for carbon, compared to silicon or hafnium, are discussed. Carbon NP passivation is attributed to structural isomerization leading to fully coordinated, fullerene-like NP surfaces, while for silicon and hafnium, passivation results from delayed formation of an oxide layer, triggered by accumulation of oxygen in the NP sub-surface region.

36 MATERIALS SCIENCE↗

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗

Response of Atmospheric River Width and Intensity to Aquaplanet Warming: A Detection Algorithm‐ and Background Moisture‐Independent Approach

The width of an atmospheric river (AR) is an important parameter when evaluating its impact. Although previous research suggests ARs will widen with global warming, a precise response has been muddled by the large sensitivity of width to a diverse set of AR detection techniques (ARDTs). Here, we propose a methodology that removes the influence of the ARDT by modeling AR‐integrated vapor transport (IVT) profiles as idealized exponential curves with free parameters given by background IVT, intensity above background IVT, and profile width. Notably, our definition for AR profile width does not include any explicit numerical thresholds, relative or absolute, for IVT. We apply our approach to a series of idealized aquaplanet experiments, first with a baseline sea surface temperature (SST) distribution, and then with +2K, +4K, and +6K uniform warming, so as to determine the contributions of each free parameter to AR width. We also apply our approach to high‐resolution atmosphere‐only models forced with SSTs modified to preserve historical variability but following projected warming over 2016–2050. Our results show that contributions to impacts‐relevant AR widening comes primarily from enhancements in background IVT and AR intensity, as opposed to from dynamic width changes.

54 ENVIRONMENTAL SCIENCES↗

Turbulent Parameters by airborne measurements over BNF in March 2025

The original data were collected during the AAF Engineering Flights (AEF2025) in the vicinity of the ARM Bankhead National Forest (BNF) Atmospheric Observatory (https://www.arm.gov/capabilities/observatories/bnf ) in northwestern Alabama in March 2025. The ARM Aerial Facility ArcticShark uncrewed aerial system (UAS, https://www.arm.gov/capabilities/observatories/aaf/uas) was based at the public-use airport of Posey Field, Alabama (FAA LID: 1M4, 34.28027778° N, 87.60055556° W, 283m MSL) from March 10 through March 24, 2025. The ArcticShark UAS performed nine flights, including eight research flights over the AMF3 (BNF Main Site) and Supplemental Facilities to measure atmospheric state, turbulence, surface IR temperature and imagery, and aerosol number concentration and size distribution. The current data set presents a collection of turbulent parameters in the atmospheric boundary layer or lower free troposphere based on airborne measurement throughout the field campaign. The primary instruments used to create the current data set were the Aircraft Integrated Meteorological Measurement System (AIMMS-30) and the fine-wire thermocouple probe.

Atmosphere↗

CHELAX-BNF: Turbulent Parameters by airborne measurements

The original data were collected on board the ARM Aerial Facility ArcticShark uncrewed aerial system (UAS; https://www.arm.gov/capabilities/observatories/aaf/uas ) during the “Characterizing HEterogeneous Land-Atmosphere eXchanges at BNF” field campaign (CHEAX-BNF; https://arm.gov/research/campaigns/aaf2025CHELAX-BNF ). The ARM Aerial Facility ArcticShark UAS was based at the public-use airport of Posey Field, AL (FAA LID: 1M4, 34.28027778° N, 87.60055556° W, 283m MSL) from May 28 through June 23, 2025. The ArcticShark UAS performed 5 flights, including 4 research flights over the BNF Main Site (ARM Mobile Facility 3, https://arm.gov/capabilities/observatories/amf ) and Supplemental Facilities to measure atmospheric state, turbulence, surface IR temperature and imagery, aerosol number concentration, and aerosol size distribution. The current data set presents a collection of turbulent parameters in the atmospheric boundary layer or lower free troposphere based on airborne measurement throughout the field campaign. The primary instruments used to create the current data set were the Aircraft Integrated Meteorological Measurement System (AIMMS-30) and the fine-wire thermocouple probe.

Aircraft Integrated Meteorological Measurement Sys↗

Enhanced joining strength in additive-manufactured polylactic-acid structures fused by embedded heated metallic meshes

Additively manufactured thermoplastic polymers, such as polylactic acid (PLA), hold significant promise for sustainable engineering structures, including wind turbine blades. Upscaling these structures beyond the limitations of 3D printer build volumes is a challenge; fusion joining presents a potential solution. This paper introduces a displacement-controlled resistance welding process for PLA, as an alternative to the typical force-controlled methods. Here, we investigated the bonding quality of resistance-welded and adhesive-bonded PLA beams through three-point bending and measured the surface deformations using digital image correlation. Different metal meshes (30%/0.11 mm Ni—Cu, 34%/0.07 mm Ni—Cu, and 36%/0.25 mm Co—Ni) served as heating elements. The process parameters were varied for the 34%/0.07 mm Ni—Cu mesh to identify an optimum set of parameters. Results showed that this optimized displacement-controlled welding achieved 94% of the original strength of monolithic samples. This indicates that the new welding process not only ensures high-quality bonding and fine surface finishing but also promotes sustainability, recyclability, and economic efficiency in various polymer and composite structural applications.

Additive manufacturing↗

High-Throughput Microstructural Characterization and Process Correlation Using Automated Electron Backscatter Diffraction

The need to optimize the processing conditions of additively manufactured (AM) metals and alloys has driven advances in throughput capabilities for material property measurements such as tensile strength or hardness. High-throughput (HT) characterization of AM metal microstructure has fallen significantly behind the pace of property measurements due to intrinsic bottlenecks associated with the artisan and labor-intensive preparation methods required to produce highly polished surfaces. This inequality in data throughput has led to a reliance on heuristics to connect process to structure or structure to properties for AM structural materials. In this study, we show a transformative approach to achieve laser powder bed fusion (LPBF) printing, HT preparation using dry electropolishing and HT electron backscatter diffraction (EBSD). This approach was used to construct a library of > 600 experimental EBSD sample sets spanning a diverse range of LPBF process conditions for AM Kovar. This vast library is far more expansive in parameter space than most state-of-the-art studies, yet it required only approximately 10 labor hours to acquire. Build geometries, surface preparation methods, and microscopy details, as well as the entire library of >600 EBSD data sets over the two sample design versions, have been shared with intent for the materials community to leverage the data and further advance the approach. Using this library, we investigated process–structure relationships and uncovered an unexpected, strong dependence of microstructure on location within the build, when varied, using otherwise identical laser parameters.

Characterization and Analytical Technique↗

Experimental investigation into coal wettability changes caused by reactions with scCO 2 -H 2 O

Geological CO 2 sequestration (GCS) can help mitigate global warming and enhance methane recovery from coal beds. However, few studies have linked the effects of CO 2 to surface chemistry changes controlling wetting behavior in deep coal beds. Contact angles (CAs) of CO 2 /N 2 -high volatile bituminous coal-water systems were measured under different temperatures and pressures. The surface chemistry and physical structure of coals were characterized to investigate changes in physicochemical properties and their relations with wettability after reactions. For N 2 treatment, the time-dependence of static and dynamic CAs were insignificant, ranging within 4°. For gaseous CO 2 treatment, the static CAs and the average advancing angles increased slightly. With supercritical (sc) CO 2 , both the static and dynamic CAs increased significantly, and θ adv changed to intermediate-wet (92°). Reactions with minerals exposed to scCO 2 resulted in greater surface roughness and heterogeneity, greater contact angle hysteresis and more surface sites occupied by scCO 2 rather than H 2 O. Increases in hydrophobic functional groups and decreases in hydrophilicity were shown by FTIR spectra, reflecting the shedding of polar oxygen-containing functional groups, reduction of hydrogen bonds, and increasing percentage of hydrocarbons. XRD patterns obtained following scCO 2 -treatment showed that crystallite growth and molecular polymerization were higher toward graphite-like. The calculated structural parameters of functional groups and crystallites both showed elevated coal rank. Changes in crystallite structure, notably higher carbon content and decreased negative surface charge, are unfavorable for water-wetting. Finally, this study contributes to understanding surface chemistry changes responsible for decreased wettability during CO 2 -enhanced coal bed methane recovery and GCS in coal reservoirs.

01 COAL, LIGNITE, AND PEAT↗

Microstructural characterization and equibiaxial flexural strength of CeO 2 and Ti-doped CeO 2

In this study, the synthesis of CeO 2 and titanium dioxide (TiO 2 ) doped CeO 2 (TDC) monoliths are investigated, and their fracture strength is assessed using an equibiaxial flexure testing technique at room temperature. Pellets were synthesized using conventional powder processing and sintering methods to produce the desired characteristics. The TiO 2 dopant concentration was optimized at 0.1 wt % TiO 2 to obtain dense, solid-solution pellets with an enhanced grain microstructure. A ball-on-ring fixture was used to obtain the TRS and Weibull parameters of over 30 pellets for CeO 2 and 0.1 wt % TDC to compare fracture behavior. The TRS of CeO 2 pellets ranged from 88 to 160 MPa and the TRS of 0.1 wt % TDC pellets ranged from 102 to 171 MPa, both being consistent with published values. Weibull parameters, such as characteristic strength and Weibull modulus, were extracted as 129 MPa and 8.5 for CeO 2 and 150 MPa and 9.3 for 0.1 wt % TDC, respectively. Although Hertzian contact damage was observed on compressive surfaces, failure initiation occurred on the tensile surfaces of both types of samples. Fracture surface analysis for CeO 2 indicated a predominantly intergranular fracture while 0.1 wt % TDC had a predominantly transgranular fracture mode. Finally, the TRS of 0.1 wt % TDC resulted in increased Weibull parameters when compared to CeO 2 , indicating sample chemistry and microstructure impact mechanical behavior for these samples.

36 MATERIALS SCIENCE↗

Interface, bulk and surface structure of heteroepitaxial altermagnetic α-MnTe films grown on GaAs(111)

Epitaxial MnTe films have recently seen a surge in research into their altermagnetic semiconducting properties. However, those properties may be extremely sensitive to structural and chemical modifications. We report a detailed investigation of the synthesis of the altermagnet α-MnTe on GaAs(111), which reveals the bulk defect structure of this material, the mechanism by which it releases strain from the underlying substrate, and the impact of oxidation on its surface. X-ray diffraction measurements show that α-MnTe layers with thicknesses spanning 45 to 640 nm acquire lattice parameters different from bulk, mostly due to thermal strain caused by the substrate rather than strain from the lattice mismatch. Through high-resolution transmission electron microscopy (TEM) measurement, we then unveil a misfit dislocation array at the interface, revealing the mechanism by which lattice strain is relaxed. TEM also reveals a stacking fault in the bulk, occurring along a glide plane parallel to the interface. The combination of TEM with polarized neutron reflectometry measurements finally reveals the impact of oxidation on the chemistry of the surface of uncapped MnTe. Furthermore, or findings highlight the subtle role of epitaxy in altering the structure of α-MnTe, providing potential opportunities to tune the altermagnetic properties of this material.

Altermagnets↗

Surface acoustic waves and lattice vibrations in two-dimensional Ti 3 ⁢C 2 ⁢𝑇 𝑥 (𝑇 =O, F) MXene films

We investigated surface acoustic wave (SAW) propagation and lattice vibrations in two-dimensional (2D) titanium carbide (Ti 3 ⁢C 2 ⁢𝑇 x ) MXene films as a function of surface termination and layer stacking, using atomistic simulations. We found that SAW propagation velocity is highly sensitive to both single-layer properties and interlayer bonding. Surface terminations significantly modulate wave behavior, with oxygen and fluorine terminations producing distinct effects on wave propagation, with oxygen-terminated monolayers exhibiting 20% higher wave speeds than fluorine counterparts due to strengthened intralayer bonds. Key observations include the transition from one to two layers causing wave speed variations, and the development of interlayer modes that generate more dispersed lattice vibrations. As the film layer thickness increases, SAW propagation becomes predominantly confined to the upper surface, with coherence of vibrational modes diminishing in multilayer structures. These findings suggest MXene terminations and layer stacking are crucial parameters for controlling SAW behavior, offering promising avenues for novel acoustic wave device applications.

2-dimensional systems↗

Toward a 2D Local Implementation of Quantum Low-Density Parity-Check Codes

Geometric locality is an important theoretical and practical factor for quantum low-density parity-check (qLDPC) codes that affects code performance and ease of physical realization. For device architectures restricted to two-dimensional (2D) local gates, naively implementing the high-rate codes suitable for low-overhead fault-tolerant quantum computing incurs prohibitive overhead. In this work, we present an error-correction protocol built on a bilayer architecture that aims to reduce operational overheads when restricted to 2D local gates by measuring some generators less frequently than others. We investigate the family of bivariate-bicycle qLDPC codes and show that they are well suited for a parallel syndrome-measurement scheme using fast routing with local operations and classical communication (LOCC). Through circuit-level simulations, we find that in some parameter regimes, bivariate-bicycle codes implemented with this protocol have logical error rates comparable to the surface code while using fewer physical qubits. Published by the American Physical Society 2025

Berthusen, Noah (ORCID:0000000275862786)↗

LBNL_DIATOMICS v1

LBNL_DIATOMICS is a software package developed by P. Zarzycki (Lawrence Berkeley National Laboratory) for computing isotopic effects in diatomic molecules using mass-dependent spectroscopic parameters obtained experimentally or from first principles. The package provides a streamlined pipeline from the potential energy surface to the equilibrium constant of an isotopic exchange reaction.

Zarzycki, Piotr [Lawrence Berkeley National Labora↗

Savannah River Site Climatology for Research Activities Conducted Onsite and in P-Area (1994-2023)

This report presents a climatology of meteorological variables at the Savannah River Site (SRS) using measurements from the P-Area and Central Climatology meteorological towers to support research activities at the Savannah River Site. Key values included in this report are wind, surface pressure, temperature, relative humidity, soil moisture, cloud fraction, cloud base, lightning, atmospheric stability, and boundary layer parameters. The Central Climatology tower is useful in that it is instrumented at four levels, allowing for an examination of height dependence on some variables.

54 ENVIRONMENTAL SCIENCES↗