SEARCH · Engineering Papers

Results for “spatial constraints”

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 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Atacama Cosmology Telescope: DR6 gravitational lensing and SDSS BOSS cross-correlation measurement and constraints on gravity with the 𝐸 𝐺 statistic

We derive new constraints on the 𝐸 𝐺 statistic as a test of gravity, combining the cosmic microwave background (CMB) lensing map estimated from Data Release 6 (DR6) of the Atacama Cosmology Telescope with Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS BOSS) CMASS and LOWZ galaxy data. We develop an analysis pipeline to measure the cross-correlation between CMB lensing maps and galaxy data, following a blinding policy and testing the approach through null and consistency checks. By testing the equivalence of the spatial and temporal gravitational potentials, the 𝐸 𝐺 statistic can distinguish Λ⁢ CDM from alternative models of gravity. We find 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢1$^{+0.06}_{−0.05}$ for Atacama Cosmology Telescope (ACT) and CMASS data at 68.28% confidence level, and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢9$^{+0.14}_{−0.11}$ for the ACT and LOWZ. Systematic errors are estimated to be 3% and 4%, respectively. Including CMB lensing information from Planck PR4 results in 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢4$^{+0.05}_{−0.05}$ with CMASS and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢3$^{+0.11}_{−0.09}$ with LOWZ. These are consistent with predictions for the Λ⁢ CDM model that best fits the Planck CMB anisotropy and SDSS BOSS baryon acoustic oscillations (BAO), where 𝐸$^{GR}_{𝐺⁡}$(𝑧 eff =0.555) =0.401 ± 0.005 for CMB lensing combined with CMASS and 𝐸$^{GR}_{𝐺}$⁡(𝑧 eff = 0.316) = 0.452 ± 0.005 combined with LOWZ. We also find 𝐸 𝐺 to be scale independent, with probability to exceed >5%, as predicted by general relativity. The methods developed in this work are also applicable to improved future analyses with upcoming spectroscopic galaxy samples and CMB lensing measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1982

Spatially resolved reflectivities of Saturn - 3000-6000 A

Absolute reflectivity reductions are presented for spectrophotometric observations of seven spatially-resolved regions on the disk of Saturn in the 3000-6000 A region. An asymmetry is found between northern and southern hemisphere data, in the equivalent width of the CH4, 5430-A band and in the continuum reflectivity. In addition, the spectra are found not to include the 4700 A reflectivity peak found by Bergstralh et al (1981), and to exhibit general agreement with the Podolack and Danielson (1977) model. The blue and visual spectral region continuum reflectivity of a planetary atmosphere provides constraints for determining aerosol particle single-scattering albedo and vertical distribution, which are in turn important to an understanding of planetary energy balance, atmospheric dynamics, and planetary molecular absorption spectrum interpretation.

Cochran, W. D.↗

NASA NTRS2013

Chemical Enrichment and Physical Conditions in IZw18*

Low-metallicity star-forming dwarf galaxies are prime targets to understand the chemical enrichment of the interstellar medium. The H I region contains the bulk of the mass in blue compact dwarfs, and it provides important constraints on the dispersal and mixing of heavy elements released by successive star-formation episodes. The metallicity of the H I region is also a critical parameter to investigate the future star-formation history, as metals provide most of the gas cooling that will facilitate and sustain star formation. Aims. Our primary objective is to study the enrichment of the H I region and the interplay between star-formation history and metallicity evolution. Our secondary obje ctive is to constrain the spatial- and time-scales over which the HI and H II regions are enriched, and the mass range of stars responsible for the heavy element production. Finally, we aim to examine the gas heating and cooling mechanisms in the H I region. Methods. We observed the most metal-poor star-forming galaxy in the Local Universe, I Zw 18, with the Cosmic Origin Spectrograph onboard Hubble. The abundances in the neutral gas are derived from far-ultraviolet absorption-lines (H I, CIII, CIIi*, N I, OI,...) and are compared to the abundances in the H II region. Models are constructed to calculate the ionization structure and the thermal processes. We investigate the gas cooling in the HIi region through physical diagnostics drawn from the fine-structure level of C+. Results. We find that H I region abundances are lower by a factor of approx 2 as compared to the H II region. There is no differential depletion on dust between the H I and H II region. Using sulfur as a metallicity tracer, we calculate a metallicity of 1/46 Z(solar) (vs. 1/31 Z(solar) in the H II region). From the study of the C/O, [O/Fe], and N/O abundance ratios, we propose that C, N, O, and Fe are mainly produced in massive stars. We argue that the H I envelope may contain pockets of pristine gas with a metallicity essentially null. Finally, we derive the physical conditions in the H I region by investigating the C II* absorption line. The cooling rate derived from C II* is consistent with collisions with H(sup 0) atoms in the diffuse neutral gas. We calculate the star-formation rate from the CII* cooling rate assuming that photoelectric effect on dust is the dominant gas heating mechanism. Our determination is in good agreement with the values in the literature if we assume a low dust-to-gas ratio (approx 2000 times lower than the Milky Way value).

