Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial Characterization”

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 1,315 records · Page 73

ECCI based characterization of dislocation shear in polycrystalline arrays during heterogeneous deformation of commercially pure titanium

A new electron channeling contrast imaging (ECCI) based approach for the characterization of heterogeneous plastic deformation across multiple grains in a polycrystal is presented. The nature of the shear accommodation at grain boundaries has been assessed by evaluating the spatial distribution of dislocations and slip system activities near the boundaries and relating these observations to standard slip transfer criteria and the global state of stress. The direction of slip band propagation, and consequently the location of dislocation slip band nucleation, was determined by characterizing the slip band widths across individual grains. This assessment leads to a better understanding of the factors involved in polycrystalline deformation, including the grain boundary shear accommodation associated with both dislocation pile-ups and grain boundary dislocation nucleation. Because ECCI allows study of dislocation structures over large surface areas, and thus larger overall volumes than typically are accessible by other techniques, the overall approach allows the dislocation distribution across individual grains to be understood in terms of the broader polycrystalline constraints.

36 MATERIALS SCIENCE↗

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Integrating very-high-resolution UAS data and airborne imaging spectroscopy to map the fractional composition of Arctic plant functional types in Western Alaska

Widespread changes in vegetation cover and composition are driving strong impacts on Arctic ecosystem functioning and global climate feedbacks. An accurate characterization of tundra vegetation composition is required to understand how the Arctic will respond to future climate change. However, quantifying tundra vegetation composition over large areas is challenging as commonly-used satellite observations are too coarse, spatially and spectrally, to differentiate low-lying tundra vegetation types. Recent airborne and spaceborne imaging spectroscopy platforms provide better data to characterize vegetation composition. Yet, our ability to characterize vegetation composition with imaging spectroscopy remains largely unexplored in the Arctic, particularly due to a lack of ground observations needed to train and test classification models. To address this problem, we collected very-high-resolution (VHR, ~5 cm) unoccupied aerial system (UAS) imagery at three low-Arctic tundra sites located on the Seward Peninsula, western Alaska. In this paper, we examine the feasibility of integrating imagery from the UAS and the hyperspectral Airborne Visible/Infrared Imaging Spectrometer, Next Generation (AVIRIS-NG) airborne instrument to map the fractional composition of 12 key Arctic plant functional types (PFTs). To this end, we first mapped the 12 PFTs from our VHR UAS imagery using random forest classification. We then used these UAS-derived PFT maps as ground truth to develop partial least squares regression (PLSR) models to predict the fractional cover (FCover) of each PFT from AVIRIS-NG imagery. Further, we evaluated the performance of our PLSR models using reserved UAS samples, as well as by mapping PFT FCover and dominant PFT for large tundra landscapes. Our results show that 1) Arctic PFTs can be effectively mapped using VHR UAS imagery, with overall accuracy between 86% and 92%, 2) when the UAS mapped PFTs were used to inform PLSR scaling models, the FCover of the 12 PFTs could be effectively estimated from AVIRIS-NG imagery with a mean absolute error (MAE) <0.13, and 3) our PLSR models outperformed traditional, fully constrained least-squares (FCLS) linear mixture analysis and produced high-quality, spatially contiguous PFT FCover and PFT maps that captured vegetation spatial patterns with similar accuracy to those developed from UAS imagery. The developed PLSR models have the potential to be broadly applied for quantifying vegetation composition with AVIRIS-NG images to help monitor tundra vegetation dynamics and improve process-based modeling of tundra ecosystems.

54 ENVIRONMENTAL SCIENCES↗

A multi-criteria CCUS screening evaluation of the Gulf of Mexico, USA

