Semiparametric space–time survival modeling of chronic wasting disease in deer

In this paper, we propose a semiparametric survival model to investigate the pattern of spatial and temporal variation in disease prevalence of chronic wasting disease (CWD) in wild deer in Wisconsin over the years 2002 and 2006. The semiparametric survival model we suggested allows to build a more flexible model than the parametric model with fewer parametric assumptions by modeling the baseline hazard using a Gamma process prior. Based on the proposed model, we investigate the geographical distribution of CWD, and assess the effect of sex on disease prevalence. We use a Bayesian hierarchical framework where latent parameters capture temporal and spatial trends in disease incidence, incorporating sex and spatially correlated random effects. We also propose bivariate baseline hazard which change over age and time simultaneously to adopt different effects of age and time on the baseline hazard. Inference is carried out by using MCMC simulation techniques in a fully Bayesian framework. Our results suggest that disease has been spreaded mainly in the disease eradication zone and male deer show a significantly higher infection probability than female deer.     

Identifying landscape scale patterns from individual scale processes

Extrapolating across scales is a critical problem in ecology. Explicit mechanistic models of ecological systems provide a bridge from measurements of processes at small and short scales to larger scales; spatial patterns at large scales can be used to test the outcomes of these models. However, it is necessary to identify patterns that are not dependent on initial conditions, because small scale initial conditions will not normally be measured at large scales. We examined one possible pattern that could meet these conditions, the relationship between mean and variance in abundance of a parasitic tick in an individual based model of a lizard tick interaction. We scaled discrepancies between the observed and simulated patterns with a transformation of the variance–covariance matrix of the observed pattern to objectively identify patterns that are “close”.     

Analyzing the spatial structure of a Sri Lankan tree species with multiple scales of clustering

Clustering at multiple critical scales may be common for plants since many different factors and processes may cause clustering. This is especially true for tropical rain forests for which theories explaining species coexistence and community structure rest heavily on spatial patterns. We used point pattern analysis to analyze the spatial structure of Shorea congestiflora, a dominant species in a 25-ha forest dynamics plot in a rain forest at Sinharaja World Heritage Site (Sri Lanka), which apparently shows clustering at several scales. We developed cluster processes incorporating two critical scales of clustering for exploring the spatial structure of S. congestiflora and interpret it in relation to factors such as competition, dispersal limitation, recruitment limitation, and Janzen-Connell effects. All size classes showed consistent large-scale clustering with a cluster radius of approximately 25 m. Inside the larger clusters, small-scale clusters with a radius of 8 m were evident for recruits and saplings, weak for intermediates, and disappeared for adults. The pattern of all trees could be divided into two independent patterns: a random pattern (nearest neighbor distance > 8 m) comprising approximately 12% of the trees and a nested double-cluster pattern. This finding suggests two independent recruitment and/or seed dispersal mechanisms. Saplings were several times as abundant as recruits and may accumulate several recruit generations. Recruits were only weakly associated with adults and occupied about half of the large-scale clusters, but saplings almost all. This is consistent with recruitment limitation. For approximately 70% (95%) of all juveniles the nearest adult was less than 26 m away (53 m), suggesting a dispersal limitation that may also be related to the critical large-scale clustering. Our example illustrates the manner in which the use of a specific and complex null hypothesis of spatial structure in point pattern analysis can help us better understand the biology of a species and generate specific hypotheses to be further investigated in the field.     

Two-dimensional wavelet analysis of spruce budworm host basal area in the Border Lakes landscape

