Climate change models for the northern Rocky Mountains predict changes in temperature and water availability that in turn will alter vegetation. Changes include timing of plant life-history events, or phenology, such as green-up, flowering and senescence, and shifts in species composition. Moreover, climate changes may favor different species, such as nonnative, annual grasses over native species. Changes in vegetation could make forage for ungulates, sage-grouse, and livestock available earlier in the growing season, but shifts in species composition and phenology may also result in earlier senescence (die-off or dormancy) and reduced overall forage production.

The information presented here provides the five-year science agenda for the North Central Climate Science Center. It is meant to be a high-level guide that describes the spatial context of the center, the primary partners and stakeholders, and the strategic framework the center will use in applying climate science to inform management.

Seasonal change is important to consider when managing conservation areas at landscape scales. The study of such patterns throughout the year is referred to as phenology. Recurring life-cycle events that are initiated and driven by environmental factors include animal migration and plant flowering. Phenological events capture public attention, such as fall color change in deciduous forests, the first flowering in spring, and for those with allergies, the start of the pollen season. These events can affect our daily lives, provide clues to help understand and manage ecosystems, and provide evidence of how climate variability can affect the natural cycle of plants and animals. Phenological observations can be gathered at a range of scales, from plots smaller than an acre to landscapes of hundreds to thousands of acres. Linking these observations to diverse disciplines such as evolutionary biology or climate sciences can help further research in species and ecosystem responses to climate change scenarios at appropriate scales. A cooperative study between the National Park Service (NPS), the U.S. Geological Survey (USGS), and the National Aeronautics and Space Administration (NASA) has been exploring how satellite information can be used to summarize phenological patterns observed at the park or landscape scale and how those summaries can be presented to both park managers and visitors. This study specifically addressed seasonal changes in plants, including the onset of growth, photosynthesis in the spring, and the senescence of deciduous vegetation in the fall. The primary objective of the work is to demonstrate that seasonality even in protected areas changes considerably across years. A major challenge is to decouple natural variability from possible trends—directional change that can lead to a permanent and radically different ecosystem state. Trends can be either a gradual degradation of the landscape (often from external influences) or steady improvement (by implementing long-term conservation plans). In either case, it is important to first grasp the magnitude of natural variation so that it is not confused with actual trends. This work used existing and freely available remote sensing data, specifically the NASA-funded 250-meter (m) spatial resolution land-surface phenology product for North America. This product is calculated from an annual record of vegetation health observed by NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. The land-surface phenology product is, in essence, a method to summarize all the observations throughout a year into a few key, ecologically relevant “metrics”.

Managers need new tools for detecting the movement and spread of nonnative, invasive species. Habitat suitability models are a popular tool for mapping the potential distribution of current invaders, but the ability of these models to prioritize monitoring efforts has not been tested in the field. We tested the utility of an iterative sampling design (i.e., models based on field observations used to guide subsequent field data collection to improve the model), hypothesizing that model performance would increase when new data were gathered from targeted sampling using criteria based on the initial model results. We also tested the ability of habitat suitability models to predict the spread of invasive species, hypothesizing that models would accurately predict occurrences in the field, and that the use of targeted sampling would detect more species with less sampling effort than a nontargeted approach. We tested these hypotheses on two species at the state scale (Centaurea stoebe and Pastinaca sativa) in Wisconsin (USA), and one genus at the regional scale (Tamarix) in the western United States. These initial data were merged with environmental data at 30-m2 resolution for Wisconsin and 1-km2 resolution for the western United States to produce our first iteration models. We stratified these initial models to target field sampling and compared our models and success at detecting our species of interest to other surveys being conducted during the same field season (i.e., nontargeted sampling). Although more data did not always improve our models based on correct classification rate (CCR), sensitivity, specificity, kappa, or area under the curve (AUC), our models generated from targeted sampling data always performed better than models generated from nontargeted data. For Wisconsin species, the model described actual locations in the field fairly well (kappa = 0.51, 0.19, P < 0.01), and targeted sampling did detect more species than nontargeted sampling with less sampling effort (χ2) = 47.42, P < 0.01). From these findings, we conclude that habitat suitability models can be highly useful tools for guiding invasive species monitoring, and we support the use of an iterative sampling design for guiding such efforts.

Buffelgrass, a highly competitive and flammable African bunchgrass, is spreading rapidly across both urban and natural areas in the Sonoran Desert of southern and central Arizona. Damages include increased fire risk, losses in biodiversity, and diminished revenues and quality of life. Feasibility of sustained and successful mitigation will depend heavily on rates of spread, treatment capacity, and cost–benefit analysis. We created a decision support model for the wildland–urban interface north of Tucson, AZ, using a spatial state-and-transition simulation modeling framework, the Tool for Exploratory Landscape Scenario Analyses. We addressed the issues of undetected invasions, identifying potentially suitable habitat and calibrating spread rates, while answering questions about how to allocate resources among inventory, treatment, and maintenance. Inputs to the model include a state-and-transition simulation model to describe the succession and control of buffelgrass, a habitat suitability model, management planning zones, spread vectors, estimated dispersal kernels for buffelgrass, and maps of current distribution. Our spatial simulations showed that without treatment, buffelgrass infestations that started with as little as 80 ha (198 ac) could grow to more than 6,000 ha by the year 2060. In contrast, applying unlimited management resources could limit 2060 infestation levels to approximately 50 ha. The application of sufficient resources toward inventory is important because undetected patches of buffelgrass will tend to grow exponentially. In our simulations, areas affected by buffelgrass may increase substantially over the next 50 yr, but a large, upfront investment in buffelgrass control could reduce the infested area and overall management costs.

