Shape selection in Landsat time series: a tool for monitoring forest dynamics.

We present a new methodology for fitting nonparametric shape-restricted regression splines to time series of Landsat imagery for the purpose of modeling, mapping, and monitoring annual forest disturbance dynamics over nearly three decades. For each pixel and spectral band or index of choice in temporal Landsat data, our method delivers a smoothed rendition of the trajectory constrained to behave in an ecologically sensible manner, reflecting one of seven possible 'shapes'. It also provides parameters summarizing the patterns of each change including year of onset, duration, magnitude, and pre- and postchange rates of growth or recovery. Through a case study featuring fire, harvest, and bark beetle outbreak, we illustrate how resultant fitted values and parameters can be fed into empirical models to map disturbance causal agent and tree canopy cover changes coincident with disturbance events through time. We provide our code in the r package ShapeSelectForest on the Comprehensive R Archival Network and describe our computational approaches for running the method over large geographic areas. We also discuss how this methodology is currently being used for forest disturbance and attribute mapping across the conterminous United States.

Surveillance of Rocky Mountain wood ticks (Dermacentor andersoni) and American dog ticks (Dermacentor variabilis) in Colorado.

Ticks pose an emerging threat of infectious pathogen transmission in the United States in part due to expanding suitable habitat ranges in the wake of climate change. Active and passive tick surveillance can inform maps of tick distributions to warn the public of their risk of exposure to ticks. In Colorado, widespread active surveillance programs have difficulty due to the state's diverse terrain. However, combining multiple citizen science techniques can create a more accurate representation of tick distribution than any passive surveillance dataset alone. Our study uses county-level tick distribution data from Northern Arizona University, the Colorado Department of Public Health and the Environment, and veterinary surveillance in addition to literature data to assess the distribution of the Rocky Mountain wood tick, Dermacentor andersoni, and the American dog tick, Dermacentor variabilis. We found that D. andersoni for the most part inhabits counties at higher elevations than D. variabilis in Colorado.

Utilizing citizen science to model the distribution of Aedes aegypti in West Africa.

In the rapidly urbanizing region of West Africa, Aedes mosquitoes pose an emerging threat of infectious disease that is compounded by limited vector surveillance. Citizen science has been proposed as a way to fill surveillance gaps by training local residents to collect and share information on disease vectors. Understanding the distribution of arbovirus vectors in West Africa can inform researchers and public health officials on where to conduct disease surveillance and focus public health interventions. We utilized citizen science data collected through NASA's GLOBE Observer mobile phone application and data from a previously published literature review on Aedes mosquito distribution to examine the contribution of citizen science to understanding the distribution of Ae. aegypti in West Africa using Maximum Entropy modeling. Combining citizen science and literature-derived observations improved the fit of the model compared to models created by each data source alone but did not alleviate location bias within the models, likely due to lack of widespread observations. Understanding Ae. aegypti distribution will require greater investment in Aedes mosquito surveillance in the region, and citizen science should be utilized as a tool in this mission to increase the reach of surveillance.

A sample design for globally consistent biomass estimation using lidar data from the Geoscience Laser Altimeter System (GLAS).

BACKGROUND: Lidar height data collected by the Geosciences Laser Altimeter System (GLAS) from 2002 to 2008 has the potential to form the basis of a globally consistent sample-based inventory of forest biomass. GLAS lidar return data were collected globally in spatially discrete full waveform "shots," which have been shown to be strongly correlated with aboveground forest biomass. Relationships observed at spatially coincident field plots may be used to model biomass at all GLAS shots, and well-established methods of model-based inference may then be used to estimate biomass and variance for specific spatial domains. However, the spatial pattern of GLAS acquisition is neither random across the surface of the earth nor is it identifiable with any particular systematic design. Undefined sample properties therefore hinder the use of GLAS in global forest sampling. RESULTS: We propose a method of identifying a subset of the GLAS data which can justifiably be treated as a simple random sample in model-based biomass estimation. The relatively uniform spatial distribution and locally arbitrary positioning of the resulting sample is similar to the design used by the US national forest inventory (NFI). We demonstrated model-based estimation using a sample of GLAS data in the US state of California, where our estimate of biomass (211 Mg/hectare) was within the 1.4% standard error of the design-based estimate supplied by the US NFI. The standard error of the GLAS-based estimate was significantly higher than the NFI estimate, although the cost of the GLAS estimate (excluding costs for the satellite itself) was almost nothing, compared to at least US$ 10.5 million for the NFI estimate. CONCLUSIONS: Global application of model-based estimation using GLAS, while demanding significant consolidation of training data, would improve inter-comparability of international biomass estimates by imposing consistent methods and a globally coherent sample frame. The methods presented here constitute a globally extensible approach for generating a simple random sample from the global GLAS dataset, enabling its use in forest inventory activities.

Evaluating kriging as a tool to improve moderate resolution maps of forest biomass.

The USDA Forest Service, Forest Inventory and Analysis program (FIA) recently produced a nationwide map of forest biomass by modeling biomass collected on forest inventory plots as nonparametric functions of moderate resolution satellite data and other environmental variables using Cubist software. Efforts are underway to develop methods to enhance this initial map. We explored the possibility of modeling spatial structure to make such improvements. Spatial structure in the field biomass data as well as in residuals from the map was investigated across 18 ecological zones in the Interior Western U.S. Exploratory tools included directional graphs of summary statistics, three dimensional maps, Moran's I correlograms, and variograms. Where spatial pattern was present, field and residual biomass were kriged, and predictions made for an independent test set were evaluated for improvement over predictions in the initial biomass map. While kriging has some potential benefit when analyzing the field data and exploring spatial structure, kriging residuals resulted in little or no improvement in the initial biomass map developed using Cubist software. Stationarity assumptions, variogram behavior, and appropriate model fitting strategies are discussed.