Increases in the extent and severity of spruce budworm (Choristoneura fumiferana Clem.) outbreaks over the last century are thought to be the result of changes in forest structure due to forest management. A corollary of this hypothesis is that manipulations of forest structure and composition can be used to reduce future forest vulnerability. However, to what extent historical forest management has influenced current spatial patterns of spruce budworm host species is unknown. To identify landscape-scale spatial legacies of forest management in patterns of spruce budworm host species (i.e., Abies balsamea and Picea spp.), we analyzed remotely sensed forest data from the Border Lakes landscape of northern Minnesota and northwestern Ontario. Our study area contains three regions with different management histories: (1) fine-scale logging patterns in Minnesota, (2) coarse-scale logging patterns in Ontario, and (3) very limited logging history in the Boundary Waters Canoe Area and adjacent Quetico Provincial Park. We analyzed forest basal-area data using wavelets and null models to identify: (1) at which scales forest basal area is structured, (2) where those scales of pattern are significantly present, and (3) whether regions of local significance correspond to regional boundaries that separate the study area. Results indicate that spatial patterns in host basal area are created by nonstationary processes and that these processes are further constrained by lakes and wetlands. Wavelet analysis combined with significance testing revealed a bimodal distribution of scale-specific wavelet variance and separate zones of host species basal area that partially correspond with regional boundaries, particularly between Minnesota and the Wilderness region. This research represents one of the first comparisons of forest spatial structure in this region across an international border and presents a novel method of two-dimensional wavelet analysis that can be used to identify significant scale-specific structure in spatial data.     

Emerging prion disease drives host selection in a wildlife population

Infectious diseases are increasingly recognized as an important force driving population dynamics, conservation biology, and natural selection in wildlife populations. Infectious agents have been implicated in the decline of small or endangered populations and may act to constrain population size, distribution, growth rates, or migration patterns. Further, diseases may provide selective pressures that shape the genetic diversity of populations or species. Thus, understanding disease dynamics and selective pressures from pathogens is crucial to understanding population processes, managing wildlife diseases, and conserving biological diversity. There is ample evidence that variation in the prion protein gene (PRNP) impacts host susceptibility to prion diseases. Still, little is known about how genetic differences might influence natural selection within wildlife populations. Here we link genetic variation with differential susceptibility of white-tailed deer to chronic wasting disease (CWD), with implications for fitness and disease-driven genetic selection. We developed a single nucleotide polymorphism (SNP) assay to efficiently genotype deer at the locus of interest (in the 96th codon of the PRNP gene). Then, using a Bayesian modeling approach, we found that the more susceptible genotype had over four times greater risk of CWD infection; and, once infected, deer with the resistant genotype survived 49% longer (8.25 more months). We used these epidemiological parameters in a multi-stage population matrix model to evaluate relative fitness based on genotype-specific population growth rates. The differences in disease infection and mortality rates allowed genetically resistant deer to achieve higher population growth and obtain a long-term fitness advantage, which translated into a selection coefficient of over 1% favoring the CWD-resistant genotype. This selective pressure suggests that the resistant allele could become dominant in the population within an evolutionarily short time frame. Our work provides a rare example of a quantifiable disease-driven selection process in a wildlife population, demonstrating the potential for infectious diseases to alter host populations. This will have direct bearing on the epidemiology, dynamics, and future trends in CWD transmission and spread. Understanding genotype-specific epidemiology will improve predictive models and inform management strategies for CWD-affected cervid populations.     

抚育间伐对栓皮栎种群空间分布格局的影响

