Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial-temporal correlation”

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 19 records

A Graph Convolutional Network for Active Distribution System Anomaly Detection Considering Measurement Spatial-Temporal Correlations

The accuracy of distribution system state estimation may be significantly impacted by the existence of bad measure-ments and unexpected topology errors. This paper proposes a data-driven Graph Convolutional Network (GCN) for anomaly detection, including bad measurements and topology change events. Compared to many existing machine learning approaches, the proposed approach embeds both spatial-temporal measure-ment correlations, which allows us to detect and distinguish different anomalies. Numerical results carried out on the IEEE 37-node system demonstrate that the proposed-based method can obtain high accuracy in detecting bad data and topology changes as compared to other approaches, even in the presence of high PV penetrations.

active distribution system↗

Comparative Study of Aerosol and Cloud Detected by CALIPSO and OMI

The Ozone Monitoring Instrument (OMI) on the Aura Satellite detects the presence of desert dust and smoke particles (also known as aerosols) in terms of a parameter known as the UV Aerosol Index (UV AI). The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission measures the vertical distribution of aerosols and clouds. Aerosols and clouds play important roles in the atmosphere and climate system. Accurately detecting their presence, altitude, and properties using satellite radiance measurements is a very important task. This paper presents a comparative analysis of the CALIPSO Version 2 Vertical Feature Mask (VFM) product with the (OMI) UV Aerosol Index (UV AI) and reflectivity datasets for a full year of 2007. The comparison is done at regional and global scales. Based on CALIPSO arid OMI observations, the vertical and horizontal extent of clouds and aerosols are determined and the effects of aerosol type selection, load, cloud fraction on aerosol identification are discussed. It was found that the spatial-temporal correlation found between CALIPSO and OMI observations, is strongly dependent on aerosol types and cloud contamination. CALIPSO is more sensitivity to cloud and often misidentifies desert dust aerosols as cloud, while some small scale aerosol layers as well as some pollution aerosols are unidentified by OMI UV AI. Large differences in aerosol distribution patterns between CALIPSO and OMI are observed, especially for the smoke and pollution aerosol dominated areas. In addition, the results found a significant correlation between CALIPSO lidar 1064 nm backscatter and the OMI UV AI over the study regions.

Chen, Zhong↗

Calibration Plans for the Global Precipitation Measurement (GPM)

The Global Precipitation Measurement (GPM) is an international effort led by the National Aeronautics and Space Administration (NASA) of the U.S.A. and the National Space Development Agency of Japan (NASDA) for the purpose of improving research into the global water and energy cycle. GPM will improve climate, weather, and hydrological forecasts through more frequent and more accurate measurement of precipitation world-wide. Comprised of U.S. domestic and international partners, GPM will incorporate and assimilate data streams from many spacecraft with varied orbital characteristics and instrument capabilities. Two of the satellites will be provided directly by GPM, the core satellite and a constellation member. The core satellite, at the heart of GPM, is scheduled for launch in November 2007. The core will carry a conical scanning microwave radiometer, the GPM Microwave Imager (GMI), and a two-frequency cross-track-scanning radar, the Dual-frequency Precipitation Radar (DPR). The passive microwave channels and the two radar frequencies of the core are carefully chosen for investigating the varying character of precipitation over ocean and land, and from the tropics to the high-latitudes. The DPR will enable microphysical characterization and three-dimensional profiling of precipitation. The GPM-provided constellation spacecraft will carry a GMI radiometer identical to that on the core spacecraft. This paper presents calibration plans for the GPM, including on-board instrument calibration, external calibration methods, and the role of ground validation. Particular emphasis is on plans for inter-satellite calibration of the GPM constellation. With its Unique instrument capabilities, the core spacecraft will serve as a calibration transfer standard to the GPM constellation. In particular the Dual-frequency Precipitation Radar aboard the core will check the accuracy of retrievals from the GMI radiometer and will enable improvement of the radiometer retrievals. Observational intersections of the core with the constellation spacecraft are essential in applying this technique to the member satellites. Information from core spacecraft retrievals during intersection events will be transferred to the constellation radiometer instruments in the form of improved calibration and, with experience, improved radiometric algorithms. In preparation for the transfer standard technique, comparisons using the Tropical Rainfall Measuring Mission (TRMM) with sun-synchronous radiometers have been conducted. Ongoing research involves study of critical variables in the inter-comparison, such as correlation with spatial-temporal separation of intersection events, frequency of intersection events, variable azimuth look angles, and variable resolution cells for the various sensors.

