Results for “Analysis of Variance” — Page 83

DOE OSTI2026

Using AI to Reproduce Neutrino Cross Section Analysis - Prototyping the Neutrino Discovery Platform

The Neutrino Discovery Platform (NDP) aims to accelerate DUNE-era science by making the neutrino program's existing datasets analyzable through fast, reproducible, and auditable workflows. We report a working version of two of its layers, data curation and agentic orchestration, built and tested end to end on MINERvA open data. The guiding lesson throughout is that a cross section is a measurement, and not just a plotted shape, only if it carries a defensible systematic-uncertainty budget, a trustworthy unfolding, and a reproducible record. Using a single medium-energy playlist pair from the MINERvA open-data release (about $2.05\times10^{17}$ protons on target of data), we first reproduced the shapes of two published charged-current inclusive $\nu_\mu$ measurements through a complete extraction ladder: selection, background subtraction, D'Agostini unfolding, efficiency correction, and flux normalization. These shape-level reproductions ran and tracked the published results, but they lacked the systematic-uncertainty machinery that defines a MINERvA cross section. To supply it, we vendored and built the MINERvA Analysis Toolkit and developed a many-universe systematic-uncertainty tool that produces a portable covariance artifact, a parallel event-loop runner, and a per-run auditability harness. Validated against a published covariance release, the toolchain reproduces the released statistical, flux, and muon-energy-scale terms and shows that they account for roughly 63\% of the total variance, with the remainder unreleased. Using this same infrastructure, we then performed a measurement of our own design, the hadronic recoil-energy distribution of low-energy ($E_\nu<2.5$~GeV) charged-current inclusive events, and found data/simulation shape agreement of $\chi^2/\mathrm{ndf}=1.26$. Together these results show that the platform supports original physics and not only reproductions.

Breaux, Auto [Tulane U. (main)]↗

NASA NTRS1993

Atmospheric moisture fields derived by satellite observations over the tropical Pacific Ocean

Values of precipitable water are retrieved over the tropical and subtropical Pacific Ocean from TOVS infrared and microwave channel brightness temperature and OLR observations by means of stepwise linear regression. The most useful temperature and moisture sensing channels are pre-selected from sensitivity tests of a radiative transfer model. Numerous models are developed and tested against collocated radiosonde observations and Nimbus-7 SMMR precipitable water estimates. For RAOB comparisons, the best estimator used 15 TOVS predictors and captured 71.1 deg percent of the variance (+0.62 g/sq cm standard error) for column precipitable water; for precipitable water of 700-500 mb bulk layer, these values were 71.7 percent and +/- 0.17 g/sq cm. Little skill of estimated precipitable water was obtained for moisture above 500 mb. Regressions were less skillful against SMMR, unless collocation parameters were tightly controlled; SMMR was less acceptable than RAOB's because of observational drift and errors. Generally, the most skillful predictors were boundary layer brightness temperatures of TOVS channels and satellite estimated stability indices. 'Moisture channels' were hardly useful except for estimating middle and upper tropospheric moisture. Additional regression models were constructed testing the sensitivity to different observational and meteorological characteristics. Models which used some in situ observations surface observations or stabilities calculated from RAOB, were the most successful. The best of these explained 87.5 percent of the variance but the regression selected almost no TOVS channels, relying instead on conventional RAOB and surface observations. A set of four regression models were developed, stratifying atmospheric characteristics on the basis of collocated OLR values. These models improved the variance explained by 5.0 percent; the model associated with the highest OLR values (275 W/sq m less than or equal to OLR; that is, no cloud) showed only marginal skill. Precipitable water fields were generated from the best TOVS-only model for seven days in January 1983 and compared with SMMR-estimated fields. OLR fields and ECMWF precipitable water analysis. The TOVS regression model compared favorably to the SMMR analysis in amplitudes and features. It revealed more evolving synoptic signal than the ECMWF analysis. In synoptically active regions, it differed with respect to the OLR analysis, primarily because of actual differences in the vertical distribution of water vapor.

Chung, Hyosang↗

NASA NTRS2023

Geometric Analysis and Visualization of Maxillofacial Anthropometry

