Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “topography”

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

Real-Time Ammonia and Humidity Monitoring with Ultra-Fast Conductometric Sensors Based on Porphyrin and Phthalocyanine Complexes

Organic semiconductors like porphyrins and phthalocyanines are attracting a wide range of researchers due to their versatile electrical properties and sensing performances in conductometric sensors. In this study, we investigate two types of π-extended porphyrins, which share the same macrocyclic structure but differ in their central metal. These porphyrins are employed as sublayers in bilayer heterojunction devices, with the lutetium bisphthalocyanine complex, LuPc 2 , serving as the common top layer. Remarkably, the central metal in the porphyrin macrocycle significantly influences the solubility of the materials and, consequently, the surface topography of the resulting bilayer heterojunction devices. This structural variation translates into distinct electrical and sensing performances. The device incorporating nickel as metal centre (AM2) demonstrates superior sensitivity towards NH3, with a relative response (RR) of ca. -7% at 90 ppm, an ultra-fast response time of about 9 s, and an impressive limit of detection (LOD) of 250 ppb, whereas, the device that has zinc as metal centre in sublayer (AM3) exhibits RR value of ca. -0.9% at 90 ppm with t 90 of ca. 120 s and LOD of 2 ppm. Both devices are evaluated under randomly varying NH 3 concentration and RH value. The results shows that the AM2-based sensor allows following NH 3 in real-time, while the AM3-based sensor delivers an average concentration over time. On the other hand, the AM2-based sensor exhibits slow kinetics under RH exposure, while the AM3-based sensor precisely mirrors the pattern of random RH changes generated by the software, demonstrating its exceptional responsiveness and accuracy in tracking humidity fluctuations. In conclusion, these findings underscore the critical role of the metal centre in tuning the electrical and sensing properties of the heterojunction devices.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

The Sensitivity of Joint Inversions of Seismic and Geodynamic Data to Mantle Viscosity

Seismic tomography has revealed the existence of large-scale velocity heterogeneities in the mantle. The interpretation of seismic velocity anomalies in terms of temperature and chemical composition is nonunique. We use geodynamic observations including gravity, plate motions, dynamic topography, and excess ellipticity of the core-mantle boundary combined with seismic observations to investigate the thermo-chemical structure of the mantle through joint inversions. An outstanding issue, however, is the physical connection between mantle density anomalies and the surface geodynamic observations, which requires knowledge of the mantle viscosity structure. Here we perform joint inversions assuming different viscosity profiles and examine the dependence of the results on the viscosity. We first assume that mantle heterogeneity is due to thermal variations, which places a constraint on the relation between seismic velocity and density, and we subsequently relax the constraint to allow for potential nonthermal effects. In all of our joint inversions, a nonthermal origin of density anomalies is required to explain the geodynamic data, though the amount varies with the assumed viscosity structure. A common observation is a high-density chemical signal in the center of the large low-shear-velocity provinces at the base of the mantle resulting in a near neutral or slightly dense overall buoyancy there. Using the derived density models and their corresponding viscosity profiles, we also calculate instantaneous mantle flow fields. The predicted flow fields derived from joint inversions are generally similar but are quite different from flow fields using density models derived from a posteriori scaling of pure seismic tomography models.

58 GEOSCIENCES↗
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 Plume Liftoff Dynamics and Surface Expressions

The water surface expression of liftoff and its dependence on discharge are examined using numerical simulations with the Regional Ocean Modeling System (ROMS). Liftoff is the process by which buoyant river water separates from the bed and flows over denser saltwater. During low-discharge conditions liftoff occurs in the river and is accompanied by a change in the surface slope. During high-discharge conditions liftoff occurs outside the mouth and generates a ridge on the water surface. The location and height of the ridge can be described by analytical equations in terms of discharge, shelf slope, and river mouth aspect ratio. The offshore distance and height of the ridge are proportional to the river discharge and vary inversely with river mouth aspect ratio. For steep shelf slopes liftoff occurs close to the river mouth and generates a large ridge. The ridge is modified, but not eliminated, by the presence of tides. The water surface slope change at the ridge peak is large enough to be detected by the upcoming Surface Water and Ocean Topography (SWOT) altimeter and can be used to identify the liftoff location during high discharge. However, during low discharge the water surface slope change at the liftoff location is too small to be detected by SWOT. These results indicate that remote measurements of the presence or absence of the ridge may be useful to distinguish between low and high flows, and remote measurements of the ridge location or height could be used to estimate freshwater discharge.

