Search NASA⌕ Search

SEARCH · Search NASA

Results for “data intensive computing”

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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Dynamics of atoms in strong laser fields I: A quasi analytical model in momentum space based on a Sturmian expansion of the interacting nonlocal Coulomb potential

In this study we present a model that we have formulated in the momentum space to describe atoms interacting with intense laser fields. As a further step, it follows our recent theoretical approach in which the kernel of the reciprocal-space time-dependent Schrödinger equation (TDSE) is replaced by a finite sum of separable potentials, each of them supporting one bound state of atomic hydrogen (Tetchou Nganso et al. 2013). The key point of the model is that the nonlocal interacting Coulomb potential is expanded in a Coulomb Sturmian basis set derived itself from a Sturmian representation of Bessel functions of the first kind in the position space. As a result, this decomposition allows a simple spectral treatment of the TDSE in the momentum space. In order to illustrate the credibility of the model, we have considered the test case of atomic hydrogen driven by a linearly polarized laser pulse, and have evaluated analytically matrix elements of the atomic Hamiltonian and dipole coupling interaction. For various regimes of the laser parameters used in computations our results are in very good agreement with data obtained from other time-dependent calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Level 1 Processing of MODIS Direct Broadcast Data From Terra

In February 2000, an effort was begun to adapt the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 1 production software to process direct broadcast data. Three Level 1 algorithms have been adapted and packaged for release: Level 1A converts raw (level 0) data into Hierarchical Data Format (HDF), unpacking packets into scans; Geolocation computes geographic information for the data points in the Level 1A; and the Level 1B computes geolocated, calibrated radiances from the Level 1A and Geolocation products. One useful aspect of adapting the production software is the ability to incorporate enhancements contributed by the MODIS Science Team. We have therefore tried to limit changes to the software. However, in order to process the data immediately on receipt, we have taken advantage of a branch in the geolocation software that reads orbit and altitude information from the packets themselves, rather than external ancillary files used in standard production. We have also verified that the algorithms can be run with smaller time increments (2.5 minutes) than the five-minute increments used in production. To make the code easier to build and run, we have simplified directories and build scripts. Also, dependencies on a commercial numerics library have been replaced by public domain software. A version of the adapted code has been released for Silicon Graphics machines running lrix. Perhaps owing to its origin in production, the software is rather CPU-intensive. Consequently, a port to Linux is underway, followed by a version to run on PC clusters, with an eventual goal of running in near-real-time (i.e., process a ten-minute pass in ten minutes).

Lynnes, Christopher↗
DOE OSTI2015

Accelerating Asymptotically Exact MCMC for Computationally Intensive Models via Local Approximations

We construct a new framework for accelerating Markov chain Monte Carlo in posterior sampling problems where standard methods are limited by the computational cost of the likelihood, or of numerical models embedded therein. Our approach introduces local approximations of these models into the Metropolis–Hastings kernel, borrowing ideas from deterministic approximation theory, optimization, and experimental design. Previous efforts at integrating approximate models into inference typically sacrifice either the sampler’s exactness or efficiency; our work seeks to address these limitations by exploiting useful convergence characteristics of local approximations. We prove the ergodicity of our approximate Markov chain, showing that it samples asymptotically from the exact posterior distribution of interest. We describe variations of the algorithm that employ either local polynomial approximations or local Gaussian process regressors. Our theoretical results reinforce the key observation underlying this article: when the likelihood has some local regularity, the number of model evaluations per Markov chain Monte Carlo (MCMC) step can be greatly reduced without biasing the Monte Carlo average. Numerical experiments demonstrate multiple order-of-magnitude reductions in the number of forward model evaluations used in representative ordinary differential equation (ODE) and partial differential equation (PDE) inference problems, with both synthetic and real data. Supplementary materials for this article are available online.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Cloud Properties Derived From GOES-7 for Spring 1994 ARM Intensive Observing Period Using Version 1.0.0 of ARM Satellite Data Analysis Program

