Coupled Downscaled Climate Models and Ecophysiological Metrics Forecast Habitat Compression for an Endangered Estuarine Fish.

Climate change is driving rapid changes in environmental conditions and affecting population and species' persistence across spatial and temporal scales. Integrating climate change assessments into biological resource management, such as conserving endangered species, is a substantial challenge, partly due to a mismatch between global climate forecasts and local or regional conservation planning. Here, we demonstrate how outputs of global climate change models can be downscaled to the watershed scale, and then coupled with ecophysiological metrics to assess climate change effects on organisms of conservation concern. We employed models to estimate future water temperatures (2010-2099) under several climate change scenarios within the large heterogeneous San Francisco Estuary. We then assessed the warming effects on the endangered, endemic Delta Smelt, Hypomesus transpacificus, by integrating localized projected water temperatures with thermal sensitivity metrics (tolerance, spawning and maturation windows, and sublethal stress thresholds) across life stages. Lethal temperatures occurred under several scenarios, but sublethal effects resulting from chronic stressful temperatures were more common across the estuary (median >60 days above threshold for >50% locations by the end of the century). Behavioral avoidance of such stressful temperatures would make a large portion of the potential range of Delta Smelt unavailable during the summer and fall. Since Delta Smelt are not likely to migrate to other estuaries, these changes are likely to result in substantial habitat compression. Additionally, the Delta Smelt maturation window was shortened by 18-85 days, revealing cumulative effects of stressful summer and fall temperatures with early initiation of spring spawning that may negatively impact fitness. Our findings highlight the value of integrating sublethal thresholds, life history, and in situ thermal heterogeneity into global change impact assessments. As downscaled climate models are becoming widely available, we conclude that similar assessments at management-relevant scales will improve the scientific basis for resource management decisions.

Correspondence of biological condition models of California streams at statewide and regional scales.

We used boosted regression trees (BRT) to model stream biological condition as measured by benthic macroinvertebrate taxonomic completeness, the ratio of observed to expected (O/E) taxa. Models were developed with and without exclusion of rare taxa at a site. BRT models are robust, requiring few assumptions compared with traditional modeling techniques such as multiple linear regression. The BRT models were constructed to provide baseline support to stressor delineation by identifying natural physiographic and human land use gradients affecting stream biological condition statewide and for eight ecological regions within the state, as part of the development of numerical biological objectives for California's wadeable streams. Regions were defined on the basis of ecological, hydrologic, and jurisdictional factors and roughly corresponded with ecoregions. Physiographic and land use variables were derived from geographic information system coverages. The model for the entire state (n = 1,386) identified a composite measure of anthropogenic disturbance (the sum of urban, agricultural, and unmanaged roadside vegetation land cover) within the local watershed as the most important variable, explaining 56% of the variance in O/E values. Models for individual regions explained between 51 and 84% of the variance in O/E values. Measures of human disturbance were important in the three coastal regions. In the South Coast and Coastal Chaparral, local watershed measures of urbanization were the most important variables related to biological condition, while in the North Coast the composite measure of human disturbance at the watershed scale was most important. In the two mountain regions, natural gradients were most important, including slope, precipitation, and temperature. The remaining three regions had relatively small sample sizes (n ≤ 75 sites) and had models that gave mixed results. Understanding the spatial scale at which land use and land cover affect taxonomic completeness is imperative for sound management. Our results suggest that invertebrate taxonomic completeness is affected by human disturbance at the statewide and regional levels, with some differences among regions in the importance of natural gradients and types of human disturbance. The construction and application of models similar to the ones presented here could be useful in the planning and prioritization of actions for protection and conservation of biodiversity in California streams.

Stream macroinvertebrate response models for bioassessment metrics: addressing the issue of spatial scale.