far-ultraviolet absorption-lines↗

DOE OSTI2024

Analysis of static Wilson line correlators from lattice QCD at finite temperature with T -matrix approach

The thermodynamic T-matrix approach is used to study Wilson line correlators (WLCs) for a static quark-antiquark pair in the quark-gluon plasma (QGP). Selfconsistent results that incorporate constraints from the QGP equation of state can approximately reproduce WLCs computed in 2+1-flavor lattice-QCD (lQCD), provided the input potential exhibits less screening than in previous studies. Utilizing the updated potential to calculate pertinent heavylight T-matrices we evaluate thermal relaxation rates of heavy quarks in the QGP. We find a more pronounced temperature dependence for low-momentum quarks than in our previous results (with larger screening), which turns into a weaker temperature dependence of the (temperature-scaled) spatial diffusion coefficient, in fair agreement with the most recent lQCD data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2002

High-Resolution Gravity and Time-Varying Gravity Field Recovery using GRACE and CHAMP

This progress report summarizes the research work conducted under NASA's Solid Earth and Natural Hazards Program 1998 (SENH98) entitled High Resolution Gravity and Time Varying Gravity Field Recovery Using GRACE (Gravity Recovery and Climate Experiment) and CHAMP (Challenging Mini-satellite Package for Geophysical Research and Applications), which included a no-cost extension time period. The investigation has conducted pilot studies to use the simulated GRACE and CHAMP data and other in situ and space geodetic observable, satellite altimeter data, and ocean mass variation data to study the dynamic processes of the Earth which affect climate change. Results from this investigation include: (1) a new method to use the energy approach for expressing gravity mission data as in situ measurements with the possibility to enhance the spatial resolution of the gravity signal; (2) the method was tested using CHAMP and validated with the development of a mean gravity field model using CHAMP data, (3) elaborate simulation to quantify errors of tides and atmosphere and to recover hydrological and oceanic signals using GRACE, results show that there are significant aliasing effect and errors being amplified in the GRACE resonant geopotential and it is not trivial to remove these errors, and (4) quantification of oceanic and ice sheet mass changes in a geophysical constraint study to assess their contributions to global sea level change, while the results improved significant over the use of previous studies using only the SLR (Satellite Laser Ranging)-determined zonal gravity change data, the constraint could be further improved with additional information on mantle rheology, PGR (Post-Glacial Rebound) and ice loading history. A list of relevant presentations and publications is attached, along with a summary of the SENH investigation generated in 2000.

Shum, C. K.↗

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↗

DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Leveraging observed soil heterotrophic respiration fluxes as a novel constraint on global-scale models

Microbially-explicit models may improve understanding and projections of carbon dynamics in response to future climate change, but their fidelity in simulating global-scale soil heterotrophic respiration (RH), a stringent test for soil biogeochemical models, has never been evaluated. We used statistical global RH products, as well as 7,821 daily site-scale RH measurements, to evaluate the spatio-temporal performance of one first-order decay model (CASA-CNP) and two microbially-explicit biogeochemical models (CORPSE and MIMICS) that were forced by two different climate datasets. CORPSE and MIMICS did not provide any measurable performance improvement; instead, the models were highly sensitive to the meteorological input data used to drive them. Spatial RH variability was generally well simulated except in the northern middle latitudes (~50°N) and arid regions; models captured the seasonal variability of RH well, but showed more divergence in tropic and arctic regions. Our results demonstrate that the next generation of biogeochemical models shows promise, but also needs to be improved for realistic spatio-temporal variability of RH. Finally, we emphasize the importance of net primary production, soil moisture, and soil temperature inputs, and that jointly evaluating soil models for their spatial (global scale) and temporal (site scale) performance provides crucial benchmarks for improving biogeochemical models.