42 ENGINEERING↗
DOE OSTI2020

Variations of Stress Parameters in the Southern California Plate Boundary Around the South Central Transverse Ranges

We examine the stress field in Southern California with a focus on the region near the South Central Transverse Ranges (SCTR), using a refined stress inversion methodology. We independently invert declustered and aftershocks focal mechanisms from the years 1981 to 2017. Comparisons between the stress field orientations and the stress ratios from these inversions provide information on local sources of the stress field. Regionally, the SHmax tends to be toward NNE and the stress ratios vary from transtensional stress regime near the Eastern California Shear Zone (ECSZ), to strike-slip faulting near the SCTR, and toward transpression near the Western Transverse Ranges. Detailed analysis of stress parameters near the SCTR indicates deviations from the regional strike-slip faulting including the following: (1) The San Bernardino Mountain area shows S Hmax direction toward NNW and transpressional stress components likely associated with the motions of the San Andreas Fault and ECSZ; (2) the Cajon Pass and San Gorgonio Pass show transpressional stress regime near the bottom of the seismogenic zone likely associated with the elevated topography; and (3) in Crafton Hills, rotation of the principal stress plunges and SHmax direction and transtensional stress regime below ~10 km, along with lower estimated apparent friction coefficient, suggest a weak fault possibly associated with deep creep. The performed multiscale analysis resolves effects of regional and local loadings. The average stress parameters in the analyzed 37 years do not show significant temporal variations of stress field near the SCTR region.

58 GEOSCIENCES↗
DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗
DOE OSTI2020

Improvement of Atmospheric Objective Analysis Over Sloping Terrain and Its Impact on Shallow-Cumulus Clouds in Large-Eddy Simulations

Surface topography strongly impacts the regional atmospheric circulation dynamically and thermodynamically. Over the Great Plains in the United States, the gently tilted slope is an important factor that impacts clouds, convection and regional circulations. This study enhances an atmospheric constrained variational analysis by using a terrain-following sigma vertical coordinate and applies to the data collected at the Southern Great Plain (SGP) site of the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program. Sensitivity studies are performed to examine the impact of the terrain effects on the derived large-scale atmospheric forcing fields and the simulated shallow-cumulus clouds in large-eddy simulation (LES) driven by the forcing. We found that the terrain impacts on the large-scale forcing fields are mainly at lower levels and are strongly controlled by up/downslope winds. The response of the derived forcing to the slope of the terrain is monotonic, but the response of the simulated shallow-cumulus clouds is more complex and depends on several factors. Overall, the terrain impact is small over SGP due to the small slope angle. However, the flat-surface assumption may cause larger biases in the large-scale forcing fields at locations with steeper terrain. The new sigma-coordinate algorithm, with its consideration of surface slope, should be more suitable to derive large-scale objective analysis over regions with steep terrain for application to force single-column models, cloud resolving models and LES models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluation of Atmospheric Boundary Layer Height From Wind Profiling Radar and Slab Models and Its Responses to Seasonality of Land Cover, Subsidence, and Advection

In this study, we evaluated the effect of land cover, atmospheric subsidence, and advection on the annual dynamics of atmospheric boundary layer (ABL) height from two contrasting sites. The first site is the Walker Branch forest, a deciduous forest of temperate climate, complex topography, and cloudy summers. The second site is the Sacramento-San Joaquin River Delta, a site of Mediterranean climate, flat terrain on a local scale, and clear summers. After testing a new algorithm to calculate ABL heights from 915 MHz radar wind profilers, we evaluated a hierarchy of three slab models to recreate the diurnal and annual patterns of ABL growth. We found that the lower ABL heights in the Delta, particularly during late summer, are driven by the combined effects of increased atmospheric subsidence and marine air advection. In both sites, the annual pattern of ABL height was strongly correlated to total daily incoming radiation, and in the Delta, the annual pattern of ABL height closely followed the seasonal patterns of sensible heat flux from a mosaic of different land covers. A land composite of latent and sensible heat fluxes obtained through a meso-network of eddy covariance measurements and the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission resulted in higher model skill, thus showing that land cover heterogeneity is an important driver of ABL growth. To conclude, model simulations show that in the Delta, restoring agricultural land to wetlands with large open water areas could result in a reduction of ABL height during those months with low subsidence and advection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluating Variable‐Resolution CESM Over China and Western United States for Use in Water‐Energy Nexus and Impacts Modeling

