Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Climate”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16Linked to original sources

Climate and hydrological changes in the northeastern United States: recent trends and implications for forested and aquatic ecosystems

We review twentieth century and projected twenty-first century changes in climatic and hydrologic conditions in the northeastern United States and the implications of these changes for forest ecosystems. Climate warming and increases in precipitation and associated changes in snow and hydrologic regimes have been observed over the last century, with the most pronounced changes occurring since 1970. Trends in specific climatic and hydrologic variables differ in their responses spatially (e.g., coastal vs. inland) and temporally (e.g., spring vs. summer). Trends can differ depending on the period of record analyzed, hinting at the role of decadal-scale climatic variation that is superimposed over the longer-term trend. Model predictions indicate that continued increases in temperature and precipitation across the northeastern United States can be expected over the next century. Ongoing increases in growing season length (earlier spring and later autumn) will most likely increase evapotranspiration and frequency of drought. In turn, an increase in the frequency of drought will likely increase the risk of fire and negatively impact forest productivity, maple syrup production, and the intensity of autumn foliage coloration. Climate and hydrologic changes could have profound effects on forest structure, composition, and ecological functioning in response to the changes discussed here and as described in related articles in this issue of the Journal.

Connecticut, Maine, Massachusetts, New Hampshire,

One-way coupling of an atmospheric and a hydrologic model in Colorado

This paper examines the accuracy of high-resolution nested mesoscale model simulations of surface climate. The nesting capabilities of the atmospheric fifth-generation Pennsylvania State University (PSU)-National Center for Atmospheric Research (NCAR) Mesoscale Model (MM5) were used to create high-resolution, 5-yr climate simulations (from 1 October 1994 through 30 September 1999), starting with a coarse nest of 20 km for the western United States. During this 5-yr period, two finer-resolution nests (5 and 1.7 km) were run over the Yampa River basin in northwestern Colorado. Raw and bias-corrected daily precipitation and maximum and minimum temperature time series from the three MM5 nests were used as input to the U.S. Geological Survey's distributed hydrologic model [the Precipitation Runoff Modeling System (PRMS)] and were compared with PRMS results using measured climate station data. The distributed capabilities of PRMS were provided by partitioning the Yampa River basin into hydrologic response units (HRUs). In addition to the classic polygon method of HRU definition, HRUs for PRMS were defined based on the three MM5 nests. This resulted in 16 datasets being tested using PRMS. The input datasets were derived using measured station data and raw and bias-corrected MM5 20-, 5-, and 1.7-km output distributed to 1) polygon HRUs and 2) 20-, 5-, and 1.7-km-gridded HRUs, respectively. Each dataset was calibrated independently, using a multiobjective, stepwise automated procedure. Final results showed a general increase in the accuracy of simulated runoff with an increase in HRU resolution. In all steps of the calibration procedure, the station-based simulations of runoff showed higher accuracy than the MM5-based simulations, although the accuracy of MM5 simulations was close to station data for the high-resolution nests. Further work is warranted in identifying the causes of the biases in MM5 local climate simulations and developing methods to remove them. ?? 2006 American Meteorological Society.

Journal of Hydrometeorology

Impact of longer-term modest climate shifts on architecture of high-frequency sequences (Cyclothems), Pennsylvanian of midcontinent U.S.A

