Regression analysis of spatial data

Many of the most interesting questions ecologists ask lead to analyses of spatial data. Yet, perhaps confused by the large number of statistical models and fitting methods available, many ecologists seem to believe this is best left to specialists. Here, we describe the issues that need consideration when analysing spatial data and illustrate these using simulation studies. Our comparative analysis involves using methods including generalized least squares, spatial filters, wavelet revised models, conditional autoregressive models and generalized additive mixed models to estimate regression coefficients from synthetic but realistic data sets, including some which violate standard regression assumptions. We assess the performance of each method using two measures and using statistical error rates for model selection. Methods that performed well included generalized least squares family of models and a Bayesian implementation of the conditional auto-regressive model. Ordinary least squares also performed adequately in the absence of model selection, but had poorly controlled Type I error rates and so did not show the improvements in performance under model selection when using the above methods. Removing large-scale spatial trends in the response led to poor performance. These are empirical results; hence extrapolation of these findings to other situations should be performed cautiously. Nevertheless, our simulation-based approach provides much stronger evidence for comparative analysis than assessments based on single or small numbers of data sets, and should be considered a necessary foundation for statements of this type in future.

Spatial and Temporal Changes in Soil Organic Carbon Assessment using an Integrated Geostatistical and Compositional Data Analysis Approach

The Irish and UK governments, along with other countries, have made a commitment to limit the concentrations of greenhouse gases in the atmosphere by reducing emissions from the burning of fossil fuels. This can be achieved (in part) through increasing the sequestration of CO2 from the atmosphere including monitoring the amount stored in vegetation and soils. A large proportion of soil carbon is held within peat due to the relatively high carbon density of peat and organic-rich soils. This is particularly important for a country such as Ireland, where some 16% of the land surface is covered by peat. For Northern Ireland, it has been estimated that the total amount of carbon stored in vegetation is 4.4Mt compared to 386Mt stored within peat and soils. As a result it has become increasingly important to measure and monitor changes in stores of carbon in soils. The conservation and restoration of peat covered areas, although ongoing for many years, has become increasingly important. This is summed up in current EU policy outlined by the European Commission (2012) which seeks to assess the relative contributions of the different inputs and outputs of organic carbon and organic matter to and from soil. Results are presented from the EU-funded Tellus Border Soil Carbon Project (2011 to 2013) which aimed to improve current estimates of carbon in soil and peat across Northern Ireland and the bordering counties of the Republic of Ireland.
Historical reports and previous surveys provide baseline data. To monitor change in peat depth and soil organic carbon, these historical data are integrated with more recently acquired airborne geophysical (radiometric) data and ground-based geochemical data generated by two surveys, the Tellus Project (2004-2007: covering Northern Ireland) and the EU-funded Tellus Border project (2011-2013) covering the six bordering counties of the Republic of Ireland, Donegal, Sligo, Leitrim, Cavan, Monaghan and Louth. The concept being applied is that saturated organic-rich soil and peat attenuate gamma-radiation from underlying soils and rocks. This research uses the degree of spatial correlation (coregionalization) between peat depth, soil organic carbon (SOC) and the attenuation of the radiometric signal to update a limited sampling regime of ground-based measurements with remotely acquired data. To comply with the compositional nature of the SOC data (perturbations of loss on ignition [LOI] data), a compositional data analysis approach is investigated. Contemporaneous ground-based measurements allow corroboration for the updated mapped outputs. This provides a methodology that can be used to improve estimates of soil carbon with minimal impact to sensitive habitats (like peat bogs), but with maximum output of data and knowledge.

Spatial allocation of material flow analysis in residential developments: a case study of Kildare County, Ireland

Studies of urban metabolism provide important insights for environmental management of cities, but are not widely used in planning practice due to a mismatch of data scale and coverage. This paper introduces the Spatial Allocation of Material Flow Analysis (SAMFA) model as a potential decision support tool aimed as a contribution to overcome some of these difficulties and describes its pilot use at the county level in the Republic of Ireland. The results suggest that SAMFA is capable of identifying hotspots of higher material and energy use to support targeted planning initiatives, while its ability to visualise different policy scenarios supports more effective multi-stakeholder engagement. The paper evaluates this pilot use and sets out how this model can act as an analytical platform for the industrial ecology–spatial planning nexus.

