Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution function”

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 739 records · Page 41

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

Energetic Metastable Oxygen and Nitrogen Atoms in the Terrestrial Atmosphere

This report summarizes our research performed under NASA Grant NAG5-11857. The three-year grant have been supported by the Geospace Sciences SR&T program. We have investigated the energetic metastable oxygen and nitrogen atoms in the terrestrial stratosphere, mesosphere and thermosphere. Hot atoms in the atmosphere are produced by solar radiation, the solar wind and various ionic reactions. Nascent hot atoms arise in ground and excited electronic states, and their translational energies are larger by two - three orders of magnitude than the thermal energies of the ambient gas. The relaxation kinetics of hot atoms determines the rate of atmospheric heating, the intensities of aeronomic reactions, and the rate of atom escape from the planet. Modeling of the non-Maxwellian energy distributions of metastable oxygen and nitrogen atoms have been focused on the determination of their impact on the energetics and chemistry of the terrestrial atmosphere between 25 and 250 km . At this altitudes, we have calculated the energy distribution functions of metastable O and N atoms and computed non-equilibrium rates of important aeronomic reactions, such as destruction of the water molecules by O(1D) atoms and production of highly excited nitric oxide molecules. In the upper atmosphere, the metastable O(lD) and N(2D) play important role in formation of the upward atomic fluxes. We have computed the upward fluxes of the metastable and ground state oxygen atoms in the upper atmosphere above 250 km. The accurate distributions of the metastable atoms have been evaluated for the day and night-time conditions.

Kharchenko, Vasili↗

Spectral Single and Double Power-law Formation by Sequential Particle Acceleration in Flare Magnetic Islands

Spectral single and double power laws are common in high-energy phenomena, such as solar flares and solar energetic particles, including ground level enhancement events. It is not clear what determines the energy breaks and spectral indexes of these power laws. Here, we will describe a first-principles model of pitch-angle and energy distribution function evolution, which produces power laws and provides physical interpretation for such spectral features (Guidoni et al. 2022, ApJ). In this model, a prescribed fraction of particles sequentially “hops” between shrinking magnetic islands (accelerators) formed by flare reconnection. Each accelerator increases particles’ energies by a modest amount, but particles must visit only a few accelerators to increase their energies by orders of magnitude. Data from global magnetohydrodynamic simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. We will also describe the fully analytic method for forming and interpreting power laws, which requires only a few constrained physical parameters of the acceleration region and is independent of the acceleration model, as well as preliminary results extending the analytical model to the formation of double power laws.

S E Guidoni↗

Calculation of the angular radiance distribution for a coupled atmosphere and canopy

The radiative transfer equations for a coupled atmosphere and canopy are solved numerically by an improved Gauss-Seidel iteration algorithm. The radiation field is decomposed into three components: unscattered sunlight, single scattering, and multiple scattering radiance for which the corresponding equations and boundary conditions are set up and their analytical or iterational solutions are explicitly derived. The classic Gauss-Seidel algorithm has been widely applied in atmospheric research. This is its first application for calculating the multiple scattering radiance of a coupled atmosphere and canopy. This algorithm enables us to obtain the internal radiation field as well as radiances at boundaries. Any form of bidirectional reflectance distribution function (BRDF) as a boundary condition can be easily incorporated into the iteration procedure. The hotspot effect of the canopy is accommodated by means of the modification of the extinction coefficients of upward single scattering radiation and unscattered sunlight using the formulation of Nilson and Kuusk. To reduce the computation for the case of large optical thickness, an improved iteration formula is derived to speed convergence. The upwelling radiances have been evaluated for different atmospheric conditions, leaf area index (LAI), leaf angle distribution (LAD), leaf size and so on. The formulation presented in this paper is also well suited to analyze the relative magnitude of multiple scattering radiance and single scattering radiance in both the visible and near infrared regions.

Liang, Shunlin↗

Langmuir Probe Measurements of Inductively Coupled CHF3/Ar and Ar/CHF3/O2 Plasmas