Pennsylvanian glacioeustatic cyclothems exposed in Kansas and adjacent areas provide a unique opportunity to test models of the impact of relative sea level and climate on stratal architecture. A succession of eight of these high-frequency sequences, traced along dip for 500 km, reveal that modest climate shifts from relatively dry-seasonal to relatively wet-seasonal with a duration of several sequences (???600,000 to 1 million years) had a dominant impact on facies, sediment dispersal patterns, and sequence architecture. The climate shifts documented herein are intermediate, both in magnitude and duration, between previously documented longer-term climate shifts throughout much of the Pennsylvanian and shorter-term shifts described within individual sequences. Climate indicators are best preserved at sequence boundaries and in incised-valley fills of the lowstand systems tracts (LST). Relatively drier climate indicators include high-chroma paleosols, typically with pedogenic carbonates, and plant assemblages that are dominated by gymnosperms, mostly xerophytic walchian conifers. The associated valleys are small (<2 km wide and <11 m deep) and are commonly filled with locally derived limestone clasts, which indicate small drainage networks. Relatively wetter climates are indicated by locally preserved paleosols exhibiting both low-chroma and high-chroma variants, common coals, and plant assemblages that are dominated by fern foliage, seed ferns, and sphenopsids. The incised-valley fills associated with these facies are large (>4 km wide and >20 m deep), and dominated by quartz sandstones derived from distant source areas, reflecting large drainage networks. Transgressive systems tracts (TST) in all eight sequences gen erally are characterized by thin, extensive limestones and thin marine shales, suggesting that the dominant control on TST facies distribution was the sequestration of siliciclastic sediment in updip positions. Highstand systems tracts (HST) were significantly impacted by the intermediate-scale climate cycle in that HSTs from relatively drier climates consist of thin marine shales overlain by extensive, thick regressive limestones, whereas HSTs from relatively wetter climates are dominated by thick marine shales. Previously documented relative sea-level changes do not track the climate cycles, indicating that climate played a role distinct from that of relative sea-level change. These intermediate-scale modest climate shifts had a dominant impact on sequence architecture. This independent measure of climate and relative sea level may allow the testing of models of climate and sediment supply based on modern systems. Copyright ?? 2005, SEPM.

Journal of Sedimentary Research

Episodic soil succession on basaltic lava fields in a cool, dry environment

Holocene- to late Pleistocene-aged lava flows at Craters of the Moon National Monument and Preserve provide an ideal setting to examine the early stages of soil formation under cool, dry conditions. Transects were used to characterize the amount and nature of soil cover on across basaltic lava flows ranging in age from 2.1 to 18.4 ka. Results indicate that on flows <13 ka, very shallow organic soils (Folists in Soil Taxonomy) are the dominant soil type, providing an areal coverage of up to &sim;25%. On flows &ge;13.9 ka, deeper mineral soils including Entisols, Aridisols, and Mollisols become dominant and the areal extent increases to &ge;95% on flows older than 18.4 ka. These data suggest there are two distinct pedogenic pathways associated with lava flows of the region. The first pathway is illustrated by the younger flows, where Folists dominate. In the absence of a major source of loess, relatively little mineral material accumulates and soils provide only minor coverage of the lava flows. Our results indicate that this pathway of soil development has not changed appreciably over the past &sim;10 ka. The second pedogenic pathway is illustrated by the flows older than 13.9 ka. These flows have been subject to deposition of large quantities of loess during and after the last regional glaciation, resulting in almost complete coverage. Subsequent pedogenesis has given rise to Aridisols and Mollisols with calcic and cambic horizons and mollic epipedons. This research highlights the importance of regional climate change on the evolution of Craters of the Moon soilscapes.

Idaho

Hydrochemical determination of source water contributions to Lake Lungo and Lake Ripasottile (central Italy)

Lake Lungo and Lake Ripasottile are two shallow (4-5 m) lakes located in the Rieti Basin, central Italy, that have been described previously as surface outcroppings of the groundwater table. In this work, the two lakes as well as springs and rivers that represent their potential source waters are characterized physio-chemically and isotopically, using a combination of environmental tracers. Temperature and pH were measured and water samples were analyzed for alkalinity, major ion concentration, and stable isotope (δ2H, δ18O, δ13C of dissolved inorganic carbon, and δ34S and δ18O of sulfate) composition. Chemical data were also investigated in terms of local meteorological data (air temperature, precipitation) to determine the sensitivity of lake parameters to changes in the surrounding environment. Groundwater represented by samples taken from Santa Susanna Spring was shown to be distinct with SO42- and Mg2+ content of 270 and 29 mg/L, respectively, and heavy sulfate isotopic composition(δ34S=15.2 ‰ and δ18O=10‰). Outflow from the Santa Susanna Spring enters Lake Ripasottile via a canal and both spring and lake water exhibits the same chemical distinctions and comparatively low seasonal variability. Major ion concentrations in Lake Lungo are similar to the Vicenna Riara Spring and are interpreted to represent the groundwater locally recharged within the plain. The δ13CDIC exhibit the same groupings as the other chemical parameters, providing supporting evidence of the source relationships. Lake Lungo exhibited exceptional ranges of δ13CDIC (±5 ‰) and δ2H, δ18O (±5 ‰ and ±7 ‰, respectively), attributed to sensitivity to seasonal changes. The hydrochemistry results, particularly major ion data, highlight how the two lakes, though geographically and morphologically similar, represent distinct hydrochemical facies. These data also show a different response in each lake to temperature and precipitation patterns in the basin that may be attributed to lake water retention time. The sensitivity of each lake to meteorological patterns can be used to understand the potential effects from long-term climate variability.