抚育间伐是一种重要的改善林木生长条件的经营措施,对林分结构和动态具有重要影响。为研究抚育间伐对林木种群空间结构与格局影响的内在机制,以间伐和未间伐的栓皮栎人工林为研究对象,通过设置2个100 m×100 m样地并进行每木定位和调查,在采用径级结构代替年龄结构方法将栓皮栎种群划分为幼树(2 cm≤DBH<5 cm)、小树(5 cm≤DBH<13 cm)和大树(DBH≥13 cm)3个生长阶段的基础上,分别采用Ripley’s K函数衍生的g(r)函数和双变量g12(r)函数,对栓皮栎种群空间分布点格局及不同生长阶段栓皮栎种群之间的关联性进行了研究。结果表明,间伐和未间伐样地栓皮栎种群空间分布点格局分别在0-16 m和0-33 m距离尺度内呈聚集分布,而分别在大于16 m和33 m距离尺度内呈随机分布;间伐和未间伐样地栓皮栎幼树、小树和大树的株数比分别为8?741?699和261?1134?683,且间伐样地栓皮栎幼树、小树和大树种群分别在0-14、1-16、0-6 m距离尺度内呈现均匀或聚集分布,而在其他距离尺度上表现为随机分布;栓皮栎幼树、小树和大树之间仅在间伐样地0-6 m距离尺度内呈现一定的相关性,而在未间伐样地更大的距离尺度内有显著的关联性,如栓皮栎幼树和大树之间在6-38 m距离尺度上呈显著正相关。因此,抚育间伐一定程度上使得栓皮栎种群在更大距离尺度上呈现出随机分布状态,并弱化了不同生长阶段的林木栓皮栎种群的关联性,这调整了栓皮栎种群空间竞争关系,有利于大径级林木个体的培育。该研究可以为开展抚育间伐对林木种群的影响的研究提供理论依据,也可以为制定科学合理的抚育技术措施提供参考。     

A methodological framework for integrating computational fluid dynamics and ecological models applied to juvenile freshwater mussel dispersal in the Upper Mississippi River

Interdisciplinary research in hydraulics and ecology for river management and restoration must integrate processes that occur over a wide range of spatial and temporal scales, which presents a challenge to ecohydraulics modelers. Computational fluid dynamics (CFD) models are being more widely used to determine flow fields for ecohydraulics applications. In the Upper Mississippi River (UMR), the mussel dynamics model was developed as a tool for management and conservation of freshwater mussels (Unionidae), which are benthic organisms, imperiled in North America, that are inextricably linked with the hydraulics of river flow. We updated the juvenile dispersal component of the mussel dynamics model by using stochastic Lagrangian particle tracking in a three dimensional flow field output from CFD models of reaches in the UMR. We developed a methodological framework to integrate hydrodynamic data with the mussel dynamics model, and we demonstrate the use of the juvenile dispersal model employed within the methodological framework in two reaches of the UMR. The method was used to test the hypothesis that impoundment affects the relationship of some hydraulic parameters with juvenile settling distribution. Simulation results were consistent with this hypothesis, and the relationships of bed shear stress and Froude number with juvenile settling were altered by impoundment most likely through effects on local hydraulics. The methodological framework is robust, integrates Eulerian and Lagrangian reference frameworks, and incorporates processes over a wide range of temporal and spatial scales, from watershed scale hydrologic processes (decades), to reach scale (km) processes that occur over hours or days, and turbulent processes on spatial scales of meter to millimeter and times scales of seconds. The methods are presently being used to assess the impacts of pre- and early post-settlement processes on mussel distributions, including the effects of bed shear stress, and the sensitivity of the location of the host fish when juveniles excyst, on juvenile settling distribution.     

Quantifying the importance of local niche-based and stochastic processes to tropical tree community assembly

Although niche-based and stochastic processes, including dispersal limitation and demographic stochasticity, can each contribute to community assembly, it is difficult to quantify the relative importance of each process in natural vegetation. Here, we extend Shipley's maxent model (Community Assembly by Trait Selection, CATS) for the prediction of relative abundances to incorporate both trait-based filtering and dispersal limitation from the larger landscape and develop a statistical decomposition of the proportions of the total information content of relative abundances in local communities that are attributable to trait-based filtering, dispersal limitation, and demographic stochasticity. We apply the method to tree communities in a mature, species-rich, tropical forest in French Guiana at 1-, 0.25- and 0.04-ha scales. Trait data consisted of species' means of 17 functional traits measured over both the entire meta-community and separately in each of nine 1-ha plots. Trait means calculated separately for each site always gave better predictions. There was clear evidence of trait-based filtering at all spatial scales. Trait-based filtering was the most important process at the 1-ha scale (34%), whereas demographic stochasticity was the most important at smaller scales (37-53%). Dispersal limitation from the meta-community was less important and approximately constant across scales (-9%), and there was also an unresolved association between site-specific traits and meta-community relative abundances. Our method allows one to quantify the relative importance of local niche-based and meta-community processes and demographic stochasticity during community assembly across spatial and temporal scales.     

