Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vegetation model”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The critical role of soil moisture in compound hazards

Soil moisture regulates the exchange of energy, water, and carbon across land–vegetation–atmosphere interfaces. Extremes in soil moisture can amplify natural hazards through interactions with diverse Earth system processes. Despite its mechanistic importance, soil moisture remains underrepresented in hazard research and predictive frameworks. Here, in this study, we review our current understanding of the role of soil moisture in the evolution and onset of diverse compound hazards by synthesizing the latest findings from observational and modelling studies. We highlight key soil moisture mechanisms, including atmospheric feedbacks that amplify drought–heatwave–wildfire events, precipitation couplings that promote clustered storms, and threshold responses that drive vegetation die-offs, trigger landslides, and induce flooding. Persistent challenges in observational data, model representation and operational implementation have limited the integration of soil moisture into hazard early-warning systems. Addressing these gaps through advances in observations, data assimilation, and physics-based and data-driven modelling will enhance hazard prediction and preparedness in a rapidly changing world.

Li, Chuxuan [University of California, Los Angeles↗
DOE OSTI2020

Turbulence-parameter estimation for current-energy converters using surrogate model optimization

Surrogate models maximize information utility by building predictive models in place of computational or experimentally expensive model runs. Marine hydrokinetic current energy converters require large-domain simulations to estimate array efficiencies and environmental impacts. Meso-scale models typically represent turbines as actuator discs that act as momentum sinks and sources of turbulence and its dissipation. An OpenFOAM model was developed where actuator disc k-ε turbulence was characterized using an approach developed for flows through vegetative canopies. Turbine-wake data from laboratory flume experiments collected at two influent turbulence intensities were used to calibrate parameters in the turbulence-source terms in the k-ε equations. Additionally, parameter influences on longitudinal wake profiles were estimated using Gaussian process regression with subsequent optimization minimizing the objective function within 3.1% of those obtained using the full model representation, but for 74% of the computational cost (far fewer model runs). This framework facilitates more efficient parameterization of the turbulence-source equations using turbine-wake data.

16 TIDAL AND WAVE POWER↗
NASA NTRS2002

Urban Surface Radiative Energy Budgets Determined Using Aircraft Scanner Data

It is estimated that by the year 2025, 80% of the world's population will live in cities. The extent of these urban areas across the world can be seen in an image of city lights from the Defense Meteorological Satellite Program. In many areas of North America and Europe, it is difficult to separate individual cities because of the dramatic growth and sprawl of urbanized areas. This conversion of the natural landscape vegetation into man-made urban structures such as roads and buildings drastically alter the regional surface energy budgets, hydrology, precipitation patterns, and meteorology. One of the earliest recognized and measured phenomena of urbanization is the urban heat island (UHI) which was reported as early as 1833 for London and 1862 for Paris. The urban heat island results from the energy that is absorbed by man-made materials during the day and is released at night resulting in the heating of the air within the urban area. The magnitude of the air temperature difference between the urban and surrounding countryside is highly dependent on the structure of the urban area, amount of solar immolation received during the day, and atmospheric conditions during the night. These night time air temperature differences can be in the range of 2 to 5 C. or greater. Although day time air temperature differences between urban areas and the countryside exists during the day, atmospheric mixing and stability reduce the magnitude. This phenomena is not limited to large urban areas, but also occurs in smaller metropolitan areas. The UHI has significant impacts on the urban air quality, meteorology, energy use, and human health. The UPI can be mitigated through increasing the amount of vegetation and modification of urban surfaces using high albedo materials for roofs and paved surfaces. To understand why the urban heat island phenomenon exists it is useful to define the surface in terms of the surface energy budget. Surface temperature and albedo is a major component of the surface energy budget. Knowledge of it is important in any attempt to describe the radiative and mass fluxes which occur at the surface. Use of energy terms in modeling surface energy budgets allows the direct comparison of various land surfaces encountered in a urban landscape, from vegetated (forest and herbaceous) to non-vegetated (bare soil, roads, and buildings). These terms are also easily measured using remote sensing from aircraft or satellite platforms allowing one to examine the spacial variability. The partitioning of energy budget terms depends on the surface type. In natural landscapes, the partitioning is dependent on canopy biomass, leaf area index, aerodynamic roughness, and moisture status, all of which are influenced by the development stage of the ecosystem. In urban landscapes, coverage by man-made materials substantially alters the surface face energy budget. The remotely sensed data obtained from aircraft and satellites, when properly calibrated allows the measurement of important terms in the radiative surface energy budget a urban landscape scale.

Luvall, Jeffrey C.↗
NASA NTRS1989

Effects of implementing the Simple Biosphere Model in a general circulation model

The Simple Biosphere Model (SiB) of Sellers et al., was designed to simulate the interactions between the earth's land surface and the atmosphere by treating the vegetation explicitly and realistically, thereby incorporating the biophysical controls on the exchanges of radiation, momentum, sensible and latent heat between the two systems. This paper describes the steps taken to implement SiB in a modified version of the National Meteorological Center's global spectral general circulation model (GCM) and explores the impact of the implementation on the simulated land surface fluxes and near-surface meteorological conditions. The coupled model (SiB-GCM) was used to produce summer and winter simulations. The same GCM was used with a conventional hydrological model (Ctl-GCM) to produce comparable 'control' summer and winter simulations for comparison. It was found that SiB-GCM produced a more realistic partitioning of energy at the land surface than Ctl-GCM. Generally, SiB-GCM, produced more sensible heat flux and less latent heat flux over vegetated land than did Ctl-GCM and this resulted in a much deeper daytime planetary boundary layer and reduced precipitation rates over the continents in SiB-GCM. In the summer simulation, the 200 mb jet stream was slightly weakened in the SiB-GCM relative to the Ctl-GCM results and analyses made from observations.

Sato, N.↗
NASA NTRS2007

Simulation Studies of Forest Structure using 3D Lidar and Radar Models

The use of lidars and radars to measure forest structure attributes such as height and biomass are being considered for future Earth Observation missions. Large footprint lidar makes a direct measurement of the heights of scatterers in the illuminated footprint and can yield information about the vertical profile of the canopy. Synthetic Aperture Radar (SAR) is known to sense the canopy volume, especially at longer wavelengths and is useful for estimating biomass. Interferometric SAR (InSAR) has been shown to yield forest canopy height information. For example, the height of scattering phase retrieved from InSAR data is considered to be correlated with the three height and the spatial structure of the forest stand. There is much interest in exploiting these technologies separately and together to get important information for carbon cycle and ecosystem science. More detailed information of the electromagnetic radiation interactions within forest canopies is needed. And backscattering models can be of much utility here. As part of a NASA funded project to explore data fusion, a three-dimensional (3D) coherent radar backscattering model and a 3D lidar backscatter models were used to investigate the use of large footprint lidar, SAR and InSAR for characterizing realistic forest scenes. For this paper, we use stem maps and other forest measurements to develop a realistic spatial structure of a spruce-hemlock forest canopy found in Maine, USA. The radar and lidar models used measurements of the 3D forest scene as input and simulated the coherent radar backscattering signature and 1064nm energy backscatter, respectively. The relationships of backscatter derived forest structure were compared with field measurements. In addition, we also had detailed airborne lidar (Laser Imaging Vegetation Sensor, LVIS) data available over the stem map sites that was used to study the accuracies of tree height derived from modeled SAR backscatter and the scattering phase center retrieved from the simulated InSAR data will be compared with the height indices, or other structure parameters derived from the lidar data. These results will address the possible synergies between lidar and radar in data in terms of forest structural information.

Sun, Guoqing↗
NASA NTRS2004

Estimation of Soil Moisture with L-band Multi-polarization Radar

Through analyses of the model simulated data-base, we developed a technique to estimate surface soil moisture under HYDROS radar sensor (L-band multi-polarizations and 40deg incidence) configuration. This technique includes two steps. First, it decomposes the total backscattering signals into two components - the surface scattering components (the bare surface backscattering signals attenuated by the overlaying vegetation layer) and the sum of the direct volume scattering components and surface-volume interaction components at different polarizations. From the model simulated data-base, our decomposition technique works quit well in estimation of the surface scattering components with RMSEs of 0.12,0.25, and 0.55 dB for VV, HH, and VH polarizations, respectively. Then, we use the decomposed surface backscattering signals to estimate the soil moisture and the combined surface roughness and vegetation attenuation correction factors with all three polarizations.

Shi, J.↗
NASA NTRS1989

Implementing the Simple Biosphere Model (SiB) in a general circulation model: Methodologies and results

The Simple Biosphere MOdel (SiB) of Sellers et al., (1986) was designed to simulate the interactions between the Earth's land surface and the atmosphere by treating the vegetation explicitly and relistically, thereby incorporating biophysical controls on the exchanges of radiation, momentum, sensible and latent heat between the two systems. The steps taken to implement SiB in a modified version of the National Meteorological Center's spectral GCM are described. The coupled model (SiB-GCM) was used with a conventional hydrological model (Ctl-GCM) to produce summer and winter simulations. The same GCM was used with a conventional hydrological model (Ctl-GCM) to produce comparable 'control' summer and winter variations. It was found that SiB-GCM produced a more realistic partitioning of energy at the land surface than Ctl-GCM. Generally, SiB-GCM produced more sensible heat flux and less latent heat flux over vegetated land than did Ctl-GCM and this resulted in the development of a much deeper daytime planetary boundary and reduced precipitation rates over the continents in SiB-GCM. In the summer simulation, the 200 mb jet stream and the wind speed at 850 mb were slightly weakened in the SiB-GCM relative to the Ctl-GCM results and equivalent analyses from observations.

Sato, N.↗
NASA NTRS1996

Passive Microwave Remote Sensing of Soil Moisture

Microwave remote sensing provides a unique capability for direct observation of soil moisture. Remote measurements from space afford the possibility of obtaining frequent, global sampling of soil moisture over a large fraction of the Earth's land surface. Microwave measurements have the benefit of being largely unaffected by cloud cover and variable surface solar illumination, but accurate soil moisture estimates are limited to regions that have either bare soil or low to moderate amounts of vegetation cover. A particular advantage of passive microwave sensors is that in the absence of significant vegetation cover soil moisture is the dominant effect on the received signal. The spatial resolutions of passive Microwave soil moisture sensors currently considered for space operation are in the range 10-20 km. The most useful frequency range for soil moisture sensing is 1-5 GHz. System design considerations include optimum choice of frequencies, polarizations, and scanning configurations, based on trade-offs between requirements for high vegetation penetration capability, freedom from electromagnetic interference, manageable antenna size and complexity, and the requirement that a sufficient number of information channels be available to correct for perturbing geophysical effects. This paper outlines the basic principles of the passive microwave technique for soil moisture sensing, and reviews briefly the status of current retrieval methods. Particularly promising are methods for optimally assimilating passive microwave data into hydrologic models. Further studies are needed to investigate the effects on microwave observations of within-footprint spatial heterogeneity of vegetation cover and subsurface soil characteristics, and to assess the limitations imposed by heterogeneity on the retrievability of large-scale soil moisture information from remote observations.

Njoku, Eni G.↗
NASA NTRS1980

Microwave remote sensing of soil moisture content over bare and vegetated fields

The author has identified the following significant results. Ground truth of soil moisture content, and ambient air and soil temperatures were acquired concurrently with measurements of soil moisture in bare fields and fields covered with grass, corn, and soybeans obtained with 1.4 GHz and 5 GHz radiometers mounted on a truck. The biomass of the vegetation was sampled about once a week. The measured brightness temperatures over the bare fields were compared with those of radiative transfer model calculations using as inputs the acquired soil moisture and temperatures data with appropriate values of dielectric constants for soil-water mixtures. A good agreement was found between the calculated and measured results over 10 deg to 70 deg incident angles. The presence of vegetation reduced the sensitivity of soil moisture sensing. At 1.4 GHz the sensitivity reduction ranged from about 20% for 10 cm tall grassland cover to over 50 to 60% for the dense soybean field. At 5 GHz corresponding reduction in sensitivity ranged from approximately 70% to approximately 90%.

Wang, J. R.↗
NASA NTRS2022

Modeling the Errors of a Time Series Algorithm for Retrieving Soil Moisture in the NISAR Mission

The National Aeronautics and Space Administration (NASA) - Indian Space Research Organization (ISRO) Synthetic Aperture Radar (NISAR) mission plan to launch a SAR operating at L- and S-band with a 12-day repeat frequency. A global soil moisture product at 200 m spatial resolution derived from 200 m NISAR radar measurements is currently under development. Although several retrieval algorithms are being investigated, this paper focuses on a “time series ratio” retrieval approach. In order to understand and assess the performance of this algorithm, an error model has been developed and is reported in this paper. The model is applied to examine errors as a function of the instrument characteristics and for a given location. Initial progress in including vegetation effects and in predicting errors as a function of spatial location is also described.

NISAR↗
NASA NTRSDate not supplied

Improved Peatland Hydrology and L-Band Microwave Radiative Transfer Modeling in Version 7 of the SMAP Level-4 Soil Moisture Data Assimilation Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface (0-5 cm) and root-zone (0-100 cm) soil moisture from April 2015 to present with a mean latency of 2.5 days from the time of observation. The L4_SM estimates are derived from the assimilation of SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment land surface model. The recently released Version 7 of L4_SM comprises two key advances. First, the Catchment model now includes the recently developed PEATCLSM hydrology module for peatlands and uses an updated global map of peatlands. Second, revised parameters are used in the L-band radiative transfer model that converts the simulated soil moisture and temperature estimates into Tb predictions for use in the radiance-based L4_SM analysis. Specifically, the L-band scattering albedo, soil roughness, and vegetation opacity climatology are obtained from the SMAP Level-2 dual-channel soil moisture retrieval product. In this presentation, we examine the improvements gained with the model and parameter changes. In peatlands, the dynamics of water table depth, surface soil moisture and evapotranspiration are considerably improved in L4_SM Version 7 over Version 6 when evaluated against in situ measurements. These improvements are also manifested in smaller Tb observation-minus-forecast residuals.

SMAP↗
NASA NTRS1979

Theoretical modelling and experimental data matching for active and passive microwave remote sensing of Earth terrain

Two theoretical models were developed to characterize terrain media: a random medium with a variance, a horizontal correlation length; and a homogeneous dielectric containing discrete scatterers. The earth terrain is modelled as layers of such scattering media bounded by air above and half-space below. Matching the theoretical results with experimental data collected from vegetation and snow-ice fields shows that: (1) for observation angles near nadir, rough surface effects are important; (2) for snow-ice field the horizontal correlation length is greater than the vertical correlation length whereas for vegetation field their relative sizes depend on the types of vegetation; (3) the vertically polarized backscattering cross-section is always larger than the horizontally polarized backscattering cross-section for half-space scattering media; (4) for snow field displaying diurnal change, a three-layer model including a thin top layer caused by sunlight illumination must be used; and (5) for a random medium with equal horizontal and vertical correlation lengths, the measured data can also be matched with a corresponding discrete scatterer model.

Kong, J. A.↗
NASA NTRS2014

Toward Improved Land Surface Initialization in Support of Regional WRF Forecasts at the Kenya Meteorological Department

Flooding and drought are two key forecasting challenges for the Kenya Meteorological Department (KMD). Atmospheric processes leading to excessive precipitation and/or prolonged drought can be quite sensitive to the state of the land surface, which interacts with the boundary layer of the atmosphere providing a source of heat and moisture. The development and evolution of precipitation systems are affected by heat and moisture fluxes from the land surface within weakly-sheared environments, such as in the tropics and sub-tropics. These heat and moisture fluxes during the day can be strongly influenced by land cover, vegetation, and soil moisture content. Therefore, it is important to represent the land surface state as accurately as possible in numerical weather prediction models. Enhanced regional modeling capabilities have the potential to improve forecast guidance in support of daily operations and high-end events over east Africa. KMD currently runs a configuration of the Weather Research and Forecasting (WRF) model in real time to support its daily forecasting operations, invoking the Nonhydrostatic Mesoscale Model (NMM) dynamical core. They make use of the National Oceanic and Atmospheric Administration / National Weather Service Science and Training Resource Center's Environmental Modeling System (EMS) to manage and produce the WRF-NMM model runs on a 7-km regional grid over eastern Africa. Two organizations at the National Aeronautics and Space Administration Marshall Space Flight Center in Huntsville, AL, SERVIR and the Short-term Prediction Research and Transition (SPoRT) Center, have established a working partnership with KMD for enhancing its regional modeling capabilities. To accomplish this goal, SPoRT and SERVIR will provide experimental land surface initialization datasets and model verification capabilities to KMD. To produce a land-surface initialization more consistent with the resolution of the KMD-WRF runs, the NASA Land Information System (LIS) will be run at a comparable resolution to provide real-time, daily soil initialization data in place of interpolated Global Forecast System soil moisture and temperature data. Additionally, real-time green vegetation fraction data from the Visible Infrared Imaging Radiometer Suite will be incorporated into the KMD-WRF runs, once it becomes publicly available from the National Environmental Satellite Data and Information Service. Finally, model verification capabilities will be transitioned to KMD using the Model Evaluation Tools (MET) package, in order to quantify possible improvements in simulated temperature, moisture and precipitation resulting from the experimental land surface initialization. The transition of these MET tools will enable KMD to monitor model forecast accuracy in near real time. This presentation will highlight preliminary verification results of WRF runs over east Africa using the LIS land surface initialization.

Case. Jonathan↗
NASA NTRS2011

Filling-In of Broad Far-Red Solar Lines by Terrestrial Fluorescence and Atmospheric Raman Scattering as Detected by SCIAMACHY Satellite Measurements

Global mapping of terrestrial vegetation fluorescence from space has recently been accomplished with high spectral resolution measurements from the Japanese Greenhouse gases Observing SATellite (GOSAT). These data are of interest because they can potentially provide global information on the functional status of vegetation including light use efficiency and global primary productivity that can be used for global carbon cycle modeling. Quantifying the impact of fluorescence on the O2-A band is important as this band is used for cloud- and aerosol-characterization for other trace-gas retrievals including CO2. Here, we demonstrate that fluorescence information can be derived from space using potentially lower-cost hyperspectral instrumentation, i.e., more than an order of magnitude less spectral resolution than GOSAT, with a relatively simple algorithm. As a demonstration, we use the filling-in of one of the few wide and deep solar Fraunhofer lines in the red and far-red chlorophylla fluorescence bands, the calcium II line near 866 nm, to retrieve fluorescence with the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) satellite instrument. Although the signal from vegetation fluorescence is extremely weak at 866 nm, our results suggest that useful information may be obtained after adjustments are made to the observed spectra to correct for instrumental artifacts. We compare fluorescence from SCIAMACHY with that retrieved at 758 and 770 nm from similarly-corrected GOSAT data as well with the Enhanced Vegetation Index (EVI) from the MODerate-resolution Imaging Spectroradiometer (MODIS). We also show that filling-in occurs at 866 nm over barren areas, possibly originating from luminescent minerals in rock and soil.

