Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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 559 records · Page 31

Aerosol Properties From Combined Oxygen A Band Radiances and Lidar

We have developed a new aerosol retrieval technique based on combing high-resolution A band spectra with lidar profiles. Our goal is the development of a technique to retrieve aerosol absorption, one of the critical parameters affecting the global radiation budget and one which is currently poorly constrained by satellite measurements. Our approach relies on two key factors: 1) the use of high spectral resolution (17,000:1) measurements which resolve the A-band line structure, and 2) the use of co-located lidar profile measurements to constrain the vertical distribution of scatterers in the forward model. The algorithm has been developed to be applied to observations from the CALIPSO and OCO-2 satellites, flying in formation as part of the A-train constellation. We describe the approach and present simulated retrievals to illustrate performance potential.

Winker, Dave↗

Have GRACE Satellites Overestimated Groundwater Depletion in the Northwest India Aquifer?

The Northwest India Aquifer (NWIA) has been shown to have the highest groundwater depletion (GWD) rate globally, threatening crop production and sustainability of groundwater resources. Gravity Recovery and Climate Experiment (GRACE) satellites have been emerging as a powerful tool to evaluate GWD with ancillary data. Accurate GWD estimation is, however, challenging because of uncertainties in GRACE data processing. We evaluated GWD rates over the NWIA using a variety of approaches, including newly developed constrained forward modeling resulting in a GWD rate of 3.1 plus or minus 0.1 centimeters per acre (or 14 plus or minus 0.4 cubic kilometers per acre) for Jan 2005-Dec 2010, consistent with the GWD rate (2.8 centimeters per acre or 12.3 cubic kilometers per acre) from groundwater-level monitoring data. Published studies (e.g., 4 plus or minus 1 centimeter per acre or 18 plus or minus 4.4 cubic kilometers per acre) may overestimate GWD over this region. This study highlights uncertainties in GWD estimates and the importance of incorporating a priori information to refine spatial patterns of GRACE signals that could be more useful in groundwater resource management and need to be paid more attention in future studies.

Long, Di↗

Soil Moisture Active Passive (SMAP) L-Band Microwave Radiometer Post-Launch Calibration