Need for multiscale planning for conservation of urban bats

For over a century there have been continual efforts to incorporate nature into urban planning. These efforts (i.e., urban reconciliation) aim to manage and create habitats that support biodiversity within cities. Given that species select habitat at different spatial scales, understanding the scale at which urban species respond to their environment is critical to the success of urban reconciliation efforts. We assessed species–habitat relationships for common bat species at 50‐m, 500‐m, and 1 km spatial scales in the Chicago (U.S.A.) metropolitan area and predicted bat activity across the greater Chicago region. Habitat characteristics across all measured scales were important predictors of silver‐haired bat (Lasionycteris noctivagans) and eastern red bat (Lasiurus borealis) activity, and big brown bat (Eptesicus fuscus) activity was significantly lower at urban sites relative to rural sites. Open vegetation had a negative effect on silver‐haired bat activity at the 50‐m scale but a positive effect at the 500‐m scale, indicating potential shifts in the relative importance of some habitat characteristics at different scales. These results demonstrate that localized effects may be constrained by broader spatial patterns. Our findings highlight the importance of considering scale in urban reconciliation efforts and our landscape predictions provide information that can help prioritize urban conservation work.     

Distribution of a marginal population of <Emphasis Type="Italic">Mytilus edulis</Emphasis>: responses to biotic and abiotic processes at different spatial scales

Physical and biological processes interact to produce pattern in nature. Pattern is scale dependent as processes generating pattern are heterogeneous in time and space. We tested some causes of variation in abundance and distribution of three marginal populations of sublittoral blue mussels, Mytilus edulis, in the non-tidal northeastern Baltic Sea. We studied the role of substrate inclination, perennial algae and siltation along local wave exposure gradients on mussel distribution over a regional salinity gradient. We found marked differences on regional scales (p < 0.001) with lower densities and biomasses of mussels with declining salinity. Along local gradients, mussel densities increased with increasing exposure (p < 0.001) and declining slope and sedimentation (p < 0.01). Site specifically, densities of blue mussels and the perennial red algae, Furcellaria lumbricalis, were positively related, results supported by a colonisation experiment. Also, young post-recruits showed significant relations to adult biomass, wave exposure, algal biomass, bottom slope and sediment cover. Findings showed that the relative importance of the determinants affecting blue mussels at the edge of their range vary with scale and are affected by the density and size structure of mussel populations. The study provides an indication of the types of factors that may be invoked as causes of spatial variation in marginal blue mussel populations and reinforces the need to consider multiple aspects when distributional patterns are assessed.     

Measurement error causes scale-dependent threshold erosion of biological signals in animal movement data.

Recent advances in telemetry technology have created a wealth of tracking data available for many animal species moving over spatial scales from tens of meters to tens of thousands of kilometers. Increasingly, such data sets are being used for quantitative movement analyses aimed at extracting fundamental biological signals such as optimal searching behavior and scale-dependent foraging decisions. We show here that the location error inherent in various tracking technologies reduces the ability to detect patterns of behavior within movements. Our analyses endeavored to set out a series of initial ground rules for ecologists to help ensure that sampling noise is not misinterpreted as a real biological signal. We simulated animal movement tracks using specialized random walks known as Lévy flights at three spatial scales of investigation: 100-km, 10-km, and 1-km maximum daily step lengths. The locations generated in the simulations were then blurred using known error distributions associated with commonly applied tracking methods: the Global Positioning System (GPS), Argos polar-orbiting satellites, and light-level geolocation. Deviations from the idealized Lévy flight pattern were assessed for each track after incrementing levels of location error were applied at each spatial scale, with additional assessments of the effect of error on scale-dependent movement patterns measured using fractal mean dimension and first-passage time (FPT) analyses. The accuracy of parameter estimation (Lévy mu, fractal mean D, and variance in FPT) declined precipitously at threshold errors relative to each spatial scale. At 100-km maximum daily step lengths, error standard deviations of > or = 10 km seriously eroded the biological patterns evident in the simulated tracks, with analogous thresholds at the 10-km and 1-km scales (error SD > or = 1.3 km and 0.07 km, respectively). Temporal subsampling of the simulated tracks maintained some elements of the biological signals depending on error level and spatial scale. Failure to account for large errors relative to the scale of movement can produce substantial biases in the interpretation of movement patterns. This study provides researchers with a framework for understanding the limitations of their data and identifies how temporal subsampling can help to reduce the influence of spatial error on their conclusions.     

