Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Soil Moisture Estimate Under Forest Using a Semi-Empirical Model at P-Band

Here we present the result of a semi-empirical inversion model for soil moisture retrieval using the three backscattering coefficients: sigma(sub HH), sigma(sub VV) and sigma(sub HV). In this paper we focus on the soil moisture estimate and use the biomass as an ancillary parameter estimated automatically from the algorithm and used as a validation parameter, We will first remind the model analytical formulation. Then we will sow some results obtained with real SAR data and compare them to ground estimates.

soil moisture↗
NASA NTRS2013

Testing the Empirical Shock Arrival Model Using Quadrature Observations

The empirical shock arrival (ESA) model was developed based on quadrature data from Helios (in situ) and P-78 (remote sensing) to predict the Sun-Earth travel time of coronal mass ejections (CMEs). The ESA model requires earthward CME speed as input, which is not directly measurable from coronagraphs along the Sun-Earth line. The Solar Terrestrial Relations Observatory (STEREO) and the Solar and Heliospheric Observatory (SOHO) were in quadrature during 20102012, so the speeds of Earth-directed CMEs were observed with minimal projection effects. We identified a set of 20 full halo CMEs in the field of view of SOHO that were also observed in quadrature by STEREO. We used the earthward speed from STEREO measurements as input to the ESA model and compared the resulting travel times with the observed ones from L1 monitors. We find that the model predicts the CME travel time within about 7.3 h, which is similar to the predictions by the ENLIL model. We also find that CME-CME and CME-coronal hole interaction can lead to large deviations from model predictions.

CMEs↗
NASA NTRS2015

Empirical Relationships Between Optical Properties and Equivalent Diameters of Fractal Soot Aggregates at 550 Nm Wavelength.

Soot aggregates (SAs)-fractal clusters of small, spherical carbonaceous monomers-modulate the incoming visible solar radiation and contribute significantly to climate forcing. Experimentalists and climate modelers typically assume a spherical morphology for SAs when computing their optical properties, causing significant errors. Here, we calculate the optical properties of freshly-generated (fractal dimension Df = 1.8) and aged (Df = 2.6) SAs at 550 nm wavelength using the numericallyexact superposition T-Matrix method. These properties were expressed as functions of equivalent aerosol diameters as measured by contemporary aerosol instruments. This work improves upon previous efforts wherein SA optical properties were computed as a function of monomer number, rendering them unusable in practical applications. Future research will address the sensitivity of variation in refractive index, fractal prefactor, and monomer overlap of SAs on the reported empirical relationships.

deuterium fluorides↗
NASA NTRS2016

Development of Semi-Empirical Damping Equation for Baffled Tank with Oblate Spheroidal Dome

Propellant slosh is a potential source of disturbance that can significantly impact the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring-mass-damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters required by the mechanical model include natural frequency of the slosh, slosh mass, slosh mass center location, and the critical damping ratio. A fundamental study has been undertaken at NASA MSFC to understand the fluid damping physics from a ring baffle in the barrel section of a propellant tank. An asymptotic damping equation and CFD blended equation have been derived by NASA MSFC team to complement the popularly used Miles equation at different flow regimes. The new development has found success in providing a nonlinear damping model for the Space Launch System. The purpose of this study is to further extend the semi-empirical damping equations into the oblate spheroidal dome section of the propellant tanks. First, previous experimental data from the spherical baffled tank are collected and analyzed. Several methods of taking the dome curvature effect, including a generalized Miles equation, area projection method, and equalized fill height method, are assessed. CFD simulation is used to shed light on the interaction of vorticity around the baffle with the locally curved wall and liquid-gas interface. The final damping equation will be validated by a recent subscale test with an oblate spheroidal dome conducted at NASA MSFC.

Yang, H. Q.↗
NASA NTRS2016

Empirical Modeling of the Storm Time Innermost Magnetosphere Using Van Allen Probes and THEMIS Data: Eastward and Banana Currents

