Search NASA⌕ Search

SEARCH · Search NASA

Results for “Predictors”

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 793 records · Page 44

AI ensemble for signal detection of higher order gravitational wave modes of quasi-circular, spinning, non-precessing binary black hole mergers

We introduce spatiotemporal-graph models that concurrently process data from the twin advanced LIGO detectors and the advanced Virgo detector. We trained these AI classifiers with 2.4 million IMRPhenomXPHM waveforms that describe quasi-circular, spinning, non-precessing binary black hole mergers with component masses m{1,2}∈[3M⊙,50M⊙], and individual spins sz{1,2}∈[−0.9,0.9]; and which include the (ℓ,|m|)={(2,2),(2,1),(3,3),(3,2),(4,4)} modes, and mode mixing effects in the ℓ=3,|m|=2 harmonics. We trained these AI classifiers within 22 hours using distributed training over 96 NVIDIA V100 GPUs in the Summit supercomputer. We then used transfer learning to create AI predictors that estimate the total mass of potential binary black holes identified by all AI classifiers in the ensemble. We used this ensemble, 3 classifiers for signal detection and 2 total mass predictors, to process a year-long test set in which we injected 300,000 signals. This year-long test set was processed within 5.19 minutes using 1024 NVIDIA A100 GPUs in the Polaris supercomputer (for AI inference) and 128 CPU nodes in the ThetaKNL supercomputer (for post-processing of noise triggers), housed at the Argonne Leadership Computing Facility. These studies indicate that our AI ensemble provides state-of-the-art signal detection accuracy, and reports 2 misclassifications for every year of searched data. This is the first AI ensemble designed to search for and find higher order gravitational wave mode signals.

Tian, Minyang↗

Data from: Biomass yield, yield components and growing season phenotypic measurements of Miscanthus

For sustainable biomass production of Miscanthus × giganteus (hereafter miscanthus), understanding the impact of stand age and nitrogen (N) fertilization on biomass yield is crucial. This study investigated the effects of varying N fertilization rates (0, 56, 112, and 168 kg N ha-1) on yield components (tiller height, density, and weight) and their correlations with end-of-season biomass yield in miscanthus. We also explored end-of-season biomass yield prediction using in-season traits (canopy height, leaf area index (LAI), and leaf chlorophyll content (LCC)). The study was conducted at two sites in Illinois: a previously unfertilized 10-year-old miscanthus research stand at Urbana and a 16-year-old commercial stand at Pesotum with a history of annual 56N application. Results from 2018-2021 in Urbana and 2020-2021 in Pesotum showed increased biomass yields with N fertilization, varying by rate, year, and location. Biomass yield in Pesotum peaked at 56N, while in Urbana, it increased significantly at 112 kg N ha-1. Biomass yield was strongly correlated with tiller height and weight measured at Urbana across N rates. Morphological traits measured every 2-3 weeks during the 2020 and 2021 growing seasons showed that canopy height was the strongest single predictor of miscanthus biomass yield, followed by LCC. Mid-August to September measurements of these traits were the best predictors of biomass yield. Multiple regressions involving the canopy height and LCC further improved yield predictions. We conclude that while N enhances biomass yields of aging miscanthus, the optimum rate depends on the site, environmental conditions, and management.

Aging↗

Modeling the contributions of global air temperature, synoptic-scale phenomena and soil moisture to near-surface static energy variability using artificial neural networks

