Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source term estimation”

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 289 records · Page 16

A Phenomenon-Based Decomposition of Model-Based Estimates of Boreal Winter ENSO Variability

Climate models are now routinely being used to simulate and predict climate variability on time scales ranging from sub-seasonal to seasonal and longer. As such, there are now long histories of such simulations and predictions spanning multiple decades and multiple ensemble members, both of which are crucial for separating climate signal from climate noise. A key focus of such runs has been the El Niño-Southern Oscillation (ENSO), spurred by recent improvements in our ability to predict such events, though questions remain as to how well climate models do beyond simply always predicting the “canonical” atmospheric response to an ENSO event—something simple statistical models already do reasonably well. This is a critical issue that needs addressing, given the importance of event-to-event differences for predicting regional impacts of ENSO teleconnections, and the need to justify the expense of running sophisticated climate models. Unfortunately, current diagnostic tools are not well suited for quantifying the different sources of variability associated with specific phenomena such as ENSO. More generally, while much effort has focused on addressing model bias, less has been done to address errors in second moment statistics—an issue whose importance is gaining increased attention particularly as we build climate prediction systems capable of taking advantage of forecasts of opportunity—a capability that requires reliable estimates of forecast uncertainty. In this report, we outline a phenomenon-based statistical decomposition of climate variance(in essence a detailed variance budget)that is specifically tailored to address the above questions by separating the variability (both the signal and noise) into that tied to the long-term average impact of a particular phenomenon(the composite mean) and the event-to-event(E2E) variability about the composite mean. In addition, we provide related decompositions of the correlations that allow us to quantify how much of the agreement with observations (the skill) comes from the composite mean as opposed to from the E2Evariability. As an example, we present the results of such a decomposition for ENSO based on simulations with the GEOS atmospheric general circulation model (AGCM), with a focus on the monthly mean impacts over North America during boreal winter(December –March). Here we take advantage of existing GEOS AGCM simulations that were produced as companion simulations to MERRA-2for the period 1980-2016. Comparisons are made throughout with MERRA-2.

ENSO↗

Summary and Annotated Bibliography of Measurement Error Corrections with Potential Application in Future Quesst Mission Community Noise Studies

This document is motivated by likely needs of the Quesst mission community response tests, which will culminate in data collection and estimation of dose-response regression relationships for consideration by domestic and international aviation regulators. Furthermore, basic research questions evaluating interactions between rates of community annoyance, dose levels, and indicators of the presence of rattle, vibration, and startle hinge on hypothesis testing in the context of regression models. For a variety of reasons, noise doses may be known only imprecisely and may not reflect the actual level experienced by responding subjects. These differences between true dose and estimated dose, be they systematic or random, constitute covariate measurement error. Available statistics literature speaks to the impacts of measurement error on regression models, both in terms of bias in estimated coefficients and predicted values, and in terms of the loss of statistical power for hypothesis testing. Given the particulars of a categorical annoyance response variable and a continuous noise dose predictor variable subject to measurement error during testing, the emphasis of this report is on findings and methods pertinent to generalized linear (and mixed) models likely to be employed during the Quesst mission community tests. We reach the following conclusions: 1. Of four reviewed methods, structural Bayesian measurement error models and simulation extrapolation (SIMEX) may be the most readily applicable to Quesst mission community noise study objectives. 2. If warranted, a linear measurement model can help model systematic sources of measurement error that the classical measurement error does not. 3. For its ready implementation and small additional input requirements, simulation extrapolation may be ideally suited for addressing secondary research questions involving interactions between annoyance, noise dose, and other factors through hypothesis testing. 4. For their flexibility and ability to propagate uncertainty, structural Bayesian hierarchical models have great appeal for mission purposes; some care may be needed in developing appropriate probability models describing actual noise exposure during testing. An annotated bibliography logs additional papers and resources that may be of value to analysts in other projects and disciplines.

Dose-Response Model↗

Anomalous heating of the polar E region by unstable plasma waves. II - Theory

