Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak finding”

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 865 records · Page 48

Quasinormal modes of charged fields in Reissner-Nordström backgrounds by Borel-Padé summation of Bender-Wu series

We extend recent work of Hatsuda on the computation of quasinormal mode frequencies via analytic continuation of bound state energies and Borel-Pade resummation of the Bender-Wu perturbation series to the case of charged fields in the background of Reissner-Nordstrom black holes. We compare the quasinormal mode frequencies obtained in this manner to calculations using Leaver's method of continued fractions, and find good agreement for damped modes (DMs) with imaginary part remaining finite in the extremal limit. Here, we also present numerical evidence that the frequencies of certain zero-damped modes (ZDMs) with imaginary part tending to zero in the extremal limit can be computed when constructing the Bender-Wu expansion about a peak of the potential inside the outer horizon of the black hole.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DESI Emission-line Galaxies: Unveiling the Diversity of [O II ] Profiles and Its Links to Star Formation and Morphology

We study the [O II ] profiles of emission-line galaxies (ELGs) from the Early Data Release of the Dark Energy Spectroscopic Instrument (DESI). To this end, we decompose and classify the shape of [O II ] profiles with the first two eigenspectra derived from principal component analysis. Our results show that DESI ELGs have diverse line profiles, which can be categorized into three main types: (1) narrow lines with a median width of ∼50 km s −1 , (2) broad lines with a median width of ∼80 km s −1 , and (3) two redshift systems with a median velocity separation of ∼150 km s −1 , i.e., double-peak galaxies. To investigate the connections between the line profiles and galaxy properties, we utilize the information from the COSMOS data set and compare the properties of ELGs, including star formation rate (SFR) and galaxy morphology, with the average properties of reference star-forming galaxies with similar stellar mass, sizes, and redshifts. Our findings show that, on average, DESI ELGs have a higher SFR and more asymmetrical/disturbed morphology than the reference galaxies. Moreover, we uncover a relationship between the line profiles, the excess SFR, and the excess asymmetry parameter, showing that DESI ELGs with broader [O II ] line profiles have more disturbed morphology and higher SFR than the reference star-forming galaxies. Finally, we discuss possible physical mechanisms giving rise to the observed relationship and the implications of our findings on the galaxy clustering measurements, including the halo occupation distribution modeling of DESI ELGs and the observed excess velocity dispersion of the satellite ELGs.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing Hawking radiation through capacity of entanglement

We consider the capacity of entanglement in models related with the gravitational phase transitions. The capacity is labeled by the replica parameter which plays a similar role to the inverse temperature in thermodynamics. In the end of the world brane model of a radiating black hole the capacity has a peak around the Page time indicating the phase transition between replica wormhole geometries of different types of topology. Similarly, in a moving mirror model describing Hawking radiation the capacity typically shows a discontinuity when the dominant saddle switches between two phases, which can be seen as a formation of island regions. In either case we find the capacity can be an invaluable diagnostic for a black hole evaporation process.

2D gravity↗

Geological study of an outburst flood event in the upper Yangtze River and risk of similar extreme events

Here, outburst flood sediments are studied on the upper course of the Yangtze River (known as the Jinsha River) where is a key region for the development of hydroelectric power. As the history and magnitude of outburst flooding in this stretch of the Jinsha River during historic times has remained largely undocumented, we have endeavored as part of this study to conduct field investigations into, and provide numerical dating (using OSL and 14 C dating techniques) of outburst flood evidence found in the First Bend area of the Yangtze River. Carbonized nutshells and charcoal from the ancient settlement yielded 14 C ages that precisely pinpointed an outburst flood event ~1200–1300 cal yr BP, during China's Tang Dynasty (618–907 CE). However, OSL ages using the minimum age model for all outburst flood samples ranged from 9.1 to 1.8 ka, which was ~0.5–7.8 ka older than the 14 C-dated outburst flood. This discrepancy suggests a potential overestimation of quartz age when dating outburst flood sediments from the Holocene period. Furthermore, using 2D hydraulic modeling, we estimated the magnitude of the aforementioned outburst flood, revealing a reconstructed peak flow discharge of 51,000–55,300 m 3 /s. This discharge was three times greater than the 10,000-year flood for the hydroelectric power station nearest to the study area. Our study underscores the significant role of combining archaeological and geological evidence in enhancing paleoflood hydrology research. Through a multidisciplinary approach, our findings emphasize the importance of paleoflood hydrology studies in comprehending both the latent and catastrophic flood risks in alpine mountain areas amid the context of global warming.