Abstract. The static energy content of the atmosphere is increasing on a global scale, but exhibits important subglobal and subregional scales of variability and is a useful parameter for integrating the net effect of changes in the partitioning of energy at the surface and for improving understanding of the causes of so-called warming holes (i.e., locations with decreasing daily maximum air temperatures ( T ) or increasing trends of lower magnitude than the global mean). Further, measures of the static energy content (herein the equivalent potential temperature, θ e ) are more strongly linked to excess human mortality and morbidity than air temperature alone, and have great relevance in understanding causes of past heat-related excess mortality and making projections of possible future events that are likely to be associated with negative human health and economic consequences. New nonlinear statistical models for summertime daily maximum and minimum θ e are developed and used to advance understanding of drivers of historical change and variability over the eastern USA. The predictor variables are an index of the daily global mean temperature, daily indices of the synoptic-scale meteorology derived from T and specific humidity ( Q ) at 850 and 500 hPa geopotential heights ( Z ), and spatiotemporally averaged soil moisture ( SM ). SM is particularly important in determining the magnitude of θ e over regions that have previously been identified as exhibiting warming holes, confirming the key importance of SM in dictating the partitioning of net radiation into sensible and latent heat and dictating trends in near-surface T and θ e . Consistent with our a priori expectations, models built using artificial neural networks (ANNs) out-perform linear models that do not permit interaction of the predictor variables (global T , synoptic-scale meteorological conditions and SM ). This is particularly marked in regions with high variability in minimum and maximum θ e , where more complex models built using ANN with multiple hidden layers are better able to capture the day-to-day variability in θ e and the occurrence of extreme maximum θ e . Over the entire domain, the ANN with three hidden layers exhibits high accuracy in predicting maximum θ e > 347 K. The median hit rate for maximum θ e >347 K is > 0.60, while the median false alarm rate is ≈ 0.08.

54 ENVIRONMENTAL SCIENCES↗

Machine-learning-based investigation of the variables affecting summertime lightning occurrence over the Southern Great Plains

Lightning is affected by many factors, many of which are not routinely measured, well understood, or accounted for in physical models. Several commonly used machine learning (ML) models have been applied to analyze the relationship between Atmospheric Radiation Measurement (ARM) data and lightning data from the Earth Networks Total Lightning Network (ENTLN) in order to identify important variables affecting lightning occurrence in the vicinity of the Southern Great Plains (SGP) ARM site during the summer months (June, July, August and September) of 2012 to 2020. Testing various ML models, we found that the random forest model is the best predictor among common classifiers. When convective clouds were detected, it predicts lightning occurrence with an accuracy of 76.9 % and an area under the curve (AUC) of 0.850. Using this model, we further ranked the variables in terms of their effectiveness in nowcasting lightning and identified geometric cloud thickness, rain rate and convective available potential energy (CAPE) as the most effective predictors. The contrast in meteorological variables between no-lightning and frequent-lightning periods was examined for hours with CAPE values conducive to thunderstorm formation. Besides the variables considered for the ML models, surface variables and mid-altitude variables (e.g., equivalent potential temperature and minimum equivalent potential temperature, respectively) have statistically significant contrasts between no-lightning and frequent-lightning hours. For example, the minimum equivalent potential temperature from 700 to 500 hPa is significantly lower during frequent-lightning hours compared with no-lightning hours. Finally, a notable positive relationship between the intracloud (IC) flash fraction and the square root of CAPE ($\sqrt{CAPE}$) was found, suggesting that stronger updrafts increase the height of the electrification zone, resulting in fewer flashes reaching the surface and consequently a greater IC flash fraction.

54 ENVIRONMENTAL SCIENCES↗

Drought-influenced mortality of tree species with different predawn leaf water dynamics in a decade-long study of a central US forest

Using decade-long continuous observations of tree mortality and predawn leaf water potential (ψ pd ) at the Missouri Ozark AmeriFlux (MOFLUX) site, we studied how the mortality of important tree species varied and how such variations may be predicted. Water stress determined inter-annual variations in tree mortality with a time delay of 1 year or more, which was correlated fairly tightly with a number of quantitative predictors formulated based on ψ pd and precipitation regimes. Predictors based on temperature and vapor pressure deficit anomalies worked reasonably well, particularly for moderate droughts. The exceptional drought of the year 2012 drastically increased the mortality of all species, including drought-tolerant oaks, in the subsequent year. The drought-influenced tree mortality was related to the species position along the spectrum of ψ pd regulation capacity with those in either ends of the spectrum being associated with elevated risk of death. Regardless of species and drought intensity, the ψ pd of all species recovered rapidly after sufficiently intense rain events in all droughts. This result, together with a lack of immediate leaf and branch desiccation, suggests an absence of catastrophic hydraulic disconnection in the xylem and that tree death was caused by significant but indirect effects. Species differences in the capacity of regulating ψ pd and its temporal integral were magnified under moderate drought intensities but diminished towards wet and dry extremes. Severe droughts may overwhelm the capacity of even drought-tolerant species to maintain differential levels of water potential as the soil becomes exhausted of available water in the rooting zone, thus rendering them more susceptible to death if predisposed by other factors such as age.