Bidwell, S. W.↗

Anomalously slow hot carrier cooling and insulator-to-metal transition in a photo-doped Mott insulator

Photo-doped Mott insulators can exhibit novel photocarrier transport and relaxation dynamics and non-equilibrium phases. However, time-resolved real-space imaging of these processes are still lacking. Here, we use scanning ultrafast electron microscopy (SUEM) to directly visualize the spatial-temporal evolution of photoexcited species in a spin-orbit assisted Mott insulator α-RuCl 3 . At low optical fluences, we observe extremely long hot photocarrier transport time over one nanosecond, almost an order of magnitude longer than any known values in conventional semiconductors. At higher optical fluences, we observe nonlinear features suggesting a photo-induced insulator-to-metal transition, which is unusual in a large-gap Mott insulator. Our results demonstrate the rich physics in a photo-doped Mott insulator that can be extracted from spatial-temporal imaging and showcase the capability of SUEM to sensitively probe photoexcitations in strongly correlated electron systems.

42 ENGINEERING↗

First Application of the Zeeman Technique to Remotely Measure Auroral Electrojet Intensity From Space

Using the O2 118 GHz spectral radiance measurements obtained by the Microwave Limb Sounder instrument on board the Aura spacecraft, we demonstrate that the Zeeman effect can be used to remotely measure the magnetic field perturbations produced by the auroral electrojet near the Hall current closure altitudes. Our derived current-induced magnetic field perturbations are found to be highly correlated with those coincidently obtained by ground magnetometers. These perturbations are also found to be linearly correlated with auroral electrojet strength. The statistically derived polar maps of our measured magnetic field perturbation reveal a spatial-temporal morphology consistent with that produced by the Hall current during substorms and storms. With today's technology, a constellation of compact, low-power, high spectral-resolution cubesats would have the capability to provide high precision and spatiotemporal magnetic field samplings needed for auroral electrojet measurements to gain insights into the spatiotemporal behavior of the auroral electrojet system.

spectral radiance measurements; Aur↗

A simple electrical-mechanical model of the heart applied to the study of electrical-mechanical alternans

Recent evidence has shown that a subtle alternation in the surface ECG (electrical alternans) may be correlated with the susceptibility to ventricular fibrillation. In the present work, the author presents evidence that a mechanical alternation in the heartbeat (mechanical alternans) generally accompanies electrical alternans. A simple finite-element computer model which emulates both the electrical and the mechanical activity of the heart is presented. A pilot animal study is also reported. The computer model and the animal study both found that (1) there exists a regime of combined electrical-mechanical alternans during the transition from a normal rhythm towards a fibrillatory rhythm, (2) the detected degree of alternation is correlated with the relative instability of the rhythm, and (3) the electrical and mechanical alternans may result from a dispersion in local electrical properties leading to a spatial-temporal alternation in the electrical conduction process.

Clancy, Edward A.↗

MemFriend: Understanding Memory Performance with Spatial-Temporal Affinity

In HPC applications, memory access behavior is one of the main factors affecting performance. Improving an application’s memory access behavior involves optimizing data layout and/or restructuring code, and requires studying spatial-temporal data locality. Existing data locality analyses focus on single-location metrics and are restricted to evaluating temporal locality. We introduce spatial-temporal affinity metrics that quantify temporal access proximity, forward access correlation, and nearby access correlation between pairs of memory locations. We describe methods for distinguishing between potential vs. realized affinity and for reasoning about affinity at multiple resolutions (3D, 2D, 1D). Finally, we construct spatial-temporal affinity signatures that classify memory behavior and that be used to reason about changes in software (data relayout, code refactoring) or hardware (caching, prefetching). We describe methods for signature visualization, interpretation, and quantitative comparison of signatures. We evaluate our methodology using applications with variants that contrast data structures, data layouts and algorithms. We show that spatial-temporal affinity analysis provides novel insights and enables predictive reasoning about application performance when contrasted with reuse distance analysis.