This document describes the initial formulation (Version 1.0.0) of the Atmospheric Radiation Measurement (ARM) program satellite data analysis procedures. Techniques are presented for calibrating geostationary satellite data with Sun synchronous satellite radiances and for converting narrowband radiances to top-of-the-atmosphere fluxes and albedos. A methodology is documented for combining geostationary visible and infrared radiances with surface-based temperature observations to derive cloud amount, optical depth, height, thickness, temperature, and albedo. The analysis is limited to two grids centered over the ARM Southern Great Plains central facility in north-central Oklahoma. Daytime data taken during 5 Apr. - 1 May 1994, were analyzed on the 0.3 deg and 0.5 deg latitude-longitude grids that cover areas of 0.9 deg x 0.9 deg and 10 deg x 14 deg, respectively. Conditions ranging from scattered low cumulus to thin cirrus and thick cumulonimbus occurred during the study period. Detailed comparisons with hourly surface observations indicate that the mean cloudiness is within a few percent of the surface-derived sky cover. Formats of the results are also provided. The data can be accessed through the World Wide Web computer network.

Minnis, Patrick↗
DOE OSTI2023

Sequence2Self: Self-supervised image sequence denoising of pixel-level spray breakup morphology

Optical imaging of fast and transient phenomena such as the turbulent breakup of liquid sprays exhibit low signal-to-noise ratios due to the limited illumination intensity relative to the short exposure time. Image denoising is required to facilitate physical studies over these data but is challenging due to the absence of clean ground-truths and the stringency of the denoising task (e.g., strong and complex noise, limited resolution, preserving physical fidelity), preventing supervised and existing un-/self-supervised deep learning methods. To this end, Sequence2Self (Seq2S) is proposed, an extension of Self2Self (S2S) to image sequences that leverages both the signal’s spatial and temporal correlation. Seq2S is demonstrated on time-resolved x-ray phase contrast imaging of liquid jet fuel sprays in a gas turbine combustor, which possesses all of challenges detailed above. Experiments are conducted across four fuels with different breakup morphology using various state-of-the-art methods. Overall, many of the methods failed and Seq2S was most successful: (1) Accurate spray structures were reconstructed with consistent evolution across frames void of artifacts. (2) The performance was robust, invariant to the hyperparameter choice. (3) Computational time is short and can be made eligible for real-time denoising. In particular, the images denoised by Seq2S showed spray droplet diameter distributions with near-zero Kullback–Leibler divergence (0.01 ± 0.01) to a cleaner reference, whereas the second best method yielded 0.06 ± 0.03. In conclusion, this suggests that Seq2S can be reliably used prior to subsequent quantitative spray analyses as it retains (if not, improves) the statistical physical properties of the data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

New Particle Formation and Growth in the Houston Atmosphere During TRACER (Final Report)

From 2020-2025, researchers from UC Irvine, UC Riverside, and Colorado State University collaborated on a Department of Energy-funded project to understand how airborne particles form and grow in urban atmospheres, conducting an intensive field campaign in Houston, Texas during summer 2022. Using advanced instruments to measure gas-phase chemicals, particle composition, and a specialized chamber to study particle growth, the team discovered that sulfur-containing compounds from industrial and power plant emissions are the dominant driver of new particle formation in Houston, with particles typically forming locally in the city and growing as air moves away in the urban plume. The research revealed an important methodological insight: measurements from fixed ground stations can be misleading when interpreting how particles actually evolve as air masses move, which has significant implications for how scientists worldwide interpret atmospheric observations. These findings improve understanding of urban air quality and help reduce uncertainties in climate models, since these particles play critical roles in cloud formation and Earth's radiation balance, while also providing detailed information about ultrafine particle composition relevant to public health. The project trained three doctoral students, developed enhanced computer models for urban particle formation, and made all data publicly available through the DOE Atmospheric Radiation Measurement data archive for use by the broader scientific community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

A patch-based pseudo-CT approach for MRI-only radiotherapy in the pelvis