54 ENVIRONMENTAL SCIENCES↗

Scaling impacts on environmental controls and spatial heterogeneity of soil organic carbon stocks

The spatial heterogeneity of land surfaces affects energy, moisture, and greenhouse gas exchanges with the atmosphere. However, representing the heterogeneity of terrestrial hydrological and biogeochemical processes in Earth system models (ESMs) remains a critical scientific challenge. We report the impact of spatial scaling on environmental controls, spatial structure, and statistical properties of soil organic carbon (SOC) stocks across the US state of Alaska. We used soil profile observations and environmental factors such as topography, climate, land cover types, and surficial geology to predict the SOC stocks at a 50 m spatial scale. These spatially heterogeneous estimates provide a data set with reasonable fidelity to the observations at a sufficiently high resolution to examine the environmental controls on the spatial structure of SOC stocks. We upscaled both the predicted SOC stocks and environmental variables from finer to coarser spatial scales ( s = 100, 200, and 500 m and 1, 2, 5, and 10 km) and generated various statistical properties of SOC stock estimates. We found different environmental factors to be statistically significant predictors at different spatial scales. Only elevation, temperature, potential evapotranspiration, and scrub land cover types were significant predictors at all scales. The strengths of control (the median value of geographically weighted regression coefficients) of these four environmental variables on SOC stocks decreased with increasing scale and were accurately represented using mathematical functions ( R 2 = 0.83–0.97). The spatial structure of SOC stocks across Alaska changed with spatial scale. Although the variance (sill) and unstructured variability (nugget) of the calculated variograms of SOC stocks decreased exponentially with scale, the correlation length (range) remained relatively constant across scale. The variance of predicted SOC stocks decreased with spatial scale over the range of 50 m to ~ 500 m, and remained constant beyond this scale. The fitted exponential function accounted for 98 % of variability in the variance of SOC stocks. We found moderately accurate linear relationships between mean and higher-order moments of predicted SOC stocks ( R 2 ∼ 0.55–0.63). Current ESMs operate at coarse spatial scales (50–100 km), and are therefore unable to represent environmental controllers and spatial heterogeneity of high-latitude SOC stocks consistent with observations. We conclude that improved understanding of the scaling behavior of environmental controls and statistical properties of SOC stocks could improve ESM land model benchmarking and perhaps allow representation of spatial heterogeneity of biogeochemistry at scales finer than those currently resolved by ESMs.

54 ENVIRONMENTAL SCIENCES↗

Shifting mineral and redox controls on carbon cycling in seasonally flooded mineral soils

Although wetland soils represent a relatively small portion of the terrestrial landscape, they account for an estimated 20%–30% of the global soil carbon (C) reservoir. C stored in wetland soils that experience seasonal flooding is likely the most vulnerable to increased severity and duration of droughts in response to climate change. Redox conditions, plant root dynamics, and the abundance of protective mineral phases are well-established controls on soil C persistence, but their relative influence in seasonally flooded mineral soils is largely unknown. To address this knowledge gap, we assessed the relative importance of environmental (temperature, soil moisture, and redox potential) and biogeochemical (mineral composition and root biomass) factors in controlling CO 2 efflux, C quantity, and organic matter composition along replicated upland–lowland transitions in seasonally flooded mineral soils. Specifically, we contrasted mineral soils under temperature deciduous forests in lowland positions that undergo seasonal flooding with adjacent upland soils that do not, considering both surface (A) and subsurface (B and C) horizons. We found the lowland soils had lower total annual CO 2 efflux than the upland soils, with monthly CO 2 efflux in lowlands most strongly correlated with redox potential (E h ). Lower CO 2 efflux as compared to the uplands corresponded to greater C content and abundance of lignin-rich, higher-molecular-weight, chemically reduced organic compounds in the lowland surface soils (A horizons). In contrast, subsurface soils in the lowland position (C g horizons) showed lower C content than the upland positions (C horizons), coinciding with lower abundance of root biomass and oxalate-extractable Fe (Fe o , a proxy for protective Fe phases). Our linear mixed-effects model showed that Fe o served as the strongest measured predictor of C content in upland soils, yet Fe o had no predictive power in lowland soils. Instead, our model showed that E h and oxalate-extractable Al (Al o , a proxy of protective Al phases) became significantly stronger predictors in the lowland soils. Combined, our results suggest that low redox potentials are the primary cause for C accumulation in seasonally flooded surface soils, likely due to selective preservation of organic compounds under anaerobic conditions. In seasonally flooded subsurface soils, however, C accumulation is limited due to lower C inputs through root biomass and the removal of reactive Fe phases under reducing conditions. Our findings demonstrate that C accrual in seasonally flooded mineral soil is primarily due to low redox potential in the surface soil and that the lack of protective metal phases leaves these C stocks highly vulnerable to climate change.