Abstract Climate models are critical tools to study earth system processes and are often further downscaled to provide refined‐resolution data sets for regional water‐energy impact analyses. In this study, the Variable‐Resolution Community Earth System Model (VR‐CESM), a new dynamical downscaling method, is employed to generate 1/8° simulations for China and the western United States over 1970–2006. Precipitation, temperature, snowpack, radiation, and wind speed are compared with two 1° global climate models, regional climate model ensemble, and reference data sets to evaluate their fidelity. Simulated precipitation skill is dependent upon the reference data sets used, a slight (<+10%) and spatially consistent wet bias in China during summer and the western United States during winter. Precipitation bias in VR‐CESM is the smallest among models and demonstrates the added value provided by resolution refinement. However, temperature biases are generally greater than precipitation due to more complicated interactions among processes. VR‐CESM simulates +2°C warm bias in the low‐elevation regions in China but −2°C cold bias in mountains. However, DJF cold bias (−2.23°C to−3.37°C) and JJA warm bias (+1.13°C to +1.87°C) exist in the western United States. The modeling skills of precipitation and temperature at various topographies are similar in the two evaluated regions. The cold bias was likely affected by downward shortwave radiation deficit (−15%) and slow bias of surface wind speeds (−16% to −29%). Despite the wet and cold bias, snow water equivalent and snow cover are underestimated in mountains. We conclude that VR‐CESM outperforms other models evaluated and provides the best estimation of downscaled climate model data for water‐energy studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Exploring Hydrologic Model Process Connectivity at the Continental Scale Through an Information Theory Approach

Exploring water fluxes between hydrological model (HM) components is essential to assess and improve model realism. Many classical metrics for HM diagnosis rely solely on streamflow and hence provide limited insights into model performance across processes. This study applies an information theory measure known as “transfer entropy” (TE) to systematically quantify the transfer of information among major HM components. To test and demonstrate the benefits of TE, we use the Framework for Understanding Structural Errors (FUSE) model to mimic and compare four commonly used HM structures, VIC, PRMS, SACRAMENTO, and TOPMODEL, across 671 catchments spanning a variety of hydrologic regimes in the conterminous United States. We explore connections between HM components and catchment landscape characteristics (e.g., climate, topography, soil, and vegetation) and characterize their nonlinear associations using distance correlation and Spearman correlation coefficients. Our results indicate that while the information transferred from precipitation to runoff is similar across model structures (likely as a result of calibration), the information transferred among other components can vary significantly from a FUSE structure to another. We find that aridity, precipitation duration and frequency, snow fraction, mean elevation, forest area, and leaf area index are often significantly associated with TE between the main HM components. Here, we propose that the presence of meaningful nonlinear associations can be used to diagnose process representation in HMs. Our results highlight the necessity to enhance the conventional streamflow-only calibration approach for a more realistic representation of water dynamics in the models.

large-sample hydrology↗
DOE OSTI2021

The First Detection of an Earthquake From a Balloon Using Its Acoustic Signature

Extreme temperature and pressure conditions on the surface of Venus present formidable technological challenges against performing ground-based seismology. Efficient coupling between the Venusian atmosphere and the solid planet theoretically allows the study of seismically generated acoustic waves using balloons in the upper atmosphere, where conditions are far more clement. However, earthquake detection from a balloon has never been demonstrated. We present the first detection of an earthquake from a balloon-borne microbarometer near Ridgecrest, CA in July 2019 and include a detailed analysis of the dependence of seismic infrasound, as measured from a balloon on earthquake source parameters, topography, and crustal and atmospheric structure. Our comprehensive analysis of seismo-acoustic phenomenology demonstrates that seismic activity is detectable from a high-altitude platform on Earth, and that Rayleigh wave-induced infrasound can be used to constrain subsurface velocities, paving the way for the detection and characterization of such signals on Venus.

58 GEOSCIENCES↗
DOE OSTI2021

Climate Signatures on Lake And Wetland Size Distributions in Arctic Deltas

