SEARCH · Engineering Papers

Results for “Spatial statistics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,729 records · Page 96

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Theorem on extensive spectral degeneracy for systems with rigid higher symmetries in general dimensions

We establish, in the spirit of the Lieb-Schultz-Mattis theorem, lower bounds on the spectral degeneracy of quantum systems with higher (gaugelike) symmetries with rather generic physical boundary conditions in an arbitrary number of spatial dimensions. Contrary to applying twists or equivalent adiabatic operations, we exploit the effects of modified boundary conditions. We find that when a general choice of boundary geometry is immaterial in determining spectral degeneracies while approaching the thermodynamic limit, systems that exhibit rigid noncommuting gaugelike symmetries, such as the orbital compass model, must have an exponential (in the size of the boundary) degeneracy of each of their spectral levels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

The Missing Satellite Problem outside of the Local Group. II. Statistical Properties of Satellites of Milky Way–like Galaxies

We present a new observation of satellite galaxies around seven Milky Way (MW)–like galaxies located outside of the Local Group (LG) using Subaru/Hyper Suprime-Cam imaging data to statistically address the missing satellite problem. We select satellite galaxy candidates using magnitude, surface brightness, Sérsic index, axial ratio, FWHM, and surface brightness fluctuation cuts, followed by visual screening of false positives such as optical ghosts of bright stars. We identify 51 secure dwarf satellite galaxies within the virial radius of nine host galaxies, two of which are drawn from the pilot observation presented in Paper I. We find that the average luminosity function of the satellite galaxies is consistent with that of the MW satellites, although the luminosity function of each host galaxy varies significantly. We observe an indication that more massive hosts tend to have a larger number of satellites. Physical properties of the satellites such as the size–luminosity relation are also consistent with the MW satellites. However, the spatial distribution is different; we find that the satellite galaxies outside of the LG show no sign of concentration or alignment, while that of the MW satellites is more concentrated around the host and exhibits a significant alignment. As we focus on relatively massive satellites with M V < –10, we do not expect that the observational incompleteness can be responsible here. This trend might represent a peculiarity of the MW satellites, and further work is needed to understand its origin.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1990

Summary of results and conclusions based on analysis of volume imaging and high spectral resolution lidar data acquired during FIRE phase 1, part 1

The collection of long term global statistics on cloud cover may be most easily accomplished with satellite based observations; however, measurements derived from passive satellite retrieval methods must be calibrated and verified by in situ or ground based remote sensor observations. Verification is not straight forward, however, because the highly variable nature of cloud altitude, morphology, and optical characteristics complicates the scaling of point measurements to satellite footprint sized areas. This is particularly evident for cirrus clouds which may be organized on horizontal scales of 10's of meters to 8 km or more, and have optical depths ranging from less than .003 to greater than 3. Cirrus clouds can strongly influence earths' radiative balance, but, because they are often transmissive, cirrus clouds are difficult to detect and characterize from satellite measurements. Because of its precise ranging capabilities, spatial resolution and sensitivity, lidar observations have played an important role in the detection, depiction, and characterization of cirrus clouds. Some of the characteristics of cirrus clouds are summarized which observed the High Spectral Resolution and Volume Imaging Lidars during the phase 1 IFO and ETO periods.

Grund, Christian J.↗

NASA NTRS1997

Foreshock Langmuir Waves for Unusually Constant Solar Wind Conditions: Data and Implications for Foreshock Structure