58 GEOSCIENCES↗

Observational benchmarks inform representation of soil organic carbon dynamics in land surface models

Abstract. Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon–climate feedbacks. Machine learning models can help identify dominant environmental controllers and establish their functional relationships with SOC stocks. The resulting knowledge can be integrated into ESMs to reduce uncertainty and improve predictions of SOC dynamics over space and time. In this study, we used a large number of SOC field observations (n=54 000), geospatial datasets of environmental factors (n=46), and two machine learning approaches (namely random forest, RF, and generalized additive modeling, GAM) to (1) identify dominant environmental controllers of global and biome-specific SOC stocks, (2) derive functional relationships between environmental controllers and SOC stocks, and (3) compare the identified environmental controllers and predictive relationships with those in models used in Phase 6 of the Coupled Model Intercomparison Project (CMIP6). Our results showed that the diurnal temperature, drought index, cation exchange capacity, and precipitation were important observed environmental predictors of global SOC stocks. While the RF model identified 14 environmental factors that describe climatic, vegetation, and edaphic conditions as important predictors of global SOC stocks (R2=0.61, RMSE = 0.46 kg m−2), current ESMs oversimplify the relationships between environmental factors and SOC, with precipitation, temperature, and net primary productivity explaining > 96 % of the variability in ESM-modeled SOC stocks. Further, our study revealed notable disparities among the functional relationships between environmental factors and SOC stocks simulated by ESMs compared with observed relationships. To improve SOC representations in ESMs, it is imperative to incorporate additional environmental controls, such as the cation exchange capacity, and refine the functional relationships to align more closely with observations.

54 ENVIRONMENTAL SCIENCES↗

Scaling impacts on environmental controls and spatial heterogeneity of soil organic carbon stocks

The spatial heterogeneity of land surfaces affects energy, moisture, and greenhouse gas exchanges with the atmosphere. However, representing heterogeneity of terrestrial hydrological and biogeochemical processes in earth system models (ESMs) remains a critical scientific challenge. We report the impact of spatial scaling on environmental controls, spatial structure, and statistical properties of soil organic carbon (SOC) stocks across the US state of Alaska. We used soil profile observations and environmental factors such as topography, climate, land cover types, and surficial geology to predict the SOC stocks at a 50 m spatial scale. These spatially heterogeneous estimates provide a dataset with reasonable fidelity to the observations at a sufficiently high resolution to examine the environmental controls on the spatial structure of SOC stocks. We upscaled both the predicted SOC stocks and environmental variables from finer to coarser spatial scales ( s = 100, 200, 500 m, 1, 2, 5, 10 km) and generated various statistical properties of SOC stock estimates. We found different environmental factors to be statistically significant predictors at different spatial scales. Only elevation, temperature, potential evapotranspiration, and scrub land cover types were significant predictors at all scales. The strengths of control (the median value of geographically weighted regression coefficients) of these four environmental variables on SOC stocks decreased with increasing scale and were accurately represented using mathematical functions ( R 2 = 0.83–0.97). The spatial structure of SOC stocks across Alaska changed with spatial scale. Although the variance (sill) and unstructured variability (nugget) of the calculated variograms of SOC stocks decreased exponentially with scale, the correlation length (range) remained relatively constant across scale. The variance of predicted SOC stocks decreased with spatial scale over the range of 50 to ~ 500 m, and remained constant beyond this scale. The fitted exponential function accounted for 98% of variability in the variance of SOC stocks. We found moderately-accurate linear relationships between mean and higher-order moments of predicted SOC stocks ( R 2 ~ 0.55–0.63). Current ESMs operate at coarse spatial scales (50–100 km), and are therefore unable to represent environmental controllers and spatial heterogeneity of high-latitude SOC stocks consistent with observations. We conclude that improved understanding of the scaling behavior of environmental controls and statistical properties of SOC stocks can improve ESM land model benchmarking and perhaps allow representation of spatial heterogeneity of biogeochemistry at scales finer than those currently resolved by ESMs.