Abstract Understanding how thermokarst lakes on arctic river deltas will respond to rapid warming is critical for projecting how carbon storage and fluxes will change in those vulnerable environments. Yet, this understanding is currently limited partly due to the complexity of disentangling significant interannual variability from the longer‐term surface water signatures on the landscape, using the short summertime window of optical spaceborne observations. Here, we rigorously separate perennial lakes from ephemeral wetlands on 12 arctic deltas and report distinct size distributions and climate trends for the two waterbodies. Namely, we find a lognormal distribution for lakes and a power‐law distribution for wetlands, consistent with a simple proportionate growth model and inundated topography, respectively. Furthermore, while no trend with temperature is found for wetlands, a statistically significant decreasing trend of mean lake size with warmer temperatures is found, attributed to colder deltas having deeper and thicker permafrost preserving larger lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Determining Spatial Scales of Soil Moisture—Cloud Coupling Pathways Using Semi-Idealized Simulations

The spatial scales of land-atmosphere interactions over heterogeneous, realistic soil moisture patterns are critical to improve our understanding of land-atmosphere interactions. However, these scales are poorly quantified. We use high-resolution numerical experiments to quantify the spatial scales for a selected day during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) campaign in Oklahoma. A semi-idealized approach (removing background winds and surface topography while retaining the real-world data for other boundary conditions) enables examination of the local influence from the surface. Simulations that apply uniform surface, variable land cover with uniform soil moisture, and variable land cover with variable soil moisture are then compared. Land cover variations, especially water and urban surfaces, create two variance peaks in the surface sensible heat flux in the scales >10 km (meso-β) and <6 km (meso-γ). Both peaks are enhanced by the observation-based soil moisture distribution that has most of its variance in the meso-β scale. Atmospheric response in the meso-β scale is hydrostatic with weak vertical motion. The land cover variations induce secondary circulations in the scales of 2.5–6 km. The additional soil moisture variability intensifies extreme circulations over dry soil patches and sharp sensible heat flux gradients, slightly extending the characteristic scales to 2.5–9 km. While their areal coverage is small, soil moisture variations on such “hotspots” lead to substantially larger convective cloud structure and domain-mean rainfall, calling for a consideration of their non-linear effect in large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Mesoscale Convective Systems Simulated by a High–Resolution Global Nonhydrostatic Model Over the United States and China

Mesoscale convective systems (MCSs) contribute a large fraction of warm-season precipitation and generate hazardous weather with substantial socio-economic impacts. Uncertainties in convection parameterizations in climate models limit our understanding of MCS characteristics and reliability of future projection. We examine the MCS simulation from the global 14-km Nonhydrostatic ICosahedral Atmospheric Model (NICAM) without cumulus parameterization against satellite observation from Global Precipitation Measurement (GPM) during 2001-2008. We focus on MCSs over the central U.S. and eastern China where MCSs prevail from March to August. A process-oriented tracking method incorporating both cloud and precipitation criteria is used to identify and track MCSs. About 140/100 MCSs initiate in the central U.S./eastern China per warm season and most of them initiate east of high mountains and in coastal regions. The frequency distribution of MCS lifetime is well captured in NICAM. But the simulated MCSs have stronger precipitation, smaller precipitation area, and larger cold cloud system than observed in both regions, which may be caused by weak entrainment as it is not well resolved at 14 km resolution. Here, the simulated MCS number is also underestimated in summer. By examining the climatological and MCS large-scale environments, the significant underestimation of MCS number in summer over the central U.S. may be attributed to the large climatological dry bias in the atmosphere. For China, mean moisture in summer is well simulated but deficiency in capturing the dynamic condition related to the coastal topography for triggering convection may have contributed to underestimation of MCS even in a sufficiently moist environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Hyporheic Exchange Due to Cobbles on Sandy Beds

Hyporheic exchange induced by riverbed topography and roughness provides important ecosystem functions. We investigated how cobble clusters embedded in a finer-grained sand bed affect the near-bed channel flow and the exchange of surface and subsurface water. We tested how the spacing and embeddedness of cobbles altered hyporheic exchange through a three-dimensional fully coupled surface-subsurface model. The 3D modeling framework captured the full physics of the exchange process locating the lateral position of upwelling zones on the side of the cobble. As the cobbles protrude more into the channel, eddies appear downstream, with smaller vortices where cobbles are closer together. In our simulations, hyporheic exchange increased with the spacing and protrusion ratio. The travel time in the hyporheic zone also increased with the protrusion ratio because the hyporheic flow paths are longer when the cobble is more protruded. With different spacings at the highest protrusion ratio, the average slope of early time breakthrough curve was steeper for large spacing settings: the travel time increases as the spacing decreases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Coastal Marsh Eco‐Geomorphologic Change on Saltwater Intrusion Under Future Sea Level Rise