Lake Lungo, Lake Ripasottile

Is April to July runoff really decreasing in the Western United States?

Global warming has been the topic of a great deal of heated discussion and debate in recent years, both in the lay press and in scientific journals. The debate is about whether we are beginning to detect signs of a buildup of greenhouse gases on a global scale. A major part of the debate concerns the possible effects on climate and on the future availability of water resources. The ongoing drought in California has added impetus to the debate, serving notice of the serious consequences of any prolonged decrease in the availability of adequate water supplies. This paper has three primary objectives: (1) To evaluate the ramifications of using fractional runoff rather than total runoff to define trends in runoff; (2) to analyze additional streamflow data for the presence and extent of trends in annual and seasonal runoff volume for the conterminous Western United States; and (3) to examine the influence of the current California drought on indicators of trend.

Conference Paper

Detecting changes in the spatial distribution of nitrate contamination in ground water

Many studies of ground water pollution in general and nitrate contamination in particular have often relied on a one-time investigation, tracking of individual wells, or aggregate summaries. Studies of changes in spatial distribution of contaminants over time are lacking. This paper presents a method to compare spatial distributions for possible changes over time. The large-scale spatial distribution at a given time can be considered as a surface over the area (a trend surface). The changes in spatial distribution from period to period can be revealed by the differences in the shape and/or height of surfaces. If such a surface is described by a polynomial function, changes in surfaces can be detected by testing statistically for differences in their corresponding polynomial functions. This method was applied to nitrate concentration in a population of wells in an agricultural drainage basin in Iowa, sampled in three different years. For the period of 1981-1992, the large-scale spatial distribution of nitrate concentration did not show significant change in the shape of spatial surfaces; while the magnitude of nitrate concentration in the basin, or height of the computed surfaces showed significant fluctuations. The change in magnitude of nitrate concentration is closely related to climatic variations, especially in precipitation. The lack of change in the shape of spatial surfaces means that either the influence of land use/nitrogen management was overshadowed by climatic influence, or the changes in land use/management occurred in a random fashion.

Journal of the American Water Resources Associatio

Succession of desert plants on debris flow terraces, Grand Canyon, Arizona, U.S.A.

Vegetation sampling on 11 debris flow terraces in Grand Canyon National Park, Arizona, U.S.A., showed that plant assemblages changed as age of surface increased. The terraces ranged in age from about 5 to about 3100 years. There were distinct differences among sites in the life history characteristics of the dominant plants. Young terraces (5-55 years) were dominated by short-lived plants that had high reproductive potential. Older surfaces were dominated by species with longer life-spans and lower reproductive potential. Density and cover of long-lived species increased with age of surface; for short-lived plants, density was inversely related to surface age. Species composition was also correlated with site age; however, location, exposure, and other factors ensured that no two debris flows supported identical mixtures of species. Succession on recent Grand Canyon debris flows is driven in part by life-history strategies, particularly life-span and seed dispersal traits, and also by climatic factors, especially those that control germination and establishment of the long-lived dominants.

Journal of Arid Environments

Does Ferocactus wislizeni (Cactaceae) have a between-year seed bank?

Field and laboratory experiments at Tumamoc Hill, Tucson, Arizona, U.S.A., demonstrated that Ferocactus wislizeni, a common perennial cactus in the northern Sonoran Desert, has a between-year seed bank. In laboratory studies, F. wislizeni seeds lost dormancy during storage at room temperature and had a light requirement for germination. Field experiments suggested that as much as 2% of the annual seed crop might escape post-dispersal predation even when unprotected; where suitable safe sites exist, a higher percentage might escape. Germination of seed recovered monthly from above- and below-ground components of an artificial seed bank showed that seeds can survive at least 18 months in and on the soil. Seed banks enable F. wislizeni to take advantage of favorable rains and temperatures throughout the growing season, thus increasing the number of opportunities for germination. Moreover, seed banks enable F. wislizeni to respond hugely when the climate seems especially favorable, thus producing the large cohorts necessary to compensate for high seedling mortality. (C) 2000 Academic Press.