54 ENVIRONMENTAL SCIENCES↗

Upscaled diurnal cycles of land–atmosphere fluxes: a new global half-hourly data product

Interactions between the biosphere and the atmosphere can be well characterized by fluxes between the two. In particular, carbon and energy fluxes play a major role in understanding biogeochemical processes on an ecosystem level or global scale. However, the fluxes can only be measured at individual sites, e.g., by eddy covariance towers, and an upscaling of these local observations is required to analyze global patterns. Previous work focused on upscaling monthly, 8-day, or daily average values, and global maps for each flux have been provided accordingly. In this paper, we raise the upscaling of carbon and energy fluxes between land and atmosphere to the next level by increasing the temporal resolution to subdaily timescales. We provide continuous half-hourly fluxes for the period from 2001 to 2014 at 0.5° spatial resolution, which allows for analyzing diurnal cycles globally. The data set contains four fluxes: gross primary production (GPP), net ecosystem exchange (NEE), latent heat (LE), and sensible heat (H). We propose two prediction approaches for the diurnal cycles based on large-scale regression models and compare them in extensive cross-validation experiments using different sets of predictor variables. We analyze the results for a set of FLUXNET tower sites showing the suitability of our approaches for this upscaling task. Finally, we have selected one approach to calculate the global half-hourly data products based on predictor variables from remote sensing and meteorology at daily resolution as well as half-hourly potential radiation. In addition, we provide a derived product that only contains monthly average diurnal cycles, which is a lightweight version in terms of data storage that still allows studying the important characteristics of diurnal patterns globally. We recommend to primarily use these monthly average diurnal cycles, because they are less affected by the impacts of day-to-day variation, observation noise, and short-term fluctuations on subdaily timescales compared to the full half-hourly flux products. The global half-hourly data products are available at https://doi.org/10.17871/BACI.224 .

54 ENVIRONMENTAL SCIENCES↗

Modeling the topographic influence on aboveground biomass using a coupled model of hillslope hydrology and ecosystem dynamics

Abstract. Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling.

54 ENVIRONMENTAL SCIENCES↗

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗

Low flow characteristics for regulated and unregulated streams in North Carolina and prediction using climate signals

In this work, low flow statistics of regulated and unregulated streams in the state of North Carolina were updated with streamflow data through 2019. About 22% of the streams considered show a significant downward trend, but considerable low-frequency variability confounding trends. The relationship between regional index time series of groundwater-depth low flow and the Atlantic Multidecadal Oscillation (AMO) was examined to determine whether low-frequency climate modes can account for the long-term pattern in low flows. Consequently, a significant correlation was found between AMO and groundwater-depth low flow, such that positive AMO is associated with lower groundwater-depth low flow and vice versa, particularly for the Piedmont region. Predictive equations for annual low flows at the ecoregion level shows that springtime average streamflow and AMO were selected as the primary predictors of low flow for coastal and Piedmont regions, whereas springtime average streamflow and the November–December–January average Oceanic Niño Index were used as the primary predictors for the mountain region. The relative root mean square error (RMSE) of the disaggregated predictions to US Geologic Survey gauge locations was <23% at 79% of the stations, between 24% and 43% at 10% of the stations, and greater than 44% at 1% of the stations. The remaining 10% of stations showed large RMSEs. This latter percentage is characterized by smaller drainage basins and intermittent flows, suggesting the prediction models are not applicable to drainage basins smaller than roughly 20 square km and intermittent streams.

54 ENVIRONMENTAL SCIENCES↗

Local and latitudinal variation in abundance: the mechanisms shaping the distribution of an ecosystem engineer