Abstract Coastal saltwater intrusion (SWI) is one key factor that affects the hydrology, ecology, and biogeochemistry of coastal ecosystems. Future climate change, especially intensified sea level rise (SLR), is expected to trigger SWI to encroach on coastal freshwater aquifers more intensively. Numerous studies have investigated decadal/century scale SWI under SLR by assuming a static coastal landscape topography. However, coastal landscapes are highly dynamic in response to SLR, and the impact of coastal landscape evolution on SWI has received very little attention. Thus, this study used a coastal marsh landscape as an example and investigated how coastal marsh evolution affects future SWI with a physically‐based coastal hydro‐eco‐geomorphologic model, Advanced Terrestrial Simulator. Our numerical experiments showed that it is very likely that the marsh elevation increases with future SLR due to sediment deposition, and a depression zone is formed due to different marsh accretion rates between the ocean boundary and the inland. We found that marsh accretion may significantly reduce the surface saltwater inflow at the ocean boundary, and the evolved topographic depression zone may prolong the residence time of surface ponded saltwater, affecting subsurface salinity distribution differently. We also predicted that marshlands may become more sensitive to upland freshwater supply under future SLR, compared with previous predictions without marsh evolution. This study demonstrates the importance of coastal evolution to coastal freshwater‐saltwater interaction. The eco‐geomorphologic effect may not be ignored when evaluating coastal SWI under SLR at decadal or century scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Topographically Scattered Infrasound Waves Observed on Microbarometer Arrays in the Lower Stratosphere

When an acoustic wave strikes a topographic feature, some of its energy is scattered. Sensors on the ground cannot capture these scattered signals when they propagate at high angles. We report observations of upwardly-scattered acoustic waves prior to refraction back to the ground, intercepting them with a set of balloon-borne infrasound microbarometers in the lower stratosphere over northern Sweden. We show that these scattered waves generate a coda whose presence can be related to topography beneath balloons and low-altitude acoustic ducts. The inclination of the coda signals changes systematically with time, as expected from waves arriving from scatterers successively closer to receivers. The codas are present when a temperature inversion channels infrasound from a set of ground chemical explosions along the ground, but are absent following the inversion's dissipation. Since scattering partitions energy away from the main arrival, these observations imply a mechanism of amplitude loss that had previously been inaccessible to measurement. As such, these data and results allow for a better comprehension of interactions between atmospheric infrasound propagation and the solid earth

58 GEOSCIENCES↗
DOE OSTI2023

Response of Extreme Rainfall to Atmospheric Warming and Wetting: Implications for Hydrologic Designs Under a Changing Climate

Understanding the processes of rainfall extremes and their response to anthropogenic climate change is pivotal for improved adaptation of unprecedented flood hazards around the world. Here we take the record-breaking 20 July 2021 storm over central China as an example. We investigate the response of this particular storm to atmospheric warming (i.e., increase in air temperature) and wetting (i.e., increase in atmospheric moisture content) based on a series of convection-permitting model simulations. Our results show non-monotonic changes of the space-time rainfall variability to either increased temperature or atmospheric moisture content. The most extreme rain rate is produced when relative humidity is increased by 20~40 % or temperature is increased by less than 2°C. The non-monotonic rainfall response is more clearly revealed at fine spatial (100-1000 km 2 ) and temporal scales (less than 6 hours) rather than over the entire domain (~10 4 km 2 ) and aggregated over the storm duration (around two days). This is mainly attributable to the distinct feedbacks from atmospheric dynamics (i.e., moisture convergence and interaction with regional topography) rather than regulated by thermodynamic changes alone. Furthermore, atmospheric warming poses notable changes in the vertical structure of storm cells, contributing to reduced areal reduction factors at small spatial scales and short durations, while atmospheric wetting additionally modifies storm evolution properties. Our modeling analyses challenge the existing practices for hydrologic designs under a changing climate, highlighting particular vulnerability for cities or small basins to short-duration rainfall extremes and the resultant flash flood hazards.

54 ENVIRONMENTAL SCIENCES↗