As part of head borne systems, the maxillofacial and mandibular components are critical for Soldier protection in extreme and dynamic conditions. Often this equipment must not only cover the lower face of the wearer, but also accommodate motion of the lower jaw and resultant deformation of the face shape. Although the application is different, head borne equipmentdesigned for astronauts has similarly critical design requirements including comfort and protection from other environmental hazards while allowing for maximum performance. To improve the design of such systems and to ensure adequate accommodation of the user population, NASA Johnson Space Center is collaborating with the U.S. Army’s Development Command Soldier Center (DEVCOM SC) to statistically analyze facial landmark and measurement data from large Soldier anthropometric databases ( ANSUR 1988 and ANSUR II). Anthropometric measurements were incorporated from 5383 males and 3594 females for analysis, including individuals aged 17-58 across a range of races. Facial landmarks were analyzed, specifically the 3D coordinates of the Menton, Prementon, Gonion, Tragion, Stomion, Zygion, Pronasale, and Subnasale. These data were statistically correlated with demographiccharacteristics, such as sex and population origin. A statistical model of the lower head-face surface geometry is currently in development to describe the variance of head and face shape across the population. The model will be able to visually show the statistical trends of the different shapes and sizes of the head, along with the synthetically represented boundary cases. The specific analytical methods and outcomes will be detailed in the final proceeding. This study will provide a more comprehensive assessment of how face shape varies across the military population compared to traditional linear measurements. The outcome of this study is expected to provide critical insight and recommendations for improved protective equipment design and development, not only with Soldiers but also astronauts in future exploration missions, as well as first responders and others that require critical equipment fitted to the maxillofacial and or mandibular region.

W Green↗

DOE OSTI2026

Toward the Neutrino Discovery Platform: An Auditable, Uncertainty-Bearing Toolchain for MINERvA Open-Data Cross-Section Analysis

The Neutrino Discovery Platform (NDP) aims to accelerate DUNE-era science by making the neutrino program's existing datasets analyzable through fast, reproducible, and auditable workflows. We report a working version of two of its layers, data curation and agentic orchestration, built and tested end to end on MINERvA open data. The guiding lesson throughout is that a cross section is a measurement, and not just a plotted shape, only if it carries a defensible systematic-uncertainty budget, a trustworthy unfolding, and a reproducible record. Using a single medium-energy playlist pair from the MINERvA open-data release (about $2.05\times10^{17}$ protons on target of data), we first reproduced the shapes of two published charged-current inclusive $\nu_\mu$ measurements through a complete extraction ladder: selection, background subtraction, D'Agostini unfolding, efficiency correction, and flux normalization. These shape-level reproductions ran and tracked the published results, but they lacked the systematic-uncertainty machinery that defines a MINERvA cross section. To supply it, we vendored and built the MINERvA Analysis Toolkit and developed a many-universe systematic-uncertainty tool that produces a portable covariance artifact, a parallel event-loop runner, and a per-run auditability harness. Validated against a published covariance release, the toolchain reproduces the released statistical, flux, and muon-energy-scale terms and shows that they account for roughly 63\% of the total variance, with the remainder unreleased. Using this same infrastructure, we then performed a measurement of our own design, the hadronic recoil-energy distribution of low-energy ($E_\nu<2.5$~GeV) charged-current inclusive events, and found data/simulation shape agreement of $\chi^2/\mathrm{ndf}=1.26$. Together these results show that the platform supports original physics and not only reproductions.

Breaux, Auto [Tulane U. (main)]↗

DOE OSTI2019

Performance Analysis and Tank Test Validation of a Hybrid Wave-Current Energy Converter with a Single Power Takeoff

Marine and hydrokinetic (MHK) energy, including ocean waves, tidal current, ocean current and river current, has been recognized as a promising power source due to its full-day availability and high energy potential. At this stage, ocean current energy, tidal energy and ocean wave energy are currently the most competitive sourves among all the categories of MHK. The state of art MHK energy harvesting technology mainly focus on harvesting either ocean wave energy or current energy, but not both. However, a significant amount of ocean waves and tidal/ ocean current coexist in many sites and traditional devices that harvest from a single form of MHK energy, cannot make full use of the coexisting ocean energy. Furthermore, MHK energy harvesting devices need to advance to be cost-effective and competitive with other energy sources. This is difficult to achieve. Ocean wave excitation is irregular, which means that ocean wave height and wave periods are unpredictable and excitation forces on energy harvesting devices can have large variance in amplitude and frequency. Such problems/ restrictions can be possibly addressed by the concept of a hybrid energy converter. In this sense, a hybrid wave-current ocean energy conveter (HWCEC) that simutaneously harvests energy from current and wave with one single power takeoff (PTO) is designed.The wave energy is extracted through relative heaving motion between a floating buoy and a submerged second body, while the current energy is extracted using a marine current turbine (MCT). Energy from both sources are integrated by a hybrid PTO whose concept is based on a mechanical motion rectifier (MMR). In this study, different working modes are investigated together with switching criteria.Simulations were conducted with hydrodynamic coefficients obtained from computational fluid dynamics analysis and boundary element method. Tank tests were conducted for a HWCEC under coexisting wave and current inputs. For comparison, separate baseline tests of a turbine and a two-body point absorber, each acting in isolation, are conducted. Experimental results validate the dynamic modeling and show that a HWCEC can increase the output power with a range between 29-87 percent over either current turbine and wave energy converter acting individually, and it can reduce by up to 70 percent the peak-to-average power ratio compared with the wave energy converter on the tested conditions.Such results demonstrate the potential of the HWCEC as an efficient and cost-effective design.

