Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vulnerabilities”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Atomically Ordered PdCu Electrocatalysts for Selective and Stable Electrochemical Nitrate Reduction

Electrochemical nitrate reduction (NO 3 RR) has attracted attention as an emerging approach to mitigate nitrate pollution in groundwater. Here, we report that a highly ordered PdCu alloy-based electrocatalyst exhibits selective (91% N 2 ), stable (480 h), and near complete (94%) removal of nitrate without loss of catalyst. In situ and ex situ XAS provide evidence that structural ordering between Pd and Cu improves long-term catalyst stability during NO 3 RR. In contrast, we also report that a disordered PdCu alloy-based electrocatalyst exhibits non-selective (44% N 2 and 49% NH 4 + ), unstable, and incomplete removal of nitrate. The copper within disordered PdCu alloy is vulnerable to accepting electrons from hydrogenated neighboring Pd atoms. This resulted in copper catalyst losses which were 10× greater than that of the ordered catalyst. The design of stable catalysts is imperative for water treatment because loss of the catalyst adds to the system cost and environmental impacts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mitigating Hydrogen-Induced Degradation of Iridium Anodes in Proton Exchange Membrane Water Electrolyzers

One promising method for reducing precious metal usage in proton exchange membrane water electrolysis is lowering iridium (Ir) loading at the anode. However, low-loading Ir catalysts often suffer from poor stability under high current densities. In this study, hydrogen (H2) crossover from the cathode to the anode is identified as a key degradation pathway. Temperature-programmed reduction confirms the reduction of IrO2 at 80 °C in a hydrogen environment, highlighting the vulnerability of IrO2-based catalysts to H2 exposure. To mitigate this effect, palladium (Pd) is introduced as an anode additive, acting as an H2 oxidation catalyst and mitigating IrO2 reduction. This protective role is verified by inductively coupled plasma optical emission spectroscopy and in situ X-ray absorption spectroscopy, showing significantly suppressed Ir dissolution at 80 °C under H2 flow when an O-covered Pd surface is present at oxygen evolution reaction potentials. Results from the current study identify a new strategy in improving activity and durability of catalysts in electrolyzers.

58 GEOSCIENCES↗
DOE OSTI2024

Porous Iron Electrodes Reduce Energy Consumption During Electrocoagulation of a Virus Surrogate: Insights into Performance Enhancements Using Three-Dimensional Neutron Computed Tomography

Electrocoagulation has attracted significant attention as an alternative to conventional chemical coagulation because it is capable of removing a wide range of contaminants and has several potential advantages. In contrast to most electrocoagulation research that has been performed with nonporous electrodes, in this study, we demonstrate energy-efficient iron electrocoagulation using porous electrodes. In batch operation, investigation of the external pore structures through optical microscopy suggested that a low porosity electrode with sparse connection between pores may lead to mechanical failure of the pore network during electrolysis, whereas a high porosity electrode is vulnerable to pore clogging. Electrodes with intermediate porosity, instead, only suffered a moderate surface deposition, leading to electrical energy savings of 21% and 36% in terms of electrocoagulant delivery and unit log virus reduction, respectively. Neutron computed tomography revealed the critical role of electrode porosity in utilizing the electrode’s internal surface for electrodissolution and effective delivery of electrocoagulant to the bulk. Energy savings of up to 88% in short-term operation were obtained with porous electrodes in a continuous flow-through system. Further investigation on the impact of current density and porosity in long-term operation is desired as well as the capital cost of porous electrodes.

42 ENGINEERING↗
DOE OSTI2024

Engineering the Novel Extremophile Alga Chlamydomonas pacifica for High Lipid and High Starch Production as a Path to Developing Commercially Relevant Strains