It is found that anomalous electron temperatures in the disturbed high-latitude E region can be quantitatively explained in terms of heating by unstable plasma waves. The electron temperatures at 110 km have been measured to be as high as 1500 K instead of the expected value of about 300 K. It is shown that by using quasi-linear theory there is an ample source of heat in the unstable waves and that the measured electron temperature profiles have a shape very similar to what is expected from plasma wave heating by the modified two-stream instability. It is found that there is even more heating going to the ion gas, but that the resulting effect on the ion temperature may be difficult to measure. The best estimate of the wave heating rates leads to the conclusion that wave heating can be as much as 50% of the Joule heating for dc electric field strengths of the order of 45 mV/m or greater.

St.-Maurice, J. P.↗

Orbital precession in Hercules X-1

It is noted that the 35-day X-ray high-low cycle of Her X-1 could be explained in terms of periodic obscuration of the source by a tilted and precessing 'slaved' accretion disk, which would require precession of the mass-shedding star and, consequently, a precessing binary orbit. Four months of Uhuru data obtained in 1972 are analyzed to determine whether the orbital precession can be detected as a Doppler modulation in the arrival times of the 1.24-sec X-ray pulsations. The precession geometry is examined, the test for precession is outlined, and a marginally significant indication of orbital precession in the Hz Her/Her X-1 system is obtained. Possible biases in the analysis are evaluated, and the spin-axis misalignment of HZ Her is estimated from the measured tilt of the orbit.

Deeter, J. E.↗

Infrared polar brightening on Jupiter. III - Spectrometry from the Voyager 1 IRIS experiment

Spectra from the Voyager 1 IRIS experiment confirm the existence of enhanced infrared emission near Jupiter's north magnetic pole in March 1979. The spectral characteristics of the enhanced emission are consistent with a Planck source function. A temperature-pressure profile is derived for the region near the north magnetic pole, from which quantitative abundance estimates of minor species are made. Some species previously detected on Jupiter, including CH3D, C2H2, and C2H6, have been observed again near the pole. Newly discovered species, not previously observed on Jupiter, include C2H4, C3H4, and C6H6. All of these species except CH3D appear to have enhanced abundances at the north polar region with respect to midlatitudes. Upper limits are determined for C4H2 and C3H8. The quantitative results are compared with model calculations based on ultraviolet results from the IUE satellite. The plausibility of the C6H6 identification is discussed in terms of the literature on C2H2 polymerization. The relation of C6H6 to cuprene is also discussed.

Kim, S. J.↗

A Forecast of Reduced Solar Activity and Its Implications for NASA

The "Solar Dynamo" method of solar activity forecasting is reviewed. Known generically as a 'precursor" method, insofar as it uses observations which precede solar activity generation, this method now uses the Solar Dynamo Amplitude (SODA) Index to estimate future long-term solar activity. The peak amplitude of the next solar cycle (#24), is estimated at roughly 124 in terms of smoothed F10.7 Radio Flux and 74 in terms of the older, more traditional smoothed international or Zurich Sunspot number (Ri or Rz). These values are significantly smaller than the amplitudes of recent solar cycles. Levels of activity stay large for about four years near the peak in smoothed activity, which is estimated to occur near the 2012 timeflame. Confidence is added to the prediction of low activity by numerous examinations of the Sun s weakened polar field. Direct measurements are obtained by the Mount Wilson Solar Observatory and the Wilcox Solar Observatory. Further support is obtained by examining the Sun s polar faculae (bright features), the shape of coronal soft X-ray "holes," and the shape of the "source surface" - a calculated coronal feature which maps the large scale structure of the Sun s field. These features do not show the characteristics of well-formed polar coronal holes associated with typical solar minima. They show stunted polar field levels, which are thought to result in stunted levels of solar activity during solar cycle #24. The reduced levels of solar activity would have concomitant effects upon the space environment in which satellites orbit. In particular, the largest influences would affect orbit determination of satellites in LEO (Low Earth Orbit), based upon the altered thermospheric and exospheric densities. A decrease in solar activity would result in smaller satellite decay rates, as well as fewer large solar events that can destroy satellite electronic functions. Other effects of reduced solar activity upon the space environment include enhanced galactic cosmic rays and more space debris at low altitudes (from the decay of old satellite parts, etc.). The reasons are well known: namely, solar activity serves to sweep the inner heliosphere of galactic cosmic rays, and lower exospheric densities result in decreased drag on LEO debris, allowing longer lifetimes.