16 TIDAL AND WAVE POWER↗

DOE OSTI2018

Modeling the He II Transverse Proximity Effect: Constraints on Quasar Lifetime and Obscuration

The He II transverse proximity effect - enhanced He II Lyα transmission in a background sightline caused by the ionizing radiation of a foreground quasar - offers a unique opportunity to probe the emission properties of quasars, in particular the emission geometry (obscuration, beaming) and the quasar lifetime. Building on the foreground quasar survey published in Schmidt et al., we present a detailed model of the He II transverse proximity effect, specifically designed to include light travel time effects, finite quasar ages, and quasar obscuration. Here, we post-process outputs from a cosmological hydrodynamical simulation with a fluctuating He ii ultraviolet background model, with the added effect of the radiation from a single bright foreground quasar. We vary the age t age and obscured sky fractions Ω obsc of the foreground quasar, and explore the resulting effect on the He ii transverse proximity effect signal. Fluctuations in intergalactic medium density and the ultraviolet background, as well as the unknown orientation of the foreground quasar, result in a large variance of the He II Lyα transmission along the background sightline. We develop a fully Bayesian statistical formalism to compare far-ultraviolet He ii Lyα transmission spectra of the background quasars to our models, and extract joint constraints on t age and Ω obsc for the six Schmidt et al. foreground quasars with the highest implied He II photoionization rates. Our analysis suggests a bimodal distribution of quasar emission properties, whereby one foreground quasar, associated with a strong He II transmission spike, is relatively old (22 Myr) and unobscured (Ω obsc < 35%), whereas three others are either younger than 10 Myr or highly obscured (Ω obsc > 70%).

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2016

Green Propellant Landing Demonstration at U.S. Range

The Green Propellant Loading Demonstration (GPLD) was conducted December 2015 at Wallops Flight Facility (WFF), leveraging work performed over recent years to bring lower toxicity hydrazine replacement green propellants to flight missions. The objective of this collaboration between NASA Goddard Space Flight Center (GSFC), WFF, the Swedish National Space Board (SNSB), and Ecological Advanced Propulsion Systems (ECAPS) was to successfully accept LMP-103S propellant at a U.S. Range, store the propellant, and perform a simulated flight vehicle propellant loading. NASA GSFC Propulsion (Code 597) managed all aspects of the operation, handling logistics, preparing the procedures, and implementing the demonstration. In addition to the partnership described above, Moog Inc. developed an LMP-103S propellant-compatible titanium rolling diaphragm flight development tank and loaned it to GSFC to act as the GPLD flight vessel. The flight development tank offered the GPLD an additional level of flight-like propellant handling process and procedures. Moog Inc. also provided a compatible latching isolation valve for remote propellant expulsion. The GPLD operation, in concert with Moog Inc. executed a flight development tank expulsion efficiency performance test using LMP-103S propellant. As part of the demonstration work, GSFC and WFF documented Range safety analyses and practices including all elements of shipping, storage, handling, operations, decontamination, and disposal. LMP-103S has not been previously handled at a U.S. Launch Range. Requisite for this activity was an LMP-103S Risk Analysis Report and Ground Safety Plan. GSFC and WFF safety offices jointly developed safety documentation for application into the GPLD operation. The GPLD along with the GSFC Propulsion historical hydrazine loading experiences offer direct comparison between handling green propellant versus safety intensive, highly toxic hydrazine propellant. These described motives initiated the GPLD operation in order to investigate the handling and process safety variances in project resources between LMP-103S and typical in-space propellants. The GPLD risk reduction operation proved successful for many reasons including handling the green propellant at a U.S. Range, loading and pressurizing a flight-like tank, expelling the propellant, measuring the tank expulsion efficiency, and most significantly, GSFC propulsion personnel's new insight into the LMP-103S propellant handling details.