14C dating↗

Holocene paleoclimate change in the western US: The importance of chronology in discerning patterns and drivers

We report sediment in lakes and meadows forms a powerful archive that can be used to reconstruct environmental change through time. Reconstructions of lake level, of chemical, biological, and hydrological conditions, and of surrounding vegetation provide detailed information about past climate conditions, both locally and regionally. Indeed, most of our current knowledge of centennial- to millennial-scale climate variability in the arid western United States, where information about past hydroclimate is particularly important, comes from such sediment-based reconstructions. The pressing need for robust, precise predictions of future conditions is a significant motivation for paleoclimate science, and current research questions frequently require Holocene reconstructions to be resolved at sub-centennial timescales. Increasingly, regional syntheses seek to identify synoptic-scale patterns similar to those defined from modern observations (seasonal, interannual, multi-decadal, etc.) or to compare with the output of climate model simulations. However, the age control on existing records, especially those more than about 20 years old, is often sufficient only for millennial-scale interpretation. Here we assess the age control for 84 published and unpublished records from lakes and meadows in the Great Basin, California, and desert southwest, and use Bayesian modeling to evaluate the 95% uncertainty ranges for the 42 best-dated records. In the Late Holocene, about half of the 42 records have <400-year mean uncertainty ranges; however, high-precision age control is especially critical for young records, used to develop an accurate understanding of a proxy’s response to known climate variations. In the Middle Holocene, records vary from 400 to >800-year mean uncertainty and records of the Early Holocene have 600- to >1400-year mean uncertainty ranges. We find that the largest control on modeled uncertainties is dating density, with at least 2 dates/kyr being optimal and suggest obtaining “range-finder” dates at the onset of a study to better predict the total number of dates needed for an adequate age model. Such a density avoids a commonly observed phenomenon of significant peaks in uncertainty arising in gaps between age control points. Analysis of the uncertainties associated with proxy shifts reveal that more than half are >400 years. Although such large uncertainties currently prevent sub-centennial interpretations in most cases, increased dating density, strategic use of limited funds (including budgeting for a 2 date/kyr minimum at the proposal stage), construction of age-depth models with Bayesian methods, and critical evaluation of chronological uncertainty will shed light on past climate variability at finer timescales, enhancing our understanding of global and regional drivers of western U.S. climate.

54 ENVIRONMENTAL SCIENCES↗

Partitioning of Ca to metastable precipitates in a Mg-rare earth alloy

The potential effect of the element Ca on the precipitation behavior was investigated in an Mg-rare earth alloy. A combination of metastable β''' and β' precipitates was observed for the peak aging condition at 200°C. Ca addition was found to have no significant effect on the precipitating phases and evolution sequence. Composition analysis showed that the Ca partitioned to both β''' and β' precipitate phases. First-principles calculations indicated that Ca partitions to the rare-earth sublattice in the precipitate phase. This finding suggests the potential of Ca to partially replace costly rare-earth elements in precipitation-hardened Mg-rare earth alloys.

36 MATERIALS SCIENCE↗

Drift kinetic effects on plasma response to resonant magnetic perturbation for EU DEMO design

Abstract A systematic investigation of the plasma response to an externally applied resonant magnetic perturbation (RMP) field, for the purpose of controlling edge localized modes, is carried out for an EU DEMO reference plasma. Particular emphasis is placed on the role of kinetic effects associated with both thermal particles and fusion-born alphas. The single fluid, resistive model predicts a large peak amplification of the n = 1 ( n is the toroidal mode number) plasma response to the target equilibrium, which is found to be close to the Troyon no-wall limit. A more advanced response model, including kinetic resonances between the RMP perturbation and drift motions of thermal and energetic particles, on the other hand, finds a strong suppression of the n = 1 field amplification. A major role is played by the precessional drift resonance of fusion-born alphas. A strong parallel sound wave damping model is found to well reproduce the full kinetic response results for the DEMO plasma, in terms of both the resonant field response amplitude and the plasma displacement. Finally, both fluid and kinetic models produce similar responses for the n = 2 and 3 RMP fields for the considered DEMO plasma, whilst kinetic effects again become important for the n = 4 RMP due to proximity of the reference plasma to the no-wall limit for the n = 4 ideal kink instability.