Ecological processes that determine the abundance of species within ecological communities vary across space and time. These scale-dependent processes are especially important when they affect key members of a community, such as ecosystem engineers that create shelter and food resources for other species. Yet, few studies have examined the suite of processes that shape the abundance of ecosystem engineers. Here, we evaluated the relative influence of temporal variation, local processes, and latitude on the abundance of an engineering insect—a rosette-galling midge, Rhopalomyia solidaginis (Diptera: Cecidomyiidae). Over a period of 3–5 years, we studied the density and size of galls across a suite of local experiments that manipulated genetic variation, soil nutrient availability, and the removal of other insects from the host plant, Solidago altissima (tall goldenrod). We also surveyed gall density within a single growing season across a 2,300 km latitudinal transect of goldenrod populations in the eastern United States. At the local scale, we found that host-plant genotypic variation was the best predictor of rosette gall density and size within a single year. We found that the removal of other insect herbivores resulted in an increase in gall density and size. The amendment of soil nutrients for four years had no effect on gall density, but galls were smaller in carbon-added plots compared to control and nitrogen additions. Finally, we observed that gall density varied several fold across years. At the biogeographic scale, we observed that the density of rosette gallers peaked at mid-latitudes. Using meta-analytic approaches, we found that the effect size of time, followed by host-plant genetic variation and latitude were the best predictors of gall density. Taken together, our study provides a unique comparison of multiple factors across different spatial and temporal scales that govern engineering insect herbivore density.

54 ENVIRONMENTAL SCIENCES↗

The Adams formulas for numerical integration of differential equations from 1st to 20th order

The Adams Bashforth predictor coefficients and the Adams Moulton corrector coefficients for the integration of differential equations are presented for methods of 1st to 20th order. The order of the method as presented refers to the highest order difference formula used in Newton's backward difference interpolation formula, on which the Adams method is based. The Adams method is a polynomial approximation method derived from Newton's backward difference interpolation formula. The Newton formula is derived and expanded to 20th order. The Adams predictor and corrector formulas are derived and expressed in terms of differences of the derivatives, as well as in terms of the derivatives themselves. All coefficients are given to 18 significant digits. For the difference formula only, the ratio coefficients are given to 10th order.

Kirkpatrick, J. C.↗

Experimental and analytical investigation of a nonaxisymmetric wedge nozzle at static conditions

An experimental investigation of a nonaxisymmetric wedge nozzle was conducted at static conditions. The resulting data, in the form of detailed pressure distributions and oil flow photographs, expand the current nonaxisymmetric nozzle data base. An analytical investigation has been conducted to evaluate a two-dimensional, inviscid, time-dependent theory as a nonaxisymmetric nozzle performance predictor. For the range of nozzle pressure ratios investigated, results indicate good agreement between theory and experiment in regions of predominately two-dimensional flow and limited agreement in regions of three-dimensional flow. For the wedge nozzle and related nozzle configurations, the two dimensional, inviscid theory may be applied as a limited performance predictor.

Carson, G. T., Jr.↗

Updating of ADAPT predictions of sunspot activity

An eigenvector analysis procedure was used to analyze and develop algorithms for predicting sunspot numbers. The predictors in these algorithms consist of sunspot numbers from the preceding two solar cycles and magnetic index data from the preceding cycle. Predictions are presented for cycles 21 and 22. The sunspot activity for cycle 21 is predicted to remain below 100 until early 1980 when it will rapidly reach a peak of approximately 120. The two sigma accuracy on these estimates is approximately 20 sunspot numbers in the region of the peak and 10 sunspot numbers early and late in the cycle. Algorithms were also developed for predicting the period of future sunspot cycles using the same predictor vector.

Hunter, H. C.↗

Estimation of the accuracy of dynamic flight-determined coefficients

This paper discusses means of assessing the accuracy of maximum likelihood parameter estimates obtained from dynamic flight data. The commonly used analytical predictors of accuracy are compared from both statistical and simplified geometric standpoints. Emphasizing practical considerations, such as modeling error, the accuracy predictions are evaluated with real and simulated data. Improved computations of the Cramer-Rao bound to correct large discrepancies caused by colored noise and modeling error are presented. This corrected Cramer-Rao bound is the best available analytical predictor of accuracy. Engineering judgement, aided by such analytical tools, is the final arbiter of accuracy estimation.

Maine, R. E.↗