Plasma parameters, such as, electron number density (ne), electron temperature (Te), y electron energy distribution function (EEDF), mean electron energy (Ee), ion number density (ni), and plasma potential (Vp), have been measured by using Langmuir probe in low-pressure (10-50 mTorr) inductively coupled CHF3/Ar and CHF3/Ar/O2 plasmas generated in the GEC cell. The measurements were made at the center of the plasma, keeping the lower electrode grounded, for various CHF3/Ar and Ar/CHF3/O2 mixtures operating at 10-50 mTorr pressures and two input RF power levels, 200 and 300 W. EEDF data show a strong Druyvesteyn distribution with relatively lower number of low energy electrons as compared to a Maxwell distribution and a large electron population with energies higher than the plasma potential. The results further show that at low CHF3 concentrations (less than 50%) the electron number density remains nearly constant with increase in pressure. At higher CHF3 concentrations, however, it decreases with increase in pressure. Plasma potential and electron temperature increase with decrease in pressure and with increase in CHF3 concentration. An analysis of the above observations and mechanisms will be presented.

Sharma, S. P.↗

The AERONET Version 3 Aerosol Retrieval Algorithm, Associated Uncertainties and Comparisons to Version 2

The Aerosol Robotic Network (AERONET) Version 3 (V3) aerosol retrieval algorithm is described, which is based on the Version 2 (V2) algorithm with numerous updates. Comparisons of V3 aerosol retrievals to those of V2 are presented, along with a new approach to estimate uncertainties in many of the retrieved aerosol parameters. Changes in the V3 aerosol retrieval algorithm include (1) a new polarized radiative transfer code (RTC), which replaced the scalar RTC of V2, (2) detailed characterization of gas absorption by adding NO2 and H2O to specify total gas absorption in the atmospheric column, specification of vertical profiles of all the atmospheric species, (3) new bidirectional reflectance distribution function (BRDF) parameters for land sites adopted from the MODIS BRDF/Albedo product, (4) a new version of the extraterrestrial solar flux spectrum, and (5) a new temperature correction procedure of both direct Sun and sky radiance measurements. The potential effect of each change in V3 on single scattering albedo (SSA) retrievals was analyzed. The operational almucantar retrievals of V2 versus V3 were compared for four AERONET sites: GSFC, Mezaira, Mongu, and Kanpur. Analysis showed very good agreement in retrieved parameters of the size distributions. Comparisons of SSA retrievals for dust aerosols (Mezaira) showed a good agreement in 440 nm SSA, while for longer wavelengths V3 SSAs are systematically higher than those of V2, with the largest mean difference at 675 nm due to cumulative effects of both extraterrestrial solar flux and BRDF changes. For non-dust aerosols, the largest SSA deviation is at 675 nm due to differences in extraterrestrial solar flux spectrums used in each version. Further, the SSA 675 nm mean differences are very different for weakly (GSFC) and strongly (Mongu) absorbing aerosols, which is explained by the lower sensitivity to a bias in aerosol scattering optical depth by less absorbing aerosols. A new hybrid (HYB) sky radiance measurement scan is introduced and discussed. The HYB combines features of scans in two different planes to maximize the range of scattering angles and achieve scan symmetry, thereby allowing for cloud screening and spatial averaging, which is an advantage over the principal plane scan that lacks robust symmetry. We show that due to an extended range of scattering angles, HYB SSA retrievals for dust aerosols exhibit smaller variability with solar zenith angles (SZAs) than those of almucantar (ALM), which allows extension of HYB SSA retrievals to SZAs less than 50∘ to as small as 25∘. The comparison of SSA retrievals from closely time-matched HYB and ALM scans in the 50 to 75∘ SZA range showed good agreement with the differences below ∼0.005. We also present an approach to estimate retrieval uncertainties which utilizes the variability in retrieved parameters generated by perturbing both measurements and auxiliary input parameters as a proxy for retrieval uncertainty. The perturbations in measurements and auxiliary inputs are assumed as estimated biases in aerosol optical depth (AOD), radiometric calibration of sky radiances combined with solar spectral irradiance, and surface reflectance. For each set of Level 2 Sun/sky radiometer observations, 27 inputs corresponding to 27 combinations of biases were produced and separately inverted to generate the following statistics of the inversion results: average, standard deviation, minimum and maximum values. From these statistics, standard deviation (labeled U27) is used as a proxy for estimated uncertainty, and a lookup table (LUT) approach was implemented to reduce the computational time. The U27 climatological LUT was generated from the entire AERONET almucantar (1993–2018) and hybrid (2014–2018) scan databases by binning U27s in AOD (440 nm), Angström exponent (AE, 440–870 nm), and SSA (440, 675, 870, 1020 nm). Using this LUT approach, the uncertainty estimates U27 for each individual V3 Level 2 retrieval can be obtained by interpolation using the corresponding measured and inverted combination of AOD, AE, and SSA.