Independent effects of fragmentation on forest songbirds: an organism-based approach.

The degree to which spatial patterns influence the dynamics and distribution of populations is a central question in ecology. This question is even more pressing in the context of rapid habitat loss and fragmentation, which threaten global biodiversity. However, the relative influence of habitat loss and landscape fragmentation, the spatial patterning of remaining habitat, remains unclear. If landscape pattern affects population size, managers may be able to design landscapes that mitigate habitat loss. We present the results of a mensurative experiment designed to test four habitat loss vs. fragmentation hypotheses. Unlike previous studies, we measured landscape structure using quantitative, spatially explicit habitat distribution models previously developed for two species: Blackburnian Warbler (Dendroica fusca) and Ovenbird (Seiurus aurocapilla). We used a stratified sampling design that reduced the confounding of habitat amount and fragmentation variables. Occurrence and reoccurrence of both species were strongly influenced by characteristics at scales greater than the individual territory, indicating little support for the random-sample hypothesis. However, the type and spatial extent of landscape influence differed. Both occurrence and reoccurrence of Blackburnian Warblers were influenced by the amount of poor-quality matrix at 300- and 2000-m spatial extents. The occurrence and reoccurrence of Ovenbirds depended on a landscape pattern variable, patch size, but only in cases when patches were isolated. These results support the hypothesis that landscape pattern is important for some species only when the amount of suitable habitat is low. Although theoretical models have predicted such an interaction between landscape fragmentation and composition, to our knowledge this is the first study to report empirical evidence of such nonlinear fragmentation effects. Defining landscapes quantitatively from an organism-based perspective may increase power to detect fragmentation effects, particularly in forest mosaics where boundaries between patches and matrix are ambiguous. Our results indicate that manipulating landscape pattern may reduce negative impacts of habitat loss for Ovenbird, but not Blackburnian Warbler. We emphasize that most variance in the occurrence of both species was explained by local scale or landscape composition variables rather than variables reflecting landscape pattern.     

Diversity pattern and spatial scale: a study of a tropical rain forest of Malaysia

Scale is emerging as one of the critical problems in ecology because our perception of most ecological variables and processes depends upon the scale at which the variables are measured. A conclusion obtained at one scale may not be valid at another scale without sufficient knowledge of the scaling effect, which is also a source of misinterpretation for many ecological problems, such as the design of reserves in conservation biology.This paper attempts to study empirically how scaling may affect the spatial patterns of diversity (tree density, richness and Shannon diversity) that we may perceive in tropical forests, using as a test-case a 50 ha forest plot in Malaysia. The effect of scale on measurements of diversity patterns, the occurrence of rare species, the fractal dimension of diversity patterns, the spatial structure and the nearest-neighbour autocorrelation of diversity are addressed. The response of a variable to scale depends on the way it is measured and the way it is distributed in space.We conclude that, in general, the effect of scaling on measures of biological diversity is non-linear; heterogeneity increases with the size of the sampling units, and fine-scale information is lost at a broad scale. Our results should lead to a better understanding of how ecological variables and processes change over scale.     