Propellant↗

NASA NTRS2016

Green Propellant Loading Demonstration at U.S. Range

The Green Propellant Loading Demonstration (GPLD) was conducted December 2015 at Wallops Flight Facility (WFF), leveraging work performed over recent years to bring lower toxicity hydrazine replacement green propellants to flight missions. The objective of this collaboration between NASA Goddard Space Flight Center (GSFC), WFF, the Swedish National Space Board (SNSB), and Ecological Advanced Propulsion Systems (ECAPS) was to successfully accept LMP-103S propellant at a U.S. Range, store the propellant, and perform a simulated flight vehicle propellant loading. NASA GSFC Propulsion (Code 597) managed all aspects of the operation, handling logistics, preparing the procedures, and implementing the demonstration. In addition to the partnership described above, Moog Inc. developed an LMP-103S propellant-compatible titanium rolling diaphragm flight development tank and loaned it to GSFC to act as the GPLD flight vessel. The flight development tank offered the GPLD an additional level of flight-like propellant handling process and procedures. Moog Inc. also provided a compatible latching isolation valve for remote propellant expulsion. The GPLD operation, in concert with Moog Inc. executed a flight development tank expulsion efficiency performance test using LMP-103S propellant. As part of the demonstration work, GSFC and WFF documented Range safety analyses and practices including all elements of shipping, storage, handling, operations, decontamination, and disposal. LMP-103S has not been previously handled at a U.S. Launch Range. Requisite for this activity was an LMP-103S Risk Analysis Report and Ground Safety Plan. GSFC and WFF safety offices jointly developed safety documentation for application into the GPLD operation. The GPLD along with the GSFC Propulsion historical hydrazine loading experiences offer direct comparison between handling green propellant versus safety intensive, highly toxic hydrazine propellant. These described motives initiated the GPLD operation in order to investigate the handling and process safety variances in project resources between LMP-103S and typical in-space propellants. The GPLD risk reduction operation proved successful for many reasons including handling the green propellant at a U.S. Range, loading and pressurizing a flight-like tank, expelling the propellant, measuring the tank expulsion efficiency, and most significantly, GSFC propulsion personnel's new insight into the LMP-103S propellant handling details.

Green Propellant↗

NASA NTRSDate not supplied

Analysis and Evaluation of NASA's MK40B Six-Component TASK Balance

The paper discusses the physical characteristics NASA's MK40B TASK balance. It is NASA's newest TASK balance that was optimized for the testing of wind tunnel models in the NASA Ames 11-ft Transonic Wind Tunnel. The calibration data of the MK40B balance was obtained in Calspan's Automatic Balance Calibration System. A machine load schedule design was used for the calibration that was developed at the NASA Ames Balance Calibration Laboratory. The calibration data was analyzed using both the Non-Iterative Method and the Iterative Method. Then, the chosen regression models of the calibration data were assessed and compared using the percent contribution, the variance inflation factor, and the calibration load residuals. Interactions of single-component loads were also quantified. In addition, the bi-directional characteristics of the six gage outputs of the balance were assessed by applying a semi-empirical test that was developed at the Ames Balance Calibration Laboratory. It was concluded that four of the six gages of the balance have bi-directional characteristics. Finally, the accuracy of the load prediction equations of the balance was tested by applying the equations to independent check load data that was obtained at the NASA Ames Balance Calibration Laboratory. The standard deviation of the check load residuals met expectations. The check load data also confirmed that the interactions of the forward and aft normal forces on the axial force gage output are very small. - This paper is to be presented in the invited special session titled: "65 Years of TASK Wind Tunnel Balances - Past, Present, and Future".

wind tunnel balance↗

DOE OSTI2015

A Calibration to Predict the Concentrations of Impurities in Plutonium Oxide by Prompt Gamma Analysis Revision 2

This report documents the new PG calibration regression equation. These calibration equations incorporate new data that have become available since revision 1 of “A Calibration to Predict the Concentrations of Impurities in Plutonium Oxide by Prompt Gamma Analysis” was issued [3] The calibration equations are based on a weighted least squares (WLS) approach for the regression. The WLS method gives each data point its proper amount of influence over the parameter estimates. This gives two big advantages, more precise parameter estimates and better and more defensible estimates of uncertainties. The WLS approach makes sense both statistically and experimentally because the variances increase with concentration, and there are physical reasons that the higher measurements are less reliable and should be less influential. The new magnesium calibration includes a correction for sodium and separate calibration equation for items with and without chlorine. These additional calibration equations allow for better predictions and smaller uncertainties for sodium in materials with and without chlorine. Chlorine and sodium have separate equations for RICH materials. Again, these equations give better predictions and smaller uncertainties chlorine and sodium for RICH materials.