The structure of storm time currents in the inner magnetosphere, including its innermost region inside 4R(sub E), is studied for the first time using a modification of the empirical geomagnetic field model TS07D and new data from Van Allen Probes and Time History of Events and Macroscale Interactions during Substorms missions. It is shown that the model, which uses basis-function expansions instead of ad hoc current modules to approximate the magnetic field, consistently improves its resolution and magnetic field reconstruction with the increase of the number of basis functions and resolves the spatial structure and evolution of the innermost eastward current. This includes a connection between the westward ring current flowing largely at R > or approx. 3R(sub E) and the eastward ring current concentrated at R < or approx. 3R(sub E) resulting in a vortex current pattern. A similar pattern coined 'banana current' was previously inferred from the pressure distributions based on the energetic neutral atom imaging and first-principles ring current simulations. The morphology of the equatorial currents is dependent on storm phase. During the main phase, it is complex, with several asymmetries forming banana currents. Near SYM-H minimum, the banana current is strongest, is localized in the evening-midnight sector, and is more structured compared to the main phase. It then weakens during the recovery phase resulting in the equatorial currents to become mostly azimuthally symmetric.

Stephens, G. K.↗
NASA NTRS2017

Security Vulnerability Profiles of Mission Critical Software: Empirical Analysis of Security Related Bug Reports

While some prior research work exists on characteristics of software faults (i.e., bugs) and failures, very little work has been published on analysis of software applications vulnerabilities. This paper aims to contribute towards filling that gap by presenting an empirical investigation of application vulnerabilities. The results are based on data extracted from issue tracking systems of two NASA missions. These data were organized in three datasets: Ground mission IVV issues, Flight mission IVV issues, and Flight mission Developers issues. In each dataset, we identified security related software bugs and classified them in specific vulnerability classes. Then, we created the security vulnerability profiles, i.e., determined where and when the security vulnerabilities were introduced and what were the dominating vulnerabilities classes. Our main findings include: (1) In IVV issues datasets the majority of vulnerabilities were code related and were introduced in the Implementation phase. (2) For all datasets, around 90 of the vulnerabilities were located in two to four subsystems. (3) Out of 21 primary classes, five dominated: Exception Management, Memory Access, Other, Risky Values, and Unused Entities. Together, they contributed from 80 to 90 of vulnerabilities in each dataset.

Goseva-Popstojanova, Katerina↗
NASA NTRS2018

Improving NASA GEOS Atmospheric CO2 Simulations by Calibrating CASA Surface Fluxes with an Empirical Sink

With the adoption of the Paris climate accord, efforts to monitor and understand both anthropogenic and natural carbon sources and sinks are increasing across the world. Given their low latency and global coverage, satellite observations of atmospheric carbon dioxide (CO2) are poised to make important contributions to this field. The combination of satellite data and high resolution global models can be used to monitor changes in carbon fluxes and to evaluate the consistency of nationally reported emissions estimates in support of multiple stakeholder communities. However, a consistent challenge to such work has been the high latency of surface carbon flux estimates, which are often not available for a year or more. This presentation describes the construction of surface carbon flux estimates meant to improve the near real time simulation of atmospheric CO2 with NASA's Goddard Earth Observing System (GEOS) general circulation model. The surface flux estimates begin with a collection of bottom-up fluxes which incorporate satellite measurements in their construction, e.g. vegetation indices in the Carnegie-Ames-Stanford Approach (CASA) and nighttime lights in the Open-source Data Inventory for Anthropogenic CO2 (ODIAC). From there, we take the additional step of using an empirical sink to calibrate terrestrial net biospheric exchange (NBE) to estimated values from atmospheric inversion systems. This approach removes a known, systematic bias in predicted atmospheric mixing ratios. Using these fluxes in a free running simulation, the model is able to reproduce in situ measurements with the same skill as when it uses gridded fluxes from a flux inversion system. Using these fluxes as a prior in an assimilation system, e.g. one incorporating retrievals of column CO2 from the Orbiting Carbon Observatory 2 (OCO-2), allows the analysis to capture variability in CO2 on scales that would be missed otherwise. This approach supports NASA's capability to forecast atmospheric CO2 up to two weeks in advance by leveraging a GEOS system used to produce quasi-operational weather analyses and forecasts, providing a valuable new tool to the carbon monitoring research and applications communities.