Schatten, Kenneth↗

Tropical Diabatic Heating and the Role of Convective Processes as Represented in Several Contemporary Climate Models

One of the most fundamental properties of the global heat balance is the net heat input into the tropical atmosphere that helps drive the planetary atmospheric circulation. Although broadly understood in terms of its gross structure and balance of source / sink terms, incorporation of the relevant processes in predictive models is still rather poor. The work reported here examines the tropical radiative and water cycle behavior as produced by four contemporary climate models. Among these are the NSIPP-2 (NASA Seasonal to Interannual Prediction Project) which uses the RAS convective parameterization; the FVCCM, a code using finite volume numerics and the CCM3.6 physics; FVCCM-MCRAS again having the finite volume numerics, but MCRAS convective parameterization and a different radiation treatment; and, finally, the NCEP GSM which uses the RAS. Using multi-decadal integrations with specified SSTs we examine the statistics of radiative / convective processes and associated energy transports, and then estimate model energy flux sensitivities to SST changes. In particular the behavior of the convective parameterizations is investigated. Additional model integrations are performed specifically to assess the importance representing convective inhibition in regulating convective cloud-top structure and moisture detrainment as well as controlling surface energy fluxes. To evaluate the results of these experiments, a number of satellite retrievals are used: TRMM retrievals of vertical reflectivity structure, rainfall rate, and inferred diabatic heating are analyzed to show both seasonal and interannual variations in vertical structure of latent heat release. Top-of-atmosphere radiative fluxes from ERBS and CERES are used to examine shortwave and longwave cloud forcing and to deduce required seasonal energy transports. Retrievals of cloud properties from ISCCP and water vapor variations from SSM/T-2 are also used to understand behavior of the humidity fields. These observations are supplemented with output form the DOE Reanalysis-2.

Robertson, Franklin R.↗

Large Interferometer for Exoplanets (Life) I. Improved Exoplanet Detection Yield Estimates for A Large Mid-Infrared Space-Interferometer Mission

Context. One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale that can spatially separate the signals from exoplanets and their host stars and thus directly scrutinize the exoplanets and their atmospheres. Aims. We seek to quantify the exoplanet detection performance of a space-based mid-infrared (MIR) nulling interferometer that measures the thermal emission of exoplanets. We study the impact of various parameters and compare the performance with that of large single-aperture mission concepts detecting exoplanets in reflected light. Methods. We have developed an instrument simulator that considers all major astrophysical noise sources and coupled it with Monte Carlo of a synthetic exoplanet population around main-sequence stars within 20 pc. This allows us to quantify the number (and types) of exoplanets that our mission concept could detect over a certain time period. Two different scenarios to distribute the observing time among the stellar targets are discussed and different apertures sizes and wavelength ranges are considered. Results. An interferometer consisting of four 2 m apertures working in the 4–18.5 μ.m wavelength range with a total instrument throughput of 5% could detect up to ≈550 exoplanets with radii between 0.5 and 6 R ⊕ with an integrated S/N ≥ 7. At least ≈160 of the detected exoplanets have radii ≤1.5 R ⊕ . Depending on the observing scenario, ≈25–45 rocky exoplanets (objects with radii between 0.5 and 1.5 R ⊕ ) orbiting within the empirical habitable zone (eHZ) of their host stars are among the detections. With four 3.5 m apertures, the total number of detections can increase to up to ≈770, including ≈60–80 rocky eHZ planets. With four times 1 m apertures, the maximum detection yield is ≈315 exoplanets, including ≤20 rocky eHZ planets. The vast majority of small, temperate exoplanets are detected around M dwarfs. The impact of changing the wavelength range to 3–20 μm or 6–17 μm on the detection yield is negligible. Conclusions. A large space-based MIR nulling interferometer will be able to directly detect hundreds of small, nearby exoplanets, tens of which would be habitable world candidates. In terms of predicted detection yield, such a mission can compete with large single-aperture reflected light missions. Given that we assumed only 2.5 years for the search phase, a significant sub-set of the detected exoplanets can be followed-up in a second mission phase to obtain high SNR thermal emission spectra, leveraging the superior diagnostic power of the MIR wavelength regime compared to shorter wavelengths.

planets and satellites: terrestrial planets↗

