Temporal and spatial variation of nitrogen transformations in nitrogen-saturated soils of a central Appalachian hardwood forest

We studied temporal and spatial patterns of soil nitrogen (N) dynamics from 1993 to 1995 in three watersheds of Fernow Experimental Forest, W.V.: WS7 (24-year-old, untreated); WS4 (mature, untreated); and WS3 (24-year-old, treated with (NH4)2SO since 1989 at the rate of 35 kg Nha–1year–1). Net nitrification was 141, 114, and115 kg Nha–1year–1, for WS3, WS4, and WS7, respectively, essentially 100% of net N mineralization for all watersheds. Temporal (seasonal) patterns of nitrification were significantly related to soil moisture and ambient temperaturein untreated watersheds only. Spatial patterns of soil water NO3–of WS4 suggest that microenvironmental variabilitylimits rates of N processing in some areas of this N-saturated watershed, in part by ericaceous species in the herbaceous layer. Spatial patterns of soil water NO3–in treated WS3 suggest that later stages of N saturation may result inhigher concentrations with less spatial variability. Spatial variability in soil N variables was lower in treated WS3 versus untreated watersheds. Nitrogen additions have altered the response of N-processing microbes to environmental factors, becoming less sensitive to seasonal changes in soil moisture and temperature. Biotic processes responsible forregulating N dynamics may be compromised in N-saturated forest ecosystems.

Factors Influencing Spatial Variability in Nitrogen Processing in Nitrogen-Saturated Soils

Nitrogen (N) saturation is an environmental concern for forests in the eastern U.S. Although several watersheds of the Fernow Experimental Forest (FEF), West Virginia exhibit symptoms of Nsaturation, many watersheds display a high degree of spatial variability in soil N processing. This study examined the effects of temperature on net N mineralization and nitrification in N-saturatedsoils from FEF, and how these effects varied between high N-processing vs. low N-processingsoils collected from two watersheds, WS3 (fertilized with [NH4]2SO4) and WS4 (untreated control). Samples of forest floor material (O2 horizon) and mineral soil (to a 5-cm depth) were taken from three subplots within each of four plots that represented the extremes of highest and lowest ratesof net N mineralization and nitrification (hereafter, high N and low N, respectively) of untreated WS4 and N-treated WS3: control/low N, control/high N, N-treated/low N, N-treated/high N. Forest floor material was analyzed for carbon (C), lignin,and N. Subsamples of mineral soil were extractedimmediately with 1 N KCl and analyzed for NH4+and NO3– to determine preincubation levels. Extracts were also analyzed for Mg, Ca, Al, and pH. To test the hypothesis that the lack of net nitrification observed in field incubations on the untreated/low N plot was the result of absence ofnitrifier populations, we characterized the bacterial community involved in N cycling by amplification of amoA genes. Remaining soil was incubated for 28 d at three temperatures (10, 20, and30°C), followed by 1 N KCl extraction and analysis for NH4+ and NO3–. Net nitrification was essentially 100% of net N mineralization for all samples combined. Nitrification rates from lab incubation sat all temperatures supported earlier observations based on field incubations. At 30°C, rates from N- t reated/high N were three times those of N-treated/low N. Highest rates were found for untreated/high N (two times greater than those of N-treated/high N), whereas untreated/low N exhibited no net nitrification. However, soils exhibitingno net nitrification tested positive for presence of nitrifying bacteria, causing us to reject our initial hypothesis. We hypothesize that nitrifier populations in such soil are being inhibited by a combination of low Ca:Al ratios in mineral soil and allelopathic interactions with mycorrhizae of ericaceous species in the herbaceous layer.

Spatial heterogeneity of soil nutrients and plant species in herb-dominated communities of contrasting land use

Recent interest in spatial pattern in terrestrial ecosystems has come from an awareness of theintimate relationship between spatial heterogeneity of soil resources and maintenance of plant species diversity. Soil and vegetation can vary spatially inresponse to several state factors of the system. In this study, we examined fine-scale spatial variability of soil nutrients and vascular plant species in contrasting herb-dominated communities (a pasture and an oldfield) to determine degree of spatial dependenceamong soil variables and plant community characteristics within these communities by sampling at 1-m intervals. Each site was divided into 25 1-m 2 plots. Mineral soil was sampled (2-cm diameter, 5-cm depth) from each of four 0.25-m2 quarters and combined into a single composite sample per plot. Soil organic matter was measured as loss-on-ignition. Extractable NH4 and NO3 were determined before and after laboratory incubation to determine potential net N mineralization and nitrification. Cations were analyzed using inductively coupled plasma emission spectrometry. Vegetation was assessed using estimated percent cover. Most soiland plant variables exhibited sharp contrasts betweenpasture and old-field sites, with the old field having significantly higher net N mineralization/nitrification, pH, Ca, Mg, Al, plant cover, and species diversity, richness, and evenness. Multiple regressions revealedthat all plant variables (species diversity, richness,evenness, and cover) were significantly related to soil characteristics (available nitrogen, organic matter,moisture, pH, Ca, and Mg) in the pasture; in the old field only cover was significantly related to soil characteristics (organic matter and moisture). Both sites contrasted sharply with respect to spatial pattern of soil variables, with the old field exhibiting a higher degree of spatial dependence. These results demonstrate that land-use practices can exert profound influence on spatial heterogeneity of both soil properties and vegetation in herb-dominated communities.

The Importance of Scale for Spatial-confounding Bias and Precision of Spatial Regression Estimators