Physics↗

Fine-scale vegetation composition and structure shape spatiotemporal variation in surface albedo across a low Arctic tundra landscape

The unprecedented rate of warming in the Arctic is driving changes in the structure and composition of tundra vegetation. Increases in deciduous tall shrub cover, height, and density are of particular concern, as these changes alter local surface albedo in ways that could amplify effects on the regional surface energy budget (SEB). Despite this importance, significant uncertainties remain in understanding the interplay between fine-scale vegetation patterns and emergent albedo dynamics across space and time. Here, we address these uncertainties by (1) quantifying spatiotemporal variation in surface shortwave albedo and (2) determining the relative influence of fine-scale vegetation composition, structure, and environmental conditions on albedo across a representative low-Arctic tundra landscape on Alaska’s Seward Peninsula. To do this, we synthesized multi-scale, multi-platform remote sensing observations, including a novel Landsat-derived albedo time series, a fine-scale map of Arctic plant functional type (PFT) fractional cover, and airborne LiDAR estimates of canopy height and topography. We show that there are substantial reductions in winter albedo for pixels dominated by tall, woody PFTs (28.13%) relative to pixels dominated by non-woody vegetation, but almost no change in summer albedo (3% increase). Further, we identified a unimodal trend in the relationship between canopy height and the timing of the springtime transition from high (snowy) to low (leafy) albedo (peak at 5.5 m), possibly because of competing ‘snow-fence’ and ‘protrusion’ snow-shrub interactions. To explore the primary drivers of albedo, we constructed a random forest model and found that canopy height and the fractional cover of woody PFTs were as- or more important predictors of winter albedo than topographic features. These findings provide strong evidence for the impacts of local vegetation characteristics on regional surface albedo, highlighting the need for better quantification of snow-shrub interactions to accurately predict the Arctic’s SEB under future environmental change.

Arctic↗

Final-state interactions and spin structure in E 1 breakup of Li 11 in halo effective field theory

In this work, we calculate the E1 breakup of the 2n halo nucleus 11 Li in halo effective field theory (Halo EFT) at leading order. In Halo EFT, 11 Li is treated as a three-body system of a 9 Li core and two neutrons. We present a detailed investigation of final-state interactions (FSIs) in the neutron-neutron (nn) and neutron-core (nc) channels. We employ Møller operators to formulate an expansion scheme that satisfies the non-energy-weighted cluster sum rule and successively includes higher-order terms in the multiple-scattering series for the FSI. Computing the E1 strength up to third order in this scheme, we observe apparent convergence and good agreement with experiment. The neutron-neutron FSI is by far the most important contribution and largely determines the maximum value of the E1 distribution. However, inclusion of nc FSI does shift the peak position to slightly lower energies. Moreover, we investigate the sensitivity of the E1 response to the spin structure of the neutron- 9 Li interaction. We contrast results for an interaction that is the same in the spin-1 and spin-2 channels with one that is only operative in the spin-2 channel, and find that good agreement with experimental data is only obtained if the interaction is present in both spin channels. The latter case is shown to be equivalent to a calculation in which the spin of 9 Li is neglected.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of an Efficient Conjugate Heat Transfer Modeling Framework to Optimize Mixing-Limited Combustion of Ethanol in a Diesel Engine