Sector analysis provides additional spatial information on the permanent prostate brachytherapy learning curve

PURPOSE: To evaluate the permanent prostate brachytherapy (PPB) learning curve using postimplant multisector dosimetric analysis and to assess the correlation between sector -specific dosimetry and patient-reported outcome measures (PROMs).

METHODS AND METHODS: First 200 patients treated with (125)I PPB monotherapy (145 Gy) at a single institution were assessed. Postimplant dosimetry (PID) using CT was evaluated for whole prostate (global) and 12 sectors, assessing minimum dose to 90% of prostate (D90) and dose to 0.1 cm(3) of rectum (D0.1cc). Global and sector PID results were evaluated to investigate changes in D90 with case number. Urinary and bowel PROMs were assessed using the International Prostate Symptom Score and the Expanded Prostate Cancer Index Composite questionnaire. The correlation between global and individual sector PID and urinary/bowel PROMs was also evaluated.

RESULTS: Linear regression confirmed a significant improvement in global D90 with case number (r(2) = 0.20; p = 0.001) at a rate of 0.11 Gy/case. Postimplant D90 of base sectors increased at a rate of 0.11-0.15 Gy/case (p = 0.0001) and matched global improvement. The regression lines of midgland and apex sectors were significantly different from global D90 (p = 0.01). Posterior midgland sectors showed a significant reduction in D90 with case number at a rate of 0.13-0.19 Gy/case (p = 0.01). Dose to posterior midgland sectors correlated with rectal D0.1cc dose but not bowel PROMs. Dose to posterior midgland sectors correlated with urinary International Prostate Symptom Score change, which was not apparent when global D90 alone was considered.

CONCLUSIONS: Sector analysis provided increased spatial information regarding the PPB learning curve. Furthermore, sector analysis correlated with urinary PROMs and rectal dose.

Design and Qualitative/Quantitative Analysis of Multi-Agent Spatial Simulation Library

Thesis (Master's)--University of Washington, 2012

Analysis of temporal and spatial variations of forest. A case of study in notheastern Armenia

The forest has a crucial ecological role and the continuous forest loss can cause colossal effects on the environment. As Armenia is one of the low forest covered countries in the world, this problem is more critical. Continuous forest disturbances mainly caused by illegal logging started from the early 1990s had a huge damage on the forest ecosystem by decreasing the forest productivity and making more areas vulnerable to erosion. Another aspect of the Armenian forest is the lack of continuous monitoring and absence of accurate estimation of the level of cuts in some years. In order to have insight about the forest and the disturbances in the long period of time we used Landsat TM/ETM + images. Google Earth Engine JavaScript API was used, which is an online tool enabling the access and analysis of a great amount of satellite imagery. To overcome the data availability problem caused by the gap in the Landsat series in 1988- 1998, extensive cloud cover in the study area and the missing scan lines, we used pixel based compositing for the temporal window of leaf on vegetation (June-late September). Subsequently, pixel based linear regression analyses were performed. Vegetation indices derived from the 10 biannual composites for the years 1984-2014 were used for trend analysis. In order to derive the disturbances only in forests, forest cover layer was aggregated and the original composites were masked. It has been found, that around 23% of forests were disturbed during the study period.

Spatial and temporal analysis of laser induced plasma from a polymer sample

Spatial and temporal analyses of the spectra of the laser induced plasma from a polytetrafluroethylene (PTFE) target obtained with the 1.06 mu m radiation from a Q-switched Nd:YAG laser have been carried out. The spatially resolved spectra of the plasma emission show that molecular bands of C2 (Swan bands) and CN are very intense in the outer regions of the plasma, whereas higher ionized states of carbon are predominant in the core region of the plasma emission. The vibrational temperature and population distribution in the different vibrational levels have been studied as a function of laser energy. From the time resolved studies, it has been observed that there exist fairly large time delays for the onset of emission from all the species in the outer region of the plasma. The molecular bands in each region exhibit much larger time delays in comparison to the ionic lines in the plasma.