Journal of Arid Environments

Potential effects of climate change on the growth of fishes from different thermal guilds in Lakes Michigan and Huron

We used a bioenergetics modeling approach to investigate potential effects of climate change on the growth of two economically important native fishes: yellow perch ( Perca flavescens ), a cool-water fish, and lake whitefish ( Coregonus clupeaformis ), a cold-water fish, in deep and oligotrophic Lakes Michigan and Huron. For assessing potential changes in fish growth, we contrasted simulated fish growth in the projected future climate regime during the period 2043-2070 under different prey availability scenarios with the simulated growth during the baseline (historical reference) period 1964-1993. Results showed that effects of climate change on the growth of these two fishes are jointly controlled by behavioral thermoregulation and prey availability. With the ability of behavioral thermoregulation, temperatures experienced by yellow perch in the projected future climate regime increased more than those experienced by lake whitefish. Thus simulated future growth decreased more for yellow perch than for lake whitefish under scenarios where prey availability remains constant into the future. Under high prey availability scenarios, simulated future growth of these two fishes both increased but yellow perch could not maintain the baseline efficiency of converting prey consumption into body weight. We contended that thermal guild should not be the only factor used to predict effects of climate change on the growth of a fish, and that ecosystem responses to climate change should be also taken into account.

Lake Huron, Lake Michigan

Controls on mangrove forest‐atmosphere carbon dioxide exchanges in western Everglades National Park

We report on net ecosystem production (NEP) and key environmental controls on net ecosystem exchange (NEE) of carbon dioxide (CO 2 ) between a mangrove forest and the atmosphere in the coastal Florida Everglades. An eddy covariance system deployed above the canopy was used to determine NEE during January 2004 through August 2005. Maximum daytime NEE ranged from −20 to −25 μ mol (CO 2 ) m −2 s −1 between March and May. Respiration (R d ) was highly variable (2.81 ± 2.41 μ mol (CO 2 ) m −2 s −1 ), reaching peak values during the summer wet season. During the winter dry season, forest CO 2 assimilation increased with the proportion of diffuse solar irradiance in response to greater radiative transfer in the forest canopy. Surface water salinity and tidal activity were also important controls on NEE. Daily light use efficiency was reduced at high (>34 parts per thousand (ppt)) compared to low (<17 ppt) salinity by 46%. Tidal inundation lowered daytime R d by ∼0.9 μ mol (CO 2 ) m −2 s −1 and nighttime R d by ∼0.5 μ mol (CO 2 ) m −2 s −1 . The forest was a sink for atmospheric CO 2 , with an annual NEP of 1170 ± 127 g C m −2 during 2004. This unusually high NEP was attributed to year‐round productivity and low ecosystem respiration which reached a maximum of only 3 g C m −2 d −1 . Tidal export of dissolved inorganic carbon derived from belowground respiration likely lowered the estimates of mangrove forest respiration. These results suggest that carbon balance in mangrove coastal systems will change in response to variable salinity and inundation patterns, possibly resulting from secular sea level rise and climate change.

Florida

Source, migration pathways, and atmospheric release of geologic methane associated with the complex permafrost regimes of the outer Mackenzie River Delta, Arctic, Canada

Sources and fluxes of methane to the atmosphere from permafrost are significant but poorly constrained in global climate models. We present data collected from the variable permafrost setting of the outer Mackenzie River Delta, including observations of aquatic methane seepage, core determinations of in situ methane occurrence and seep gas isotope geochemistry. The sources and locations of in situ geologic methane occurrence and aquatic and atmospheric gas release appear to be controlled by the regional geology and permafrost conditions. Where permafrost is >250 m thick, thermogenic gas deposits at depth are isolated by laterally continuous, low permeability ice-bearing sediments with few through-going thawed taliks. Thus, the observed in situ methane and aquatic gas seepage appears to be dominated by microbial methane. In contrast, where permafrost is <80 m thick, taliks are more likely to be through-going, providing permeable conduits from depth and migration pathways for both thermogenic and biogenic gas. Continuous annual fluid sampling of two lakes and a river channel documents aquatic methane flux from microbial sources, more deeply buried thermogenic sources, and mixtures of both. Using estimates of in situ methane concentration from deep core samples and observations of in situ free gas occurrences, we conclude that the reservoir of in situ geologic methane within ice bonded permafrost is substantial and that this methane is presently migrating with ongoing atmospheric release. It is our assessment that the permafrost setting, and processes described are sensitive to future climate change as the permafrost warms.