Purpose: In radiotherapy based only on magnetic resonance imaging (MRI), knowledge about tissue electron densities must be derived from the MRI. This can be achieved by converting the MRI scan to the so-called pseudo-computed tomography (pCT). An obstacle is that the voxel intensities in conventional MRI scans are not uniquely related to electron density. The authors previously demonstrated that a patch-based method could produce accurate pCTs of the brain using conventional T{sub 1}-weighted MRI scans. The method was driven mainly by local patch similarities and relied on simple affine registrations between an atlas database of the co-registered MRI/CT scan pairs and the MRI scan to be converted. In this study, the authors investigate the applicability of the patch-based approach in the pelvis. This region is challenging for a method based on local similarities due to the greater inter-patient variation. The authors benchmark the method against a baseline pCT strategy where all voxels inside the body contour are assigned a water-equivalent bulk density. Furthermore, the authors implement a parallelized approximate patch search strategy to speed up the pCT generation time to a more clinically relevant level. Methods: The data consisted of CT and T{sub 1}-weighted MRI scans of 10 prostate patients. pCTs were generated using an approximate patch search algorithm in a leave-one-out fashion and compared with the CT using frequently described metrics such as the voxel-wise mean absolute error (MAE{sub vox}) and the deviation in water-equivalent path lengths. Furthermore, the dosimetric accuracy was tested for a volumetric modulated arc therapy plan using dose–volume histogram (DVH) point deviations and γ-index analysis. Results: The patch-based approach had an average MAE{sub vox} of 54 HU; median deviations of less than 0.4% in relevant DVH points and a γ-index pass rate of 0.97 using a 1%/1 mm criterion. The patch-based approach showed a significantly better performance than the baseline water pCT in almost all metrics. The approximate patch search strategy was 70x faster than a brute-force search, with an average prediction time of 20.8 min. Conclusions: The authors showed that a patch-based method based on affine registrations and T{sub 1}-weighted MRI could generate accurate pCTs of the pelvis. The main source of differences between pCT and CT was positional changes of air pockets and body outline.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A method for volumetric imaging in radiotherapy using single x-ray projection

Purpose: It is an intriguing problem to generate an instantaneous volumetric image based on the corresponding x-ray projection. The purpose of this study is to develop a new method to achieve this goal via a sparse learning approach. Methods: To extract motion information hidden in projection images, the authors partitioned a projection image into small rectangular patches. The authors utilized a sparse learning method to automatically select patches that have a high correlation with principal component analysis (PCA) coefficients of a lung motion model. A model that maps the patch intensity to the PCA coefficients was built along with the patch selection process. Based on this model, a measured projection can be used to predict the PCA coefficients, which are then further used to generate a motion vector field and hence a volumetric image. The authors have also proposed an intensity baseline correction method based on the partitioned projection, in which the first and the second moments of pixel intensities at a patch in a simulated projection image are matched with those in a measured one via a linear transformation. The proposed method has been validated in both simulated data and real phantom data. Results: The algorithm is able to identify patches that contain relevant motion information such as the diaphragm region. It is found that an intensity baseline correction step is important to remove the systematic error in the motion prediction. For the simulation case, the sparse learning model reduced the prediction error for the first PCA coefficient to 5%, compared to the 10% error when sparse learning was not used, and the 95th percentile error for the predicted motion vector was reduced from 2.40 to 0.92 mm. In the phantom case with a regular tumor motion, the predicted tumor trajectory was successfully reconstructed with a 0.82 mm error for tumor center localization compared to a 1.66 mm error without using the sparse learning method. When the tumor motion was driven by a real patient breathing signal with irregular periods and amplitudes, the average tumor center error was 0.6 mm. The algorithm robustness with respect to sparsity level, patch size, and presence or absence of diaphragm, as well as computation time, has also been studied. Conclusions: The authors have developed a new method that automatically identifies motion information from an x-ray projection, based on which a volumetric image is generated.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Angle-dependent strong-field molecular ionization rates with tuned range-separated time-dependent density functional theory