Continued research into reservoir characterization along with offshore carbon dioxide (CO 2 ) transportation and infrastructure assets is needed to facilitate development of safe and successful carbon capture, utilization, and storage (CCUS) projects. This paper outlines a multi-criteria evaluation methodology that incorporates disparate sets of quantitative, spatially variable data into a decision-making framework for screening the Gulf of Mexico (GOM) outer continental shelf (OCS) for potentially viable CO 2 storage and enhanced oil recovery (EOR) sites. Criteria categories include favorable geologic characteristics, logistics, and potential risks. Data compiled for 14 criteria from several publicly available geographic information system (GIS) layers was aggregated over 2559 spatially balanced points across the study area using the National Energy Technology Laboratory (NETL)-developed Cumulative Spatial Impact Layers™ (CSIL) GIS tool. Criteria are weighted by qualitative expert opinion relative to their perceived importance to given scenarios— the output of combined criteria values and weights enables regional CO 2 storage suitability differentiation. The methodology considers both technical and non-technical factors impacting CCUS decision-making. The flexible methodology enables a systematic approach to regional ranking at high spatial resolution over a large study domain. Additionally, the framework enables high-grading of priority sites that warrant further characterization and follow-on analysis. Areas along the Louisiana coast and Mississippi River Delta consistently rank high for all scenarios largely a result of the favorable geology with the potential for stacked storage, as well as the density of existing pipelines and platforms, and proximity to several onshore CO 2 sources. High-graded regions for the CO 2 EOR-related scenarios are typically located further offshore towards the middle and edge of the OCS compared to higher priority regions for the geologic storage scenarios which fall closer to the Louisiana coastline.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity analysis for the dynamic response of viscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.

Kulkarni, Makarand↗

Heliostat Consortium: Gap Analysis on Wind Load for Achieving a Fully Competitive Heliostat Industry

The Heliostat Consortium (HelioCon) Wind Load Subtask was initiated with the aim of bringing research work pertaining to wind load measurement, characterization, and prediction taking place across several tasks, including Advanced Manufacturing, Components and Controls, and Field Deployment. The cross-cutting wind load subtopic in the HelioCon roadmap report [1] highlighted standardized methods and tools that are needed for a more detailed understanding of the static and dynamic loads on a heliostat. This will enable cost reduction of wind-dependent heliostat components to avoid unnecessarily conservative, overly constrained designs and increase field efficiency/reliability, to reduce the risk of component failures due to high-wind events (>15 m/s). Gaps related to heliostat wind load include site characterization for wind measurements, critical load cases for heliostat design, turbulence impacts on heliostat tracking error, testing of heliostat array configurations, understanding spatial variation of maximum loads across the solar field, and heliostat field layout and operating strategies. The recommended highest-priority pathway as first steps taken by HelioCon to address these gaps are to develop site characterization guidelines for heliostat design and field load measurements.

Emes, Matthew (ORCID:0000000341474387)↗

Potential for Small Unmanned Aircraft Systems Applications for Identifying Groundwater-Surface Water Exchange in a Meandering River Reach: sUAS in Groundwater-Surface Water Study

The exchange of groundwater and surface water (GW-SW), including dissolved constituents and energy, represents a critical yet challenging characterization problem for hydrogeologists and stream ecologists. Here we describe the use of a suite of high spatial resolution remote sensing techniques, collected using a small unmanned aircraft system (sUAS), to provide novel and complementary data to analyze GW-SW exchange. sUAS provided centimeter-scale resolution topography and water surface elevations, which are often drivers of exchange along the river corridor. Additionally, sUAS-based vegetation imagery, vegetation-top elevation, and normalized difference vegetation index mapping indicated GW-SW exchange patterns that are difficult to characterize from the land surface and may not be resolved from coarser satellite-based imagery. We combined these data with estimates of sediment hydraulic conductivity to provide a direct estimate of GW “shortcutting” through meander necks, which was corroborated by temperature data at the riverbed interface.

54 ENVIRONMENTAL SCIENCES↗

Optical emission spectroscopy and imaging of low-pressure N2 plasmas generated by intense fast-pulsed electron beams

An optical emission spectroscopic (OES) and imaging characterization is conducted on N2 plasmas generated by a 100 keV fast-pulsed electron beam. The electron beams are injected into an N2 gas filled volume with a current of 4.5 kA (300 A/cm2) and a 100 ns pulse width. The characterization is conducted at the pressures, 1 Torr and 0.1 Torr, corresponding to two distinct regimes that exhibit significantly different plasma dynamics. Beam impact ionization is shown to be a primary mechanism for producing low temperature plasmas at 1 Torr during beam output. After beam termination, ionization by an inductive electric field becomes the primary mechanism for plasma formation later in time for both pressures. OES and plasma imaging are used in this work as a diagnostic tool to track the distribution of electronic, vibrational, and rotational state transitions and ionized species. This is achieved with the use of a multi-resolution suite of spectrometers capable of acquiring time-resolved spectra. Vibrational and rotational bands of the N2 second positive system (C3Πu→B3Πg) and the N2+ first negative system (B2Σu+→X2Σg+) are identified in both regimes as well as N+ states exclusively in the lower pressure regime. Vibrational and rotational spectra are shown to track the evolution of the time-varying plasma current. Plasma imaging also reveals spatially nonuniform plasma emission at 0.1 Torr. A few hundred shots are recorded to fully characterize the emissions, and results are shown to be highly reproducible (≤±1% with 2σ confidence).