Aliaksandr Siniuk↗

On the use of the noncentral chi-square density function for the distribution of helicopter spectral estimates

A probability density function for the variability of ensemble averaged spectral estimates from helicopter acoustic signals in Gaussian background noise was evaluated. Numerical methods for calculating the density function and for determining confidence limits were explored. Density functions were predicted for both synthesized and experimental data and compared with observed spectral estimate variability.

Garber, Donald P.↗

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗

Application of neural networks to prediction of advanced composite structures mechanical response and behavior

Two types of neural networks were used to evaluate acousto-ultrasonic (AU) data for material characterization and mechanical reponse prediction. The neural networks included a simple feedforward network (backpropagation) and a radial basis functions network. Comparisons of results in terms of accuracy and training time are given. Acousto-ultrasonic (AU) measurements were performed on a series of tensile specimens composed of eight laminated layers of continuous, SiC fiber reinforced Ti-15-3 matrix. The frequency spectrum was dominated by frequencies of longitudinal wave resonance through the thickness of the specimen at the sending transducer. The magnitude of the frequency spectrum of the AU signal was used for calculating a stress-wave factor based on integrating the spectral distribution function and used for comparison with neural networks results.

Cios, K. J.↗

Determination of a Limited Scope Network's Lightning Detection Efficiency

This paper outlines a modeling technique to map lightning detection efficiency variations over a region surveyed by a sparse array of ground based detectors. A reliable flash peak current distribution (PCD) for the region serves as the technique's base. This distribution is recast as an event probability distribution function. The technique then uses the PCD together with information regarding: site signal detection thresholds, type of solution algorithm used, and range attenuation; to formulate the probability that a flash at a specified location will yield a solution. Applying this technique to the full region produces detection efficiency contour maps specific to the parameters employed. These contours facilitate a comparative analysis of each parameter's effect on the network's detection efficiency. In an alternate application, this modeling technique gives an estimate of the number, strength, and distribution of events going undetected. This approach leads to a variety of event density contour maps. This application is also illustrated. The technique's base PCD can be empirical or analytical. A process for formulating an empirical PCD specific to the region and network being studied is presented. A new method for producing an analytical representation of the empirical PCD is also introduced.

Rompala, John T.↗

Systems engingeering for the Kepler Mission : a search for terrestrial planets

The Kepler mission will launch in 2007 and determine the distribution of earth-size planets (0.5 to 10 earth masses) in the habitable zones (HZs) of solar-like stars. The mission will monitor > 100,000 dwarf stars simultaneously for at least 4 years. Precision differential photometry will be used to detect the periodic signals of transiting planets. Kepler will also support asteroseismology by measuring the pressure-mode (p-mode) oscillations of selected stars. Key mission elements include a spacecraft bus and 0.95 meter, wide-field, CCD-based photometer injected into an earth-trailing heliocentric orbit by a 3-stage Delta II launch vehicle as well as a distributed Ground Segment and Follow-up Observing Program. The project is currently preparing for Preliminary Design Review (October 2004) and is proceeding with detailed design and procurement of long-lead components. In order to meet the unprecedented photometric precision requirement and to ensure a statistically significant result, the Kepler mission involves technical challenges in the areas of photometric noise and systematic error reduction, stability, and false-positive rejection. Programmatic and logistical challenges include the collaborative design, modeling, integration, test, and operation of a geographically and functionally distributed project. A very rigorous systems engineering program has evolved to address these challenges. This paper provides an overview of the Kepler systems engineering program, including some examples of our processes and techniques in areas such as requirements synthesis, validation & verification, system robustness design, and end-to-end performance modeling.

systems engineering↗

Effects of a parallel electric field and the geomagnetic field in the topside ionosphere on auroral and photoelectron energy distributions