We assessed the ability of climatic, environmental, and anthropogenic variables to predict areas of high-risk for plant invasion and consider the relative importance and contribution of these predictor variables by considering two spatial scales in a region of rapidly changing climate. We created predictive distribution models, using Maxent, for three highly invasive plant species (Canada thistle, white sweetclover, and reed canarygrass) in Alaska at both a regional scale and a local scale. Regional scale models encompassed southern coastal Alaska and were developed from topographic and climatic data at a 2 km (1.2 mi) spatial resolution. Models were applied to future climate (2030). Local scale models were spatially nested within the regional area; these models incorporated physiographic and anthropogenic variables at a 30 m (98.4 ft) resolution. Regional and local models performed well (AUC values > 0.7), with the exception of one species at each spatial scale. Regional models predict an increase in area of suitable habitat for all species by 2030 with a general shift to higher elevation areas; however, the distribution of each species was driven by different climate and topographical variables. In contrast local models indicate that distance to right-of-ways and elevation are associated with habitat suitability for all three species at this spatial level. Combining results from regional models, capturing long-term distribution, and local models, capturing near-term establishment and distribution, offers a new and effective tool for highlighting at-risk areas and provides insight on how variables acting at different scales contribute to suitability predictions. The combinations also provides easy comparison, highlighting agreement between the two scales, where long-term distribution factors predict suitability while near-term do not and vice versa.

The development of different energy balance models has allowed users to choose a model based on its suitability in a region. We compared four commonly used models—Mapping EvapoTranspiration at high Resolution with Internalized Calibration (METRIC) model, Surface Energy Balance Algorithm for Land (SEBAL) model, Surface Energy Balance System (SEBS) model, and the Operational Simplified Surface Energy Balance (SSEBop) model—using Landsat images to estimate evapotranspiration (ET) in the Midwestern United States. Our models validation using three AmeriFlux cropland sites at Mead, Nebraska, showed that all four models captured the spatial and temporal variation of ET reasonably well with an R2 of more than 0.81. Both the METRIC and SSEBop models showed a low root mean square error (<0.93 mm·day−1) and a high Nash–Sutcliffe coefficient of efficiency (>0.80), whereas the SEBAL and SEBS models resulted in relatively higher bias for estimating daily ET. The empirical equation of daily average net radiation used in the SEBAL and SEBS models for upscaling instantaneous ET to daily ET resulted in underestimation of daily ET, particularly when the daily average net radiation was more than 100 W·m−2. Estimated daily ET for both cropland and grassland had some degree of linearity with METRIC, SEBAL, and SEBS, but linearity was stronger for evaporative fraction. Thus, these ET models have strengths and limitations for applications in water resource management.

Correlative species distribution models are becoming commonplace in the scientific literature and public outreach products, displaying locations, abundance, or suitable environmental conditions for harmful invasive species, threatened and endangered species, or species of special concern. Accurate species distribution models are useful for efficient and adaptive management and conservation, research, and ecological forecasting. Yet, these models are often presented without fully examining or explaining the caveats for their proper use and interpretation and are often implemented without understanding the limitations and assumptions of the model being used. We describe common pitfalls, assumptions, and caveats of correlative species distribution models to help novice users and end users better interpret these models. Four primary caveats corresponding to different phases of the modeling process, each with supporting documentation and examples, include: (1) all sampling data are incomplete and potentially biased; (2) predictor variables must capture distribution constraints; (3) no single model works best for all species, in all areas, at all spatial scales, and over time; and (4) the results of species distribution models should be treated like a hypothesis to be tested and validated with additional sampling and modeling in an iterative process.

Evapotranspiration (ET) mapping at the Landsat spatial resolution (100 m) is essential to fully understand water use and water availability at the field scale. Water use estimates in the Colorado River Basin (CRB), which has diverse ecosystems and complex hydro-climatic regions, will be helpful to water planners and managers. Availability of Landsat 8 images, starting in 2013, provides the opportunity to map ET in the CRB to assess spatial distribution and patterns of water use. The Operational Simplified Surface Energy Balance (SSEBop) model was used with 528 Landsat 8 images to create seamless monthly and annual ET estimates at the inherent 100 m thermal band resolution. Annual ET values were summarized by land use/land cover classes. Croplands were the largest consumer of “blue” water while shrublands consumed the most “green” water. Validation using eddy covariance (EC) flux towers and water balance approaches showed good accuracy levels with R2 ranging from 0.74 to 0.95 and the Nash–Sutcliffe model efficiency coefficient ranging from 0.66 to 0.91. The root mean square error (and percent bias) ranged from 0.48 mm (13%) to 0.60 mm (22%) for daily (days of satellite overpass) ET and from 7.75 mm (2%) to 13.04 mm (35%) for monthly ET. The spatial and temporal distribution of ET indicates the utility of Landsat 8 for providing important information about ET dynamics across the landscape. Annual crop water use was estimated for five selected irrigation districts in the Lower CRB where annual ET per district ranged between 681 mm to 772 mm. Annual ET by crop type over the Maricopa Stanfield irrigation district ranged from a low of 384 mm for durum wheat to a high of 990 mm for alfalfa fields. A rainfall analysis over the five districts suggested that, on average, 69% of the annual ET was met by irrigation. Although the enhanced cloud-masking capability of Landsat 8 based on the cirrus band and utilization of the Fmask algorithm improved the removal of contaminated pixels, the ability to reliably estimate ET over clouded areas remains an important challenge. Overall, the performance of Landsat 8 based ET compared to available EC datasets and water balance estimates for a complex basin such as the CRB demonstrates the potential of using Landsat 8 for annual water use estimation at a national scale. Future efforts will focus on (a) use of consistent methodology across years, (b) integration of multiple sensors to maximize images used, and (c) employing cloud-computing platforms for large scale processing capabilities.