Spatial variance in abundance and occupancy of corals across broad geographic scales

Species assemblages vary in structure due to a wide variety of processes operating at ecological and much broader biogeographical scales. Cross-scale studies of assemblage structure are necessary to fully understand this variability. Here, we evaluate the abundance and occupancy patterns of hierarchically sampled coral assemblages in three habitats (reef flat, crest, and slope) and five regions (Indonesia, Papua New Guinea, the Solomon Islands, American Samoa, and the Society Islands) across the west-central Pacific Ocean. Specifically, we compare two alternative models that unify spatial variance and occupancy via the negative binomial distribution. The first assumes a power-law scaling between the mean and variance of abundance; the second assumes a quadratic variance-mean relationship and a constant abundance-invariant aggregation parameter. Surprisingly, the well-established power-law model performs worse than the model assuming abundance-invariant aggregation, for both variance-mean and occupancy-abundance relationships. We also find strong evidence for regional and habitat variation in these relationships and in the levels of aggregation estimated by the abundance-invariant aggregation model. Among habitats, corals on reef flats exhibited lower occupancy and higher levels of aggregation compared to reef crests and slopes. Among regions, low occupancy and high aggregation were most pronounced across all habitats in American Samoa. These patterns may be related to habitat and regional differences in disturbance and recovery processes. Our results suggest that the spatial scaling of abundance and occupancy is sensitive to processes operating among these habitats and at regional scales. However, the consistency of these relationships across species within assemblages suggests that a theoretical unification of spatial variance and occupancy patterns is indeed possible.     

In search of quorum effects in metacommunity structure: species co-occurrence analyses

Alternative models of community assembly emphasize regional, stochastic, dispersal-based processes or local, deterministic, niche-based processes. Community ecology's historical focus on local processes implicitly assumes that local processes surpass regional processes over time or across space to derive nonrandom metacommunity structure (i.e., a quorum effect). Quorum effects are expected late in succession among nearby sites, whereas quorum effects are not expected early in succession among distant sites. I conducted a meta-analysis of zooplankton data sets encompassing time scales of one to thousands of years and spatial scales of <1 m to thousands of kilometers. Species co-occurrence analyses statistically evaluated presence/absence patterns relative to random patterns obtained with Monte Carlo null models. A series of weighted analyses was conducted and alternative randomization algorithms and null models were evaluated. Most zooplankton metacommunities were randomly structured in unweighted analyses, and the distribution of significant structure did not follow quorum effect predictions. Weighted analyses (e.g., by habitat area) revealed significant, nonrandom structure in most zooplankton metacommunities, but the distribution of significant structure still did not adhere well to quorum effect predictions. Finally, additional weighting for study scale (number of sites) nullified most significant area-weighted structure, and again, the distribution of significant structure did not follow quorum effect predictions. Overall, a quorum effect was not supported, perhaps related to zooplankton life histories and energetics and/or the quorum effect itself. Results at the presence/absence level of resolution indicated that local processes did not generally override regional processes over time or across space to drive community structure. A full integration of dispersal- and niche-based concepts in metacommunity dynamics will be most fruitful for unraveling community assembly. Species co-occurrence analyses were scale dependent (habitat area and study size). Future analyses should use weights for important factors (e.g., habitat area), and meta-analyses should include study scale as an additional factor contributing to apparent patterns.     

Spatial scale affects community concordance among fishes, benthic macroinvertebrates, and bryophytes in streams.