Microalgae offer a compelling platform for the production of commodity products, due to their superior photosynthetic efficiency, adaptability to nonarable lands and nonpotable water, and their capacity to produce a versatile array of bioproducts, including biofuels and biomaterials. However, the scalability of microalgae as a bioresource has been hindered by challenges such as costly biomass production related to vulnerability to pond crashes during large-scale cultivation. This study presents a pipeline for the genetic engineering and pilot-scale production of biodiesel and thermoplastic polyurethane precursors in the extremophile species Chlamydomonas pacifica. This extremophile microalga exhibits exceptional resilience to high pH (>11.5), high salinity (up to 2% NaCl), and elevated temperatures (up to 42 °C). Initially, we evolved this strain to also have a high tolerance to high light intensity (>2000 μE/m 2 /s) through mutagenesis, breeding, and selection. We subsequently genetically engineered C. pacifica to significantly enhance lipid production by 28% and starch accumulation by 27%, all without affecting its growth rate. We demonstrated the scalability of these engineered strains by cultivating them in pilot-scale raceway ponds and converting the resulting biomass into biodiesel and thermoplastic polyurethanes. This study showcases the complete cycle of transforming a newly discovered species into a commercially relevant commodity production strain. This research underscores the potential of extremophile algae, including C. pacifica, as a key species for the burgeoning sustainable bioeconomy, offering a viable path forward in mitigating environmental challenges and supporting global bioproduct demands.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Can N95 Respirators Be Reused after Disinfection? How Many Times?

The coronavirus disease 2019 (COVID-19) pandemic has led to a major shortage of N95 respirators, which are essential for protecting healthcare professionals and the general public who may come into contact with the virus. Thus, it is essential to determine how we can reuse respirators and other personal protective equipment in these urgent times. We investigated multiple commonly used disinfection schemes on media with particle filtration efficiency of 95%. Heating was recently found to inactivate the virus in solution within 5 min at 70 °C and is among the most scalable, user-friendly methods for viral disinfection. We found that heat (≤85 °C) under various humidities (≤100% relative humidity, RH) was the most promising, nondestructive method for the preservation of filtration properties in meltblown fabrics as well as N95-grade respirators. At 85 °C, 30% RH, we were able to perform 50 cycles of heat treatment without significant changes in the filtration efficiency. At low humidity or dry conditions, temperatures up to 100 °C were not found to alter the filtration efficiency significantly within 20 cycles of treatment. Ultraviolet (UV) irradiation was a secondary choice, which was able to withstand 10 cycles of treatment and showed small degradation by 20 cycles. However, UV can potentially impact the material strength and subsequent sealing of respirators. Finally, treatments involving liquids and vapors require caution, as steam, alcohol, and household bleach all may lead to degradation of the filtration efficiency, leaving the user vulnerable to the viral aerosols.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

From Filamentary Failure to Durable Halide Perovskite Memristors

Halide perovskites have emerged as promising materials for memristive devices. While their pronounced electrochemical reactivity and fast ionic mobility enable numerous advantages including low-voltage operation and fast switching, the same features also render perovskite-based memristors vulnerable to metallic shunts and poor endurance, limiting their practical applications. Here, we elucidate both the resistive switching and failure mechanisms in perovskite memristors with a fluorine-doped tin oxide (FTO)/methylammonium lead triiodide (MAPbI3)/Ag structure and demonstrate a strategy to substantially enhance device durability. As opposed to commonly invoked filamentary mechanisms, electrical, structural, and spectroscopic analyses reveal that resistive switching arises from interfacial barrier modulation by reversible Ag redox reactions that drive electrochemical doping/dedoping within the perovskite. Device failure, however, originates from metallic Ag0 filamentation that ultimately forms permanent conductive pathways. Introducing an ultrathin Al2O3 interlayer at the inert-electrode interface improves device endurance by more than 30-fold, exceeding 15,000 switching cycles, without compromising other performance metrics. Interfacial characterization indicates that the Al2O3 layer modifies wettability of Ag deposits, promoting planar island growth rather than through-film filamentation. These findings establish a clear link between interfacial electrochemistry, metal precipitation behavior, and memristor reliability, highlighting inert-electrode interfacial engineering as an effective pathway toward durable perovskite-based memristors.

14 SOLAR ENERGY↗
DOE OSTI2026

An Activity-Based Sensing Approach to Monitor Nanomaterial-Promoted Changes in Labile Metal Pools in Living Systems