The relationship of sensor parameters to applications data analysis

A stochastic model for the data acquisition system in a multispectral scanner system, like the one utilized by the LANDSAT satellites, is presented. A list of noise sources which are known or presumed to have a significant effect in the information extraction process was constructed. Since the shot noise introduced by the photodetectors in the sensor system is signal level dependent, an atmospheric model was adopted which could adequately describe the amount of radiation that gets into the sensors based on the atmospheric transmittance. An analysis was carried out to find the output spectral statistics in terms of the input signal statistics and the system parameters. This was integrated into a set of FORTRAN programs that when supplied with, the class statistics, the noise levels introduced by the sensor system, the atmospheric transmittance, and the atmospheric path radiance, can be used to estimate the classification performance. In order to show the beneficts of this model a series of runs were performed in which the Thematic Mapper multispectral scanner was the system under consideration.

Landgrebe, D. A.↗

En route noise levels from propfan test assessment airplane

The en route noise test was designed to characterize propagation of propfan noise from cruise altitudes to the ground. In-flight measurements of propfan source levels and directional patterns were made by a chase plane flying in formation with the propfan test assessment (PTA) airplane. Ground noise measurements were taken during repeated flights over a distributed microphone array. The microphone array on the ground was used to provide ensemble-averaged estimates of mean flyover noise levels, establish confidence limits for those means, and measure propagation-induced noise variability. Even for identical nominal cruise conditions, peak sound levels for individual overflights varied substantially about the average, particularly when overflights were performed on different days. Large day-to-day variations in peak level measurements appeared to be caused by large day-to-day differences in propagation conditions and tended to obscure small variations arising from operating conditions. A parametric evaluation of the sensitivity of this prediction method to weather measurement and source level uncertainties was also performed. In general, predictions showed good agreement with measurements. However, the method was unable to predict short-term variability of ensemble-averaged data within individual overflights. Although variations in absorption appear to be the dominant factor in variations of peak sound levels recorded on the ground, accurate predictions of those levels require that a complete description of operational conditions be taken into account. The comprehensive and integrated methods presented in this paper have adequately predicted ground-measured sound levels. On average, peak sound levels were predicted within 3 dB for each of the three different cruise conditions.

Garber, Donald P.↗

Ciliated protozoa of two antarctic lakes: analysis by quantitative protargol staining and examination of artificial substrates

Planktonic and artificial substrate-associated ciliates have been identified in two perennially ice-covered antarctic lakes of the McMurdo Dry Valleys. Abundances estimated by quantitative protargol staining ranged from < 5 to 31690 cells l-1, levels that are comparable to those previously obtained using other methods. Nineteen ciliate taxa were identified from these lakes, with the most frequently encountered genera being Plagiocampa, Askenasia, Monodinium, Sphaerophrya and Vorticella. The taxonomic findings compare favorably with those of previous investigators; however four previously unreported genera were observed in both Lakes Fryxell and Hoare. The variability in the depth distributions of ciliates in Lake Fryxell is explained in terms of lake physicochemical properties and ciliate prey distributions, while factors related to temporal succession in the Lake Hoare assemblage remain unexplained. Local marine or temperate zone freshwater habitats are a more likely source than the surrounding dry valleys soils for present ciliate colonists in these lakes. Although the taxonomic uncertainties require further examination, our results suggest that ciliate populations in these antarctic lakes undergo significant fluctuations and are more diverse than was previously recognized.

NASA Discipline Exobiology↗

Estimation of Model Error Variances During Data Assimilation

Data assimilation is all about understanding the error characteristics of the data and models that are used in the assimilation process. Reliable error estimates are needed to implement observational quality control, bias correction of observations and model fields, and intelligent data selection. Meaningful covariance specifications are obviously required for the analysis as well, since the impact of any single observation strongly depends on the assumed structure of the background errors. Operational atmospheric data assimilation systems still rely primarily on climatological background error covariances. To obtain error estimates that reflect both the character of the flow and the current state of the observing system, it is necessary to solve three problems: (1) how to account for the short-term evolution of errors in the initial conditions; (2) how to estimate the additional component of error caused by model defects; and (3) how to compute the error reduction in the analysis due to observational information. Various approaches are now available that provide approximate solutions to the first and third of these problems. However, the useful accuracy of these solutions very much depends on the size and character of the model errors and the ability to account for them. Model errors represent the real-world forcing of the error evolution in a data assimilation system. Clearly, meaningful model error estimates and/or statistics must be based on information external to the model itself. The most obvious information source is observational, and since the volume of available geophysical data is growing rapidly, there is some hope that a purely statistical approach to model error estimation can be viable. This requires that the observation errors themselves are well understood and quantifiable. We will discuss some of these challenges and present a new sequential scheme for estimating model error variances from observations in the context of an atmospheric data assimilation system.