Owing to the lack of information about the distribution patterns of many taxonomic groups, biodiversity conservation strategies commonly rely on a surrogate taxa approach for identifying areas of maximum conservation potential. Macroinvertebrates or fish are the most likely candidates for such a role in many freshwater systems. The usefulness of the surrogate taxa depends largely on community concordance, i.e., the degree of similarity in community patterns among taxonomic groups across a set of sites. We examined the effect of the spatial scale of a. study on the strength of community concordance among macroinvertebrates, bryophytes, and fish by comparing the concordance between ordinations of these groups in 101 boreal stream sites. We specifically asked if communities spanning several drainages are more concordant than those originating from a single drainage system. Our results indicate that community concordance is affected by spatial extent, being variable and generally weak at the scale of individual drainages, but strong across multiple drainage systems and ecoregions. We attribute this finding to different taxonomic groups responding to similar environmental factors and sharing a similar latitudinal gradient of community structure when viewed across large spatial scales. We also identified a "gradient of concordance," with sites contributing disproportionately to community concordance being in relatively large streams with high microhabitat variability. Overall, our results suggest that the degree of community concordance among freshwater organism groups depends critically on the spatial extent of the study, and surrogate groups at the scale of single river systems should be used with caution.     

Simulating Historical Variability in the Amount of Old Forests in the Oregon Coast Range

Abstract: We developed the landscape age-class demographics simulator ( LADS) to model historical variability in the amount of old-growth and late-successional forest in the Oregon Coast Range over the past 3,000 years. The model simulated temporal and spatial patterns of forest fires along with the resulting fluctuations in the distribution of forest age classes across the landscape. Parameters describing historical fire regimes were derived from data from a number of existing dendroecological and paleoecological studies. Our results indicated that the historical age-class distribution was highly variable and that variability increased with decreasing landscape size. Simulated old-growth percentages were generally between 25% and 75% at the province scale (2,250,000 ha) and never fell below 5%. In comparison, old-growth percentages varied from 0 to 100% at the late-successional reserve scale (40,000 ha). Province-scale estimates of current old-growth (5%) and late-successional forest (11%) in the Oregon Coast Range were lower than expected under the simulated historical fire regime, even when potential errors in our parameter estimates were considered. These uncertainties do, however, limit our ability to precisely define ranges of historical variability. Our results suggest that in areas where historical disturbance regimes were characterized by large, infrequent fires, management of forest age classes based on a range of historical variability may be feasible only at relatively large spatial scales. Comprehensive landscape management strategies will need to consider other factors besides the percentage of old forests on the landscape, including the spatial pattern of stands and the rates and pathways of landscape change.     

Evidence, exaggeration, and error in historical accounts of chaparral wildfires in California.

For more than half a century, ecologists and historians have been integrating the contemporary study of ecosystems with data gathered from historical sources to evaluate change over broad temporal and spatial scales. This approach is especially useful where ecosystems were altered before formal study as a result of natural resources management, land development, environmental pollution, and climate change. Yet, in many places, historical documents do not provide precise information, and pre-historical evidence is unavailable or has ambiguous interpretation. There are similar challenges in evaluating how the fire regime of chaparral in California has changed as a result of fire suppression management initiated at the beginning of the 20th century. Although the firestorm of October 2003 was the largest officially recorded in California (approximately 300,000 ha), historical accounts of pre-suppression wildfires have been cited as evidence that such a scale of burning was not unprecedented, suggesting the fire regime and patch mosaic in chaparral have not substantially changed. We find that the data do not support pre-suppression megafires, and that the impression of large historical wildfires is a result of imprecision and inaccuracy in the original reports, as well as a parlance that is beset with hyperbole. We underscore themes of importance for critically analyzing historical documents to evaluate ecological change. A putative 100 mile long by 10 mile wide (160 x 16 km) wildfire reported in 1889 was reconstructed to an area of chaparral approximately 40 times smaller by linking local accounts to property tax records, voter registration rolls, claimed insurance, and place names mapped with a geographical information system (GIS) which includes data from historical vegetation surveys. We also show that historical sources cited as evidence of other large chaparral wildfires are either demonstrably inaccurate or provide anecdotal information that is immaterial in the appraisal of pre-suppression fire size. Since historical evidence is inadequate for reconstructing a statistical distribution of pre-suppression fire sizes to compare with post-suppression data, other more propitious methods of evaluating change are discussed.