Plasma wave data are compared with ISEE 1's position in the electron foreshock for an interval with unusually constant (but otherwise typical) solar wind magnetic field and plasma characteristics. For this period, temporal variations in the wave characteristics can be confidently separated from sweeping of the spatially varying foreshock back and forth across the spacecraft. The spacecraft's location, particularly the coordinate D(sub f) downstream from the foreshock boundary (often termed DIFF), is calculated by using three shock models and the observed solar wind magnetometer and plasma data. Scatterplots of the wave field versus D(sub f) are used to constrain viable shock models, to investigate the observed scatter in the wave fields at constant D(sub f), and to test the theoretical predictions of linear instability theory. The scatterplots confirm the abrupt onset of the foreshock waves near the upstream boundary, the narrow width in D(sub f) of the region with high fields, and the relatively slow falloff of the fields at large D(sub f), as seen in earlier studies, but with much smaller statistical scatter. The plots also show an offset of the high-field region from the foreshock boundary. It is shown that an adaptive, time-varying shock model with no free parameters, determined by the observed solar wind data and published shock crossings, is viable but that two alternative models are not. Foreshock wave studies can therefore remotely constrain the bow shock's location. The observed scatter in wave field at constant D(sub f) is shown to be real and to correspond to real temporal variations, not to unresolved changes in D(sub f). By comparing the wave data with a linear instability theory based on a published model for the electron beam it is found that the theory can account qualitatively and semiquantitatively for the abrupt onset of the waves near D(sub f) = 0, for the narrow width and offset of the high-field region, and for the decrease in wave intensity with increasing D(sub f). Quantitative differences between observations and theory remain, including large overprediction of the wave fields and the slower than predicted falloff at large D(sub f) of the wave fields. These differences, as well as the unresolved issue of the electron beam speed in the high-field region of the foreshock, are discussed. The intrinsic temporal variability of the wave fields, as well as their overprediction based on homogeneous plasma theory, are indicative of stochastic growth physics, which causes wave growth to be random and varying in sign, rather than secular.

Cairns, Iver H.↗

NASA NTRS2003

Application of a data-mining method based on Bayesian networks to lesion-deficit analysis

Although lesion-deficit analysis (LDA) has provided extensive information about structure-function associations in the human brain, LDA has suffered from the difficulties inherent to the analysis of spatial data, i.e., there are many more variables than subjects, and data may be difficult to model using standard distributions, such as the normal distribution. We herein describe a Bayesian method for LDA; this method is based on data-mining techniques that employ Bayesian networks to represent structure-function associations. These methods are computationally tractable, and can represent complex, nonlinear structure-function associations. When applied to the evaluation of data obtained from a study of the psychiatric sequelae of traumatic brain injury in children, this method generates a Bayesian network that demonstrates complex, nonlinear associations among lesions in the left caudate, right globus pallidus, right side of the corpus callosum, right caudate, and left thalamus, and subsequent development of attention-deficit hyperactivity disorder, confirming and extending our previous statistical analysis of these data. Furthermore, analysis of simulated data indicates that methods based on Bayesian networks may be more sensitive and specific for detecting associations among categorical variables than methods based on chi-square and Fisher exact statistics.

NASA Discipline Neuroscience↗

NASA NTRS2012

Global CO2 Distributions over Land from the Greenhouse Gases Observing Satellite (GOSAT)

January 2009 saw the successful launch of the first space-based mission specifically designed for measuring greenhouse gases, the Japanese Greenhouse gases Observing SATellite (GOSAT). We present global land maps (Level 3 data) of column-averaged CO2 concentrations (X(sub CO2)) derived using observations from the GOSAT ACOS retrieval algorithm, for July through December 2009. The applied geostatistical mapping approach makes it possible to generate maps at high spatial and temporal resolutions that include uncertainty measures and that are derived directly from the Level 2 observations, without invoking an atmospheric transport model or estimates of CO2 uptake and emissions. As such, they are particularly well suited for comparison studies. Results show that the Level 3 maps for July to December 2009 on a lO x 1.250 grid, at six-day resolution capture much of the synoptic scale and regional variability of X(sub CO2), in addition to its overall seasonality. The uncertainty estimates, which reflect local data coverage, X(sub CO2) variability, and retrieval errors, indicate that the Southern latitudes are relatively well-constrained, while the Sahara Desert and the high Northern latitudes are weakly-constrained. A probabilistic comparison to the PCTM/GEOS-5/CASA-GFED model reveals that the most statistically significant discrepancies occur in South America in July and August, and central Asia in September to December. While still preliminary, these results illustrate the usefulness of a high spatiotemporal resolution, data-driven Level 3 data product for direct interpretation and comparison of satellite observations of highly dynamic parameters such as atmospheric CO2.

Hammerling, Dorit M.↗

DOE OSTI2026

GPM IMERG V07B and V06B: Evaluation Using Ground-Based Radar Observations and Application in Global Mesoscale Convective System Tracking