Jian, Jinshi↗

DOE OSTI2024

ELM-Wet: Inclusion of a wetland landunit with sub-grid representation of eco-hydrological patches and hydrological forcing in E3SM Land Model (ELM)

Wetlands emit the most biogenic methane (CH4) and present the greatest uncertainty in the global CH4 budget. Modeling these emissions is challenging due to the temporal and spatial variability in wetland structure and CH4 flux rates, along with complex interactions among hydrological, ecological, meteorological, and microbial processes that govern CH4 dynamics. To address these issues, we aim to enhance the accuracy of wetland representation in the U.S. Department of Energy’s Exascale Earth System Model (E3SM) Land Model, ELM. This effort led to the development of ELM-Wet, which incorporates a dedicated wetland landunit with subgrid representation of eco-hydrological patch types. We implement wetland-specific hydrology by imposing site-specific constraints on surface water levels, thereby allowing different patches to sustain varying depths of inundation. Additionally, we refined the calculation of aerenchyma transport diffusivity based on observed conductance across different vegetation types. We validated these enhancements through site-specific simulations of a coastal freshwater wetland, the Salvador WMA Freshwater Marsh (Ameriflux, site ID US-LA2), located in the coast of Louisiana (29.85N,90.29W) at an elevation of 0 m. The site was simulated with ELM-Wet and the default version ELMv1. We use Bayesian Optimization to parameterize CO2 and CH4 fluxes. Eddy covariance observations of CO2 and CH4 fluxes from 2012-2013 were used to train the model and data from 2021 were used for validation. In this repository, we include all the output of all the simulations performed using different versions of ELMv1 and ELM-Wet, all input required to the model, and the Matlab scripts we used to processed the output data. ELM-Wet_flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1995

Laser Doppler Velocimetry and full-field soot volume fraction

Since its introduction in the mid-sixties, Laser Doppler Velocimetry (LDV) has become one of the most widely used methods for the measurement of flows. Its remote and essentially non-intrusive nature provides an invaluable tool for a variety of difficult measurement situations which would be otherwise inaccessible. The high spatial resolution and rapid temporal response afforded by this technique are well suited to the determination of spatial and temporal details of flow fields, as well as characterization of turbulence. Advances in the understanding of the properties of LDV signals, accompanied by technological advances in coherent laser sources, detectors of high sensitivity and low noise, optical fabrication techniques and high-speed digital signal processing architectures have resulted in systems of increased accuracy and flexibility. As will be shown, recent progress in solid-state lasers and photo-detectors has been beneficial insofar as the compatibility of this method with the unique and severe constraints inherent in microgravity combustion science experiments.

Greenberg, Paul S.↗

NASA NTRS1986

Halos around planetary nebulae

Preliminary results of a CCD survey designed to detect and investigate faint halos around planetary nebulae are reported. A TI 800 x 800 pixel CCD was used to take deep exposures of 44 planetary nebulae. The exposures were obtained through an H-alpha filter at the Cassegrain focus of the Palomar 1.5 m telescope. Spatial resolutions of 1 to 2 arcsec were obtained across 400 arcsec wide fields. The images, which are in many cases considerably deeper than any previously taken, reveal numerous planetary nebula halos. About two-thirds of the studied nebulae possess extensive outer halos, here defined as any extended emission beyond the 10 percent isophote. Ionized sulphur electron density measurements show that in some nebulae, the mass in the halo is comparable to the mass contained in the primary H II region. The data have been used to place constraints on the mode of origin of the halos. It is likely that the halos originate either by dynamical separation of a single ejected shell of gas or by the ejection of two or more such shells from the central star. It is possible but less likely that the halos are caused by excitation of the preplanetary stellar wind and improbable that the halos represent reflection nebulae.

Jewitt, D. C.↗

NASA NTRS1986

A comparison of the bounded derivative and the normal-mode initialization methods using real data