Mixing controlled combustion of alcohol fuels has been identified as a promising technology based on their low propensity for particulate and NOx production, but the higher heats of vaporization and auto-ignition temperatures of these fuels make their direct use in diesel engine architectures a challenge. To realize the potential of alcohol-fueled combustion, a computational fluid dynamics (CFD) modeling framework is developed, validated, and exercised to identify designs that maximize engine thermal efficiency. To evaluate the use of thermal barrier coatings (TBCs), a simplified one-dimensional (1D) conjugate heat transfer (CHT) modeling framework is employed. The addition of the 1D CHT model only increases the computational expense by 15% relative to traditional approaches, yet offers more accurate heat transfer predictions over constant temperature boundary conditions. The validated model is then used to explore a range of injector orientations and piston bowl geometries. Using a design of experiments (DoE) approach, several designs were identified that improved fuel-air mixing, shortened the combustion duration, and increased thermal efficiency. The most promising design was fabricated and tested in a Caterpillar 1Y3700 single-cylinder oil test engine (SCOTE). Engine testing confirmed the findings from the CFD simulations and found that the co-optimized injector and piston bowl design yielded over 2-percentage point increase in thermal efficiency at the same equivalence ratio (0.96) and over 6-percentage point increase at the same engine load (10.1 bar indicated mean effective pressure (IMEP)), while satisfying design constraints for peak pressure and maximum pressure rise rate.

Magnotti, Gina M.↗

ASSESSING THE EFFECTIVENESS OF ULTRASONIC IMPACT TREATMENT ON RESIDUAL STRESS PROFILES IN DISSIMILAR WELDED JOINTS

Residual stresses (RS) induced during welding processes are a critical concern in materials engineering, as they can significantly impair the mechanical performance of components by reducing fatigue strength and tensile load capacity. This challenge is especially pronounced in dissimilar metal welds (DMWs), where variations in thermal expansion properties between the joined alloys exacerbate the formation of tensile RS. Conventional post-weld heat treatments, though effective for homogeneous materials, often require substantial energy, specialized equipment, and extensive processing time, making them less practical for DMW applications. Thus, there is a clear need for innovative, energy-efficient techniques to mitigate these detrimental stresses. This study investigates ultrasonic impact treatment (UIT) as a possible alternative for mitigating tensile RS in both similar and dissimilar metal welds. To evaluate UIT’s effectiveness, neutron diffraction (ND) was employed as a nondestructive technique to quantify RS in three orthogonal directions—longitudinal, transverse, and normal. The results showed that UIT significantly reduced peak tensile RS, particularly in the longitudinal direction, by up to 180 MPa in similar welds and up to 150 MPa in dissimilar welds. Given the limited literature on UIT application in DMWs, this work contributes valuable data on stress redistribution mechanisms and highlights UIT’s potential as a practical stress-relief method. The findings lay the groundwork for further investigations aimed at optimizing process parameters and understanding long-term performance in welded joints.