Metal-based nanoparticles are a promising class of materials for diagnosis and treatment of cancer and other diseases. However, mechanisms of action of these nanomedicines remain insufficiently understood due in large part to our limited understanding of the dynamic equilibria between solid metal nanoparticles and labile metal ions generated from these nanoparticles within complex biological milieus. Here, we apply activitybased sensing to directly identify and investigate the fate of labile copper pools with metal and oxidation state-specificity generated by anticancer copper nanomedicines. We found that treatment of cells with copper-releasing nanoparticles alter labile Cu(I)/Cu(II) ratios through an increase in labile Cu(II), while overall labile copper levels decrease. Labile copper release triggers compensatory responses in two major antioxidant pathways, glutathione (GSH) and nuclear factor erythroid 2-related factor 2 (NRF2), as well as in metal homeostasis to limit copper availability via regulation of copper export (ATP7B) and copper import (CTR1) proteins. These findings establish the value of activity-based sensing as a generalizable approach for labile metal imaging to help decipher molecular mechanisms of bioactive metal nanoparticles and guide the development of more effective nanomedicine diagnostics and therapies to target metal-dependent disease vulnerabilities.

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

Molecular Interaction Regulates the Performance and Longevity of Defect Passivation for Metal Halide Perovskite Solar Cells

Defect passivation constitutes one of the most commonly used strategies to fabricate highly efficient perovskite solar cells (PSCs). However, the durability of the passivation effects under harsh operational conditions has not been extensively studied regardless of the weak and vulnerable secondary bonding between the molecular passivation agents and perovskite crystals. Here, we incorporated strategically designed passivating agents to investigate the effect of their interaction energies on the perovskite crystals and correlated these with the performance and longevity of the passivation effects. We unraveled that the passivation agents with a stronger interaction energy are advantageous not only for effective defect passivation but also to suppress defect migration. The prototypical PSCs treated with the optimal passivation agent exhibited superior performance and operational stability, retaining 81.9 and 85.3% of their initial performance under continuous illumination or nitrogen at 85 °C after 1008 h, respectively, while the reference device completely degraded during that time. Furthermore, this work provides important insights into designing operationally durable defect passivation agents for perovskite optoelectronic devices.

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

Origins of Acid-Gas Stability Behavior in Zeolitic Imidazolate Frameworks: The Unique High Stability of ZIF-71

Zeolitic imidazolate frameworks (ZIFs) are promising materials for industrial process separations, but recent literature reports have highlighted their vulnerability to acid gases (e.g., SO 2 , CO 2 , NO 2 , H 2 S), often present in practical applications. While previous work has documented the widely varying stability behavior of many ZIFs under varying (humid and dry) acid gas environments, efforts to explain or correlate these experimental observations via empirical descriptors have not succeeded. A key observation is that ZIF-71 (RHO topology) is an extraordinarily stable ZIF material, retaining both structure and porosity under prolonged humid SO 2 exposure whereas many other well-known ZIFs with different linkers and topologies (such as ZIF-8) were shown to degrade. Through a combination of hybrid quantum mechanics/molecular mechanics (QM/MM) based methods and statistical mechanical models, we successfully explain this important experimental observation via atomistic investigations of the reaction mechanism. Our holistic approach reveals an ~9 times lower average defect formation rate in ZIF-71 RHO compared to ZIF-8 SOD, leading to the conclusion that the observed experimental stability of this material rises from kinetic effects. Moreover, our analysis reveals that differing stability of the two materials is determined by the distributions of acid gas molecules, which is difficult to capture using empirical descriptors. Furthermore, our results suggest wider applicability of the present approach, toward identifying tuned functional groups and topologies that move the acid gas distributions away from more reactive sites and thus allow enhanced kinetic stability.

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

Enriched Molecular-Level View of Saline Wetland Soil Carbon by Sensitivity-Enhanced Solid-State NMR

Soil organic matter (SOM) plays a major role in mitigating greenhouse gas emission and regulating earth’s climate, carbon cycle, and biodiversity. Wetland soils account for one-third of all SOM; however, globally, coastal wetland soils are eroding faster due to increasing sea-level rise. Our understanding of carbon sequestration dynamics in wetlands lags behind that of upland soils. Here, we employ solid-state nuclear magnetic resonance (ssNMR) to investigate the molecular-level structure of biopolymers in wetland soils spanning 11 centuries. High-resolution multidimensional spectra, enabled by dynamic nuclear polarization (DNP), demonstrate enduring preservation of molecular structures within herbaceous plant cores, notably condensing aromatic motifs and carbohydrates, even over a millennium, with the preserved cores constituting a decreasing minority among molecules from decomposition and repolymerization with depth and age. Such preserved cores occur alongside molecules from the decomposition of loosely packed parent biopolymers. These findings emphasize the relative vulnerability of coastal wetland SOM when exposed to oxygenated water due to geological and anthropogenic changes.

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