Application of the bounded-derivative and normal-mode methods to a simple linear barotropic model at a typical middle latitude shows that the two methods lead to identical constraints up to a certain degree of approximation. Beyond this accuracy the two methods may differ from each other. When applied to a global nonlinear barotropic model using real data, again the two methods lead to similar balanced initial states. The gravity oscillations in the unbalanced height field, which have amplitudes of up to 60 m with a dominant periodicity of about 5 to 6 h, are practically eliminated by both initialization methods. The rotational wind component is smooth even for the unbalanced initial state. The small-scale spatial features of the irrotational wind component are drastically reduced by initialization. Both the nonlinear normal-mode and the bounded-derivative initialization methods yield similar divergence fields centered around the areas of highest orography. The comparison shows that there is no significant loss of information in the mass and momentum fields, despite the fact that the bounded-derivative method employs only the original, rotational wind component to construct a balanced initial state compared to the normal-mode method, which, in addition, makes use of the unbalanced divergent wind and height fields.

Semazzi, F. H. M.↗

NASA NTRS2018

Seasonal and Interannual Assessment of Cloud Cover and Atmospheric Constituents Across the Amazon (2000 - 2015): Insights for Remote Sensing and Climate Analysis

The quantitative assessment of cloud cover and atmospheric constituents improves our ability to exploit the climate feedback into the Amazon basin. In the 21st century, three droughts have already occurred in the Amazonia (e.g. 2005, 2010, 2015), inducing regional changes in the seasonal patterns of atmospheric constituents. In addition to climate, the atmospheric dynamic and attenuation properties are long-term challenges for satellite-based remote sensing of this ecosystem: high cloudiness, abundant water vapor content and biomass burning season. Therefore, while climatology analysis supports the understanding of atmospheric variability and trends, it also offers valuable insights for remote sensing applications. In this study, we evaluate the seasonal and interannual variability of cloud cover and atmospheric constituents (aerosol loading, water vapor and ozone content) over the Amazon basin, with focus on both climate analysis and remote sensing implications. We take the advantage of new atmosphere daily products at 1 km resolution derived from Multi-Angle Implementation for Atmospheric Correction (MAIAC) algorithm developed for Moderate Resolution Imaging Spectroradiometer (MODIS) data. An intercomparison of Aerosol Robotic Network (AERONET) and MAIAC aerosol optical depth (AOD) and columnar water vapor (CWV) showed quantitative information with a correlation coefficient higher than 0.81. Our results show distinct regional patterns of cloud cover across the Amazon basin: northwestern region presets a persistent cloud cover (> 80%) throughout the year, while low cloud cover (0-20%) occurs in the southern Amazon during the dry season. The cloud-free period in the southern Amazon is followed by an increase in the atmospheric burden due to fire emissions. Our results reveal that AOD records are changing in terms of area and intensity. During the 2005 and 2010 droughts, the positive AOD anomalies (δ > 0.1) occurred over 39.03% (240.3 million ha) and 27.14% (165.99 million ha) of total basin in the SON season, respectively. In contrast, the recent 2015 drought occurred towards the end of year (October through December) and these anomalies were observed over 23.72% (145 million ha) affecting areas in the central and eastern Amazon unlike previous droughts. The water vapor presents high concentration values (4.0-5.0 g/sq.cm) in the wet season (DJF), while we observed a strong spatial gradient from northwestern to southeastern of the basin during the dry season. In addition, we also found a positive trend of water vapor content ( ∼0.3 g/sq.cm) between 2000 and 2015. The total ozone typically varies between 220 and 270 DU, and it has a seasonal change of ∼25-35 DU from wet season to dry season caused by large emissions of ozone precursors and long-range transport. Finally, while this study contributes to climatological analysis of atmospheric constituents, the remote sensing users can also understand the regional constraints caused by atmospheric attenuation, such as high aerosol loading and cloud obstacles for surface observations.

ozone↗

NASA NTRS1992

A principled approach to the measurement of situation awareness in commercial aviation

The issue of how to support situation awareness among crews of modern commercial aircraft is becoming especially important with the introduction of automation in the form of sophisticated flight management computers and expert systems designed to assist the crew. In this paper, cognitive theories are discussed that have relevance for the definition and measurement of situation awareness. These theories suggest that comprehension of the flow of events is an active process that is limited by the modularity of attention and memory constraints, but can be enhanced by expert knowledge and strategies. Three implications of this perspective for assessing and improving situation awareness are considered: (1) Scenario variations are proposed that tax awareness by placing demands on attention; (2) Experimental tasks and probes are described for assessing the cognitive processes that underlie situation awareness; and (3) The use of computer-based human performance models to augment the measures of situation awareness derived from performance data is explored. Finally, two potential example applications of the proposed assessment techniques are described, one concerning spatial awareness using wide field of view displays and the other emphasizing fault management in aircraft systems.