Mackenzie River Delta

Impact of climate variability on runoff in the north-central United States

Large changes in runoff in the north-central United States have occurred during the past century, with larger floods and increases in runoff tending to occur from the 1970s to the present. The attribution of these changes is a subject of much interest. Long-term precipitation, temperature, and streamflow records were used to compare changes in precipitation and potential evapotranspiration (PET) to changes in runoff within 25 stream basins. The basins studied were organized into four groups, each one representing basins similar in topography, climate, and historic patterns of runoff. Precipitation, PET, and runoff data were adjusted for near-decadal scale variability to examine longer-term changes. A nonlinear water-balance analysis shows that changes in precipitation and PET explain the majority of multidecadal spatial/temporal variability of runoff and flood magnitudes, with precipitation being the dominant driver. Historical changes in climate and runoff in the region appear to be more consistent with complex transient shifts in seasonal climatic conditions than with gradual climate change. A portion of the unexplained variability likely stems from land-use change.

Iowa, Kansas, Minnesota, Missouri, Nebraska, North

Changes in agricultural cropland areas between a water-surplus year and a water-deficit year impacting food security, determined using MODIS 250 m time-series data and spectral matching techniques, in the Krishna river basin (India)

The objective of this study was to investigate the changes in cropland areas as a result of water availability using Moderate Resolution Imaging Spectroradiometer (MODIS) 250 m time-series data and spectral matching techniques (SMTs). The study was conducted in the Krishna River basin in India, a very large river basin with an area of 265 752 km 2 (26 575 200 ha), comparing a water-surplus year (2000–2001) and a water-deficit year (2002–2003). The MODIS 250 m time-series data and SMTs were found ideal for agricultural cropland change detection over large areas and provided fuzzy classification accuracies of 61–100% for various land‐use classes and 61–81% for the rain-fed and irrigated classes. The most mixing change occurred between rain-fed cropland areas and informally irrigated (e.g. groundwater and small reservoir) areas. Hence separation of these two classes was the most difficult. The MODIS 250 m-derived irrigated cropland areas for the districts were highly correlated with the Indian Bureau of Statistics data, with R 2 -values between 0.82 and 0.86. The change in the net area irrigated was modest, with an irrigated area of 8 669 881 ha during the water-surplus year, as compared with 7 718 900 ha during the water-deficit year. However, this is quite misleading as most of the major changes occurred in cropping intensity, such as changing from higher intensity to lower intensity (e.g. from double crop to single crop). The changes in cropping intensity of the agricultural cropland areas that took place in the water-deficit year (2002–2003) when compared with the water-surplus year (2000–2001) in the Krishna basin were: (a) 1 078 564 ha changed from double crop to single crop, (b) 1 461 177 ha changed from continuous crop to single crop, (c) 704 172 ha changed from irrigated single crop to fallow and (d) 1 314 522 ha changed from minor irrigation (e.g. tanks, small reservoirs) to rain-fed. These are highly significant changes that will have strong impact on food security. Such changes may be expected all over the world in a changing climate.

International Journal of Remote Sensing

Problems of snowmelt runoff modelling for a variety of physiographic and climatic conditions

Problems include: a) definition of the spatial and temporal distribution of model input; b) measurement or estimation of snow accumulation, snowmelt, and runoff process parameters for a range of applications and scales; and c) development of accurate short term and long term snowmelt runoff forecasts. Procedures being investigated to solve these problems include: a) integrating conventional and remote-sensing data to improve estimates of input data; b) developing snowmelt process algorithms which have parameters that are closely related to measurable basin and climatic characteristics; and c) updating model paramters and components using measured data or knowledge of past uncertainty. -from Author

Hydrological Sciences Journal

Movements of four native Hawaiian birds across a naturally fragmented landscape