Strong-field ionization and the resulting electronic dynamics are important for a range of processes such as high harmonic generation, photodamage, charge resonance enhanced ionization, and ionization-triggered charge migration. Modeling ionization dynamics in molecular systems from first-principles can be challenging due to the large spatial extent of the wavefunction which stresses the accuracy of basis sets, and the intense fields which require non-perturbative time-dependent electronic structure methods. In this paper, we develop a time-dependent density functional theory approach which uses a Gaussian-type orbital (GTO) basis set to capture strong-field ionization rates and dynamics in atoms and small molecules. This involves propagating the electronic density matrix in time with a time-dependent laser potential and a spatial non-Hermitian complex absorbing potential which is projected onto an atom-centered basis set to remove ionized charge from the simulation. For the density functional theory (DFT) functional we use a tuned range-separated functional LC-PBE*, which has the correct asymptotic 1/r form of the potential and a reduced delocalization error compared to traditional DFT functionals. Ionization rates are computed for hydrogen, molecular nitrogen, and iodoacetylene under various field frequencies, intensities, and polarizations (angle-dependent ionization), and the results are shown to quantitatively agree with time-dependent Schrödinger equation and strong-field approximation calculations. This tuned DFT with GTO method opens the door to predictive all-electron time-dependent density functional theory simulations of ionization and ionization-triggered dynamics in molecular systems using tuned range-separated hybrid functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A comparison of deep-learning-based inpainting techniques for experimental X-ray scattering

The implementation is proposed of image inpainting techniques for the reconstruction of gaps in experimental X-ray scattering data. The proposed methods use deep learning neural network architectures, such as convolutional autoencoders, tunable U-Nets, partial convolution neural networks and mixed-scale dense networks, to reconstruct the missing information in experimental scattering images. In particular, the recovered pixel intensities are evaluated against their corresponding ground-truth values using the mean absolute error and the correlation coefficient metrics. The results demonstrate that the proposed methods achieve better performance than traditional inpainting algorithms such as biharmonic functions. Overall, tunable U-Net and mixed-scale dense network architectures achieved the best reconstruction performance among all the tested algorithms, with correlation coefficient scores greater than 0.9980.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

External beam planning module of Eclipse for external beam radiation therapy

Eclipse is a 3-dimensional (3D) treatment planning system for radiation therapy offered by Varian Medical Systems, Inc. The system has the network connectivity for the electronic transfer of image datasets and digital data communication among different equipment. The scope of this project for this special issue of Medical Dosimetry on 3D treatment planning systems is the assessment of planning tools in the external beam planning module of Eclipse to generate optimized treatment plans for patients undergoing external beam radiation therapy. This treatment planning system is relatively mature to be able to generate (1) simple treatment plans, (2) conformal radiation therapy plans, (3) static intensity-modulated radiation therapy (IMRT) plans, (4) volumetric-modulated arc therapy (VMAT) plans, and (5) treatment plans for electron beam therapy. The treatment planning tools are relatively plentiful to assist in the radiation therapy treatment planning. Some new features have been incorporated in the latest version and are helpful for making high-quality treatment plans. However, the location of the tools is not intuitive, and hence, familiarity with the user interface is essential to the efficient use of the treatment planning system. In addition, there are a number of dose algorithms available for the computation of dose distributions. The understanding of each dose computation algorithm is essential for the optimal use of this treatment planning system.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2015

An Investigation of Transonic Resonance in a Mach 2.2 Round Convergent-Divergent Nozzle

Hot-wire and acoustic measurements were taken for a round convergent nozzle and a round convergent-divergent (C-D) nozzle at a jet Mach number of 0.61. The C-D nozzle had a design Mach number of 2.2. Compared to the convergent nozzle jet flow, the Mach 2.2 nozzle jet flow produced excess broadband noise (EBBN). It also produced a transonic resonance tone at 1200 Herz. Computational simulations were performed for both nozzle flows. A steady Reynolds-Averaged Navier-Stokes simulation was performed for the convergent nozzle jet flow. For the Mach 2.2 nozzle flow, a steady RANS simulation, an unsteady RANS (URANS) simulation, and an unsteady Detached Eddy Simulation (DES) were performed. The RANS simulation of the convergent nozzle showed good agreement with the hot-wire velocity and turbulence measurements, though the decay of the potential core was over-predicted. The RANS simulation of the Mach 2.2 nozzle showed poor agreement with the experimental data, and more closely resembled an ideally-expanded jet. The URANS simulation also showed qualitative agreement with the hot-wire data, but predicted a transonic resonance at 1145 Herz. The DES showed good agreement with the hot-wire velocity and turbulence data. The DES also produced a transonic tone at 1135 Herz. The DES solution showed that the destabilization of the shock-induced separation region inside the nozzle produced increased levels of turbulence intensity. This is likely the source of the EBBN.