Suriyakumar, Yasodhadevi↗

The Spectral Information Based Angular Correction Methodology for Satellite Intercalibration Applications

Satellite inter-calibration often requires collocated observations with minimized discrepancies in sun-view angles, observation times, and sensor characteristics. The collocation criteria directly impact achievable inter-calibration accuracy. Addressing potential angular mismatches in inter-calibration samples is critical but not as fully recognized and addressed as spatial-temporal mismatches in many studies. To achieve high-accuracy corrections for errors due to mismatched sun-view geometry angles, an angular correction algorithm has been developed for the Climate Absolute Radiance and Refractivity Observatory Pathfinder (CPF) mission. This algorithm uses spectral correlation relationships to estimate differences in spectral radiances measured at different angles. This methodology can be extended for inter-calibrations between sensors measuring band radiances across a broad spectral region. We demonstrate its application in reducing angular mismatch errors between collocated measurements of multi-spectral imaging sensors, using the inter-calibration between the Moderate Resolution Imaging Spectrometer (MODIS) and the Visible Infrared Imaging Radiometer Suite (VIIRS) as an example. The angular correction allows for more relaxed collocation criteria so that more satellite-based inter-calibration samples can be utilized. Furthermore, implementing the angular correction algorithm improves inter-calibration accuracy in applications where angular mismatch errors have not been explicitly addressed previously.

Wan Wu↗

Human Plague Risk: Spatial-Temporal Models

This chpater reviews the use of spatial-temporal models in identifying potential risks of plague outbreaks into the human population. Using earth observations by satellites remote sensing there has been a systematic analysis and mapping of the close coupling between the vectors of the disease and climate variability. The overall result is that incidence of plague is correlated to positive El Nino/Southem Oscillation (ENSO).

Pinzon, Jorge E.↗

LIS Validation at The KSC-ER

Each year, thousands of lightning electric field disturbances are recorded and archived by the ground-based field mill (FM) network at the NASA Kennedy Space Center (KSC) and USAF Eastern Range (ER). The FM network has a range of several tens of kilometers, and a digital accuracy of 4 V/m. It has provided years of continuous lightning warning surveillance to KSC-ER space vehicle launch operations, and has undergone one major hardware upgrade since its inception in the early 1970s. Additional KSC lightning warning data is derived from a multistation radio time-of-arrival system called Lightning Detection and Ranging (LDAR). This system provides the location and space-time mapping of individual lightning channels (for both cloud and ground flashes). Additional lightning information for the KSC region is available from the National Lightning Detection Network (NLDN) and a 5-station local magnetic direction finder network. In this study, all of the above mentioned data are used to ground-validate data derived from the Lightning Imaging Sensor (LIS) onboard the Tropical Rainfall Measuring Mission (TRMM). The FM network can be used to retrieve the charges deposited in a lightning flash, provided the flash is within a few kilometers of the FM Network. Although it is rare to obtain a TRMM overpass of thunderstorms hat occur this close to the FM network, seven such storms have been found and examined in this study. We compare the times and locations of LIS optical pulses with the spatial-temporal character of the FM, LDAR, and magnetic direction finder data. We also inter-compare LIS optical pulse amplitude data with FM-derived charge magnitudes, number of LDAR radio sources, and peak current values from magnetic direction finder data. Generally speaking, LIS lightning locations and times agree favorably with the KSC ground-based systems for most cases, but little correlation appears to exist between optical pulse amplitude and any of charge, # LDAR sources, peak current), owing possibly to the effects of source complexity and/or cloud multiple scattering.

Koshak, William J.↗

Continental Spatio-Temporal Data Analysis with Linear Spectral Mixture Model Using FOSS