We developed independent predictive disturbance models for a full regional data set and four individual ecoregions (Full Region vs. Individual Ecoregion models) to evaluate effects of spatial scale on the assessment of human landscape modification, on predicted response of stream biota, and the effect of other possible confounding factors, such as watershed size and elevation, on model performance. We selected macroinvertebrate sampling sites for model development (n = 591) and validation (n = 467) that met strict screening criteria from four proximal ecoregions in the northeastern U.S.: North Central Appalachians, Ridge and Valley, Northeastern Highlands, and Northern Piedmont. Models were developed using boosted regression tree (BRT) techniques for four macroinvertebrate metrics; results were compared among ecoregions and metrics. Comparing within a region but across the four macroinvertebrate metrics, the average richness of tolerant taxa (RichTOL) had the highest R(2) for BRT models. Across the four metrics, final BRT models had between four and seven explanatory variables and always included a variable related to urbanization (e.g., population density, percent urban, or percent manmade channels), and either a measure of hydrologic runoff (e.g., minimum April, average December, or maximum monthly runoff) and(or) a natural landscape factor (e.g., riparian slope, precipitation, and elevation), or a measure of riparian disturbance. Contrary to our expectations, Full Region models explained nearly as much variance in the macroinvertebrate data as Individual Ecoregion models, and taking into account watershed size or elevation did not appear to improve model performance. As a result, it may be advantageous for bioassessment programs to develop large regional models as a preliminary assessment of overall disturbance conditions as long as the range in natural landscape variability is not excessive.

Responses of benthic macroinvertebrates to environmental changes associated with urbanization in nine metropolitan areas.

Responses of benthic macroinvertebrates along gradients of urban intensity were investigated in nine metropolitan areas across the United States. Invertebrate assemblages in metropolitan areas where forests or shrublands were being converted to urban land were strongly related to urban intensity. In metropolitan areas where agriculture and grazing lands were being converted to urban land, invertebrate assemblages showed much weaker or nonsignificant relations with urban intensity because sites with low urban intensity were already degraded by agriculture. Ordination scores, the number of EPT taxa, and the mean pollution-tolerance value of organisms at a site were the best indicators of changes in assemblage condition. Diversity indices, functional groups, behavior, and dominance metrics were not good indicators of urbanization. Richness metrics were better indicators of urban effects than were abundance metrics, and qualitative samples collected from multiple habitats gave similar results to those of single habitat quantitative samples (riffles or woody snags) in all metropolitan areas. Changes in urban intensity were strongly correlated with a set of landscape variables that was consistent across all metropolitan areas. In contrast, the instream environmental variables that were strongly correlated with urbanization and invertebrate responses varied among metropolitan areas. The natural environmental setting determined the biological, chemical, and physical instream conditions upon which urbanization acts and dictated the differences in responses to urbanization among metropolitan areas. Threshold analysis showed little evidence for an initial period of resistance to urbanization. Instead, assemblages were degraded at very low levels of urbanization, and response rates were either similar across the gradient or higher at low levels of urbanization. Levels of impervious cover that have been suggested as protective of streams (5-10%) were associated with significant assemblage degradation and were not protective.

A quantitative tool for assessing the integrity of southern coastal California streams.

We developed a benthic macro-invertebrate index of biological integrity (B-IBI) for the semi-arid and populous southern California coastal region. Potential reference sites were screened from a pool of 275 sites, first with quantitative GIS landscape analysis at several spatial scales and then with local condition assessments (in-stream and riparian) that quantified stressors acting on study reaches. We screened 61 candidate metrics for inclusion in the B-IBI based on three criteria: sufficient range for scoring, responsiveness to watershed and reach-scale disturbance gradients, and minimal correlation with other responsive metrics. Final metrics included: percent collector-gatherer + collector-filterer individuals, percent non-insect taxa, percent tolerant taxa, Coleoptera richness, predator richness, percent intolerant individuals, and EPT richness. Three metrics had lower scores in chaparral reference sites than in mountain reference sites and were scored on separate scales in the B-IBI. Metrics were scored and assembled into a composite B-IBI, which was then divided into five roughly equal condition categories. PCA analysis was used to demonstrate that the B-IBI was sensitive to composite stressor gradients; we also confirmed that the B-IBI scores were not correlated with elevation, season, or watershed area. Application of the B-IBI to an independent validation dataset (69 sites) produced results congruent with the development dataset and a separate repeatability study at four sites in the region confirmed that the B-IBI scoring is precise. The SoCal B-IBI is an effective tool with strong performance characteristics and provides a practical means of evaluating biotic condition of streams in southern coastal California.