DES↗
NASA NTRS2020

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗
DOE OSTI2024

A robust approach to Gaussian process implementation

Abstract. Gaussian process (GP) regression is a flexible modeling technique used to predict outputs and to capture uncertainty in the predictions. However, the GP regression process becomes computationally intensive when the training spatial dataset has a large number of observations. To address this challenge, we introduce a scalable GP algorithm, termed MuyGPs, which incorporates nearest-neighbor and leave-one-out cross-validation during training. This approach enables the evaluation of large spatial datasets with state-of-the-art accuracy and speed in certain spatial problems. Despite these advantages, conventional quadratic loss functions used in the MuyGPs optimization, such as root mean squared error (RMSE), are highly influenced by outliers. We explore the behavior of MuyGPs in cases involving outlying observations and, subsequently, develop a robust approach to handle and mitigate their impact. Specifically, we introduce a novel leave-one-out loss function based on the pseudo-Huber function (LOOPH) that effectively accounts for outliers in large spatial datasets within the MuyGPs framework. Our simulation study shows that the LOOPH loss method maintains accuracy despite outlying observations, establishing MuyGPs as a powerful tool for mitigating unusual observation impacts in the large data regime. In the analysis of US ozone data, MuyGPs provides accurate predictions and uncertainty quantification, demonstrating its utility in managing data anomalies. Through these efforts, we advance the understanding of GP regression in spatial contexts.

Mukangango, Juliette↗
NASA NTRS2001

The Numerical Analysis of a Turbulent Compressible Jet

A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).

DeBonis, James R.↗
NASA NTRS2005

Modeling of the EUV and X-Ray Emission Spectra Induced by the Solar Winds Ions in the Heliosphere

We have carried out investigation of the EUV and X-ray emission spectra induced in interaction between the Solar Wind (SW) and interstellar neutral gas. The spectra of most important SW ions have been computed for the charge-exchange mechanism of X-ray emission using new accurate spectroscopic data from recent laboratory measurements and theoretical calculations. Total spectra have been constructed as a sum of spectra induced in the charge-exchange collisions by individual O(exp q+), C(exp q+), N(exp q+), Ne(exp q+), Mg (exp q+) and Fe(exp q+) ions. Calculations have been performed for X-ray emission from the heliospheric hydrogen and helium gas. X-ray maps of the heliosphere have been computed. The power density of X-ray sources in the heliospheric ecliptic plane is shown for the H gas and for the He gas. Distances from the Sun (0,0) are given in AU. The helium cone is clear seen in the X-ray map of the charge-exchange emission induced by the solar wind. X-ray emission spectra detected by the Chandra X-ray telescope from the "dark" side of Moon has been identified as a X-ray background emission induced by the solar wind from the geocorona. Spectra and intensities of this charge-exchange X-rays have been compared with the heliospheric component of the X-ray background. Observations and modeling of the SW spectra induced from the geocorona indicate a strong presence of emission lines of highly charged oxygen ions. Anisotropy in distribution of heliospheric X-rays has been predicted and calculated for the regions of the fast and slow solar winds.

Kharchenko, Vasili↗
DOE OSTI2020

Climate-Related Trends of Within-Storm Intensities Using Dimensionless Temporal-Storm Distributions