This study evaluates the latest Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals for GPM (IMERG V07B) against its predecessor V06B, for studying mesoscale convective systems (MCSs). Both versions are compared using ground-based radar and rain gauge data from five meteorologically diverse regions: the contiguous United States (including eastern coastlines), Amazon rainforest, central Argentina mountains, equatorial Indian Ocean, and northern Australia across multiple temporal (0.5–6 hours) and spatial scales (0.1°–0.25°). An updated global MCS tracking dataset is developed by integrating satellite-observed infrared brightness temperature with IMERG V07B. Comparation of IMERG against radar observations reveals that IMERG demonstrates better performance in capturing the probability distribution and quantitative contributions of rainfall (from no-rain to intense-rain conditions) over tropical oceans than over land, with marked improvements in IMERG V07B for heavy-to-intense rain (> 10 mm h-1). Over land, systematic biases persist: IMERG tends to overestimate light-to-moderate rain (1–10 mm h-1) while underestimating heavy-to-intense rain. Additionally, aggregating IMERG to coarser resolutions (3-hourly or 0.25°) improves consistency with radar observations, outperforming the 1-hourly/0.1° resolution. The new IMERG V07B-based global MCS dataset exhibits consistent statistical characteristics with the V06B-based dataset, despite lower mean rain rates and reduced heavy precipitation contributions. These findings offer valuable insights for utilizing IMERG V07B in global precipitation studies, MCS characterization, and model evaluation.

Zhang, Sihan↗

DOE OSTI2024

Data and scripts associated with a manuscript on residence time distribution simulation in two 10-kilometer long river sections

This data package is associated with the publication “On the Transferability of Residence Time Distributions in Two 10-km Long River Sections with Similar Hydromorphic Units” submitted to the Journal of Hydrology (Bao et al. 2024).Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface, along with their residence time distributions (RTDs) in the subsurface, is crucial for managing water quality and ecosystem health in dynamic river corridors. However, directly simulating high-spatial resolution HEFs and RTDs can be a time-consuming process, particularly for watershed-scale modeling. Efficient surrogate models that link RTDs to hydromorphic units (HUs) may serve as alternatives for simulating RTDs in large-scale models. One common concern with these surrogate models, however, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this, we evaluated the HEFs and the resulting RTD-HU relationships for two 10-kilometer-long river corridors along the Columbia River, using a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework that we previously developed. Applying this framework to the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. This data package includes the model inputs files and the simulation results data. This data package contains 10 folders. The modeling simulation results data are in the folders 100H_pt_data and 300area_pt_data, for the study domain Hanford 100H and 300 area respectively. The remaining eight folders contain the scripts and data to generate the manuscript figures. The file-level metadata file (Bao_2024_Residence_Time_Distribution _flmd.csv) includes a list of all files contained in this data package and descriptions for each. The data dictionary file (Bao_2024_Residence_Time_Distribution _dd.csv) includes column header definitions and units of all tabular files.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Analysis of Polarized Dust Emission Using Data from the First Flight of SPIDER

Using data from the first flight of Spider and from the Planck High Frequency Instrument, we probe the properties of polarized emission from interstellar dust in the Spider observing region. Component-separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses of diffuse Galactic dust emission spanning the full Spider region demonstrate (i) a spectral energy distribution that is broadly consistent with a modified-blackbody (MBB) model with a spectral index of β d = 1.45 ± 0.05 (1.47 ± 0.06) for E (B)-mode polarization, slightly lower than that reported by Planck for the full sky; (ii) an angular power spectrum broadly consistent with a power law; and (iii) no significant detection of line-of-sight polarization decorrelation. Tests of several modeling uncertainties find only a modest impact (~10% in σ r ) on Spider's sensitivity to the cosmological tensor-to-scalar ratio. The size of the Spider region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature T d = 19.6 K, an analysis of two independent subregions of that field results in inferred values of β d = 1.52 ± 0.06 and β d = 1.09 ± 0.09, which are inconsistent at the 3.9σ level. Furthermore, a joint analysis of Spider and Planck 217 and 353 GHz data within one subregion is inconsistent with a simple MBB at more than 3σ, assuming a common morphology of polarized dust emission over the full range of frequencies. This evidence of variation may inform the component-separation approaches of future cosmic microwave background polarization experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1983

Soil temperature investigations using satellite acquired thermal-infrared data in semi-arid regions