The SMAP microwave radiometer is a fully-polarimetric L-band radiometer flown on the SMAP satellite in a 6 AM/ 6 PM sun-synchronous orbit at 685 km altitude. Since April, 2015, the radiometer is under calibration and validation to assess the quality of the radiometer L1B data product. Calibration methods including the SMAP L1B TA2TB (from Antenna Temperature (TA) to the Earth's surface Brightness Temperature (TB)) algorithm and TA forward models are outlined, and validation approaches to calibration stability/quality are described in this paper including future work. Results show that the current radiometer L1B data satisfies its requirements.

simulation↗

Soil Moisture ActivePassive (SMAP) L-Band Microwave Radiometer Post-Launch Calibration

The SMAP microwave radiometer is a fully-polarimetric L-band radiometer flown on the SMAP satellite in a 6 AM/ 6 PM sun-synchronous orbit at 685 km altitude. Since April, 2015, the radiometer is under calibration and validation to assess the quality of the radiometer L1B data product. Calibration methods including the SMAP L1B TA2TB (from Antenna Temperature (TA) to the Earth’s surface Brightness Temperature (TB)) algorithm and TA forward models are outlined, and validation approaches to calibration stability/quality are described in this paper including future work. Results show that the current radiometer L1B data satisfies its requirements.

radiometers↗

A Comparative Analysis of the Magnetic Field Signals over Impact Structures on the Earth, Mars and the Moon

An improved description of magnetic fields of terrestrial bodies has been obtained from recent space missions, leading to a better characterization of the internal fields including those of crustal origin. One of the striking differences in their crustal magnetic field is the signature of large impact craters. A comparative analysis of the magnetic characteristics of these structures can shed light on the history of their respective planetary-scale magnetic dynamos. This has motivated us to identify impact craters and basins, first by their quasi-circular features from the most recent and detailed topographic maps and then from available global magnetic field maps. We have examined the magnetic field observed above 27 complex craters on the Earth, 34 impact basins on Mars and 37 impact basins on the Moon. For the first time, systematic trends in the amplitude and frequency of the magnetic patterns, inside and outside of these structures are observed for all three bodies. The demagnetization effects due to the impact shock wave and excavation processes have been evaluated applying the Equivalent Source Dipole forward modeling approach. The main characteristics of the selected impact craters are shown. The trends in their magnetic signatures are indicated, which are related to the presence or absence of a planetary-scale dynamo at the time of their formation and to impact processes. The low magnetic field intensity at center can be accepted as the prime characteristic of a hypervelocity impact and strongly associated with the mechanics of impact crater formation. In the presence of an active internal field, the process of demagnetization due to the shock impact is associated with post-impact remagnetization processes, generating a more complex magnetic signature.

Earth↗

Mass Evolution of Mediterranean, Black, Red, and Caspian Seas from GRACE and Altimetry: Accuracy Assessment and Solution Calibration

We present new measurements of mass evolution for the Mediterranean, Black, Red, and Caspian Seas as determined by the NASA Goddard Space Flight Center (GSFC) GRACE time-variable global gravity mascon solutions. These new solutions are compared to sea surface altimetry measurements of sea level anomalies with steric corrections applied. To assess their accuracy, the GRACE and altimetry-derived solutions are applied to the set of forward models used by GSFC for processing the GRACE Level-1B datasets, with the resulting inter-satellite range acceleration residuals providing a useful metric for analyzing solution quality.

Loomis, B. D.↗

Spatial and Temporal Antarctic Ice Sheet Mass Trends, Glacio-Isostatic Adjustment, and Surface Processes from a Joint Inversion of Satellite Altimeter, Gravity, and GPS Data

We present spatiotemporal mass balance trends for the Antarctic Ice Sheet from a statistical inversion of satellite altimetry, gravimetry, and elastic-corrected GPS data for the period 2003-2013. Our method simultaneously determines annual trends in ice dynamics, surface mass balance anomalies, and a time-invariant solution for glacio-isostatic adjustment while remaining largely independent of forward models. We establish that over the period 2003-2013, Antarctica has been losing mass at a rateof -84 +/- 22 Gt per yr, with a sustained negative mean trend of dynamic imbalance of -111 +/- 13 Gt per yr. West Antarctica is the largest contributor with -112 +/- 10 Gt per yr, mainly triggered by high thinning rates of glaciers draining into the Amundsen Sea Embayment. The Antarctic Peninsula has experienced a dramatic increase in mass loss in the last decade, with a mean rate of -28 +/- 7 Gt per yr and significantly higher values for the most recent years following the destabilization of the Southern Antarctic Peninsula around 2010. The total mass loss is partly compensated by a significant mass gain of 56 +/- 18 Gt per yr in East Antarctica due to a positive trend of surface mass balance anomalies.

Martin-Espanol, Alba↗

Survey of Galileo Plasma Observations in Jupiter's Plasma Sheet

The plasma science (PLS) Instrument on the Galileo spacecraft (orbiting Jupiter from December 1995 to September 2003) measured properties of the ions that were trapped in the magnetic field. The PLS data provide a survey of the plasma properties between approx. 5 and 30 Jupiter radii [R(sub J)] in the equatorial region. We present plasma properties derived via two analysis methods: numerical moments and forward modeling. We find that the density decreases with radial distance by nearly 5 orders of magnitude from approx. 2 to 3000 cm(exp.-3) at 6R(sub j) to approx. 0.05cm(sub -3) at 30 R(sub j). The density profile did not show major changes from orbit to orbit, suggesting that the plasma production and transport remained constant within about a factor of 2. The radial profile of ion temperature increased with distance which implied that contrary to the concept of adiabatic cooling on expansion, the plasma heats up as it expands out from Io's orbit (where TI is approx.60-80 eV) at approx. 6R(sub j) to a few keV at 30R(sub j).There does not seem to be a long-term, systematic variation in ion temperature with either local time or longitude. This latter finding differs from earlier analysis of Galileo PLS data from a selection of orbits. Further examination of all data from all Galileo orbits suggests that System Ill variations are transitory on timescales of weeks, consistent with the modeling of Cassini Ultraviolet Imaging Spectrograph observations. The plasma flow is dominated by azimuthal flow that is between 80% and 100% of corotation out to 25 R(sub j).

Bagenal, Fran↗

Challenges in Discerning Atmospheric Composition in Directly Imaged Planets

One of the justifications motivating efforts to detect and characterize young extrasolar giant planets has been to measure atmospheric composition for comparison with that of the primary star. If the enhancement of heavy elements in the atmospheres of extrasolar giant planets, like it is for their solar system analogs, is inversely proportional to mass, then it is likely that these worlds formed by core accretion. However in practice it has been very difficult to constrain metallicity because of the complex effect of clouds. Cloud opacity varies both vertically and, in some cases, horizontally through the atmosphere. Particle size and composition, both of which impact opacity, are difficult challenges both for forward modeling and retrieval studies. In my presentation I will discuss systematic efforts to improve cloud studies to enable more reliable determinations of atmospheric composition. These efforts are relevant both to discerning composition of directly imaged young planets from ground based telescopes and future space based missions, such as WFIRST and LUVOIR.

young extrasolar giant planets↗

Comparison of the Coronal Mass Ejection Shock Acceleration of Three Widespread SEP Events During Solar Cycle 24

We studied three solar energetic particle (SEP) events observed on 14 August 2010, 3 November 2011, and 5 March 2013 by Solar Terrestrial Relations Observatory (STEREO) A, B, and near-Earth (L1) spacecraft with a longitudinal distribution of particles greater than 90 degrees. Using a forward modeling method combined with extreme ultraviolet and white-light images, we determined the angular extent of the shock, the time and location (cobpoint) of the shock intersection with the magnetic field line connecting to each spacecraft, and compute the shock speed at the cobpoint of each spacecraft. We then examine whether the observations of SEPs at each spacecraft were accelerated and injected by the spatially extended shocks or whether another mechanism such as cross-field transport is required for an alternative explanation. Our analyses results indicate that the SEPs observed at the three spacecraft on 3 November, STEREO B (STB) and L1 on 14 August, and the 5 March SEP event at STEREO A (STA) can be explained by the direct shock acceleration. This is consistent with the observed significant anisotropies, short time delays between particle release times and magnetic connection times, and sharp rises in the SEP time profiles. Cross-field diffusion is the likely cause for the 14 August SEP event observed by STA and the 5 March SEPs observed by STB and L1 spacecraft, as particle observations featured weak electron anisotropies and slow rising intensity profiles. Otherwise, the wide longitudinal spread of these SEP increases would require an existence of a circumsolar shock, which may not be a correct assumption in the corona and heliosphere.

Xie, H.↗

A Sun-to-Earth Analysis of Magnetic Helicity of the 2013 March 17-18 Interplanetary Coronal Mass Ejection

We compare the magnetic helicity in the 2013 March 17-18 interplanetary coronal mass ejection (ICME) flux rope at 1 au and in its solar counterpart. The progenitor coronal mass ejection (CME) erupted on 2013 March 15 from NOAA active region 11692 and is associated with an M1.1 flare. We derive the source region reconnection flux using the post-eruption arcade (PEA) method that uses the photospheric magnetogram and the area under the PEA. The geometrical properties of the near-Sun flux rope is obtained by forward-modeling of white-light CME observations. Combining the geometrical properties and the reconnection flux, we extract the magnetic properties of the CME flux rope. We derive the magnetic helicity of the flux rope using its magnetic and geometric properties obtained near the Sun and at 1 au. We use a constant-α force-free cylindrical flux rope model fit to the in situ observations in order to derive the magnetic and geometric information of the 1 au ICME. We find a good correspondence in both amplitude and sign of the helicity between the ICME and the CME, assuming a semicircular (half torus) ICME flux rope with a length of π au. We find that about 83% of the total flux rope helicity at 1 au is injected by the magnetic reconnection in the low corona. We discuss the effect of assuming flux rope length in the derived value of the magnetic helicity. This study connecting the helicity of magnetic flux ropes through the Sun-Earth system has important implications for the origin of helicity in the interplanetary medium and the topology of ICME flux ropes at 1 au and hence their space weather consequences.

Pal, Sanchita↗

Remote Sensing, Uncertainty Quantification, and a Theory of Data Systems; Workshop Report

The purpose of the workshop was to invite statisticians, applied mathematicians, computer scientists, data system architects, experts in remote sensing technology, and Climate and Earth System scientists to review, discuss, and plan research on issues related to large-scale, efficient analysis of distributed data using spatial statistical methods. Our motivation in organizing this event was to catalyze interchange among experts on the fast-emerging problem of analysis of distributed data. As part of SAMSI's 2017-2018 Program on Mathematical and Statistical Methods for Climate and the Earth System, a Working Group on Remote Sensing was established to address statistical and mathematical research problems in the analysis of remote sensing data. The Working Group has five subgroups: 1) Spatial Retrieval Methodology (the so-called \Spatial-X" subgroup); 2) Spatial Analysis for Hyperspectral Data (the so-called \Spatial-Y" subgroup); 3) Emulators for Complex Forward Models; 4) Optimization for Remote Sensing Retrievals; and 5) Theory of Data Systems (ToDS). The ToDS subgroup spent the first half of this academic year formulating a framework in which to consider the joint problem of a) optimizing statistical methods for environments where data are distributed and too large to move to a central location, and b) the design of data system infrastructures within which to implement those statistical methods. To x ideas, the Workshop focused on spatial statistical methods. To date there are many new spatial statistical methods designed with massive data sets in mind, in the literature. However, very few have been implemented for remote sensing data, and none have been implemented in operational settings like those used by NASA and NOAA. A major impediment to their use in these cases is that the data are not only massive, but are stored in different physical locations. These data must be brought together in some way in order to estimate spatial covariance functions, but moving data to a central location for analysis is tedious at best and impossible at worst. Some remote data reduction is almost certainly necessary, but how much? What are the consequences for inference? The fundamental issue underlying these questions is how to navigate the trade-space between costs and uncertainty in the estimates or inferences that are ultimately produced.

Braverman, Amy↗

Physics-Based Retrieval of Surface Roughness Parameters for Bare Soils from Combined Active-Passive Microwave Signatures

In the past the effect of soil roughness was often considered secondary within the determination of soil moisture from remote sensing data. Several studies showed that accurate determination of soil roughness leads to an improved estimation of soil moisture. Two default parameters to describe the surface roughness are the standard deviation of the surface height variation 𝑠 and the surface correlation length 𝑙 with its corresponding autocorrelation function. Both parameters (𝑠,𝑙) affect the emissivity measured by radiometers as well as the backscattering observed by radars. In this study, we develop a physics-based approach to retrieve 𝑠 and 𝑙 by combining both microwave signals based on active-passive microwave covariation. To test the approach, containing a forward model and a retrieval algorithm, we used active/passive microwave data measured with the ComRAD truck-based SMAP simulator at L-band. Results and validations with corresponding field measurements on ground show that 𝑠 and 𝑙 can be estimated simultaneously when using this approach. The physics-based retrieval algorithm works robustly for two investigated test fields having an RMS-Error of 0.68 cm and 0.69 cm between the microwave-based and field-measured 𝑠-values, and of 3.13 cm and 3.04 cm for 𝑙-values. The first validation of the results reveals that the influence of the autocorrelation function, needed within the retrieval, is distinct.

correlation length↗

GRAIL-Identified Gravity Anomalies in Oceanus Procellarum: Insight into Subsurface Impact and Magmatic Structures on the Moon

Four, quasi-circular, positive Bouguer gravity anomalies (PBGAs) that are similar in diameter (~90-190 km) and gravitational amplitude (> 140 mGal contrast) are identified within the central Oceanus Procellarum region of the Moon. These spatially associated PBGAs are located south of Aristarchus Plateau, north of Flamsteed crater, and two are within the Marius Hills volcanic complex (north and south). Each is characterized by distinct surface geologic features suggestive of ancient impact craters and/or volcanic/plutonic activity. Here, we combine geologic analyses with forward modeling of high-resolution gravity data from the Gravity Recovery and Interior Laboratory (GRAIL) mission in order to constrain the subsurface structures that contribute to these four PBGAs. The GRAIL data presented here, at spherical harmonic degrees 6–660, permit higher resolution analyses of these anomalies than previously reported, and reveal new information about subsurface structures. Specifically, we find that the amplitudes of the four PBGAs cannot be explained solely by mare-flooded craters, as suggested in previous work; an additional density contrast is required to explain the high-amplitude of the PBGAs. For Northern Flamsteed (190 km diameter), the additional density contrast may be provided by impact-related mantle uplift. If the local crust has a density ~2800 kg/cu.m, then ~7 km of uplift is required for this anomaly, although less uplift is required if the local crust has a lower mean density of ~2500 kg/cu.m. For the Northern and Southern Marius Hills anomalies, the additional density contrast is consistent with the presence of a crustal complex of vertical dikes that occupies up to ~50% of the regionally thin crust. The structure of Southern Aristarchus Plateau (90 km diameter), an anomaly with crater-related topographic structures, remains ambiguous. Based on the relatively small size of the anomaly, we do not favor mantle uplift; however, understanding mantle response in a region of especially thin crust needs to be better resolved. It is more likely that this anomaly is due to subsurface magmatic material given the abundance of volcanic material in the surrounding region. Overall, the four PBGAs analyzed here are important in understanding the impact and volcanic/plutonic history of the Moon, specifically in a region of thin crust and elevated temperatures characteristic of the Procellarum KREEP Terrane.

Deutsch, Ariel N.↗

Improved Earth Oblateness Rate Reveals Increased Ice Sheet Losses and Mass-Driven Sea Level Rise

Satellite laser ranging (SLR) observations are routinely applied toward the estimation of dynamic oblateness, C(sub 20), which is the largest globally integrated component of Earth's time-variable gravity field. Since 2002, GRACE and GRACE Follow-On have revolutionized the recovery of higher spatial resolution features of global time-variable gravity, with SLR continuing to provide the most reliable estimates of C (sub 20).We quantify the effect of various SLR processing strategies on estimating C(sub 20) and demonstrate better signal recovery with the inclusion of GRACE-derived low-degree gravity information in the forward model. This improved SLR product modifies the Antarctic and Greenland Ice Sheet mass trends by -15.4 and -3.5 Gt/year, respectively, as compared to CSR TN11, and improves global mean sea level budget closure by modifying sea level rise by +0.08 mm/year. We recommend that this new C(sub 20) product be applied to RL06 GRACE data products for enhanced accuracy and scientific interpretation.

Loomis, B. D.↗

Statistical Study on Multispacecraft Widespread Solar Energetic Particle Events During Solar Cycle 24

We conduct a statistical study on the large three-spacecraft widespread solar energetic particle (SEP) events. Longitudinal distributions of the peak intensities, onset delays, and relation between the SEP intensity, coronal mass ejection (CME) shock speed, width, and the kinetic energy of the CME have been investigated. We apply a Gaussian fit to obtain the SEP intensity Io and distribution width and a forward-modeling fit to determine the true shock speed and true CME width. We found a good correction between and connection angel to the flare site and Io and the kinetic energy of the CME. By including the true chock speed and true CME widths, we reduce root-mean-square errors on the predicted SEP intensity by ~41% for protons compared to Richards et al.'s (2014, https://doi.org/10.1007/s11207-014-0524-8) prediction. The improved correction between the CME kinetic energy and SEP intensity provides strong evidence for the CME-shock acceleration theory of SEPs. In addition, we found that electron and proton release time delays (DTs) relative to Type II radio bursts increase with connection angles. The average electron (proton) DT is ~14 (32) min for strongly anisotropic events and ~2.5 (4.4) hr for weakly anisotropic events. Poor magnetic connectivity and large scattering effects are two main reasons to cause large delays.

Xie, H.↗

Informing Algorithms from Ground Validation: The GPM Combined Algorithm

Ground validation (GV) for the Global Precipitation Measurement (GPM) mission encompasses in-situ (e.g., gauge, disdrometer) measurements, ground radar products, and comprehensive datasets from dedicated airborne field campaigns. These datasets are used for direct validation of the precipitation products from GPM as well as to inform assumptions used by the algorithms that produce these products. This presentation will focus on current and potential such uses of GV data in the GPM combined radar-radiometer algorithm.The GPM combined algorithm, by virtue of using data from the Dual-frequency Precipitation Radar (DPR) and GPM Microwave Imager (GMI), is the most well-constrained instantaneous precipitation product from GPM. It also plays an important in role the passive microwave algorithms as a basis for the construction of brightness temperature-precipitation profile databases and is a calibrator for the multi- satellite gridded product IMERG. Since both radar and radiometer data are used as observational inputs, and even with these data the retrievals are underconstrained, microphysical properties of the hydrometeor profile that are relevant over the range of GMI and DPR wavelengths and incidence angles need to be assumed by the algorithm forward models. Such properties include the normalized intercept parameter (N w ) and shape parameter () of the particle size distribution, ice size-density relationship, and particle size-aspect ratio relationships. The sub-beam variability must also be prescribed in order to accurately simulate the observed radar reflectivity profiles without introducing significant biases. The sensitivity of the algorithm to these parameters, along with their mean values and variability, will be discussed. Examples of future directions include refinement with new data (particularly the assumptions related to ice), and seeking relationships between assumptions and observable or environmental data.

Munchak, Joe.↗

Advantages of Measuring the Q Stokes Parameter in Addition to the Total Radiance I in the Detection of Absorbing Aerosols: a Proof of Concept Using Vector Radiative Transfer, Neural Networks, and Optimal Estimation

A simple but novel study was conducted to investigate whether an imager-type spectroradiometer instrument like MODIS, currently flying on board the Aqua and Terra satellites, or MERIS, which flew on board Envisat, could detect absorbing aerosols if they could measure the Q Stokes parameter in addition to the total radiance I, that is if they could also measure the linear polarization of the light. Accurate radiative transfer calculations were used to train a fast neural network forward model, which together with a simple statistical optimal estimation scheme was used to retrieve three aerosol parameters: aerosol optical depth, fraction of absorbing aerosols, and aerosol location. The aerosols were assumed to be located either between 0 and 2 km or between 2 and 4 km in the Earth’s atmosphere. From simulated data it was found that by itself the total radiance I was generally insufficient to retrieve all three aerosol parameters, but that together with the Q Stokes component it was possible to reliably retrieve values of aerosol optical depth, fraction, and location.

Stamnes, Snorre↗