Dee, Dick↗

Effects of solar radiation on the orbits of small particles

A modification of the Robertson (1937) equations of particle motion in the presence of solar radiation is developed which allows for partial reflection of sunlight as a result of rapid and varying particle rotations caused by interaction with the solar wind. The coefficients and forces in earlier forms of the equations are compared with those in the present equations, and secular rates of change of particle orbital elements are determined. Orbital dimensions are calculated in terms of time, probable sizes and densities of meteoric and cometary particles are estimated, and times of infall to the sun are computed for a particle moving in an almost circular orbit and a particle moving in an elliptical orbit of high eccentricity. Changes in orbital elements are also determined for particles from a long-period sun-grazing comet. The results show that the time of infall to the sun from a highly eccentric orbit is substantially shorter than from a circular orbit with a radius equal to the mean distance in the eccentric orbit. The possibility is considered that the free orbital kinetic energy of particles drawn into the sun may be the energy source for the solar corona.

Lyttleton, R. A.↗

Atmospheric Airborne Pressure Measurements Using the Oxygen A Band for the ASCENDS Mission

We report on airborne atmospheric pressure measurements using new fiber-based laser technology and the oxygen A-band at 765 nm. Remote measurements of atmospheric temperature and pressure are required for a number of NASA Earth science missions and specifically for the Active Sensing of CO2 Emissions Over Nights, Days, and Seasons (ASCENDS) mission. Accurate measurements of tropospheric CO2 on a global scale are very important in order to better understand its sources and sinks and to improve predictions on any future climate change. The ultimate goal of a CO2 remote sensing mission, such as ASCENDS, is to derive the CO2 concentration in the atmosphere in terms of mole fraction in unit of parts-per-million (ppmv) with regard to dry air. Therefore, both CO2 and the dry air number of molecules in the atmosphere are needed in deriving this quantity. O2 is a stable molecule and uniformly mixed in the atmosphere. Measuring the O2 absorption in the atmosphere can thus be used to infer the dry air number of molecules and then used to calculate CO2 concentration. With the knowledge of atmospheric water vapor, we can then estimate the total surface pressure needed for CO2 retrievals. Our work, funded by the ESTO IIP program, uses fiber optic technology and non-linear optics to generate 765 nm laser radiation coincident with the Oxygen A-band. Our pulsed, time gated technique uses several on- and off-line wavelengths tuned to the O2 absorption line. The choice of wavelengths allows us to measure the pressure by using two adjacent O2 absorptions in the Oxygen A-band. Our retrieval algorithm fits the O2 lineshapes and derives the pressure. Our measurements compare favorably with a local weather monitor mounted outside our laboratory and a local weather station.

Riris, Haris↗

Atmospheric Airborne Pressure Measurements Using the Oxygen A Band for the ASCENDS Mission

We report on airborne atmospheric pressure measurements using new fiber-based laser technology and the oxygen A-band at 765 nm. Remote measurements of atmospheric temperature and pressure are required for a number of NASA Earth science missions and specifically for the Active Sensing of CO2 Emissions Over Nights, Days, and Seasons (ASCENDS) mission. Accurate measurements of tropospheric CO2 on a global scale are very important in order to better understand its sources and sinks and to improve predictions on any future climate change. The ultimate goal of a CO2 remote sensing mission, such as ASCENDS, is to derive the CO2 concentration in the atmosphere in terms of mole fraction in unit of parts-per-million (ppmv) with regard to dry air. Therefore, both CO2 and the dry air number of molecules in the atmosphere are needed in deriving this quantity. O2 is a stable molecule and uniformly mixed in the atmosphere. Measuring the O2 absorption in the atmosphere can thus be used to infer the dry air number of molecules and then used to calculate CO2 concentration. With the knowledge of atmospheric water vapor, we can then estimate the total surface pressure needed for CO2 retrievals. Our work, funded by the ESTO IIP program, uses fiber optic technology and non-linear optics to generate 765 nm laser radiation coincident with the Oxygen A-band. Our pulsed, time gated technique uses several on- and off-line wavelengths tuned to the O2 absorption line. The choice of wavelengths allows us to measure the pressure by using two adjacent O2 absorptions in the Oxygen A-band. Our retrieval algorithm fits the O2 lineshapes and derives the pressure. Our measurements compare favorably with a local weather monitor mounted outside our laboratory and a local weather station.

