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,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS2009

Kepler Data Validation and Follow-up Programs

The approach that the Kepler Mission uses to remove false positive events and to validate the discoveries consists of two parts; data validation (DPI) and follow up observations program (FOP). DV consists of several methods of examining the data from the spacecraft observations. First, to rule out statistical fluctuations in the data, accept only signals that show 3 or more transits and that have a total signal-to-noise ratio that exceeds 7 sigma. Second. to identify small stellar companions to the target star, we check for secondary eclipses and determine if the transit characteristics are appropriate for a planetary companion. Third, check for background binaries that are in the target aperture. Here we measure the movement of the image centroid before, during, and after the transit. If the target is producing the signal, a dimming wi11 move the image centroid in a known direction and magnitude. If the signal comes from a nearby star, the amplitude and direction of the motion wi11 be different, This test is expected to rule out the hundreds of binary signals expected from background stars. The precision of the measurement depends on the stellar fluxes and positions but can be better than 0.01 pixel; i.e., 0.04". Those candidates that pass these tests are examined using ground-based telescopes and radial velocity spectrometers. First medium precision RV is used to rule out any remaining stellar companions. Then high spatial resolution imaging is used to check for nearby stars that are in the aperture- (The Kepler apertures depend on magnitude but are of order 36 sq are sec in area.) If no stars are present that quid generate the observed signal, then the candidate goes to a large telescope such as Keck, HET, or Wi1lilam Herschel for high precision observations to get the planet mass or an upper limit to it, if there are some stars in the aperture, then the photometric observations are employed to look for the transit by cane of the confounding stars. If none are seen, then the candidate is sent to a large telescope such as Keck

William, Borucki J.↗

NASA NTRS2002

A geostatistical analysis of small-scale spatial variability in bacterial abundance and community structure in salt marsh creek bank sediments

Small-scale variations in bacterial abundance and community structure were examined in salt marsh sediments from Virginia's eastern shore. Samples were collected at 5 cm intervals (horizontally) along a 50 cm elevation gradient, over a 215 cm horizontal transect. For each sample, bacterial abundance was determined using acridine orange direct counts and community structure was analyzed using randomly amplified polymorphic DNA fingerprinting of whole-community DNA extracts. A geostatistical analysis was used to determine the degree of spatial autocorrelation among the samples, for each variable and each direction (horizontal and vertical). The proportion of variance in bacterial abundance that could be accounted for by the spatial model was quite high (vertical: 60%, horizontal: 73%); significant autocorrelation was found among samples separated by 25 cm in the vertical direction and up to 115 cm horizontally. In contrast, most of the variability in community structure was not accounted for by simply considering the spatial separation of samples (vertical: 11%, horizontal: 22%), and must reflect variability from other parameters (e.g., variation at other spatial scales, experimental error, or environmental heterogeneity). Microbial community patch size based upon overall similarity in community structure varied between 17 cm (vertical) and 35 cm (horizontal). Overall, variability due to horizontal position (distance from the creek bank) was much smaller than that due to vertical position (elevation) for both community properties assayed. This suggests that processes more correlated with elevation (e.g., drainage and redox potential) vary at a smaller scale (therefore producing smaller patch sizes) than processes controlled by distance from the creek bank. c2002 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved.

Non-NASA Center↗

NASA NTRS1999

A Global Data Assimilation System for Atmospheric Aerosol

We will give an overview of an aerosol data assimilation system which combines advances in remote sensing of atmospheric aerosols, aerosol modeling and data assimilation methodology to produce high spatial and temporal resolution 3D aerosol fields. Initially, the Goddard Aerosol Assimilation System (GAAS) will assimilate TOMS, AVHRR and AERONET observations; later we will include MODIS and MISR. This data assimilation capability will allows us to integrate complementing aerosol observations from these platforms, enabling the development of an assimilated aerosol climatology as well as a global aerosol forecasting system in support of field campaigns. Furthermore, this system provides an interactive retrieval framework for each aerosol observing satellites, in particular TOMS and AVHRR. The Goddard Aerosol Assimilation System (GAAS) takes advantage of recent advances in constituent data assimilation at DAO, including flow dependent parameterizations of error covariances and the proper consideration of model bias. For its prognostic transport model, GAAS will utilize the Goddard Ozone, Chemistry, Aerosol, Radiation and Transport (GOCART) model developed at NASA/GSFC Codes 916 and 910.3. GOCART includes the Lin-Rood flux-form, semi-Langrangian transport model with parameterized aerosol chemistry and physical processes for absorbing (dust and black carbon) and non-absorbing aerosols (sulfate and organic carbon). Observations and model fields are combined using a constituent version of DAO's Physical-space Statistical Analysis System (PSAS), including its adaptive quality control system. In this talk we describe the main components of this assimilation system and present preliminary results obtained by assimilating TOMS data.