EisaZadeh, Hamid [Western Carolina University, Cul↗

Estimation of Surface Sensible Heat Flux due to Precipitation over CONUS and Its Impact on Urban Extreme Precipitation Modeling

The surface sensible heat flux induced by precipitation (Q P ) is a consequence of the temperature difference between the surface and the rain droplets. Despite its seemingly negligible nature, Q P is frequently omitted from both meteorological and climatological models. Nevertheless, it is important to acknowledge the numerous occasions in which the instantaneous values of Q P can be significant, particularly during extreme precipitation events. This study undertakes a comprehensive assessment of Q P across the contiguous United States (CONUS) utilizing high-resolution reanalysis, observational data, and numerical modeling to examine the influence of Q P on precipitation and the surface energy budget. The findings indicate that the spatial distribution of Q P climatology is analogous to that of precipitation, with magnitudes ranging from 2 to 3 W m -2 predominantly over the Midwest and Southeast regions. A seasonal analysis of QP reveals that the highest values occurring during the June–August (JJA) period, averaging 3.18 W m -2 . Peak Q P values of approximately 4 W m -2 are observed during JJA over the Great Plains region. We hypothesize that the Q P during an extreme precipitation event would be nonnegligible and have a significant impact on the local weather. To test this conjecture, we perform high-resolution simulations with and without Q P during an extreme precipitation event over the Chicago Metropolitan Area (CMA). The results show that the Q P may be a dominant factor compared to other components of surface heat flux during the zenith of precipitation hours. Also, Q P has the potential to not only diminish precipitation but also alter and reconfigure the remaining surface energy budget components.

54 ENVIRONMENTAL SCIENCES↗

Ultimate Precision on Sin^2 Theta_23 With the Numi Experiments

There were two major goals to this proposal: the study of the discrepancy between data and MC in the NuMI Medium Energy neutrino beam and the incorporation of the MINOS and NOVA neutrino oscillation datasets. The first topic was approached via development of a simulation of the NuMI beam monitoring instrumentation suite. Initially, the Horn-Off current data sets were used to measure a set of hadron production parameters which improved the data/MC agreement in a way which was independent of the focusing parameters. The peak position discrepancy was only explained by a large (10kA) reduction in the horn current. A full simulation of the muon monitors was developed to attempt to make an independent measurement of the horn current by means of differential data taken by those monitors. While only at the 5kA level, the result corroborated the former findings. The second goal has not yet been completed but progress has focused on understanding how to modify the MINOS MC to take into account the MEC contribution which is missing from the underlying NEUGEN MC. A number of improvements to the MINOS reconstruction have been necessary. The framework to run the fits on the combined data has been completed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Early-time Light Curves of Type Ia Supernovae Observed with TESS

We present the early-time light curves of Type Ia supernovae (SNe Ia) observed in the first six sectors of Transiting Exoplanet Survey Satellite (TESS) data. Ten of these SNe were discovered by ASAS-SN, seven by ATLAS, six by ZTF, and one by Gaia. For nine of these objects with sufficient dynamic range (>3.0 mag from detection to peak), we fit power-law models and searched for signatures of companion stars. We found a diversity of early-time light-curve shapes, although most of our sources are consistent with fireball models where the flux increases as ∝t 2 . Three SNe displayed a flatter rise with flux ∝t. We did not find any obvious evidence for additional structures, such as multiple power-law components, in the early rising light curves. For assumptions about the SN properties and the observer viewing angle (ejecta mass of 1.4 M⊙, expansion velocity of 10 4 km s –1 , opacity of 0.2 cm 2 g –1 , and viewing angle of 45°) and a further assumption that any companion stars would be in Roche lobe overflow, it is possible to place upper limits on the radii of any companion stars. Six of the nine SNe had complete coverage of the early-time light curves, and we placed upper limits on the radii of companion stars of ≲32 R ⊙ for these SNe, ≲20 R ⊙ for five of the six, and ≲4 R ⊙ for two of the six. The small sample size did not allow us to put limits on the occurrence rate of companion stars in the progenitors of SNe Ia. However, we expect that TESS observed enough SNe in its two-year primary mission (26 sectors) to either detect the signature of a large companion (R > 20 R ⊙ ) or constrain the occurrence rate of such systems, at least for the fiducial SN properties adopted here. Here, we also show that TESS is capable of detecting emission from a 1 R ⊙ companion for an SN Ia within 50 Mpc and has a reasonable chance of doing so after about six years.

79 ASTRONOMY AND ASTROPHYSICS↗

Demand-side solutions in the US building sector could achieve deep emissions reductions and avoid over $100 billion in power sector costs

Buildings are energy-intensive and a primary source of US end-use sector carbon emissions. Although building emissions today are 25% below their 2005 peak, far deeper reductions are needed to reach the US 2050 net-zero emissions goal. However, plausible decarbonization pathways that consider both buildings and their interactions with the power grid remain poorly understood. Here, we couple detailed modeling of building energy use and the grid to quantify building decarbonization potential and associated grid impacts. We find up to a 91% reduction in building CO 2 emissions from 2005 levels by 2050 using a portfolio of building efficiency, demand flexibility, and electrification measures alongside rapid grid decarbonization. Building efficiency and flexibility could generate up to $107 billion in annual power system cost savings by 2050, offsetting over a third of the incremental cost of full grid decarbonization. Our results underscore multiple benefits of demand-side solutions for deep decarbonization of US buildings.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Denoising atomic resolution 4D scanning transmission electron microscopy data with tensor singular value decomposition

Tensor singular value decomposition (SVD) is a method to find a low-dimensional representation of data with meaningful structure in three or more dimensions. Here, tensor SVD has been applied to denoise atomic-resolution 4D scanning transmission electron microscopy (4D STEM) data. On data simulated from a SrTiO 3 [100] perfect crystal and a Si [110] edge dislocation, tensor SVD achieved an average peak signal-to-noise ratio (PSNR) of ~40 dB, which matches or exceeds the performance of other denoising methods, with processing times at least 100 times shorter. On experimental data from SrTiO 3 [100] and LiZnSb [11 2 ¯ 0]/GaSb [110] samples, tensor SVD denoises multiple GB 4D STEM data sets in ten minutes on a typical personal computer. Denoising with tensor SVD improves both convergent beam electron diffraction patterns and virtual-aperture annular dark field images.

36 MATERIALS SCIENCE↗

Disparities in the air quality monitoring stations and PM₂.₅ in Chicago’s air quality landscape

Fine particulate matter (PM₂.₅) poses significant public and environmental health risks in urban areas. Chicago’s dense industry and traffic create variable air quality, yet monitoring is unevenly distributed, resulting in undersampling of air quality data in some city areas. This study applied a hybrid approach using GIS-based kernel density mapping, interpolation modeling (IDW, Spline, Kriging) of USEPA monitoring data, multi-scale temporal trend analyses (hourly to annual), and ESDA. Accordingly, the density surface showed that monitors are concentrated in the affluent north, northwest, and southwest sides of Chicago (up to ~ 0.07 stations per sq mile), while the south and southeast regions, with predominantly minority communities, have virtually no coverage. Overall, citywide coverage is minimal (~ 4–5 monitors total; ~0.02 per sq mile; ≈1 per 600,000 residents). Temporal analyses showed that the city’s mean annual PM₂.₅ (~ 10.8 µg/m³) exceeds USEPA/WHO standards (9 µg/m³), with summer means (~ 17.1 µg/m³) significantly higher than other seasons. Diurnally, a clear pattern was observed, with PM₂.₅ concentrations peaking overnight (00:00–03:00) and during the morning rush hours, and dipping during midday to late afternoon. Spatial distribution of PM₂.₅ identified hotspots near O’Hare Airport, the downtown Loop area, and south-side neighborhoods, contrasting with lower concentrations on the north side, revealing Chicago’s socioeconomic divides and resulting environmental inequities. The findings underscore the need for expanded monitoring and targeted interventions in under-monitored, high-pollution communities to advance equitable community health.

54 ENVIRONMENTAL SCIENCES↗

Large Eddy Simulation of the Diurnal Cycle of Shallow Convection in the Central Amazon

Climate models often face challenges in accurately simulating the daily precipitation cycle over tropical land areas, particularly in the Amazon. One contributing factor may be the incomplete representation of the diurnal evolution of shallow cumulus (ShCu) clouds. This study aimed to enhance the understanding of the diurnal cycles of ShCu clouds—from formation to maturation and dissipation—over the Central Amazon (CAMZ). Using observational data from the Green Ocean Amazon 2014 (GoAmazon) campaign and large eddy simulation (LES) modeling, we analyzed the diurnal cycles of six selected pure ShCu cases and their composite behavior. Our results revealed a well-defined cycle, with cloud formation occurring between 10 and 11 local time (LT), maturity from 13 to 15 LT, and dissipation by 17–18 LT. The vertical extent of the liquid water mixing ratio and the intensity of the updraft mass flux were closely associated with increases in turbulent kinetic energy (TKE), enhanced buoyancy flux within the cloud layer, and reduced large-scale subsidence. We further analyzed the diurnal cycles of the convective available potential energy (CAPE), the convective inhibition (CIN), the Bowen ratio (BR), and the vertically integrated TKE in the mixed layer (ITKE-ML), exploring their relationships with the cloud base mass flux (Mb) and cloud depth across the six ShCu cases. ITKE-ML and Mb exhibited similar diurnal trends, peaking at approximately 14–15 LT. However, no consistent relationships were found between CAPE (or BR) and Mb. Similarly, comparisons of the cloud depth with CAPE, BR, ITKE-ML, CIN, and Mb revealed no clear relationships. Smaller ShCu clouds were sometimes linked to higher CAPE and lower CIN. It is important to emphasize that these findings are preliminary and based on a limited sample of ShCu cases. Further research involving an expanded dataset and more detailed analyses of the TKE budget and synoptic conditions is necessary. Such efforts would yield a more comprehensive understanding of the factors influencing ShCu clouds’ vertical development.

54 ENVIRONMENTAL SCIENCES↗