Weir, B.↗
NASA NTRS2019

The CO2-Broadened H2O Continuum in the 100-1500 cm(sup -1) Region: Measurements, Predictions and Empirical Model

Transmission spectra of H2O+CO2 mixtures have been recorded, at 296, 325 and 366 K, for various pressures and mixture compositions using two experimental setups. Their analysis enables to retrieve values of the “continuum” absorption by the CO2-broadened H2O line wings between 100 and 1500 cm(exp -1). The results are in good agreement with those, around 1300 cm(exp -1), of the single previous experimental study available. Comparisons are also made with direct predictions based on line-shape correction factors X calculated, almost thirty years ago, using a quasistatic approach and an input H2O-CO2 intermolecular potential. They show that this model quite nicely predicts, with slightly overestimated values, the continuum over a spectral range where it varies by more than three orders of magnitude. An empirical correction is proposed, based on the experimental data, which should be useful for radiative transfer and climate studies in CO2 rich planetary atmospheres.

CO2↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗
NASA NTRS2019

Empirical Absolute Calibration Model for Multiple Pseudo-Invariant Calibration Sites

This work extends an empirical absolute calibration model initially developed for the Libya 4 Pseudo-Invariant Calibration Site (PICS) to five additional Saharan Desert PICS (Egypt 1, Libya 1, Niger 1, Niger 2, and Sudan 1), and demonstrates the efficacy of the resulting models at predicting sensor top-of-atmosphere (TOA) reflectance. It attempts to generate absolute calibration models for these PICS that have an accuracy and precision comparable to or better than the current Libya 4 model, with the intent of providing additional opportunities for sensor calibration. In addition, this work attempts to validate the general applicability of the model to other sites. The method uses Terra Moderate Resolution Imaging Spectroradiometer (MODIS) as the reference radiometer and Earth Observing-1 (EO-1) Hyperion image data to provide a representative hyperspectral reflectance profile of the PICS. Data from a region of interest (ROI) in an “optimal region” of 3% temporal, spatial, and spectral stability within the PICS are used for developing the model. The developed models were used to simulate observations of the Landsat 7 (L7) Enhanced Thematic Mapper Plus (ETM+), Landsat 8 (L8) Operational Land Imager (OLI), Sentinel 2A (S2A) MultiSpectral Instrument (MSI) and Sentinel 2B (S2B) MultiSpectral Instrument (MSI) from their respective launch date through 2018. The models developed for the Egypt 1, Libya 1 and Sudan 1 PICS have an estimated accuracy of approximately 3% and precision of approximately 2% for the sensors used in the study, comparable to the current Libya 4 model. The models developed for the Niger 1 and Niger 2 sites are significantly less accurate with similar precision.

Raut, Bipin↗
NASA NTRS2021

Positional Dependence of SNPP VIIRS Solar Diffuser BRDF Change Factor: An Empirical Approach

The Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership satellite regularly performs on-orbit radiometric calibration of its reflective solar bands (RSBs), primarily through observations of an onboard sunlit solar diffuser (SD). The on-orbit change of the SD bidirectional reflectance distribution function (BRDF) value, quantified by a numerical factor called the H-factor, is determined by the onboard SD stability monitor. Our previous study showed that the H-factor is solar angle and view direction dependent. In this study, we determine the dependence of the H-factor on the detector SD view footprint location. We fit an empirical model to the NASA Collection 1 SNPP VIIRS Level 1B (L1B) spectral reflectance difference across the detectors in an RSB over uniform Earth scenes of the Libya 4 desert and deep convective clouds. We apply the model predicted SD positional dependent H-factor to calibrate the RSBs. Under this new calibration scheme, the original unreal striping is removed from the homogeneous Libya 4 desert and the deep convective cloud images, as well as the original unreal striping from the Dunhuang desert image. The SD positional dependent H-factor has been used to calculate the SNPP VIIRS RSB radiometric correction factor for the NASA Collection 2.0 SNPP VIIRS L1B products.