daSilva, Arlindo↗

NASA NTRS2017

Cloud and Sun-Glint Statistics Derived from GOES and MODIS Observations Over the Intra-Americas Sea for GEO-CAPE Mission Planning

Knowledge of cloud cover, frequency, and duration is not only important to study cloud dynamics, but also critical in determining when and where to take ocean measurements from geostationary orbits such as the Geostationary Coastal and Air Pollution Events (GEO-CAPE) mission due to the challenges in achieving complete hemispheric coverage of coastal oceans, estuaries, and inland waters at hourly frequency. Using GOES hourly measurements at 4 km nadir resolution between 2006 and 2011, the number of cloud-free hourly observations per day (N(sub cf)) for solar zenith angle Theta(sub 0) less than 80 degrees was estimated for each 0.1 degree location of the Intra-Americas Sea. The number of Sun-glint-affected hourly observations per day [Ns(sub sg)] was also calculated based on the planned GEO-CAPE observation geometry and realistic wind speed. High-latitude and equatorial oceans showed the lowest N(sub cf) (less than 2.4) in all climatological months, and highest N(sub cf) was observed in the Gulf of Mexico (GoM) and Caribbean (greater than 4.5). Different regions showed differences in seasonality of cloud-free conditions and also showed differences in the hour of a day at which the satellite observations would have the maximal cloud-free and glint-free probability (Temperature maximum). Cloud cover from Moderate Resolution Imaging Spectroradiometer (MODIS) 1 km measurements are greater than 10 degrees higher than those from the MODIS 250 m measurements, supporting ocean color missions at subkilometer resolutions to enhance both spatial coverage and temporal frequency. These findings provide valuable information for GEO-CAPE mission planning to maximize its science value through minimizing the impacts of clouds and Sun glint.

Feng, Lian↗

NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗

NASA NTRSDate not supplied

Exoplanet Yield Sensitivity for the Hybrid Lyot Coronagraph from end-to-end modeling for LUVOIR-A

NASA's future direct imaging flagship mission will require telescope mirror segmentation to survey a statistically significant population of nearby stellar systems by increasing aperture size beyond that achieved by traditional space-borne monolithic apertures thus far. Segmentation is a path for the size expansion, but it entails a challenge due to the introduction of temporal pupil discontinuities in terms of both phase and amplitude. Modeling of the integrated performance of the coronagraph with the expected dynamical profile of the telescope is an important step toward realistic estimates of the mission's science yields. We present results of a study using the hybrid Lyot coronagraph with a 15-m on-axis, segmented optical telescope assembly (OTA), akin to that of the Large UV / Optical /Infrared (LUVOIR) mission. We have developed a high-fidelity, end-to-end model as part of the Segmented Coronagraph Design and Analysis (SCDA) study that includes diffractive optical propagations at the 1e-10 contrast level across the optical train including the OTA and intermediate re-imaging optics inside the coronagraph instrument. We perform a sensitivity analysis to assess the contrast sensitivity to wavefront perturbations for a set of key telescope terms including segment-level low-order and mid-spatial frequency aberrations and evaluate expected exoplanet yields obtained in the presence of these sets of aberrations. Industry-provided segment-level wavefront aberration datasets involving the following mechanisms are input to the end-to-end model: (1) thermal gradients and accelerations, and (2) telescope dynamics including line-of-sight and segment jitter. Dynamically averaged PSFs in the presence of aberrations are fed into AYO (Altruistic Yield Optimization), a mission design reference simulator, that then assesses the scientific yield sensitivities.

high-contrast imaging↗

NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗

NASA NTRS2018

A Comparative Analysis of Micrometeorological Determinants of Evapotranspiration Rates Within a Heterogeneous Urban Environment