Persistence and Plasticity in Conifer Water–Use Strategies

The selective use of seasonal precipitation by vegetation is critical to understanding the residence time and flow path of water in watersheds, yet there are limited datasets to test how climate alters these dynamics. Here, we use measurements of the seasonal cycle of tree ring δ 18 O for two widespread conifer species in the Rocky Mountains of North America to provide a multi-decadal depiction of the seasonal origins of forest water use. The results show that while the conifer tree stands had a dominant preference for use of snowmelt, there were multi-annual periods over the last four decades when use of summer precipitation was preferential. Utilization of summer rain emerged during years with increased snowfall and tree growth, suggesting that summer rain enhanced the transpiration stream only during the periods of highest water use. We hypothesize this could be explained through shallowing of the root profile during wetter periods and/or through the influence of changing water table depths on the residence time of summer precipitation in the root zone. We suggest the tree ring proxy approach used here could be applied in other watersheds to provide critical insight into the temporal dynamics of plant water use that could not be inferred from short measurement campaigns. Furthermore, these data on the seasonal origins of forest water are critical for understanding forest vulnerability to drought, the processes that affect precipitation pathways and residence time in watersheds and the interpretation of tree ring proxy data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Rapid Viscoelastic Deformation Slows Marine Ice Sheet Instability at Pine Island Glacier

The ice sheets of the Amundsen Sea Embayment (ASE) are vulnerable to the marine ice sheet instability (MISI), which could cause irreversible collapse and raise sea levels by over a meter. The uncertain timing and scale of this collapse depend on the complex interaction between ice, ocean, and bedrock dynamics. The mantle beneath the ASE is likely less viscous (~10 18 Pa s) than the Earth's average mantle (~10 21 Pa s). In this paper we show that an effective equilibrium between Pine Island Glacier's retreat and the response of a weak viscoelastic mantle can reduce ice mass lost by almost 30% over 150 years. Other components of solid Earth response—purely elastic deformations and geoid perturbations—provide less stability than the viscoelastic response alone. Uncertainties in mantle rheology, topography, and basal melt affect how much stability we expect, if any. Our study indicates the importance of considering viscoelastic uplift during the rapid retreat associated with MISI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

River Regulation Alleviates the Impacts of Climate Change on US Thermoelectricity Production

Climate change and anthropogenic activities can both influence stream temperature, which is a critical indicator of thermoelectric power generation potential. Previous studies emphasized the vulnerability of power-generation systems to climate warming (ref 1, 2). Despite the fact that two-thirds of the global rivers are regulated by dams (ref 3), the relative impacts of water management and climate change on stream temperature and thermoelectric production have not yet been evaluated. Using an integrated modeling system applied to the conterminous United States, simulations confirm the dominant role of climate change in shaping the future trends of stream temperature. However, water management, reservoir regulation in particular, can mitigate the stream temperature warming and reduce the frequency and magnitude of extreme high stream temperatures through redistribution of the seasonal streamflow. Moreover, water management can alleviate the loss of power plant usable capacity from climate change to the level of reduced warming from emission mitigation. By regulating the flow and temperature regimes, water management is a promising adaptation measure to alleviate climate change induced effects in the water-energy nexus.

Zhang, Xiao↗
DOE OSTI2021

Climatology and Evolution of the Antarctic Peninsula Föhn Wind-Induced Melt Regime From 1979–2018