Tenney, Yvette J.↗

NASA NTRS1993

Infrared studies of circumstellar matter around Herbig Ae/Be and related stars

The study presents high spatial resolution data at 50 and 100 microns for seven young intermediate-mass stars with flat or rising spectral energy distributions. Five stars, namely, LkH-alpha 198, R Mon, CD -42 deg 11721, R CrA, and V645 Cyg, are resolved at 100 microns. The emission comes from spatially extended envelopes, with sizes varying between 0.025 and 0.4 pc. The wavelengths are optically thin at mid-IR and FIR wavelengths but become optically thick in the visual where the extinction amounts to 5-10 mag. In two cases, R CrA and V 645 Cyg, dust density distributions with steep density gradients are obtained, alpha is approximately equal to 2, while for LkH-alpha, R Mon, and CD -42 deg 11721, alpha is approximately equal to 0.5. Model fitting to the spectral energy distributions, combined with the constraints on the size of the emitting regions, indicates that inside the circumstellar envelopes a large fraction of the stellar radiation is already degraded to longer wavelengths.

Natta, Antonella↗

DOE OSTI2025

Global terrestrial nitrogen fixation and its modification by agriculture

Biological nitrogen fixation (BNF) is the largest natural source of new nitrogen (N) that supports terrestrial productivity, yet estimates of global terrestrial BNF remain highly uncertain. Here, in this study, we show that this uncertainty is partly because of sampling bias, as field BNF measurements in natural terrestrial ecosystems occur where N fixers are 17 times more prevalent than their mean abundances worldwide. To correct this bias, we develop new estimates of global terrestrial BNF by upscaling field BNF measurements using spatially explicit abundances of all major biogeochemical N-fixing niches. We find that natural biomes sustain lower BNF, 65 (52–77) Tg N yr −1 , than previous empirical bottom-up estimates, with most BNF occurring in tropical forests and drylands. We also find high agricultural BNF in croplands and cultivated pastures, 56 (54–58) Tg N yr −1 . Agricultural BNF has increased terrestrial BNF by 64% and total terrestrial N inputs from all sources by 60% over pre-industrial levels. Our results indicate that BNF may impose stronger constraints on the carbon sink in natural terrestrial biomes and represent a larger source of agricultural N than is generally considered in analyses of the global N cycle, with implications for proposed safe operating limits for N use.

Reis Ely, Carla R. [Oregon State Univ., Corvallis,↗

DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

On the temperature and velocity through the photosphere of a sunspot penumbra

We investigate the structure in depth of a sunspot penumbra by means of the inversion code of the radiative transfer equation proposed by Ruiz Cobo & del Toro Iniesta (1992), applied to a set of filtergrams of a sunspot, scanning the Fe I line at 5576.1 A, with a sampling interval of 30 mA, from -120 to 120 mA from line center (data previously analyzed by Title et al. 1993). The temperature structure of this penumbra is obtained for each of the 801 pixels selected (0.32 sec x 0.32 sec). On the average, the temperatures seem to decrease as we move inward, but the differences are of the order of the rms values (approximately equal 100-200 K) at a given distance to sunspot center. The outer parts of the penumbra have also a bigger curvature in the T versus log tau(sub 5) relation than the inner parts. We realize, however, that these differences might be influenced by possible stray light effects. Compared to the quiet Sun, penumbral temperatures are cooler at deep layers and hotter at high layers. A mean penumbral model atmosphere is presented. The asymmetries observed in the intensity profile (the line is magnetically insensitive) are deduced to be produced by strong gradients of the line-of-sight velocity that sharply vary spatially along slices of almost constant distance to sunspot center. These variations suggest that such gradients are not only needed to explain the broadband circular polarization observed in sunspots (see Sanchez Almeida & Lites 1992) but are a main characteristic of the fine-scale penumbra. The results are compatible with an Evershed flow present everywhere, but its gradient with depth turns out to vary so that the flow seems to be mainly concentrated in some penumbral fibrils when studied through Dopplergrams. Finally, as by-products of this study, we put constraints to the practical usefulness of the Eddington-Barbier relation, and we explain the values of the Fourier Dopplergrams to be carrying information of layers around the centroid of the generalized response function of Dopplergrams to velocity fluctuations.

Del Toro Iniesta, J. C.↗