Variability in micrometeorological conditions and their influence on estimated reference evapotranspiration (RET) rates were evaluated across a heterogeneous urban environment. Micrometeorological data sets (incoming solar radiation, air temperature, relative humidity and wind speed) were collected over a one-year period at six weather stations in New York City, NY (USA). Weather stations are located at four new urban green space monitoring sites and two airports. Reference evapotranspiration (RET) rates were estimated from the micrometeorological data sets for a short reference surface at a daily time-step using the ASCE Standardized Reference Evapotranspiration Equation, a Penman-Monteith based combination equation. Nonparametric comparative statistical analyses (Kruskal-Wallis) revealed statistically significant differences (at significance level α = 0.05) in micrometeorological conditions and estimated RET rates between the six sites. On a cumulative annual basis, estimated RET varied by up to 40 percent between the sites. A new technique for adjusting weather data collected at one location (e.g. regional airports) for use at another location (e.g. interior engineered urban green spaces) was evaluated. The study highlights the importance, for accurate estimation of ET, of onsite micrometeorological data sets, but concludes that additional research is needed to more thoroughly characterize micrometeorological variability across heterogeneous urban environments, and also to evaluate the influence of non-meteorological determinants, e.g. vegetation type, soil/media type, media moisture conditions and anthropogenic heat fluxes, on urban ET.

Urban environment↗

NASA NTRS1994

Progress in modeling atmospheric propagation of sonic booms

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by humidity and turbulence. A realistic atmosphere is invariably somewhat turbulent, and may be characterized by an ambient fluid velocity v and sound speed c that vary from point to point. The absolute humidity will also vary from point to point, although possibly not as irregularly. What is ideally desired is a relatively simple scheme for predicting the probable spreads in key sonic boom signature parameters. Such parameters could be peak amplitudes, rise times, or gross quantities obtainable by signal processing that correlate well with annoyance or damage potential. The practical desire for the prediction scheme is that it require a relatively small amount of knowledge, possibly of a statistical nature, concerning the atmosphere along, the propagation path from the aircraft to the ground. The impact of such a scheme, if developed, implemented, and verified, would be that it would give the persons who make planning decisions a tool for assessing the magnitude of environmental problems that might result from any given overflight or sequence of overflights. The technical approach that has been followed by the author and some of his colleagues is to formulate a hierarchy of simple approximate models based on fundamental physical principles and then to test these models against existing data. For propagation of sonic booms and of other types of acoustic pulses in nonturbulent model atmospheres, there exists a basic overall theoretical model that has evolved as an outgrowth of geometrical acoustics. This theoretical model depicts the sound as propagating within ray tubes in a manner analogous to sound in a waveguide of slowly varying cross-section. Propagation along the ray tube is quasi-one-dimensional, and a wave equation for unidirectional wave propagation is used. A nonlinear term is added to this equation to account for nonlinear steepening, and the formulation has been carried through to allow for spatially varying sound speed, ambient density, and ambient wind velocities. The model intrinsically neglects diffraction, so it cannot take into account what has previously been mentioned in the literature as possibly important mechanisms for turbulence-related distortion. The model as originally developed could predict an idealized N-waveform which often agrees with data in terms of peak amplitude and overall positive phase duration. It is possible, moreover, to develop simple methods based on the physics of relaxation processes for incorporating molecular relaxation into the quasi-one-dimensional model of nonlinear propagation along ray tubes.

Pierce, Allan D.↗

NASA NTRS2016

The Viability of Trajectory Analysis for Diagnosing Dynamical and Chemical Influences on Ozone Concentrations in the UTLS