Kaiser, E. R. (ORCID:0000000336493938)↗

Spatial Coverage Planning for Exploration Robots

A report discusses an algorithm for an onboard planning and execution technology to support the exploration and characterization of geological features by autonomous rovers. A rover that is capable of deciding which observations are more important relieves the engineering team from much of the burden of attempting to make accurate predictions of what the available rover resources will be in the future. Instead, the science and engineering teams can uplink a set of observation requests that may potentially oversubscribe resources and let the rover use observation priorities and its current assessment of available resources to make decisions about which observations to perform and when to perform them. The algorithm gives the rover the ability to model spatial coverage quality based on data from different scientific instruments, to assess the impact of terrain on coverage quality, to incorporate user-defined priorities among subregions of the terrain to be covered, and to update coverage quality rankings of observations when terrain knowledge changes. When the rover is exploring large geographical features such as craters, channels, or boundaries between two different regions, an important factor in assessing the quality of a mission plan is how the set of chosen observations spatially cover the area of interest. The algorithm allows the rover to evaluate which observation to perform and to what extent the candidate observation will increase the spatial coverage of the plan.

Gaines, Daniel↗

Characterizing and Improving the Data Reduction Pipeline for the Keck OSIRIS Integral Field Spectrograph

OSIRIS is a near-infrared (1.0–2.4 μm) integral field spectrograph operating behind the adaptive optics system at Keck Observatory and one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP, including characterization of the flux assignment algorithm, exploration of spatial rippling in the reduced data cubes, and improvements to several calibration files, including the rectification matrix, bad-pixel mask, and wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved Band-to-Band Registration Characterization for VIIRS Reflective Solar Bands Based on Lunar Observations

Spectral bands of the Visible Infrared Imaging Radiometer Suite (VIIRS) instrumentaboard the Suomi National Polar-orbiting Partnership (S-NPP) satellite are spatially co-registered.The accuracy of the band-to-band registration (BBR) is one of the key spatial parameters that must becharacterized. Unlike its predecessor, the Moderate Resolution Imaging Spectroradiometer (MODIS), VIIRS has no on-board calibrator specifically designed to perform on-orbit BBR characterization.To circumvent this problem, a BBR characterization method for VIIRS reflective solar bands (RSB) based on regularly-acquired lunar images has been developed. While its results can satisfactorily demonstrate that the long-term stability of the BBR is well within +/- 0.1 moderate resolution bandpixels, undesired seasonal oscillations have been observed in the trending. The oscillations are most obvious between the visiblenear-infrared bands and short-middle wave infrared bands. This paper investigates the oscillations and identifies their cause as the band spectral dependence of the centroid position and the seasonal rotation of the lunar images over calibration events. Accordingly, an improved algorithm is proposed to quantify the rotation and compensate for its impact. After the correction, the seasonal oscillation in the resulting BBR is reduced from up to 0.05 moderate resolution band pixels to around 0.01 moderate resolution band pixels. After removing this spurious seasonal oscillation, the BBR, as well as its long-term drift are well determined.

Wang, Zhipeng↗

Binary pseudo-random array for calibration of interferometers with transmission spheres and cylinders

Binary pseudo-random array (BPRA) “white noise” artifacts are highly effective for characterizing the instrument transfer function (ITF) of surface topography metrology tools and wavefront measurement instruments. These BPRA artifacts feature all spatial frequencies within the instrument bandpass equally, resulting in a power spectral density with a white-noise-like character. This characteristic allows for direct ITF determination with uniform sensitivity across the entire spatial frequency range. We have developed a novel BPRA calibration standard that combines the diffractive pattern of a reflection computer-generated hologram (CGH) with the white noise generating BPRA pattern. By integrating these technologies using the same lithographic techniques, the resulting calibration sample enables ITF characterization of a Fizeau interferometer with a transmission sphere, or any nulling optic.