Huff curves are probabilistic time distributions of rainfall expressed as dimensionless cumulative percentages of storm depth and duration. Previous studies have documented development factors, spatial robustness, and the utility of Huff curves in practical applications. However, the effects of trending rainfall on Huff curve intensity patterns have not yet been studied. As such, the goal of this paper is to fill this gap by studying Huff curve patterns in a watershed with demonstrated increasing trends of temperature and precipitation, with the intention that it can be generalized to other areas in the US and the world. To achieve this goal, the high temporal resolution precipitation data collected from a high spatial density, 72-year precipitation-gauge network on the 4.25-km 2 North Appalachian Experimental Watershed in east-central Ohio were used. Seasonal storm pattern trends from 1939 to 2010 were investigated using dimensionless depth (with the frequency of 50%, $d_{50}$) and the curve variability $(V = d_{80} - d_{20})$ at three dimensionless within-storms time periods (three verticals). The Spearman rank correlation procedure (correlation coefficient, ρ and significance probability, p) was used to statistically determine trends over time using 8 periods of 4-season sets of Huff curves over the 72 years. Two averages of ρ and p were computed: (1) by averaging the individual ρ and p obtained from 10 gauges (AvgI); and (2) by grouped averaging of all individual gauge values of $d_{50}$ and V and then computing ρ and p (AvgG). The test results of individual gauges showed that 4 cases for $d_{50}$ and 23 cases for V were significant for all seasons and verticals (total of 120 cases for each variable). The test results of AvgI for $d_{50}$ and V and AvgG for $d_{50}$ showed no significant trends in all seasons and verticals. Only the AvgG for V led to a significant trend for V in spring and fall at different times within storm patterns. The data do not provide sufficient evidence at the p=0.05 significance level to reject the null hypothesis of unchanging position of the dimensionless depth of the 50% Huff curves for individual or averages for all seasons and verticals. Also, there is insufficient evidence to reject the null hypothesis of unchanging variability, V, using the average ρ of individual gauges (AvgI for V) for all seasons and verticals. AvgG results showed a significant trend in V; however, this analysis may be affected by the nonindependence of storm data. The results suggest that it is likely that there is little if any effect of trending climate over approximately 70 years on Huff curve patterns. The results of this study add to the robustness characteristics of Huff curves and to their potential use in hydrological practice as design storms, as the foundation of stochastic storm generation, and for storm disaggregation, and they deserve further investigation. These different forms of inputs to watershed models have the potential to improve runoff estimation. The results suggest that, if verified in other studies, they have applicability to provide useful stationary precipitation patterns across the US and other areas of the world in areas of nonstationary climate. Also, individual rain gauge data may not be representative of trends even over small areas, and seasonal differences were noticeable as found in other studies. Recommendations are provided.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Backus Effect on a Perpendicular Errors in Harmonic Models of Real vs. Synthetic Data

Measurements of geomagnetic scalar intensity on a thin spherical shell alone are not enough to separate internal from external source fields; moreover, such scalar data are not enough for accurate modeling of the vector field from internal sources because of unmodeled fields and small data errors. Spherical harmonic models of the geomagnetic potential fitted to scalar data alone therefore suffer from well-understood Backus effect and perpendicular errors. Curiously, errors in some models of simulated 'data' are very much less than those in models of real data. We analyze select Magsat vector and scalar measurements separately to illustrate Backus effect and perpendicular errors in models of real scalar data. By using a model to synthesize 'data' at the observation points, and by adding various types of 'noise', we illustrate such errors in models of synthetic 'data'. Perpendicular errors prove quite sensitive to the maximum degree in the spherical harmonic expansion of the potential field model fitted to the scalar data. Small errors in models of synthetic 'data' are found to be an artifact of matched truncation levels. For example, consider scalar synthetic 'data' computed from a degree 14 model. A degree 14 model fitted to such synthetic 'data' yields negligible error, but amplifies 4 nT (rmss) added noise into a 60 nT error (rmss); however, a degree 12 model fitted to the noisy 'data' suffers a 492 nT error (rmms through degree 12). Geomagnetic measurements remain unaware of model truncation, so the small errors indicated by some simulations cannot be realized in practice. Errors in models fitted to scalar data alone approach 1000 nT (rmss) and several thousand nT (maximum).

Voorhies, C. V.↗