The viability of trajectory analysis for diagnosing the interplay between chemistry and dynamics is investigated by comparing ozone mixing ratios modelled using air-parcel pathways to values observed along flight tracks during ATTREX (Airborne Tropical TRopopause EXperiment). Trajectories are initiated at the locations of ozone observations and tracked backward in time to their sources at termini of backward trajectories. The modelled values of ozone utilize 3-dimensional analysis fields from WACCM (Whole Atmosphere Community Climate Model) (a chemical-climate model with dynamical fields nudged towards MERRA (Modern-Era Retrospective Analysis and Research Applications) reanalysis) and ERA-interim (product of ECMWF - the European Centre for Medium-Range Weather Forecasts) to determine source mixing ratios with chemical production and loss terms derived from the ozone chemistry used in WACCM. A statistical base of modelled ozone is constructed with 6 trajectory platforms (adiabatic, diabatic, and kinematic forced by ERA-interim and MERRA), two chemical models (WACCM chemistry and no chemistry), and 4 trajectory lengths (5, 10, 20, and 30 days). Linear regression is employed to separate systematic errors from random errors and to characterize the impact of source mixing ratios, path length, vertical motion, and chemistry on modelled ozone errors. Errors in the analysis ozone fields are large, if not dominant, contributors to model error. Random errors are particularly large for point-by-point comparisons, however averaging over 800 km (75 minutes) flight segments substantially reduces random error and exposes systematic errors. Of the two analysis ozone data sets, WACCM, which incorporates detailed chemistry, provides the smaller systematic errors while ERA-interim, which has crude chemistry but assimilates observational data, has the smaller random errors. Of the different trajectory platforms, adiabatic calculations produce the smaller random errors (irrespective of the use of chemistry) but both vertical motion and chemistry are required to optimally reduce systematic errors. These results suggest that meaningful analysis of dynamical and chemical interactions that control ozone mixing ratios are viable on spatial scales larger than a few reanalysis grid spaces, that errors in the analyzed ozone data sets are large but not prohibitively so, and that vertical velocities and heating rates from reanalysis data, while problematic, contain useful information [on the ozone concentrations in the UTLS (Upper Troposphere/Lower Stratosphere)].

ozone↗

NASA NTRS2000

A Common Calibration Source Framework for Fully-Polarimetric and Interferometric Radiometers

Two types of microwave radiometry--synthetic thinned array radiometry (STAR) and fully-polarimetric (FP) radiometry--have received increasing attention during the last several years. STAR radiometers offer a technological solution to achieving high spatial resolution imaging from orbit without requiring a filled aperture or a moving antenna, and FP radiometers measure extra polarization state information upon which entirely new or more robust geophysical retrieval algorithms can be based. Radiometer configurations used for both STAR and FP instruments share one fundamental feature that distinguishes them from more 'standard' radiometers, namely, they measure correlations between pairs of microwave signals. The calibration requirements for correlation radiometers are broader than those for standard radiometers. Quantities of interest include total powers, complex correlation coefficients, various offsets, and possible nonlinearities. A candidate for an ideal calibration source would be one that injects test signals with precisely controllable correlation coefficients and absolute powers simultaneously into a pair of receivers, permitting all of these calibration quantities to be measured. The complex nature of correlation radiometer calibration, coupled with certain inherent similarities between STAR and FP instruments, suggests significant leverage in addressing both problems together. Recognizing this, a project was recently begun at NASA Goddard Space Flight Center to develop a compact low-power subsystem for spaceflight STAR or FP receiver calibration. We present a common theoretical framework for the design of signals for a controlled correlation calibration source. A statistical model is described, along with temporal and spectral constraints on such signals. Finally, a method for realizing these signals is demonstrated using a Matlab-based implementation.

Kim, Edward J.↗

NASA NTRS2009

Winter QPF Sensitivities to Snow Parameterizations and Comparisons to NASA CloudSat Observations

Steady increases in computing power have allowed for numerical weather prediction models to be initialized and run at high spatial resolution, permitting a transition from larger scale parameterizations of the effects of clouds and precipitation to the simulation of specific microphysical processes and hydrometeor size distributions. Although still relatively coarse in comparison to true cloud resolving models, these high resolution forecasts (on the order of 4 km or less) have demonstrated value in the prediction of severe storm mode and evolution and are being explored for use in winter weather events . Several single-moment bulk water microphysics schemes are available within the latest release of the Weather Research and Forecast (WRF) model suite, including the NASA Goddard Cumulus Ensemble, which incorporate some assumptions in the size distribution of a small number of hydrometeor classes in order to predict their evolution, advection and precipitation within the forecast domain. Although many of these schemes produce similar forecasts of events on the synoptic scale, there are often significant details regarding precipitation and cloud cover, as well as the distribution of water mass among the constituent hydrometeor classes. Unfortunately, validating data for cloud resolving model simulations are sparse. Field campaigns require in-cloud measurements of hydrometeors from aircraft in coordination with extensive and coincident ground based measurements. Radar remote sensing is utilized to detect the spatial coverage and structure of precipitation. Here, two radar systems characterize the structure of winter precipitation for comparison to equivalent features within a forecast model: a 3 GHz, Weather Surveillance Radar-1988 Doppler (WSR-88D) based in Omaha, Nebraska, and the 94 GHz NASA CloudSat Cloud Profiling Radar, a spaceborne instrument and member of the afternoon or "A-Train" of polar orbiting satellites tasked with cataloguing global cloud characteristics. Each system provides a unique perspective. The WSR-88D operates in a surveillance mode, sampling cloud volumes of Rayleigh scatterers where reflectivity is proportional to the sixth moment of the size distribution of equivalent spheres. The CloudSat radar provides enhanced sensitivity to smaller cloud ice crystals aloft, as well as consistent vertical profiles along each orbit. However, CloudSat reflectivity signatures are complicated somewhat by resonant Mie scattering effects and significant attenuation in the presence of cloud or rain water. Here, both radar systems are applied to a case of light to moderate snowfall within the warm frontal zone of a cold season, synoptic scale storm. Radars allow for an evaluation of the accuracy of a single-moment scheme in replicating precipitation structures, based on the bulk statistical properties of precipitation as suggested by reflectivity signatures.