Increasingly, regression models are used when residuals are spatially correlated. Prominent examples include studies in environmental epidemiology to understand the chronic health effects of pollutants. I consider the effects of residual spatial structure on the bias and precision of regression coefficients, developing a simple framework in which to understand the key issues and derive informative analytic results. When the spatial residual is induced by an unmeasured confounder, regression models with spatial random effects and closely-related models such as kriging and penalized splines are biased, even when the residual variance components are known. Analytic and simulation results show how the bias depends on the spatial scales of the covariate and the residual; bias is reduced only when there is variation in the covariate at a scale smaller than the scale of the unmeasured confounding. I also discuss how the scales of the residual and the covariate affect efficiency and uncertainty estimation when the residuals can be considered independent of the covariate. In an application on the association between black carbon particulate matter air pollution and birth weight, controlling for large-scale spatial variation appears to reduce bias from unmeasured confounders, while increasing uncertainty in the estimated pollution effect.

Computational Techniques for Spatial Logistic Regression with Large Datasets

In epidemiological work, outcomes are frequently non-normal, sample sizes may be large, and effects are often small. To relate health outcomes to geographic risk factors, fast and powerful methods for fitting spatial models, particularly for non-normal data, are required. We focus on binary outcomes, with the risk surface a smooth function of space. We compare penalized likelihood models, including the penalized quasi-likelihood (PQL) approach, and Bayesian models based on fit, speed, and ease of implementation. A Bayesian model using a spectral basis representation of the spatial surface provides the best tradeoff of sensitivity and specificity in simulations, detecting real spatial features while limiting overfitting and being more efficient computationally than other Bayesian approaches. One of the contributions of this work is further development of this underused representation. The spectral basis model outperforms the penalized likelihood methods, which are prone to overfitting, but is slower to fit and not as easily implemented. Conclusions based on a real dataset of cancer cases in Taiwan are similar albeit less conclusive with respect to comparing the approaches. The success of the spectral basis with binary data and similar results with count data suggest that it may be generally useful in spatial models and more complicated hierarchical models.

Limitations of Remotely-sensed Aerosol as a Spatial Proxy for Fine Particulate Matter

Recent research highlights the promise of remotely-sensed aerosol optical depth (AOD) as a proxy for ground-level PM2.5. Particular interest lies in the information on spatial heterogeneity potentially provided by AOD, with important application to estimating and monitoring pollution exposure for public health purposes. Given the temporal and spatio-temporal correlations reported between AOD and PM2.5 , it is tempting to interpret the spatial patterns in AOD as reflecting patterns in PM2.5 . Here we find only limited spatial associations of AOD from three satellite retrievals with PM2.5 over the eastern U.S. at the daily and yearly levels in 2004. We then use statistical modeling to show that the patterns in monthly average AOD poorly reflect patterns in PM2.5 because of systematic, spatially-correlated error in AOD as a proxy for PM2.5 . Furthermore, when we include AOD as a predictor of monthly PM2.5 in a statistical prediction model, AOD provides little additional information to improve predictions of PM2.5 when included in a model that already accounts for land use, emission sources, meteorology and regional variability. These results suggest caution in using spatial variation in AOD to stand in for spatial variation in ground-level PM2.5 in epidemiological analyses and indicate that when PM2.5 monitoring is available, careful statistical modeling outperforms the use of AOD.

A BAYESIAN HIERARCHICAL FRAMEWORK FOR SPATIAL MODELING OF fMRI DATA

Functional neuroimaging techniques enable investigations into the neural basis of human cognition, emotions, and behaviors. In practice, applications of functional magnetic resonance imaging (fMRI) have provided novel insights into the neuropathophysiology of major psychiatric,neurological, and substance abuse disorders, as well as into the neural responses to their treatments. Modern activation studies often compare localized task-induced changes in brain activity between experimental groups. One may also extend voxel-level analyses by simultaneously considering the ensemble of voxels constituting an anatomically defined region of interest (ROI) or by considering means or quantiles of the ROI. In this work we present a Bayesian extension of voxel-level analyses that offers several notable benefits. First, it combines whole-brain voxel-by-voxel modeling and ROI analyses within a unified framework. Secondly, an unstructured variance/covariance for regional mean parameters allows for the study of inter-regional functional connectivity, provided enough subjects are available to allow for accurate estimation. Finally, an exchangeable correlation structure within regions allows for the consideration of intra-regional functional connectivity. We perform estimation for our model using Markov Chain Monte Carlo (MCMC) techniques implemented via Gibbs sampling which, despite the high throughput nature of the data, can be executed quickly (less than 30 minutes). We apply our Bayesian hierarchical model to two novel fMRI data sets: one considering inhibitory control in cocaine-dependent men and the second considering verbal memory in subjects at high risk for Alzheimer’s disease. The unifying hierarchical model presented in this manuscript is shown to enhance the interpretation content of these data sets.

BIVARIATE BINOMIAL SPATIAL MODELLING LOA loa PREVALENCE IN TROPICAL AFRICA

We present a state-of-the-art application of smoothing for dependent bivariate binomial spatial data to Loa loa prevalence mapping in West Africa. This application is special because it starts with the non-spatial calibration of survey instruments, continues with the spatial model building and assessment and ends with robust, tested software that will be used by the field scientists of the World Health Organization for online prevalence map updating. From a statistical perspective several important methodological issues were addressed: (a) building spatial models that are complex enough to capture the structure of the data but remain computationally usable; (b)reducing the computational burden in the handling of very large covariate data sets; (c) devising methods for comparing spatial prediction methods for a given exceedance policy threshold.