Joiner, J.↗
DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗
NASA NTRS2021

Assessing impacts of canopy 3D structure on chlorophyll fluorescence radiance and radiative budget of deciduous forest stands using DART

Although remote sensing (RS) of solar-induced chlorophyll fluorescence (SIF) is increasingly used as a valuable source of information about vegetation photosynthetic activity, the RS SIF observations are significantly influenced by canopy-specific structural features (i.e., canopy architecture including leaf area index and presence of woody components), atmospheric conditions during their acquisition (e.g., proportion of direct and diffuse irradiance) and observational geometric configurations (e.g., sun and viewing directions). Radiative transfer (RT) models have the potential to provide a better understanding of the canopy structural effects on the SIF emission and RS signals. Here, we used the DART model to assess the daily influence, from morning to evening, of forest 3D architecture on SIF nadir radiance, emission, escape factor and nadir yield of eight 100 m × 100 m forest study plots established in a temperate deciduous forest of the Smithsonian Environmental Research Center (Edgewater, MD, USA). The 3D architecture of each plot was derived from airborne LiDAR. DART simulations of these 3D forest plots and their 1D (i.e., vertical profile of sun-adapted and shade-adapted leaves) and 0D (i.e., homogeneous layer of sun-adapted leaves above an homogeneous layer of shade-adapted leaves) abstractions were compared to assess the relative errors (ε1D−3D and ε0D−3D) associated with horizontal and vertical structural heterogeneity, respectively. Forest 3D structure, especially horizontal heterogeneity, had a great influence on forest nadir SIF radiance, resulting in ε1D−3D up to 55% at 8:00 and 18:00 (i.e., for oblique sun directions). The key indicators of this impact, in the descending order of importance, were the SIF escape factor (ε1D−3D up to 40%), the attenuation of incident photosynthetically active radiation (ε1D−3D less than 5%), and the SIF emission yield (ε1D−3D less than 2%). The influence of forest architecture on the nadir SIF escape factor and SIF yield (ε1D−3D up to 40%) varied over time, with differences in forest stand structure, and per spectral domain, being always larger between 640 and 700 nm than between 700 and 850 nm. In addition, woody elements demonstrated a large influence on forest SIF radiance due to their “shading” effect (ε up to 17%) and their “blocking” effect (ε ≈ 10%), both of them higher for far-red than for red SIF. These results underline the importance of 3D forest canopy architecture, especially 2D heterogeneity, and inclusion of woody elements in RT modeling used for interpretation of the RS SIF signal, and subsequently for the estimation of gross primary production and detection of vegetation stress.

Omar Regaieg↗
DOE OSTI2022

Attention-based convolutional capsules for evapotranspiration estimation at scale

Evapotranspiration (ET) measures the amount of water lost from the Earth's surface to the atmosphere and is an integral metric for both agricultural and environmental sciences. Understanding and quantifying ET is critical for achieving effective management of freshwater and irrigation systems. However, current ET estimation models suffer from a trade-off between accuracy and spatial coverage. In this study, we introduce our model Quench, a neural network architecture that achieves highly-accurate ET estimates over large continuous spatial extents. Quench uses our novel Attention-Based Convolutional Capsule for its neural network layers to identify areas of focus and efficiently extract ET information from satellite imagery. Benchmarks that profile our model's performance show substantive improvements in accuracy, with up to 128% increase in accuracy compared to traditional convolutional-based and process-based models. Finally, Quench also demonstrates consistent model performance over high geospatial variability and a diverse array of regions, seasons, climates, and vegetations.

54 ENVIRONMENTAL SCIENCES↗