Molthan, Andrew↗

NASA NTRS2014

Global-Scale Comparison of Passive (SMOS) and Active (ASCAT) Satellite Based Microwave Soil Moisture Retrievals with Soil Moisture Simulations (MERRA-Land)

Global surface soil moisture (SSM) datasets are being produced based on active and passive microwave satellite observations and simulations from land surface models (LSM). This study investigates the consistency of two global satellite-based SSM datasets based on microwave remote sensing observations from the passive Soil Moisture and Ocean Salinity (SMOS;SMOSL3 version 2.5) and the active Advanced Scatterometer (ASCAT; version TUWien- WARP 5.5) with respect to LSM SSM from the MERRA-Land data product. The relationship between the global-scale SSM products was studied during the 2010-2012 period using (1) a time series statistics (considering both original SSM data and anomalies), (2) a space-time analysis using Hovmoller diagrams, and (3) a triple collocation error model. The SMOSL3 and ASCAT retrievals are consistent with the temporal dynamics of modeled SSM (correlation R (is) greater than 0.70 for original SSM) in the transition zones between wet and dry climates, including the Sahel, the Indian subcontinent, the Great Plains of North America, eastern Australia, and southeastern Brazil. Over relatively dense vegetation covers, a better consistency with MERRA-Land was obtained with ASCAT than with SMOSL3. However, it was found that ASCAT retrievals exhibit negative correlation versus MERRA-Land in some arid regions (e.g., the Sahara and the Arabian Peninsula). In terms of anomalies, SMOSL3 better captures the short term SSM variability of the reference dataset (MERRA-Land) than ASCAT over regions with limited radio frequency interference (RFI) effects (e.g., North America, South America, and Australia). The seasonal and latitudinal variations of SSM are relatively similar for the three products, although the MERRALand SSM values are generally higher and their seasonal amplitude is much lower than for SMOSL3 and ASCAT. Both SMOSL3 and ASCAT have relatively comparable triple collocation errors with similar spatial error patterns: (i) lowest errors in arid regions (e.g., Sahara and Arabian Peninsula), due to the very low natural variability of soil moisture in these areas, and Central America, and (ii) highest errors over most of the vegetated regions (e.g., northern Australia, India, central Asia, and South America). However, the ASCAT SSM product is prone to larger random errors in some regions (e.g., north-western Africa, Iran, and southern South Africa). Vegetation density was found to be a key factor to interpret the consistency with MERRA-Land between the two remotely sensed products (SMOSL3 and ASCAT) which provides complementary information on SSM. This study shows that both SMOS and ASCAT have thus a potential for data fusion into long-term data records.

ASCAT↗

NASA NTRS2012

Seismic Rate Changes Associated with Seasonal, Annual, and Decadal Changes in the Cryosphere