SNPP VIIRS↗
NASA NTRS2021

Empirical Model of 10 – 130 MeV Solar Energetic Particle Spectra at 1 AU Based on Coronal Mass Ejection Speed and Direction

We present a new empirical model to predict solar energetic particle (SEP) event integrated and peak intensity spectra between 10 and 130 MeV at 1 AU, based on multi-point spacecraft measurements from the Solar TErrestrial RElations Observatory (STEREO), the Geostationary Operational Environmental Satellites (GOES), and the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) satellite experiment. The analyzed data sample includes 32 SEP events occurring between 2010 and 2014, with a statistically significant proton signal at energies in excess of a few tens of MeV, unambiguously recorded at three spacecraft locations. The spatial distributions of SEP intensities are reconstructed by assuming an energy-dependent 2D Gaussian functional form, and accounting for the correlation between the intensity and the speed of the parent coronal mass ejection (CME), and the magnetic-field-line connection angle. The CME measurements used are from the Space Weather Database Of Notifications, Knowledge, Information (DONKI). The model performance, including its extrapolations to lower/higher energies, is tested by comparing with the spectra of 20 SEP events not used to derive the model parameters. Despite the simplicity of the model, the observed and predicted event-integrated and peak intensities at Earth and at the STEREO spacecraft for these events show remarkable agreement, both in the spectral shapes and their absolute values.

Alessandro Bruno↗
NASA NTRS2021

A Semi-Empirical Scheme for Bathymetric Mapping in Shallow Water by ICESat-2 and Sentinel-2: A Case Study in the South China Sea

To derive shallow water bathymetry for coastal areas, a common approach is to deploy a scanning airborne bathymetric light detection and ranging (LiDAR) system or a shipborne echosounder for ground surveys. However, recent advancements in satellite remote sensing, including the Ice, Cloud and land Elevation Satellite-2 (ICESat-2) offer new tools for generating satellite derived bathymetry (SDB). The key payload onboard ICESat-2 is the Advanced Topographic Laser Altimeter System (ATLAS), a micro-pulse, photon-counting LiDAR system, simultaneously emitting six separate 532 nm beams at 10 kHz pulse rate. However, despite its high resolution, the major limitation for bathymetry is that ICESat-2 only provides along-track height profiles, leaving observation gaps between the parallel ground tracks. Merging ICESat-2 observations with optical multispectral imagery, as demonstrated herein, provides an effective solution for deriving a full scene of water depth in light of the spectral attenuation behavior. This study aims to combine ICESat-2 and Sentinel-2 optical data to derive shallow water bathymetry (depth <20 m) at six islands and reefs in the South China Sea. ICESat-2 ATL03 point clouds of georeferenced photons are first filtered to determine the seafloor elevation along the ground track. Results indicate a root-mean-square error (RMSE) of 0.26–0.61 m as compared with independent observations from an airborne LiDAR campaign. Next, three semi-empirical functions, namely the Modified Linear/Polynomial/Exponential Ratio Models with its kernel formed by the log ratio between Sentinel-2′s green and blue bands, are used to fit the spectral data with ICESat-2 height profiles. After water depth mapping using the trained model, independent ICESat-2 point clouds are used to validate the Sentinel-2 derived bathymetry. The RMSE values of the three models using the weighted average of multiple images for these six islands are within 0.50–0.90 m in 0–15 m deep. We also demonstrate that a synthesis of satellite laser altimetry and optical remote sensing can produce SDB results that potentially meet the requirement of category C in Zones of Confidence (ZOC) of the Electronic Navigational Chart (ENC) in 0–8 m deep. It is foreseen that ICESat-2 will be a helpful tool for mapping coastal and shallow waters around the world especially where bathymetric data are unavailable.