Instrumentation Related to Nuclear Science & Techn↗

NASA NTRS1973

Bayesian estimation of life parameters in the Weibull distribution.

Development of a Bayesian analysis of the scale and shape parameters in the Weibull distribution and the corresponding reliability function with respect to the usual life-testing procedures. For the scale parameter theta, Bayesian estimates of theta and reliability are obtained for the uniform, exponential, and inverted gamma prior probability densities. Bhattacharya's results (1967) for the one-parameter exponential life-testing distribution are reduced to a special case of these results. A fully Bayesian analysis of both the scale and shape parameters is developed by assuming independent prior distributions; since in the latter case, analytical tractability is not possible, Bayesian estimates are obtained through a conjunction of Monte Carlo simulation and numerical-integration techniques. In both cases, a computer simulation is carried out, and a comparison is made between the Bayesian and the corresponding minimum-variance unbiased, or maximum likelihood, estimates. As expected, the Bayesian estimates are superior.

Canavos, G. C.↗

NASA NTRS1984

Using the solar elevation angle and radiance conversion to normalize forest spectral signatures

A number of spectral signals associated with forest cover were selected in order to test the usefulness of adjustments to Landsat MSS digital count values. The signatures were taken from 26 independent regional land cover inventories in the northeastern US. A regression analysis was conducted which revealed a significant relationship between digital count values and variation in the solar elevation angles in MSS bands four, six, and seven for deciduous forests and in all MSS bands for conifers. When signatures were adjusted for solar elevation angles and MSS sensor calibration differences, the dependence on solar elevation angle was reduced. The adjustments also significantly reduced the variance with Level II forest categories among scenes.

Middleton, E. M.↗

NASA NTRS1992

Detector noise statistics in the non-linear regime

The statistical behavior of an idealized linear detector in the presence of threshold and saturation levels is examined. It is assumed that the noise is governed by the statistical fluctuations in the number of photons emitted by the source during an exposure. Since physical detectors cannot have infinite dynamic range, our model illustrates that all devices have non-linear regimes, particularly at high count rates. The primary effect is a decrease in the statistical variance about the mean signal due to a portion of the expected noise distribution being removed via clipping. Higher order statistical moments are also examined, in particular, skewness and kurtosis. In principle, the expected distortion in the detector noise characteristics can be calibrated using flatfield observations with count rates matched to the observations. For this purpose, some basic statistical methods that utilize Fourier analysis techniques are described.

Shopbell, P. L.↗

DOE OSTI2016

Analysis of dpa rates in the HFIR reactor vessel using a hybrid Monte Carlo/deterministic method

The Oak Ridge High Flux Isotope Reactor (HFIR), which began full-power operation in 1966, provides one of the highest steady-state neutron flux levels of any research reactor in the world. An ongoing vessel integrity analysis program to assess radiation-induced embrittlement of the HFIR reactor vessel requires the calculation of neutron and gamma displacements per atom (dpa), particularly at locations near the beam tube nozzles, where radiation streaming effects are most pronounced. In this study we apply the Forward-Weighted Consistent Adjoint Driven Importance Sampling (FW-CADIS) technique in the ADVANTG code to develop variance reduction parameters for use in the MCNP radiation transport code. We initially evaluated dpa rates for dosimetry capsule locations, regions in the vicinity of the HB-2 beamline, and the vessel beltline region. We then extended the study to provide dpa rate maps using three-dimensional cylindrical mesh tallies that extend from approximately 12 below to approximately 12 above the axial extent of the core. The mesh tally structures contain over 15,000 mesh cells, providing a detailed spatial map of neutron and photon dpa rates at all locations of interest. Relative errors in the mesh tally cells are typically less than 1%.

Blakeman, Edward↗

DOE OSTI2021

Transformational challenge reactor analysis to inform preconceptual core design decisions: Sensitivity study of transient analysis in a hydride-moderated microreactor