Riris, Haris↗

Energy yields in the prebiotic synthesis of hydrogen cyanide and formaldehyde

Prebiotic experiments are usually reported in terms of carbon yields, i.e., the yield of product based on the total carbon in the system. These experiments usually involve a large input of energy and are designed to maximize the yields of product. However, large inputs of energy result in multiple activation of the reactants and products. A more realistic prebiotic experiment is to remove the products of the activation step so they are not exposed a second time to the energy source. This is equivalent to transporting the products synthesized in the primitive atmosphere to the ocean, and thereby protecting them from destruction by atmospheric energy sources. Experiments of this type, using lower inputs of energy, give energy yields (moles of products/joule) which can be used to estimate the relative importance of the different energy sources on the primitive earth. Simulated prebiotic atmospheres containing either CH4, CO or CO2 with N2, H2O and variable amounts of H2 were subjected to a high frequency Tesla coil. Samples of the aqueous phase were taken at various time intervals from 1 hr to 7 days, and the energy yields were obtained by extrapolation to zero time. The samples were analyzed for HCN with the cyanide electrode and for H2CO by chromotropic acid. The spark energy was estimated by calorimetry. The temperature rise in an insulated discharge flask was compared with the temperature rise from a resistance heater in the same flask. These results will be compared with calculated production rates of HCN and H2CO from lightning and a number of photochemical processes on the primitive Earth.

Stribling, R.↗

Efficient 1.6 Micron Laser Source for Methane DIAL

Methane is a potent greenhouse gas and on a per molecule basis has a warming influence 72 times that of carbon dioxide over a 20 year horizon. Therefore, it is important to look at near term radiative effects due to methane to develop mitigation strategies to counteract global warming trends via ground and airborne based measurements systems. These systems require the development of a time-resolved DIAL capability using a narrow-line laser source allowing observation of atmospheric methane on local, regional and global scales. In this work, a demonstrated and efficient nonlinear conversion scheme meeting the performance requirements of a deployable methane DIAL system is presented. By combining a single frequency 1064 nm pump source and a seeded KTP OPO more than 5 mJ of 1.6 μm pulse energy is generated with conversion efficiencies in excess of 20%. Even without active cavity control instrument limited linewidths (50 pm) were achieved with an estimated spectral purity of ~95%. Tunable operation over 400 pm (limited by the tuning range of the seed laser) was also demonstrated. This source demonstrated the critical needs for a methane DIAL system motivating additional development of the technology.

Shuman, Timothy↗

Improving Drought Monitoring for the Middle East and North Africa Region via Data Assimilation Using the NASA Land Information System (LIS)

Food and water security become an increasingly pressing concern for countries within the Middle East and North Africa (MENA) region under climate change. Droughts are among the most important issues facing this region in both economic and social terms, which have become much more severe in recent decades. An effective drought monitoring and early warning system is, therefore, critical to support drought impact assessment and risk management. This presentation describes the effort to improve the simulation of hydrological conditions for such semi-arid environment by Noah Multi-parameterization Land Surface Model within the NASA Land Information System (LIS) via assimilating remotely sensed leaf area index and soil moisture in the presence of irrigation. Multi-source satellite evapotranspiration products are used to evaluate the model performance under the impact of the choice of met-forcings, the presence of irrigation and data assimilation. The improved model configuration is aimed to deliver a better soil moisture estimation for the composite drought indicator (CDI) for drought monitoring and to provide a better initial condition in support for hydrological drought forecasting.

Wanshu Nie↗