Warm and dry föhn winds on the Antarctic Peninsula (AP) cause surface melt that can destabilize vulnerable ice shelves. Topographic funneling of these downslope winds through mountain passes and canyons can produce localized wind-induced melt that is difficult to quantify without direct measurements. Our Föhn Detection Algorithm (FöhnDA) identifies the surface föhn signature that causes melt from measurement by 12 Automatic Weather Stations on the AP, that train a machine learning model to detect föhn in 5 km Regional Atmospheric Climate Model 2 (RACMO2.3p2) simulations and in the ERA5 reanalysis model. We estimate the fraction of AP surface melt attributed to föhn and possibly katabatic winds and identify the drivers of melt, temporal variability, and long-term trends and evolution from 1979–2018. We find that föhn wind-induced melt accounts for 3.1% of the total melt on the AP and can be as high at 18% close to the mountains where the winds funnel through mountain canyons. Föhn-induced surface melt does not significantly increase from 1979–2018, despite a warmer atmosphere and more positive Southern Annular Mode. However, a significant increase (+0.1 Gt y-1) and subsequent decrease/stabilization occur in 1979–1998 and 1999–2018, consistent with the AP warming and cooling trends during the same time periods. Föhn occurrence, more than föhn strength, drives the annual variability in föhn-induced melt. Long-term föhn-induced melt trends and evolution are attributable to seasonal changes in föhn occurrence, with increased occurrence in summer, and decreased occurrence in fall, winter, and early spring over the past 20 years.

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

Sequential Abiotic-Biotic Processes Drive Organic Carbon Transformation in Peat Bogs

Peatlands, which store one-third of the terrestrial carbon (C), are subject to large disturbances under a changing climate. It is crucial to understand how microbial and physiochemical factors affect the vulnerability of these large C stores to predict climate induced greenhouse gas fluxes. Here, we used a combination of mass spectrometry and spectroscopy techniques, to understand sequential biotic and abiotic degradation pathways of Sphagnum fallax leachate in an anaerobic incubation experiment, in the presence and absence of microorganisms. Removal of microorganisms was carried out by passing aqueous samples through 0.2 µm filters. In this work, our results revealed that Sphagnum fallax leachate degradation by abiotic reactions is a significant contributor to CO 2 production. Further, abiotic factors, such as low pH, are responsible for partial dissolved organic carbon (DOC) degradation that produces bioavailable compounds that shift microbial metabolic pathways and stimulate respiration in peat bogs. Acid-catalyzed hydrolysis of Sphagnum-produced glycosides can provide the microbial communities with glucose and stimulate microbial respiration of DOC to CO 2 . These results, while unique to peatlands, demonstrate the importance and underscore the complexity of sequential abiotic and biotic processes on C cycling in peat bogs. It is therefore crucial to incorporate abiotic degradation and sequential below-ground biotic and abiotic interactions into climate models for a better prediction of the influence of climate change on DOC stability in peatlands. These findings might not be representative of other ecosystems with different environmental conditions including mineral-rich peatlands and plant matter in surface peat horizons that comprise discrete microbial populations, and DOC composition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improving CLM5.0 Biomass and Carbon Exchange Across the Western United States Using a Data Assimilation System

The Western United States is dominated by natural lands that play a critical role for carbon balance, water quality, and timber reserves. This region is also particularly vulnerable to forest mortality from drought, insect attack, and wildfires, thus requiring constant monitoring to assess ecosystem health. Carbon monitoring techniques are challenged by the complex mountainous terrain, thus there is an opportunity for data assimilation systems that combine land surface models and satellite-derived observations to provide improved carbon monitoring. Here, we use the Data Assimilation Research Testbed to adjust the Community Land Model (CLM5.0) with remotely sensed observations of leaf area and above-ground biomass. The adjusted simulation significantly reduced the above-ground biomass and leaf area, leading to a reduction in both photosynthesis and respiration fluxes. The reduction in the carbon fluxes mostly offset, thus both the adjusted and free simulation projected a weak carbon sink to the land. This result differed from a separate observation-constrained model (FLUXCOM) that projected strong carbon uptake to the land. Simulation diagnostics suggested water limitation had an important influence upon the magnitude and spatial pattern of carbon uptake through photosynthesis. We recommend that additional observations important for water cycling (e.g., snow water equivalent, land surface temperature) be included to improve the veracity of the spatial pattern in carbon uptake. Furthermore, the assimilation system should be enhanced to maximize the number of the simulated state variables that are adjusted, especially those related to the recommended observed quantities including water cycling and soil carbon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