Thermal-infrared data from the Heat Capacity Mapping Mission satellite were used to map the spatial distribution of diurnal surface temperatures and to estimate mean annual soil temperatures (MAST) and annual surface temperature amplitudes (AMP) in semi-arid east central Utah. Diurnal data with minimal snow and cloud cover were selected for five dates throughout a yearly period and geometrically co-registered. Rubber-sheet stretching was aided by the WARP program which allowed preview of image transformations. Daytime maximum and nighttime minimum temperatures were averaged to generation average daily temperature (ADT) data set for each of the five dates. Five ADT values for each pixel were used to fit a sine curve describing the theoretical annual surface temperature response as defined by a solution of a one-dimensinal heat flow equation. Linearization of the equation produced estimates of MAST and AMP plus associated confidence statistics. MAST values were grouped into classes and displayed on a color video screen. Diurnal surface temperatures and MAST were primarily correlated with elevation.

Day, R. L.↗

NASA NTRS1999

Sea Surface Slope Statistics for Intermediate and Shore Scale Ocean Waves Measured Using a Low-Altitude Aircraft

Ocean surface remote sensing techniques often rely on scattering or emission linked to shorter- scale gravity-capillary ocean wavelets. However, it is increasingly apparent that slightly longer wavelengths of O(10 to 500 cm) are vital components in the robust sea surface description needed to link varied global remote sensing data sets. This paper describes a sensor suite developed to examine sea surface slope variations in the field using an aircraft flying at very low altitude (below 30 m) and will also provide preliminary measurements detailing changes in slope characteristics versus sea state and friction velocity. Two-dimensional surface slope is measured using simultaneous range measurements from three compact short-range laser altimeters mounted in an equilateral triangle arrangement with spacing of about 1 m. In addition, all three lasers provide independent wave elevation profiles after GPS-aided correction for aircraft altitude. Laser range precision is 1 cm rms while vertical motion correction is 15 cm rms. The measurements are made along-track at approximately 1 m intervals setting the spatial scale of the measurement to cover waves of intermediate to long scale. Products available for this array then include surface elevation, two-dimensional slope distribution, and the cross- and along-track 1-D slope distributions. To complement the laser, a down-looking mm-wave radar scatterometer is centered within the laser array to measure radar backscatter simultaneously with the laser slope. The radar's footprint is nominally 1 m in diameter. Near-vertical radar backscatter is inversely proportional to the small-scale surface slope variance and to the tilt of the underlying (laser-measured) surface facet. Together the laser and radar data provide information on wave roughness from the longest scales down to about 1 cm. These measurements are complemented by aircraft turbulence probe data that provides robust surface flux information.

Vandemack, Douglas↗

NASA NTRS2019

Evaluation of MODIS and VIIRS Cloud-Gap-Filled Snow-Cover Products for Production of an Earth Science Data Record

MODerate resolution Imaging Spectroradiometer (MODIS) cryosphere products have been available since 2000 – following the 1999 launch of the Terra MODIS and the 2002 launch of the Aqua MODIS – and include global snow-cover extent (SCE) (swath, daily, and 8 d composites) at 500 m and ∼5 km spatial resolutions. These products are used extensively in hydrological modeling and climate studies. Reprocessing of the complete snow-cover data record, from Collection 5 (C5) to Collection 6 (C6) and Collection 6.1 (C6.1), has provided improvements in the MODIS product suite. Suomi National Polar-orbiting Partnership (S-NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) Collection 1 (C1) snow-cover products at a 375 m spatial resolution have been available since 2011 and are currently being reprocessed for Collection 2 (C2). Both the MODIS C6.1 and the VIIRS C2 products will be available for download from the National Snow and Ice Data Center beginning in early 2020 with the complete time series available in 2020. To address the need for a cloud-reduced or cloud-free daily SCE product for both MODIS and VIIRS, a daily cloud-gap-filled (CGF) snow-cover algorithm was developed for MODIS C6.1 and VIIRS C2 processing. MOD10A1F (Terra) and MYD10A1F (Aqua) are daily, 500 m resolution CGF SCE map products from MODIS. VNP10A1F is the daily, 375 m resolution CGF SCE map product from VIIRS. These CGF products include quality-assurance data such as cloud-persistence statistics showing the age of the observation in each pixel. The objective of this paper is to introduce the new MODIS and VIIRS standard CGF daily SCE products and to provide a preliminary evaluation of uncertainties in the gap-filling methodology so that the products can be used as the basis for a moderate-resolution Earth science data record (ESDR) of SCE. Time series of the MODIS and VIIRS CGF products have been developed and evaluated at selected study sites in the US and southern Canada. Observed differences, although small, are largely attributed to cloud masking and differences in the time of day of image acquisition. A nearly 3-month time-series comparison of Terra MODIS and S-NPP VIIRS CGF snow-cover maps for a large study area covering all or parts of 11 states in the western US and part of southwestern Canada reveals excellent correspondence between the Terra MODIS and S-NPP VIIRS products, with a mean difference of 11 070 sqkm, which is ∼0.45 % of the study area. According to our preliminary validation of the Terra and Aqua MODIS CGF SCE products in the western US study area, we found higher accuracy of the Terra product compared with the Aqua product. The MODIS CGF SCE data record beginning in 2000 has been extended into the VIIRS era, which should last at least through the early 2030s.