Coastal Bathymetry↗
NASA NTRS2022

An Empirical Algorithm to Map Perennial Firn Aquifers and Ice Slabs Within the Greenland Ice Sheet Using Satellite L-band Microwave Radiometry

Perennial firn aquifers are subsurface meltwater reservoirs consisting of a meters-thick water-saturated firn layer that can form on spatial scales as large as tens of kilometers. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting and high snow accumulation. Widespread perennial firn aquifers have been identified within the Greenland Ice Sheet (GrIS) via field expeditions, airborne ice-penetrating radar surveys, and satellite microwave sensors. In contrast, ice slabs are nearly continuous ice layers that can also form on spatial scales as large as tens of kilometers as a result of surface and subsurface water-saturated snow and firn layers sequentially refreezing following multiple melting seasons. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting but in areas where snow accumulation is at least 25 % lower as compared to perennial firn aquifer areas. Widespread ice slabs have recently been identified within the GrIS via field expeditions and airborne ice-penetrating radar surveys, specifically in areas where perennial firn aquifers typically do not form. However, ice slabs have yet to be identified from space. Together, these two ice sheet features represent distinct, but related, sub-facies within the broader percolation facies of the GrIS that can be defined primarily by differences in snow accumulation, which influences the englacial hydrology and thermal characteristics of firn layers at depth. Here, for the first time, we use enhanced-resolution vertically polarized L-band brightness temperature () imagery (2015–2019) generated using observations collected over the GrIS by NASA's Soil Moisture Active Passive (SMAP) satellite to map perennial firn aquifer and ice slab areas together as a continuous englacial hydrological system. We use an empirical algorithm previously developed to map the extent of Greenland's perennial firn aquifers via fitting exponentially decreasing temporal L-band signatures to a set of sigmoidal curves. This algorithm is recalibrated to also map the extent of ice slab areas using airborne ice-penetrating radar surveys collected by NASA's Operation IceBridge (OIB) campaigns (2010–2017). Our SMAP-derived maps show that between 2015 and 2019, perennial firn aquifer areas extended over 64 000 km2, and ice slab areas extended over 76 000 km2. Combined together, these sub-facies are the equivalent of 24 % of the percolation facies of the GrIS. As Greenland's climate continues to warm, seasonal surface melting will increase in extent, intensity, and duration. Quantifying the possible rapid expansion of these sub-facies using satellite L-band microwave radiometry has significant implications for understanding ice-sheet-wide variability in englacial hydrology that may drive meltwater-induced hydrofracturing and accelerated ice flow as well as high-elevation meltwater runoff that can impact the mass balance and stability of the GrIS.

Julie Z Miller↗
NASA NTRSDate not supplied

Photometric Correction of Hayabusa2's NIRS3 Spectra of Asteroid 162173 Ryugu Using Empirical Photometric Models

In 2018, Japanese Aerospace Exploration Agency’s spacecraft Hayabusa2 began a near infrared spectroscopic imaging survey of near-Earth asteroid 162173 Ryugu. Hayabusa2 is a successful sample-return mission with an overarching goal to provide a better understanding of the origin and evolution of our solar system. The target, Ryugu, is a low-albedo carbonaceous asteroid (Cb-type) that is linked to carbonaceous chondrite meteorites. It is thought to have originated in the main asteroid belt and migrated inward to become a near-Earth asteroid, and therefore reachable by spacecraft. Preliminary findings from processed NIRS3data have shown that hydroxyl-bearing minerals are present on the surface of Ryugu, and it is likely the result of impact fragments from an aqueously altered parent body.In this project, we have used newly calibrated and processed NIRS3 data using updated shape models of Ryugu. NIRS3 spectra have both a thermal and reflectance component. The thermal component (beyond 2.5 microns) was modeled and removed from all NIRS3 spectra. Ryugu spectra were taken at different viewing geometries, and a photometric model needed to be developed to normalize all of the spectra at a common geometry for each mission phase. We have used three empirical models: Minnaert, Lommel-Seeliger, and ROLO (RObotic LunarOrbiter). These models were chosen for their ability to relate the surface reflectance to the viewing geometry, as well as their compatibility with the asteroid’s albedo range. The three models provide the global light scattering properties of Ryugu’s surface and subsequently enable us to calculate the geometric albedo, phase integral, spherical bond albedo, and the average surface normal albedo for Ryugu.