Near the Bering Glacier Global Fiducial site in southern Alaska large cryospheric fluctuations occur in a region of upper crustal faulting and folding associated with collision and accretion of the Yakutat terrane. In this study we report constraints on seasonal, annual and decadal cryospheric changes estimated over the last decade from field, aircraft and satellite measurements, and we evaluate the influence of cryospheric changes on the background seismic rate. Multi-year images from the Bering Glacier global fiducial site are available since mid-2003 to constrain changes in extent of the Bering Glacier and to discern feature changes in the glacial surface. Starting around the same time, satellite gravimetric measurements from the Gravity Recovery and Climate experiment (GRACE) commenced. Large spatial-scale mass change calculated from the GRACE 1deg x 1deg mascon solution of Luthcke et al. [2012] indicate a general trend of annual ice mass loss for southern Alaska but with large, variable seasonal mass fluctuations. Since 2007, the station position of a continuous GPS site near Cape Yakataga (Alaska EarthScope PBO site, AB35) has been available as well. In addition to changes in the geodetic position due to tectonic motion, this GPS station shows large seasonal excursions in the detrended vertical and horizontal position components consistent with snow loading in the fall and winter and melt onset/mass decrease in the spring/summer. To better understand the timing of processes responsible for the onset of cryospheric mass loss documented in the GRACE data, we examined changes in the snow cover extent and the onset of melt in the spring. We calculated the surface displacements of the solid Earth and theoretical earthquake failure criteria associated with these annual and seasonal ice and snow changes using layered elastic half-space. Additionally, we compared the seismic rate (M>1.8) from a reference background time period against other time periods with variable ice or tectonic change characteristics to test the significance of seismic rate changes. Our earlier results suggest statistically significant changes in the background seismic rate associated with large seasonal mass changes. INDEX

Sauber-Rosenberg, Jeanne↗

NASA NTRS2023

Uncertainty Analysis of Slug Calorimeters in the HyMETS Arc-Jet Facility

The objective of this work is to perform an uncertainty analysis of the deduced stagnation heat flux environment on a slug calorimeter for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System arc-jet facility located at NASA Langley Research Center. Analytical solutions are developed for boundary-value problems on the slug element accounting for non-ideal effects, including spatial variation in the slug heat flux, multi-dimensional thermal conduction, and back-face losses, which departs from the state-of-the-art method derived from the American Society of Testing and Materials. Boundary-value problem definitions are informed by preliminary finite element thermal analysis of the slug calorimeter assembly (including both slug and housing) and just the slug element. The analytical solutions are presented in a general sense and in a truncated form from error analysis. Results are shown in optimizing and validating the analytical models against available slug back-face thermal data. The optimization results indicate that the appropriate epistemic uncertainty of the deduced stagnation heat flux on the slug calorimeter is at most±2.5% for both a high-and low-enthalpy test condition. In addition, a numerical approach is used to determine the aleatory (probabilistic) uncertainty component in the slug stagnation heat flux by applying a marching least-squares slope routine through the steady-state portion of the slug back-face thermal response. Results indicate a compromise between the number of samples and the filter frequency of slug back-face thermal data points when evaluating the standard deviation of the deduced stagnation heat flux statistics. When combining the mixed uncertainty, both aleatory and epistemic, the interval of uncertainty in the deduced stagnation heat flux is determined to be up to ±4%, which is at least a 60% reduction from the standard uncertainty used in the state-of-the-art method.

uncertainty↗

NASA NTRS2014

ENSO Effect on East Asian Tropical Cyclone Landfall via Changes in Tracks and Genesis in a Statistical Model

Improvements on a statistical tropical cyclone (TC) track model in the western North Pacific Ocean are described. The goal of the model is to study the effect of El Nino-Southern Oscillation (ENSO) on East Asian TC landfall. The model is based on the International Best-Track Archive for Climate Stewardship (IBTrACS) database of TC observations for 1945-2007 and employs local regression of TC formation rates and track increments on the Nino-3.4 index and seasonally varying climate parameters. The main improvements are the inclusion of ENSO dependence in the track propagation and accounting for seasonality in both genesis and tracks. A comparison of simulations of the 1945-2007 period with observations concludes that the model updates improve the skill of this model in simulating TCs. Changes in TC genesis and tracks are analyzed separately and cumulatively in simulations of stationary extreme ENSO states. ENSO effects on regional (100-km scale) landfall are attributed to changes in genesis and tracks. The effect of ENSO on genesis is predominantly a shift in genesis location from the southeast in El Nino years to the northwest in La Nina years, resulting in higher landfall rates for the East Asian coast during La Nina. The effect of ENSO on track propagation varies seasonally and spatially. In the peak activity season (July-October), there are significant changes in mean tracks with ENSO. Landfall-rate changes from genesis- and track-ENSO effects in the Philippines cancel out, while coastal segments of Vietnam, China, the Korean Peninsula, and Japan show enhanced La Nina-year increases.