This work demonstrates the development and implementation of a Fully Constrained Least Squares (FCLS) unmixing model developed in C++ programming language with OpenCV package and boost C++ libraries in the NASA Earth Exchange (NEX). Visualization of the results is supported by GRASS GIS and statistical analysis is carried in R in a Linux system environment. FCLS was first tested on computer simulated data with Gaussian noise of various signal-to-noise ratio, and Landsat data of an agricultural scenario and an urban environment using a set of global end members of substrate (soils, sediments, rocks, and non-photosynthetic vegetation), vegetation that includes green photosynthetic plants and dark objects which encompasses absorptive substrate materials, clear water, deep shadows, etc. For the agricultural scenario, a spectrally diverse collection of 11 scenes of Level 1 terrain corrected, cloud free Landsat-5 TM data of Fresno, California, USA were unmixed and the results were validated with the corresponding ground data. To study an urbanized landscape, a clear sky Landsat-5 TM data were unmixed and validated with coincident World View-2 abundance maps (of 2 m spatial resolution) for an area of San Francisco, California, USA. The results were evaluated using descriptive statistics, correlation coefficient, RMSE, probability of success, boxplot and bivariate distribution function. Finally, FCLS was used for sub-pixel land cover analysis of the monthly WELD (Wen-enabled Landsat data) repository from 2008 to 2011 of North America. The abundance maps in conjunction with DMSP-OLS nighttime lights data were used to extract the urban land cover features and analyze their spatial-temporal growth.

Landsat Satellites↗

Interpretable Models for Workflow Differentiation in High-Performance Scientific Networks

Scientific workflows in high-performance networks spawn hundreds of interdependent flows that must be managed collectively—yet existing network classifiers treat each flow in isolation, leading to fragmented QoS decisions and missed interflow patterns. We present a novel traffic classification solution that operates at the workflow level, distinguishing entire filetransfer operations from streaming analytics by capturing how concurrent flows interact and burst together. We introduce a workflow identification window (WIW) that ingests raw packet headers from parallel flows into unified tensors, preserving the spatial-temporal patterns that differentiate scientific workflows. This approach achieves 98.7% accuracy using CNN, LSTM, and hybrid architectures, while maintaining 84% accuracy on production traffic collected a week later—demonstrating robustness to temporal drift. By integrating SHAP and GradCAM explainability, we reveal that early-packet timing patterns and cross-flow correlations drive classification decisions, providing operators with interpretable insights. Our system enables coherent workflow-level QoS enforcement and dynamic bandwidth allocation in scientific networks, eliminating manual per-flow configuration while maintaining classification latency at millisecond level.

Giannakou, Anna [LBL, Berkeley]↗

Magnetohydrodynamic Convection in the Outer Core and its Geodynamic Consequences

The Earth's fluid outer core is in vigorous convection through much of the Earth's history. In addition to generating and maintaining Earth s time-varying magnetic field (geodynamo), the core convection also generates mass redistribution in the core and a dynamical pressure field on the core-mantle boundary (CMB). All these shall result in various core-mantle interactions, and contribute to surface geodynamic observables. For example, electromagnetic core-mantle coupling arises from finite electrically conducting lower mantle; gravitational interaction occurs between the cores and the heterogeneous mantle; mechanical coupling may also occur when the CMB topography is aspherical. Besides changing the mantle rotation via the coupling torques, the mass-redistribution in the core shall produce a spatial-temporal gravity anomaly. Numerical modeling of the core dynamical processes contributes in several geophysical disciplines. It helps explain the physical causes of surface geodynamic observables via space geodetic techniques and other means, e.g. Earth's rotation variation on decadal time scales, and secular time-variable gravity. Conversely, identification of the sources of the observables can provide additional insights on the dynamics of the fluid core, leading to better constraints on the physics in the numerical modeling. In the past few years, our core dynamics modeling efforts, with respect to our MoSST model, have made significant progress in understanding individual geophysical consequences. However, integrated studies are desirable, not only because of more mature numerical core dynamics models, but also because of inter-correlation among the geophysical phenomena, e.g. mass redistribution in the outer core produces not only time-variable gravity, but also gravitational core-mantle coupling and thus the Earth's rotation variation. They are expected to further facilitate multidisciplinary studies of core dynamics and interactions of the core with other components of the Earth.

Kuang, Weijia↗

The relation of motion sickness to the spatial-temporal properties of velocity storage