Munechika, K↗

Assessment of the feedstock supply for siting single- and multiple-feedstock biorefineries in the USA and identification of prevalent feedstocks

An integrated multi-feedstock bioenergy (i.e., biofuel, biopower, or bioproduct) supply system has potential to reduce biomass supply system uncertainties and costs. This study identifies optimal configurations of multi-feedstock biomass-to-biorefinery supply chains and pertinent feedstock combinations based on spatial distribution of feedstock and lowest delivered cost to the biorefinery. We used the Supply Characterization Model (SCM) to allocate feedstock supplies to candidate biorefinery facilities. Model runs were performed for herbaceous energy crops, agriculture residue, and woody biomass available in 2017, 2022, 2025, and 2030 as estimated by the Policy Analysis System (POLYSYS) and Forest Sustainable and Economic Analysis Model (ForSEAM) models. Three feedstock supply scenarios were compared: (a) an herbaceous scenario: switchgrass, miscanthus, biosorghum, and corn stover; (b) a woody scenario: coppice wood, noncoppice wood, whole trees, and forestry residues, and (c) a mixed scenario: a combination of all feedstocks in herbaceous and woody scenarios. By 2030 the analyses predicted that 323, 168, and 473 biorefineries were sited in the herbaceous, woody, and mixed scenario, respectively, in the conterminous USA. Feedstock mixes supplied to the biorefineries were mostly dominated by a single feedstock. The most prominent feedstock mixes identified were: (1) switchgrass and miscanthus; (2) coppice and noncoppice wood; and (3) coppice wood, noncoppice wood, switchgrass and miscanthus. Biorefineries using multi-feedstock would be beneficial for growth of bioeconomy, however flexible and cost-effective conversion platforms should be developed to efficiently utilize multiple feedstocks. This analysis identifies biorefinery locations and feedstock supply mixes while minimizing delivered feedstock costs based on spatial and temporal feedstock availability.

09 BIOMASS FUELS↗

Visualization of the surface distribution of photocatalytic activity at the anatase/rutile TiO2 interface using X-ray photoelectron emission microscopy

We applied X-ray photoelectron emission microscopy (XPEEM) to visualize the surface distribution of photocatalytic activity at the anatase/rutile (A/R) interface of titanium dioxide. Acetic acid was adsorbed on the surface, and its decomposition and desorption reactions under ultraviolet light irradiation were monitored by observing C 1s spectra, enabling direct observation of the spatial distribution of photocatalytic activity. Complementary crystal structure information was obtained using spatially resolved low-energy electron diffraction, low-energy electron microscopy, X-ray absorption spectroscopy, and micro-focused Raman spectroscopy. Comparison with these structural data allowed for the evaluation of the photocatalytic activity as a function of the distance from the A/R interface. The activity is enhanced in the vicinity of the interface and gradually decreases toward both the anatase and rutile regions. These results demonstrate that XPEEM is an effective technique for spatially resolved mapping of photocatalytic activity, which has\\r\\npreviously been inaccessible using conventional characterization techniques.

25 ENERGY STORAGE↗

Effect of time and spatial domains on monolayer 2D material interface thermal conductance measurement using ns ET-Raman

Raman-based approaches to characterizing thermal properties of nanoscale and microscale materials have proven to be powerful experimental methods. Here, this work explores the potential of energy-transport state-resolved Raman (ET-Raman) for the precise determination of interfacial thermal conductance (G) between monolayer 2D materials and their substrates. ET-Raman employs a continuous wave laser and nanosecond (ns) pulsed laser to induce both steady state and transient heating while simultaneously collecting 2D material temperature information from the Raman signal. A 3D numerical model built using the finite volume method simulates the two differentiated thermal transport processes. The simulated temperature rises are compared with experimental Raman temperature information to extract G. In this study, we investigate the temporal effects of the pulse frequency and pulse width of the ns laser on precise characterization of G. By running multiple simulations for a range of pulsed laser frequencies and pulse widths, we determine that ET-Raman is most sensitive to precise interfacial thermal conductance measurements with short pulse widths (~below 20 ns) and with pulsing frequencies that ensure sufficient cooling time between consecutive pulses (~below 500 kHz). We also conclude that substrate thickness and laser spot size play critical roles in heat distribution and subsequent sensitivity to meaningful G measurements. As both laser spot size and SiO 2 substrate thickness decrease, detection of interfacial thermal transport becomes more difficult. These results offer essential guidance in the design of future Raman-based approaches to thermal measurements using pulsed lasers.