The consequences of electric field acceleration and an inhomogeneous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one-dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogeneous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of O(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function is investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogeneous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗

Effects of a Parallel Electric Field and the Geomagnetic Field in the Topside Ionosphere on Auroral and Photoelectron Energy Distributions

The consequences of electric field acceleration and an inhomogencous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one- dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogencous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of 0(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function in investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogencous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗

Electrical characterization of Hg(1-x)Cd(x)Te alloys

Theoretical models are described for calculations of charge-carrier concentrations, Fermi energy, and conduction electron mobility as functions of x, temperature, and ionized and neutral defect concentrations of Hg(1-x)Cd(x)Te alloys. Measurements are reported of electron concentration and electron mobility from 5-300 K for alloys with x values of between 0.17 and 0.30. The electrical data are in reasonable agreement with theory, and were analyzed to obtain estimates of donor and acceptor state concentrations. The electron mobilities are calculated in terms of a microscopic theory of electrical conduction derived from the solution of the Boltzmann equation for the perturbed steady-state electron distribution function, and they show that longitudinal optical-phonon and charged and neutral defect scattering are the dominant mobility-limiting mechanisms.

Lehoczky, S. L.↗

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Ott, L.↗

Estimating Bidirectional Reflectance and Monitoring Stability of SNPP-VIIRS Reflective Solar Bands Using A Deep Neural Network

The NASA Clouds and the Earth's Radiant Energy System project provides the scientific community with observed top-of-atmosphere shortwave and longwave fluxes for climate monitoring and climate model validation. To provide consistent VIIRS cloud retrievals, the CERES Imager and Geostationary Calibration Group (IGCG) must understand and quantify the stability of the VIIRS instruments. To achieve this, the IGCG utilizes tropical deep convective clouds (DCCs) as invariant targets. Proper seasonal characterization of the DCC bidirectional reflectance distribution function (BRDF) is key to the success of DCC-based calibration methods, particularly for shortwave infrared (SWIR) bands. This article proposes the use of a deep neural network (DNN) to characterize VIIRS solar reflective band BRDF reflectance, with which individual channel trends are isolated by manipulating the DNN time input. Initial results show that the DNN method can extract statistically significant SNPP-VIIRS band trends, using only SNPP-VIIRS inputs, that are correlative to and match the magnitude of significant trends determined using methods that rely on an external angular distribution model. It may be possible to apply this approach to actively monitor the stability of new instruments without the need for predetermined seasonal BRDF corrections.

Benjamin Scarino↗

The joint far-infrared-optical luminosity function for spiral galaxies and data for the Abell 400 and Cancer clusters

Visual and IRAS data for an optically selected sample of 183 late-type galaxies are compiled in tables and graphs and analyzed in detail to determine the joint FIR-optical luminosity function Psi from the FIR/blue luminosity ratio, r = L(FIR)/L(B). It is found that Psi can be approximated by a function of a single variable psi(r-prime), where r-prime is defined as r times L(B)/L(asterisk) exp -delta, with L(asterisk) a constant and delta = about 0.08. A lognormal curve peaking at r-prime = 0.35 and with dispersion of 0.28 is shown to give a good fit to psi(r-prime). From a lack of galaxies with very low r-prime in the present sample it is inferred that there are few spiral galaxies with low interstellar-dust abundances. Also included are data on the distribution function of r-prime for the more distant clusters Abell 400 and Cancer.

Corbelli, Edvige↗

Local Scale (3-M) Soil Moisture Mapping Using SMAP and Planet Superdove

A capability for mapping meter-level resolution soil moisture with frequent temporal sampling over large regions is essential for quantifying local-scale environmental heterogeneity and eco-hydrologic behavior. However, available surface soil moisture (SSM) products generally involve much coarser grain sizes ranging from 30 m to several 10s of kilometers. Hence a new method is proposed to estimate 3-m resolution SSM using a combination of multi-sensor fusion, machine- learning (ML) and Cumulative Distribution Function (CDF) matching approaches. This method established favorable SSM correspondence between 3-m pixels and overlying 9-km grid cells from overlapping Planet SuperDove (PSD) observations and NASA Soil Moisture Active-Passive (SMAP) mission products. The resulting 3-m SSM predictions showed improved accuracy by reducing ab- solute bias and RMSE by ~0.01 cm3/cm3 over the original SMAP data in relation to in-situ soil moisture measurements for the Australian Yanco region, while preserving the high sampling frequency (1-3 day global revisit) and sensitivity to surface wetness (R 0.865) from SMAP. Heterogeneous soil moisture distributions varying with vegetation biomass gradients and irrigation regimes were generally captured within a selected study area. Further algorithm refinement and implementation for regional applications will allow for improvement in water resources management, precision agriculture, and disaster forecasts and responses.

soil moisture↗