Tilting the head in roll to or from the upright while rotating at a constant velocity (roll while rotating, RWR) alters the position of the semicircular canals relative to the axis of rotation. This produces vertical and horizontal nystagmus, disorientation, vertigo, and nausea. With recurrent exposure, subjects habituate and can make more head movements before experiencing overpowering motion sickness. We questioned whether promethazine lessened the vertigo or delayed the habituation, whether habituation of the vertigo was related to the central vestibular time constant, i.e., to the time constant of velocity storage, and whether the severity of the motion sickness was related to deviation of the axis of eye velocity from gravity. Sixteen subjects received promethazine and placebo in a double-blind, crossover study in two consecutive 4-day test series 1 month apart, termed series I and II. Horizontal and vertical eye movements were recorded with video-oculography while subjects performed roll head movements of approx. 45 degrees over 2 s to and from the upright position while being rotated at 138 degrees /s around a vertical axis. Motion sickness was scaled from 1 (no sickness) to an endpoint of 20, at which time the subject was too sick to continue or was about to vomit. Habituation was determined by the number of head movements that subjects made before reaching the maximum motion sickness score of 20. Head movements increased steadily in each session with repeated testing, and there was no difference between the number of head movements made by the promethazine and placebo groups. Horizontal and vertical angular vestibulo-ocular reflex (aVOR) time constants declined in each test, with the declines being closely correlated to the increase in the number of head movements. The strength of vertiginous sensation was associated with the amount of deviation of the axis of eye velocity from gravity; the larger the deviation of the eye velocity axis from gravity, the more severe the motion sickness. Thus, promethazine neither reduced the nausea associated with RWR, nor retarded or hastened habituation. The inverse relationship between the aVOR time constants and number of head movements to motion sickness, and the association of the severity of motion sickness with the extent, strength, and time of deviation of eye velocity from gravity supports the postulate that the spatiotemporal properties of velocity storage, which are processed between the nodulus and uvula of the vestibulocerebellum and the vestibular nuclei, are likely to represent the source of the conflict responsible for producing motion sickness.

NASA Discipline Neuroscience↗

Vertical Moist Thermodynamic Structure and Spatial-Temporal Evolution of the MJO in AIRS Observations

The atmospheric moisture and temperature profiles from the Atmospheric Infrared Sounder (AIRS)/Advanced Microwave Sounding Unit on the NASA Aqua mission, in combination with the precipitation from the Tropical Rainfall Measuring Mission (TRMM), are employed to study the vertical moist thermodynamic structure and spatial-temporal evolution of the Madden-Julian oscillation (MJO). The AIRS data indicate that, in the Indian Ocean and western Pacific, the temperature anomaly exhibits a trimodal vertical structure: a warm (cold) anomaly in the free troposphere (800-250 hPa) and a cold (warm) anomaly near the tropopause (above 250 hPa) and in the lower troposphere (below 800 hPa) associated with enhanced (suppressed) convection. The AIRS moisture anomaly also shows markedly different vertical structures as a function of longitude and the strength of convection anomaly. Most significantly, the AIRS data demonstrate that, over the Indian Ocean and western Pacific, the enhanced (suppressed) convection is generally preceded in both time and space by a low-level warm and moist (cold and dry) anomaly and followed by a low-level cold and dry (warm and moist) anomaly. The MJO vertical moist thermodynamic structure from the AIRS data is in general agreement, particularly in the free troposphere, with previous studies based on global reanalysis and limited radiosonde data. However, major differences in the lower-troposphere moisture and temperature structure between the AIRS observations and the NCEP reanalysis are found over the Indian and Pacific Oceans, where there are very few conventional data to constrain the reanalysis. Specifically, the anomalous lower-troposphere temperature structure is much less well defined in NCEP than in AIRS for the western Pacific, and even has the opposite sign anomalies compared to AIRS relative to the wet/dry phase of the MJO in the Indian Ocean. Moreover, there are well-defined eastward-tilting variations of moisture with height in AIRS over the central and eastern Pacific that are less well defined, and in some cases absent, in NCEP. In addition, the correlation between MJO-related mid-tropospheric water vapor anomalies and TRMM precipitation anomalies is considerably more robust in AIRS than in NCEP, especially over the Indian Ocean. Overall, the AIRS results are quite consistent with those predicted by the frictional Kelvin-Rossby wave/conditional instability of the second kind (CISK) theory for the MJO.

troposhere↗

Detecting outbreaks using a spatial latent field