Hall, Dorothy K.↗

NASA NTRSDate not supplied

A SATELLITE AGNOSTIC APPROACH TO QUANTIFYING HAIL DAMAGE SWATHS ACROSS THE CENTRAL UNITED STATES AND OTHER AGRICULTURAL REGIONS

Intense thunderstorms can bring damaging winds and large hail to agricultural regions during the prime growing season. In certain cases, large swaths of damage from the wind and hail are left behind and visible to satellite remote sensing instruments. Often times, Earth observing optical remote sensing from low to high spatial resolutions are able to view these damaged swaths. With the large number of moderate to high-resolution instruments in orbit, these damaged areas have potential to be viewed daily. However, during the prime growing season, clouds frequently block the viewing of the land surface by these optical instruments. Space-borne synthetic aperture radar (SAR) instruments allow for the viewing of the land surface in most weather conditions, but instead measure backscatter as opposed to optical sensors measuring reflected or emitted radiation. Additionally, the number of SAR instruments with free and open data lags behind the number of optical sensors. This paper describes the development of a methodology that attempts to characterize hail damaged swaths, through independent use of multiple optical and SAR platforms. This satellite-agnostic approach will focus statistical analysis by comparing undamaged areas to suspected damaged areas by using commonly derived indices from optical instruments and SAR backscatter from multiple polarizations.

Jordan R Bell↗

NASA NTRSDate not supplied

Evaluating Retrieval Algorithm Climate Stability: Estimating 3D Optical Thickness Bias Distributions by Cloud Type

Detecting climate trends on large spatiotemporal scales requires accurate, stable measurements and stable retrieval algorithms. We strive to estimate how time-variant retrieval algorithm biases may impact trend detection. Here we focus on the 3D cloud optical thickness (τc) bias, which is among the largest in passive cloud retrieval algorithms. If this bias is time dependent, a possibility with potential decadal changes in cloud morphology, it may obscure genuine trends in τc. Although previous studies have evaluated the cloud- and sun-view geometry-dependent 3D τc bias on small spatial scales, before our current study none have evaluated the stability of this well-known bias on climate-relevant large spatiotemporal scales. These studies must estimate large scale distributions of the 3D τc bias by cloud type and estimate how cloud type amount may change between two climate states. We employ a novel approach to estimate large scale distributions of 3D τc using a proxy of the bias that quantifies the departure of clouds from satisfying the 1D radiative transfer assumption used in passive τc retrievals. This existing globally-distributed proxy is an angular consistency metric that was developed using fused Moderate-Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) measurements. Calculating the 3D τc bias and the proxy, for known cloud fields enables us to establish statistical relationships between these two quantities, which can be used to calculate large-scale distributions of the 3D τc bias. This approach limits the number of 3D radiative transfer simulations required to only those needed to estimate a statistical relationship between the 3D τc bias for known cloud fields and a proxy of the bias. It is likely that future studies will be needed to evaluate retrieval algorithm bias stability for other geophysical variables as the community develops climate data records from satellite observations and their retrievals. This must be done in addition to monitoring and correcting measurement errors and uncertainties and understanding their impact on retrieved essential climate variables.

Yolanda Shea↗

NASA NTRS1988

Development of subminiature multi-sensor hot-wire probes