simulation↗

NASA NTRS1993

Visualization techniques to aid in the analysis of multi-spectral astrophysical data sets

The goal of this project was to support the scientific analysis of multi-spectral astrophysical data by means of scientific visualization. Scientific visualization offers its greatest value if it is not used as a method separate or alternative to other data analysis methods but rather in addition to these methods. Together with quantitative analysis of data, such as offered by statistical analysis, image or signal processing, visualization attempts to explore all information inherent in astrophysical data in the most effective way. Data visualization is one aspect of data analysis. Our taxonomy as developed in Section 2 includes identification and access to existing information, preprocessing and quantitative analysis of data, visual representation and the user interface as major components to the software environment of astrophysical data analysis. In pursuing our goal to provide methods and tools for scientific visualization of multi-spectral astrophysical data, we therefore looked at scientific data analysis as one whole process, adding visualization tools to an already existing environment and integrating the various components that define a scientific data analysis environment. As long as the software development process of each component is separate from all other components, users of data analysis software are constantly interrupted in their scientific work in order to convert from one data format to another, or to move from one storage medium to another, or to switch from one user interface to another. We also took an in-depth look at scientific visualization and its underlying concepts, current visualization systems, their contributions, and their shortcomings. The role of data visualization is to stimulate mental processes different from quantitative data analysis, such as the perception of spatial relationships or the discovery of patterns or anomalies while browsing through large data sets. Visualization often leads to an intuitive understanding of the meaning of data values and their relationships by sacrificing accuracy in interpreting the data values. In order to be accurate in the interpretation, data values need to be measured, computed on, and compared to theoretical or empirical models (quantitative analysis). If visualization software hampers quantitative analysis (which happens with some commercial visualization products), its use is greatly diminished for astrophysical data analysis. The software system STAR (Scientific Toolkit for Astrophysical Research) was developed as a prototype during the course of the project to better understand the pragmatic concerns raised in the project. STAR led to a better understanding on the importance of collaboration between astrophysicists and computer scientists.

Brugel, Edward W.↗

NASA NTRS1993

Visualization techniques to aid in the analysis of multispectral astrophysical data sets

The goal of this project was to support the scientific analysis of multi-spectral astrophysical data by means of scientific visualization. Scientific visualization offers its greatest value if it is not used as a method separate or alternative to other data analysis methods but rather in addition to these methods. Together with quantitative analysis of data, such as offered by statistical analysis, image or signal processing, visualization attempts to explore all information inherent in astrophysical data in the most effective way. Data visualization is one aspect of data analysis. Our taxonomy as developed in Section 2 includes identification and access to existing information, preprocessing and quantitative analysis of data, visual representation and the user interface as major components to the software environment of astrophysical data analysis. In pursuing our goal to provide methods and tools for scientific visualization of multi-spectral astrophysical data, we therefore looked at scientific data analysis as one whole process, adding visualization tools to an already existing environment and integrating the various components that define a scientific data analysis environment. As long as the software development process of each component is separate from all other components, users of data analysis software are constantly interrupted in their scientific work in order to convert from one data format to another, or to move from one storage medium to another, or to switch from one user interface to another. We also took an in-depth look at scientific visualization and its underlying concepts, current visualization systems, their contributions and their shortcomings. The role of data visualization is to stimulate mental processes different from quantitative data analysis, such as the perception of spatial relationships or the discovery of patterns or anomalies while browsing through large data sets. Visualization often leads to an intuitive understanding of the meaning of data values and their relationships by sacrificing accuracy in interpreting the data values. In order to be accurate in the interpretation, data values need to be measured, computed on, and compared to theoretical or empirical models (quantitative analysis). If visualization software hampers quantitative analysis (which happens with some commercial visualization products), its use is greatly diminished for astrophysical data analysis. The software system STAR (Scientific Toolkit for Astrophysical Research) was developed as a prototype during the course of the project to better understand the pragmatic concerns raised in the project. STAR led to a better understanding on the importance of collaboration between astrophysicists and computer scientists. Twenty-one examples of the use of visualization for astrophysical data are included with this report. Sixteen publications related to efforts performed during or initiated through work on this project are listed at the end of this report.

Brugel, E. W.↗