The Transformational Challenge Reactor (TCR) program aims to demonstrate a revolutionary design approach enabled by advanced manufacturing and data analytics in the design of nuclear reactors. This article discusses scoping analyses of preconceptual designs to inform TCR design decisions and the evaluation of sensitivities and uncertainties on postulated transient scenarios. The applicability of the systems codes TRACE and RELAP5-3D to TCR transient analysis are examined, and RELAP5-3D models are used to examine the transient response of two candidate core designs at multiple power levels. Then, the uncertainty quantification code RAVEN is used to quantify the effect of several design parameters on reactivity-initiated accident (RIA) progression at hot zero power (HZP) and hot full power (HFP) as well as to assess the impact of uncertainties in transient parameters for the pressurized loss of forced cooling (PLOFC).When results were compared, TRACE and RELAP5-3D showed good agreement in their ability to predict system behavior, but RELAP5-3D calculations were closer to analytical predictions for the RIA. Furthermore, models for a PLOFC accident in two designs (a UO2 and TRISO core) at multiple power levels showed greater temperature margins for the TRISO core at all power levels. Using this information, along with other scoping analyses and constraints, the TCR design team selected a power level of 3 MWth and a TRISO-based core design. For this design RAVEN was applied to vary RIA parameters in RELAP5-3D models at HZP and HFP to understand the effect on figures of merit such as peak power, fuel and coolant temperature, and energy deposition. This sensitivity study found that the inserted reactivity worth was the most important parameter controlling all figures of merit, but for insertion up to 1.5$ no failure of TCR fuel is anticipated. For constant reactivity insertions, the magnitude of the fuel temperature coefficient was found to have the greatest effect on all figures of merit under most circumstances. These results not only demonstrate the anticipated robust safety of the proposed TCR fuel form but also provide a reference for future metal-hydride moderated systems to understand RIA behavior. In the PLOFC, the impact of heat transfer enhancement due to wavy flow channel effects dominated the variance in peak temperatures, and variations in heat exchanger elevation provided the greatest control on natural circulation flow rate. Additionally, no fuel particle failure is anticipated in the PLOFC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Optimization of Binning Parameters in the Feynman Variance-to-Mean Method for Delayed Neutron Reinterrogation [Poster]

The purpose of this work is to determine the optimal input parameters for analyzing active interrogation data that utilizes the physics of the delayed neutron re-interrogation technique. This work compares the signature from an active neutron measurement of a highly enriched uranium (HEU) system to an active background measurement (i.e., without the HEU item in place). In addition, varying analyses of the NoMAD list-mode detector data are performed to determine whether it is preferential to use a particular analysis bin size for producing Feynman histograms, selection of pulses to prevent discrepancies from buildup, or window of time during the pulsed neutron interrogation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Removal of Woody Riparian Vegetation Substantially Altered a Stream Ecosystem in an Otherwise Undisturbed Grassland Watershed

Riparian zones are key interfaces between stream and terrestrial ecosystems. Yet, we know of no whole-watershed experiments that cut only woody vegetation in the riparian zone in an otherwise intact watershed to isolate the role of riparian zones on stream ecology. We removed all of the woody riparian vegetation (from 10- and 30-m-wide buffers in headwaters and main channels, respectively) for 5 km of stream in a single watershed while leaving the remainder of the grassland watershed un-impacted. We assessed water chemistry changes 3 years before and 3 years after riparian wood removal and in two neighboring control watersheds with a before–after, control-impact design and analysis. Riparian woody removal caused 10–100-fold increases in mean stream water nitrate concentrations and pulses of high nitrate for 3 years thereafter. Other nutrients and total suspended solids increased 2–25 times for the 3 years of post-removal. In-stream rates of gross primary production, ecosystem respiration, and net ecosystem production had large treatment effect sizes but also high variance among samples. Past studies of whole-watershed deforestations showed similar water quality responses to our riparian deforestation. Riparian zones of grassland streams are sensitive to disturbance and likely impart relatively greater influence on stream structure and function than the upslope of the watershed. Finally, our results further emphasize the role of riparian zones in biogeochemically linking aquatic and terrestrial habitats.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Statistical Approaches to Aerosol Dynamics for Climate Simulation

In this work, we introduce two general non-parametric regression analysis methods for errors-in-variable (EIV) models: the compound regression, and the constrained regression. It is shown that these approaches are equivalent to each other and, to the general parametric structural modeling approach. The advantages of these methods lie in their intuitive geometric representations, their distribution free nature, and their ability to offer a practical solution when the ratio of the error variances is unknown. Each includes the classic non-parametric regression methods of ordinary least squares, geometric mean regression, and orthogonal regression as special cases. Both methods can be readily generalized to multiple linear regression with two or more random regressors.

Compound regression↗