Limitations on the spatial resolution of multisensor hot wire probes have precluded accurate measurements of Reynolds stresses very near solid surfaces in wind tunnels and in many practical aerodynamic flows. The fabrication, calibration and qualification testing of very small single horizontal and X-array hot-wire probes which are intended to be used near solid boundaries in turbulent flows where length scales are particularly small, is described. Details of the sensor fabrication procedure are reported, along with information needed to successfully operate the probes. As compared with conventional probes, manufacture of the subminiature probes is more complex, requiring special equipment and careful handling. The subminiature probes tested were more fragile and shorter lived than conventional probes; they obeyed the same calibration laws but with slightly larger experimental uncertainty. In spite of these disadvantages, measurements of mean statistical quantities and spectra demonstrate the ability of the subminiature sensors to provide the measurements in the near wall region of turbulent boundary layers that are more accurate than conventional sized probes.

Westphal, Russell V.↗

NASA NTRS1996

Surface Deformation and Coherence Measurements of Kilauea Volcano, Hawaii, from SIR-C Radar Interferometry

The shuttle imaging radar C/X synthetic aperture radar (SIR-C/X-SAR) radar on board the space shuttle Endeavor imaged Kilauea Volcano, Hawaii, in April and October 1994 for the purpose of measuring active surface deformation by the methods of repeat-pass differential radar interferometry. Observations at 24 cm (L band) and 5.6 cm (C band) wavelengths were reduced to interferograms showing apparent surface deformation over the 6-month interval and over a succession of 1-day intervals in October. A statistically significant local phase signature in the 6-month interferogram is coincident with the Pu'u O'o lava vent. Interpreted as deformation, the signal implies centimeter-scale deflation in an area several kilometers wide surrounding the vent. Peak deflation is roughly 14 cm if the deformation is purely vertical, centered southward of the Pu'u O'o caldera. Delays in the radar signal phase induced by atmospheric refractivity anomalies introduce spurious apparent deformation signatures, at the level of 12 cm peak-to-peak in the radar line-of-sight direction. Though the phase observations are suggestive of the wide-area deformation measured by Global Positioning System (GPS) methods, the atmospheric effects are large enough to limit the interpretation of the result. It is difficult to characterize centimeter-scale deformations spatially distributed over tens of kilometers using differential interferometry without supporting simultaneous, spatially distributed measurements of reactivity along the radar line of sight. Studies of the interferometric correlation of images acquired at different times show that L band is far superior to C band in the vegetated areas, even when the observations are separated by only 1 day. These results imply longer wavelength instruments are more appropriate for studying surfaces by repeat-pass observations.

Rosen, P. A.↗

NASA NTRS1997

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time- dependent Navier-Stokes equations are integrated numerically using spectral, high-order finite-difference, and three-stage low-storage Runge-Kutta methods. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parameterize the subgrid-scale stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability lambda-vortex structures are in agreement with the experiments, and the streak-like-structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier-Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, without ad hoc adjustments.

Huai, Xiao-Li↗

NASA NTRS2017

Early Season Large-Area Winter Crop Mapping Using MODIS NDVI Data, Growing Degree Days Information and a Gaussian Mixture Model

Knowledge on geographical location and distribution of crops at global, national and regional scales is an extremely valuable source of information applications. Traditional approaches to crop mapping using remote sensing data rely heavily on reference or ground truth data in order to train/calibrate classification models. As a rule, such models are only applicable to a single vegetation season and should be recalibrated to be applicable for other seasons. This paper addresses the problem of early season large-area winter crop mapping using Moderate Resolution Imaging Spectroradiometer (MODIS) derived Normalized Difference Vegetation Index (NDVI) time-series and growing degree days (GDD) information derived from the Modern-Era Retrospective analysis for Research and Applications (MERRA-2) product. The model is based on the assumption that winter crops have developed biomass during early spring while other crops (spring and summer) have no biomass. As winter crop development is temporally and spatially non-uniform due to the presence of different agro-climatic zones, we use GDD to account for such discrepancies. A Gaussian mixture model (GMM) is applied to discriminate winter crops from other crops (spring and summer). The proposed method has the following advantages: low input data requirements, robustness, applicability to global scale application and can provide winter crop maps 1.5-2 months before harvest. The model is applied to two study regions, the State of Kansas in the US and Ukraine, and for multiple seasons (2001-2014). Validation using the US Department of Agriculture (USDA) Crop Data Layer (CDL) for Kansas and ground measurements for Ukraine shows that accuracies of greater than 90% can be achieved in mapping winter crops 1.5-2 months before harvest. Results also show good correspondence to official statistics with average coefficients of determination R(exp. 2) greater than 0.85.

mixture model↗