2D materials↗

Multiscale models and stochastic simulation methods for computing rare but key binding events in cell biology

The main difficulty in simulating diffusion processes at a molecular level in cell microdomains is due to the multiple scales involving nano- to micrometers. Few to many particles have to be simulated and simultaneously tracked while there are exploring a large portion of the space for binding small targets, such as buffers or active sites. Bridging the small and large spatial scales is achieved by rare events representing Brownian particles finding small targets and characterized by long-time distribution. These rare events are the bottleneck of numerical simulations. A naive stochastic simulation requires running many Brownian particles together, which is computationally greedy and inefficient. Solving the associated partial differential equations is also difficult due to the time dependent boundary conditions, narrow passages and mixed boundary conditions at small windows. We present here two reduced modeling approaches for a fast computation of diffusing fluxes in microdomains. The first approach is based on a Markov mass-action law equations coupled to a Markov chain. The second is a Gillespie's method based on the narrow escape theory for coarse-graining the geometry of the domain into Poissonian rates. The main application concerns diffusion in cellular biology, where we compute as an example the distribution of arrival times of calcium ions to small hidden targets to trigger vesicular release.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Identification and correction of temporal and spatial distortions in scanning transmission electron microscopy

Scanning transmission electron microscopy (STEM) has become the technique of choice for quantitative characterization of atomic structure of materials, where the minute displacements of atomic columns from high-symmetry positions can be used to map strain, polarization, octahedra tilts, and other physical and chemical order parameter fields. The latter can be used as inputs into mesoscopic and atomistic models, providing insight into the correlative relationships and generative physics of materials on the atomic level. However, these quantitative applications of STEM necessitate understanding the microscope induced image distortions and developing the pathways to compensate them both as part of a rapid calibration procedure for in situ imaging, and the post-experimental data analysis stage. Here, we explore the spatiotemporal structure of the microscopic distortions in STEM using multivariate analysis of the atomic trajectories in the image stacks. Based on the behavior of principal component analysis (PCA), we develop the Gaussian process (GP)-based regression method for quantification of the distortion function. The limitations of such an approach and possible strategies for implementation as a part of in-line data acquisition in STEM are discussed. Here, the analysis workflow is summarized in a Jupyter notebook that can be used to retrace the analysis and analyze the reader's data.

36 MATERIALS SCIENCE↗

Time of Emergence and Large Ensemble Intercomparison for Ocean Biogeochemical Trends

Anthropogenically forced changes in ocean biogeochemistry are underway and critical for the ocean carbon sink and marine habitat. Detecting such changes in ocean biogeochemistry will require quantification of the magnitude of the change (anthropogenic signal) and the natural variability inherent to the climate system (noise). Here we use Large Ensemble (LE) experiments from four Earth system models (ESMs) with multiple emissions scenarios to estimate Time of Emergence (ToE) and partition projection uncertainty for anthropogenic signals in five biogeochemically important upper-ocean variables. We find ToEs are robust across ESMs for sea surface temperature and the invasion of anthropogenic carbon; emergence time scales are 20–30 yr. For the biological carbon pump, and sea surface chlorophyll and salinity, emergence time scales are longer (50+ yr), less robust across the ESMs, and more sensitive to the forcing scenario considered. We find internal variability uncertainty, and model differences in the internal variability uncertainty, can be consequential sources of uncertainty for projecting regional changes in ocean biogeochemistry over the coming decades. In combining structural, scenario, and internal variability uncertainty, this study represents the most comprehensive characterization of biogeochemical emergence time scales and uncertainty to date. Our findings delineate critical spatial and duration requirements for marine observing systems to robustly detect anthropogenic change.

54 ENVIRONMENTAL SCIENCES↗