Animals often increase their fitness by moving across space in response to temporal variation in habitat quality and resource availability, and as a result of intra and inter-specific interactions. The long-term persistence of populations and even whole species depends on the collective patterns of individual movements, yet animal movements have been poorly studied at the landscape level. We quantified movement behavior within four native species of Hawaiian forest birds in a complex lava-fragmented landscape: Hawai‛i ‘amakihi Chlorodrepanis virens , ‘oma‘o Myadestes obscurus , ‘apapane Himatione sanguinea , and ‘i‘iwi Drepanis coccinea . We evaluated the relative importance of six potential intrinsic and extrinsic drivers of movement behavior and patch fidelity: 1) forest fragment size, 2) the presence or absence of invasive rats ( Rattus sp.), 3) season, 4) species, 5) age, and 6) sex. The study was conducted across a landscape of 34 forest fragments varying in size from 0.07 to 12.37 ha, of which 16 had rats removed using a treatment-control design. We found the largest movements in the nectivorous ‘apapane and ‘i‘iwi, intermediate levels in the generalist Hawai‛i ‘amakihi, and shortest average movement for the ‘oma‘o, a frugivore. We found evidence for larger patch sizes increasing patch fidelity only in the ‘oma‘o, and an effect of rat-removal increasing patch fidelity of Hawai‛i ‘amakihi only after two years of rat-removal. Greater movement during the non-breeding season was observed in all species, and season was an important factor in explaining higher patch fidelity in the breeding season for ‘apapane and ‘i‘iwi. Sex was important in explaining patch fidelity in ‘oma‘o only, with males showing higher patch fidelity. Our results provide new insights into how these native Hawaiian species will respond to a changing environment, including habitat fragmentation and changing distribution of threats from climate change.

Journal of Avian Biology

Walleye recruitment success is less resilient to warming water temperatures in lakes with abundant largemouth bass populations

Lakes respond heterogeneously to climate, with implications for fisheries management. We analyzed walleye ( Sander vitreus ) recruitment to age-0 in 359 lakes in Wisconsin, USA, to ( i ) quantify the relationship between annual water temperature degree days (DD) and walleye recruitment success and ( ii ) identify the influence of lake characteristics — area, conductivity, largemouth bass ( Micropterus salmoides ) catch rates, and mean DD — on this relationship. The relationship between walleye recruitment and annual DD varied among lakes and was not distinguishable from zero overall (posterior mean = −0.11, 90% CI = −0.34, 0.15). DD effects on recruitment were negative in 198 lakes (55%) and positive in 161 (45%). The effect of annual DD was most negative in lakes with high largemouth bass densities, and, on average, the probability of recruitment was highest in large lakes with low largemouth bass densities. Conductivity and mean DD influenced neither recruitment nor the effect of annual DD. Walleye recruitment was most resilient to warming in lakes with few largemouth bass, suggesting that the effects of climate change depend on lake-specific food-web and habitat contexts.

Wisconsin

Simulating the dynamics of linear forests in Great Plains agroecosystems under changing climates

Most forest growth models are not suitable for the highly fragmented, linear (or linearly shaped) forests in the Great Plains agroecosystems (e.g., windbreaks, riparian forest buffers), where such forests are a minor but ecologically important component of the land mosaics. This study used SEEDSCAPE, a recently modified gap model designed for cultivated land mosaics in the Great Plains, to simulate the effects of climate change on the dynamics of such linear forests. We simulated the dynamics of windbreaks with different initial planting species richness and widths (light changes as the selected resulting factor) using current climate data and nested regional circulation models (RegCMs). Results indicated that (i) it took 70-80 simulation years for the linear forests to reach a steady state under both normal (present-day) and warming climates; (ii) warming climates would reduce total aboveground tree biomass and the spatial variation in biomass, but increase dominance in the linear forests, especially in the upland forests; (iii) linear forests with higher planting species richness and smaller width produced higher aboveground tree biomass per unit area; and (iv) the same species performed very differently with different climate scenarios, initial planting diversity, and forest widths. Although the model still needs further improvements (e.g., the effects of understory species should be included), the model can serve as a useful tool in modeling the succession of linear forests in human-dominated land mosaics under changing climates and may also have significant practical implications in other systems. ?? 2004 NRC Canada.

Canadian Journal of Forest Research