Lucille Grace Williamson↗
NASA NTRS2023

An Evaluation of A Physics-Based Firn Model and A Semi-Empirical Firn Model Across the Greenland Ice Sheet (1980–2020)

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (−66.6 km 3 yr −1 ) in SNOWPACK and 1.5 % (−17.4 km 3 yr −1 ) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage.

Megan Thompson-Munson↗
NASA NTRSDate not supplied

Machine Learning Emulators and Empirical Models Combining Climate and Global Crop Models for Seasonal Agricultural Production

We present results from several connected efforts to apply machine learning methods to estimates of seasonal agricultural production anomalies around the world. First, we apply the XGBoost Random Forest method to fit emulators that mimic global crop models participating in the Agricultural Model Intercomparison and Improvement Project (AgMIP) Global Gridded Crop Model Intercomparison (GGCMI). These are the same models used in the agricultural sector simulations of the Inter-Sectoral Impacts Model Intercomparison Project (ISIMIP). These emulators use 8 climate variables split across 5 sub-seasonal representations of the growing season for each ½ degree grid cell around the world for maize, wheat, rice and soybeans. Emulators are useful for estimating conditions that have not already been simulated by GGCMI (e.g., in a seasonal prediction model) and also to diagnose model differences and capabilities. For example, emulators of the pDSSAT maize model tend to be more reliant on mean temperatures than the LPJmL model, and few models have strong responses to cold extremes. Second, we use a similar XGBoost approach to fit empirical models for national production data for the top 20 producing countries according to the United Nations Food and Agricultural Organization (FAO). Models utilize both climate observations and the GGCM models as predictors, resulting in skillful models for many (but not all) top producing-countries. The patterns of climate and crop model features selected indicate regions and systems that are better or worse simulated by the GGCMs. For example, information in cold extreme predictors is often combined with GGCM output predictors to provide sensitivity that models may underrepresent.

machine learning↗
NASA NTRS2023

Comparative Analysis of Empirical and Machine Learning Models for Chla Extraction Using Sentinel-2 and Landsat OLI Data: Opportunities, Limitations, and Challenges

Remote retrieval of near-surface chlorophyll-a (Chla) concentration in small inland waters is challenging due to substantial optical interferences of various water constituents and uncertainties in the atmospheric correction (AC) process. Although various algorithms have been developed to estimate Chla from moderate-resolution terrestrial missions (∼10–60 m), the production of both accurate distribution maps and time series of Chla has proven challenging, limiting the use of remote analyses for lake monitoring. Here, we develop a support vector regression (SVR) model, which uses satellite-derived remote-sensing reflectance spectra () from Sentinel-2 and Landsat-8 images as input for Chla retrieval in a representative eutrophic prairie lake, Buffalo Pound Lake (BPL), Saskatchewan, Canada. Validated against in situ Chla from seven ice-free seasons (N ∼ 200; 2014–2020), the SVR model outperformed both locally tuned, -fed empirical models (Normalized Difference Chlorophyll Index, 2- and 3-band, and OC3) and Mixture Density Networks (MDNs) by 15–65%, while exhibiting comparable performance to a locally trained MDN, with an error of ∼35%. Comparison of Chla retrieval models, AC processors (iCOR, ACOLITE), and radiometric products (Rayleigh-corrected, surface, and top-of-atmosphere reflectance) showed that the best Chla maps and optimal time series (up to 100 mg m−3) were produced using a coupled SVR-iCOR system.

algal blooms↗