In this paper, we present a method for estimating the infection-rate of a disease as a spatial-temporal field. Our data comprises time-series case-counts of symptomatic patients in various areal units of a region. We extend an epidemiological model, originally designed for a single areal unit, to accommodate multiple units. The field estimation is framed within a Bayesian context, utilizing a parameterized Gaussian random field as a spatial prior. We apply an adaptive Markov chain Monte Carlo method to sample the posterior distribution of the model parameters condition on COVID-19 case-count data from three adjacent counties in New Mexico, USA. Our results suggest that the correlation between epidemiological dynamics in neighboring regions helps regularize estimations in areas with high variance (i.e., poor quality) data. Using the calibrated epidemic model, we forecast the infection-rate over each areal unit and develop a simple anomaly detector to signal new epidemic waves. Our findings show that anomaly detector based on estimated infection-rates outperforms a conventional algorithm that relies solely on case-counts.

Safta, Cosmin [Sandia National Laboratories (SNL-C↗

Phase Synchronization and Desynchronization of Structural Response Induced by Turbulent and External Sound

Acoustic and turbulent boundary layer flow loadings over a flexible structure are used to study the spatial-temporal dynamics of the response of the structure. The stability of the spatial synchronization and desynchronization by an active external force is investigated with an array of coupled transducers on the structure. In the synchronous state, the structural phase is locked, which leads to the formation of spatial patterns while the amplitude peaks exhibit chaotic behaviors. Large amplitude, spatially symmetric loading is superimposed on broadband, but in the desynchronized state, the spectrum broadens and the phase space is lost. The resulting pattern bears a striking resemblance to phase turbulence. The transition is achieved by using a low power external actuator to trigger broadband behaviors from the knowledge of the external acoustic load inducing synchronization. The changes are made favorably and efficiently to alter the frequency distribution of power, not the total power level. Before synchronization effects are seen, the panel response to the turbulent boundary layer loading is discontinuously spatio-temporally correlated. The stability develops from different competing wavelengths; the spatial scale is significantly shorter than when forced with the superimposed external sound. When the external sound level decreases and the synchronized phases are lost, changes in the character of the spectra can be linked to the occurrence of spatial phase transition. These changes can develop broadband response. Synchronized responses of fuselage structure panels have been observed in subsonic and supersonic aircraft; results from two flights tests are discussed.

Maestrello, Lucio↗

Arctic Tundra Vegetation Functional Types Based on Photosynthetic Physiology and Optical Properties

Non-vascular plants (lichens and mosses) are significant components of tundra landscapes and may respond to climate change differently from vascular plants affecting ecosystem carbon balance. Remote sensing provides critical tools for monitoring plant cover types, as optical signals provide a way to scale from plot measurements to regional estimates of biophysical properties, for which spatial-temporal patterns may be analyzed. Gas exchange measurements were collected for pure patches of key vegetation functional types (lichens, mosses, and vascular plants) in sedge tundra at Barrow, AK. These functional types were found to have three significantly different values of light use efficiency (LUE) with values of 0.013 plus or minus 0.0002, 0.0018 plus or minus 0.0002, and 0.0012 plus or minus 0.0001 mol C mol (exp -1) absorbed quanta for vascular plants, mosses and lichens, respectively. Discriminant analysis of the spectra reflectance of these patches identified five spectral bands that separated each of these vegetation functional types as well as nongreen material (bare soil, standing water, and dead leaves). These results were tested along a 100 m transect where midsummer spectral reflectance and vegetation coverage were measured at one meter intervals. Along the transect, area-averaged canopy LUE estimated from coverage fractions of the three functional types varied widely, even over short distances. The patch-level statistical discriminant functions applied to in situ hyperspectral reflectance data collected along the transect successfully unmixed cover fractions of the vegetation functional types. The unmixing functions, developed from the transect data, were applied to 30 m spatial resolution Earth Observing-1 Hyperion imaging spectrometer data to examine variability in distribution of the vegetation functional types for an area near Barrow, AK. Spatial variability of LUE was derived from the observed functional type distributions. Across this landscape, a fivefold variation in tundra LUE was observed. LUE calculated from the functional type cover fractions was also correlated to a spectral vegetation index developed to detect vegetation chlorophyll content. The concurrence of these alternate methods suggest that hyperspectral remote sensing can distinguish functionally distinct vegetation types and can be used to develop regional estimates of photosynthetic LUE in tundra landscapes.

VEGETATION↗