The Prairie Ecozone is the only major region of Canada where drought is a landscape hazard; aridity is linked to soil erosion. Management of prairie ecosystems and soil landscapes therefore requires an understanding of past and future trends and variability in regional aridity. We used instrumental and paleoclimatic records to define a regional baseline for prairie aridity, to evaluate the utility of modern climate normals (i.e. 1961-1990) as a benchmark for future climatic change, and to provide a historical context for a range of General Circulation Model (GCM) forecasts of regional aridity. A warm-dry scenario derived from the Canadian GCM projects a significant increase in the area of subhumid and semiarid climate. Tree rings and diatom-inferred lake salinity record prolonged arid events and show that the climate normal period of 1961-1990 may have been the most benign climate of the past 750 years. The climate of the 20th century was anomalous in terms of the absence of sustained drought. Because both lake and tree-ring analyses recorded an abrupt amelioration of climatic conditions near the start of the instrumental record, we suggest that the immediate impacts of future global warming may be to return the prairies to past conditions in which persistent aridity was recorded for intervals of decades or longer.
L’écozone des Prairies est la seule grande région du Canada où la sécheresse constitue un risque environnemental ; l’aridité est liée à l’érosion du sol. Une compréhension des tendances et des variations de l’aridité passées et à venir des écosystèmes et des sols des Prairies est donc nécessaire pour en rendre possible l’aménagement. Les données paléoclimatiques et climatiques récentes consignées ont servi à définir un niveau d’aridité de base pour la région des Prairies, à évaluer l’utilité des normales climatiques modernes (de 1961 à 1990) comme points de référence pour les changements climatiques à venir et à fournir une base historique à divers scénarios de prévision de l’aridité à l’échelle régionale à partir de modèles de circulation générale. Le scénario du changement vers un climat chaud et sec découlant du modèle canadien prévoit une expansion importante des zones subhumide et semi-aride. Les relevés dendrochronologiques et les données de salinité lacustre obtenues à partir des diatomées signalent l’existence de périodes de sécheresse prolongées et montrent que les normales climatiques enregistrées de 1961 à 1990 pourraient avoir été les plus clémentes des 750 dernières années. Le climat du XXe siècle a été anormal en ce qu’il a été caractérisé par l’absence d’une période de sécheresse prolongée. Les données dendrochronologiques et de salinité lacustre montrent une amélioration soudaine des conditions climatiques vers le début de la période de l’enregistrement des données climatiques ; on croit, de ce fait, que le réchauffement climatique prévu pourrait avoir comme effet immédiat le retour du climat qui prévalait dans le passé dans les Prairies, alors que les périodes d’aridité persistaient pendant des décennies, parfois plus.
Aridez de la praderas canadienses
La ecozona de las praderas es la única región de Canadá donde la sequía representa un riesgo ambiental ya que está relacionada con la erosión del suelo. La gestión del ecosistema de las praderas y del suelo requiere la comprensión de las tendencias presentes y futuras de la variación regional de la aridez. En este estudio se usaron registros instrumentales y paleoclimáticos para definir un umbral regional de la aridez de las praderas; para evaluar la utilidad de valores climáticos normales actuales (i.e. 1961-1990), como un punto de referencia para futuros cambios climáticos; y para proporcionar un contexto histórico para el rango del pronostico de la aridez regional con base al Modelo de Circulación General. Un escenario de calentamiento seco derivado del modelo canadiense pronostica un incremento significativo del área de zonas de clima subhúmedo y semiárido. Los datos obtenidos a partir de anillos de crecimiento y de salinidad lacustre basados en diatomeas dan cuenta de periodos prolongados de sequía y muestran que el clima normal en el periodo 1961-1990 pudo ser el clima más clemente de los últimos 750 años. El clima del siglo XX fue anómalo en términos de la ausencia de una sequía persistente ya que tanto los datos proporcionados por anillos de creciendo de los árboles como los obtenidos del lago sugieren un mejoramiento abrupto de las condiciones climáticas ocurrido próximo al inicio del registro instrumental. Se sugiere que el impacto inmediato del futuro calentamiento global podría reconvertir las praderas a su condición pasada durante la cual una aridez persistente fue registrada por intervalos de décadas o aun por periodos de tiempo más prolongados.
Corps de l’article
Transient climate change simulations from the Canadian Centre for Climate Modelling and Analysis suggest that the largest CO2-induced rise in mean surface temperature in southern Canada will occur in the Interior Plains (Boer et al., 2000). Most general circulation models (GCM) also forecast an increase in winter precipitation for this region, coupled with decreased net soil moisture and water resources in summer, and more frequent extreme departures from mean conditions, including severe drought (Hengeveld, 2000). Unfortunately, long-term forecasting of precipitation is problematic for dry environments (Hostetler and Giorgi, 1993; Herrington et al., 1997), where precipitation is variable and occasionally extreme (Hunter et al., 2002). Subhumid ecosystems and soil landscapes are sensitive to fluctuations in the surface and soil water balances (Lemmen and Vance, 1999). Sustained periods of low precipitation and soil moisture lower the resistance to disturbance such that future climatic variability may exceed thresholds for landscape degradation. The recovery of natural systems to a previous or new equilibrium can take decades or centuries (Wolfe et al., 2001).
Both aridity and drought refer to a negative soil and surface water balance. Aridity is a long-term condition and thus is characteristic of a regional climate, while drought is a short-term (seasons to years) departure from mean conditions and thus represents climatic variability. GCM forecasts of increased aridity and more frequent drought have major implications for rates of erosion and sediment yield in the southern Canadian plains. In general, biophysical systems react to short-term climate variability and to extreme events before they respond to gradual changes in mean conditions (Knox, 1984; Hulme et al., 1999). Extreme and anomalous climate can exceed natural and engineering thresholds beyond which the impacts of climate are more severe. The link between aridity and erosion is well established from paleoenvironmental records (Wolfe et al., 2001; Last and Vance, 2002) and from the monitoring of geomorphic processes and regional sediment yields (Knox, 1984). Less protection of the soil surface from wind and rain splash is generally given or implied as the cause of higher rates of erosion in semiarid landscapes, however, plants also reduce runoff erosion through the transpiration of soil water and the positive influence of stems, roots and organic matter on the infiltration or rain and snowmelt water (Thornes, 1985).
Thus aridity lowers resistance to wind and water erosion (Wheaton, 1984; Wolfe and Nickling, 1997; Campbell, 1998). Wheaton (1990) established an exact rank correlation between the number of days with dust and the Palmer Drought Severity Index for the Prairie Province during 1977-88. Wolfe et al. (1995) described the control of soil moisture on the stability of dunes, with the onset of dune reactivation closely corresponding to the incidence of drought during the past 1000 years. Optical dates defined an interval of minimal dune activity between 560 and 220 BP, implying that drought severity did not exceed a threshold for demise of vegetation cover. In contrast, sand dunes were active between about 1000 and 900 BP and since about AD 1800, following droughts of the 1790s known from tree-ring records (Sauchyn and Skinner, 2001; Wolfe et al., 2001). Current sand dune activity in the driest area of the southern Canadian Prairies may serve as an analogue of the potential response of other eolian landscapes in the region to future global warming (Thorpe et al., 2001).
Long-term changes in net precipitation characteristics may also impact fluvial geomorphic systems in the Canadian Prairies (Campbell, 1998; Ashmore and Church, 2001). The climate sensitivity of prairie fluvial geomorphic systems is recorded in late Holocene sedimentary and morphologic records (O’Hara and Campbell, 1993; Rains et al., 1994) with episodes of channel aggradation and incision roughly coinciding with climate fluctuations recorded in lacustrine records (Lemmen and Vance, 1999). Streams that are fed wholly from sources within the plains may be expected to be more responsive to changes in precipitation than the rivers that are derived from glaciers and snowmelt in the Rocky Mountains (Campbell, 1998; Ashmore and Church, 2001). Because the stream channels of the Interior Plains are incised into Quaternary sediments and poorly consolidated sedimentary strata, slope and bank failure have been integral processes in the Holocene history of glacial meltwater channels and tributary valleys (Sauchyn, 1998). Although drier climate and lesser porewater favour slope stability, slope failure during extreme events may occur with greater frequency.
Because aridity is closely linked to soil erosion and degradation, future management of prairie ecosystems and soil landscapes will require an improved understanding of past and future trends and variability in regional aridity. Land use and management practices have more immediate impacts on rates of surface processes than climate change (Jones, 1993), and thus they have the potential to significantly mitigate the impacts of climate; forest and farm management practices enable adaptation to climate change. An expansion of agriculture into the currently forested margins of the Prairie Ecozone will require assessment of the sensitivity of these soil landscapes to both climate change and an altered vegetation cover. To address these issues, we use both instrumental and paleoclimatic records to define regional baseline aridity for the Prairies, and to evaluate the utility of modern climate normals (i.e. 1961-1990) as a metric of environmental change. The proxy and instrumental climate data are compared to a range of GCM forecasts spanning the next 80 years to evaluate the usefulness of historical climate reconstructions as analogs for future climate change.
Data and Methods
Aridity is simply a lack of moisture expressed in terms of mean annual precipitation and potential evapotranspiration as a ratio (P/PET) or deficit (P-PET). P/PET, the Aridity Index (AI), is the basis for classifying drylands as hyperarid (< 0.05), arid (0.05 ≤ AI < 0.2), semiarid (0.2 ≤ AI < 0.5) and dry subhumid (0.5 ≤ AI < 0.65) and for identifying land at risk of desertification, “land degradation in arid, semiarid and dry subhumid areas resulting mainly from adverse human impact” (Middleton and Thomas, 1992). Wolfe (1997) demonstrated that the distribution of sand dunes on the Canadian plains could be modeled with P/PET. Hogg (1994, 1997) was able to map the distribution of vegetation zones in the Canadian interior using the Climate Moisture Index, P-PET. Because our objective was to describe and compare past and future changes in aridity, rather than evaluate different methods of AI calculation, all estimates of PET were based on the widely used Thornthwaite formula (Dunne and Leopold, 1998):
where PET = potential evapotranspiration (cm month-1), T = mean monthly air temperature (°C) and
Annual heat index, I, was estimated as
With Thornthwaite PET, temperature is the sole measure of the energy available for evapotranspiration. Monthly PET calculated using equation (1) must be adjusted using a correction factor for daylength according to latitude and month. Net radiation and wind data are available for few stations in the Prairie Provinces and thus an energy balance approach to calculating PET would not enable the mapping of aridity. We use Thorntwaite PET for this reason and because it works best in mid-latitude continental climates where air temperature if strongly correlated with net radiation.
To map recent and future aridity we first loaded a gridded database of historical climate observations into a GIS. This database of mean monthly temperature and total precipitation was derived from all of the available monthly climate data in the Canadian Climate Archive for the Prairie Provinces and the surrounding climate regions. The FORTRAN program GRANAL2 used an inverse square distance weighting procedure to interpolate from the station locations to a 50 km grid on a polar stereographic secant projection true at 60o N. The values at the grid points should be viewed as representative of the region within 25 km of the grid point, although the gridded analysis cannot reproduce the high spatial variability associated with precipitation events. Further information on the database and interpolation procedure is available at http://www.cics.uvic.ca/climate/data.htm. In the GIS, a more dense 5 km grid was constructed to smooth the raster boundaries between map units created by classifying the climate data. The aridity index was computed for each original grid point and interpolated to the 5 km grid by Kriging.
Scenarios of future climate change were applied to the gridded temperature and precipitation data for 1961-1990 to produce maps of future aridity centered on the 2020s (2010-2039), 2050s (2040-2069) and 2080s (2070-2099). The IPCC Task Group on Scenarios for Climate Impact Assessment (IPCC-TGCIA, 1999) strongly recommended the use of more than one GCM for any assessment of climate impacts and that the selected GCMs show a range of changes in the key climate variables for the study region. In a recent study of climate change impacts on the island forests of the northern Great Plains, Henderson et al. (2002) thoroughly analyzed the available GCM scenarios and chose three that captured the range of annual temperature and precipitation changes projected for the region. Given the very similar scenario requirements of this very recent study and the forecasting of future aridity attempted here, we accepted the conclusions reached by Henderson et al. (2002) and used the same climate change experiments listed in Table I. HadCM3 B21 forecasts the least increase in temperature and greatest increase in precipitation. CGCM2 A21 forecasts the smallest change in precipitation and the largest increase in temperature. CSIROMk2b B11 is a mid-range scenario. The names of the models refer to the national modeling centres (Table I), the version of the coupled atmosphere-ocean model, and the emission scenario (the last three digits). The greenhouse gas emission scenarios are those most recently recommended by the IPCC from the Special Report on Emission Scenarios (Nakicenovic and Swart, 2000). The scenario data were downloaded from the web site of the Canadian Climate Impacts and Scenarios project (http://www.cics.uvic.ca/scenarios/).
High-resolution records of paleo-aridity were derived from analyses of lake sediments and tree rings. The tree ring chronology is from Pinus contorta (lodgepole pine) in the Cypress Hills of southwestern Saskatchewan and southeastern Alberta (Fig. 1). Standard procedures were used to process the tree cores and disks and measure ring width to within 0.001 mm (Cook and Kairiukstis, 1990). Cross dating of all tree-ring series enabled the detection of anomalous series and missing or false rings and ensured that the correct calendar year was assigned to each tree ring. There are 53 ring-width series with an inter-correlation of 0.481 and a mean sensitivity of 0.231. The raw ring-width data were converted to index chronologies (averaged standardized ring widths) using the program ARSTAN (Cook and Holmes, 1999). The raw data were detrended with a best-fit negative exponential curve or linear regression and filtered to remove first-order autocorrelation. A biweight robust estimate of the mean ring-width index for each year comprised the tree-ring index chronology. The reconstruction of June-July P/PET from standardized residual ring width was validated using a split sample approach; two regression models were calibrated using different halves of the instrumental record and then cross-validated using the remainder of the record.
Chauvin Lake (Fig. 1) is ideally suited for climatic research because it is located within a very small effective drainage. Analysis of aerial photographs between 1940 and the present demonstrate that lake levels have varied in accordance with historical changes in precipitation. Lake sediments cores were collected using standard piston-coring procedures and sediments were cut into 2.5 mm intervals equivalent to an average of 1-3 years of lake history. Fossil diatoms were isolated, identified to species and quantified for alternate sediment levels, and changes in species composition were used in published inference models (Wilson et al., 1996) to reconstruct past changes in lake water salinity at intervals of 3-5 years. Analysis of past changes in chemistry show that lake water salinity accurately records the occurrence of major drought events (Laird et al., 1996), with greatly elevated salinity representing the most arid conditions. Sediment age was determined by analysis of naturally-occurring radioisotopes in recent bulk sediments (210Pb) and individual plant macrofossils (AMS 14C).
Diatom-inferred changes in lake salinity were compared to instrumental climate records from 1910 to 1990 to evaluate the ability of fossil reconstructions to infer changes in aridity. PET data were filtered using a three-year running mean to match the temporal resolution of the fossil time series. All time series were detrended by linear regression, and residual curves compared using correlation analysis. Lags of 1-5 years were applied to the fossil data to evaluate the possibility of delayed hydrologic response to climatic forcing functions.
Results and Discussion
A map of “normal” (i.e. 1961-1990) aridity for the Prairie Provinces (Fig. 2) illustrates that, at 50 km resolution, the “normal” climate of the Prairie Provinces cannot be described as semiarid anywhere (P/PET < 0.5), even though this terminology is often used to describe the southern Prairies. The driest region, the dry subhumid (0.5 < P/PET < 0.65) landscape of southwestern Saskatchewan and southeastern Alberta, corresponds to the mixed grass prairie ecoregion. The moist subhumid Cypress Hills are near the geographic center of this driest region. Analysis of the longest instrumental record (Medicine Hat, Alberta) showed that annual P/PET varied from 0.34 (1943) to 1.24 (1927), with 26 years of semiarid conditions (Fig. 3). On average, these dry events occurred approximately every third year, although only eight semiarid years were seen after 1950.
This large inter-annual variation in the aridity index translates into varying geographic extent of annual drought. Each documented drought had a unique geographic extent and intensity, although semiarid conditions were consistently recorded in southeastern Alberta and southwestern Saskatchewan (Fig. 4). The “normal” dry subhumid climate of this region reflects the frequency of low annual precipitation and P/PET. The most extensive drought was in 1961.
Scenarios of future aridity centered on the 2020s, 2050s and 2080s are mapped for the Prairie Provinces in Figures 5-7 for the CGCM2 A21 (warm-dry), CSIROMk2b B11 (mid-range), HadCM3 B21 (cool-wet) climate change scenarios. Under the CGCM2 A21 scenario, a small area of continuously semiarid climate (P/PET < 0.5) appears in western Saskatchewan by the 2020s and expands significantly by the 2050s. In contrast, forecasts using the CSIROMk2b B11 mid-range scenario (Fig. 6) suggest that perpetually semiarid conditions will be evident only by 2080, but that a dry subhumid climate will extend to southwestern Manitoba. The HadCM3 B21 scenario indicates the least climate change, with a marginal increase in the area of subhumid climate and maximum aridity in mid century. Taken together, these forecasts suggest a general increase in the geographic extent of dry conditions, with possible risks of desertification in the more extreme scenarios (CGCM2 A21).
Analysis of tree-ring indices of past climate suggest that very dry conditions have been common on the Prairies during the past 300 years (Fig. 8). The residual index tree-ring widths for Pinus contorta from the Cypress Hills are significantly correlated (p < 0.01) with June, July, annual and plant year (August to July) precipitation for Medicine Hat, Alberta for the period 1900-1999 (Table II). The correlation of ring width with potential evapotranspiration (PET) is weak (p > 0.05) because Thornthwaite PET is a function of temperature and heat has both positive and negative effects on tree growth depending on the ambient air temperature and antecedent moisture conditions. The negative influence would be confined to warm dry days. Thus the aridity (P/PET) signal in the tree rings mostly reflects large interannual variability in annual, and especially summer, precipitation. Not only is the response to precipitation direct and unambiguous, but amounts of rain and snow are much more variable from year to year than mean temperature and in turn evapotranspiration.
The reconstruction of summer (June-August) aridity for 1723-2000 is plotted in Figure 8. The similar standard errors of estimation (SEe) and validation (RMSEv) in Table III suggest an acceptable regression model. Events of high aridity occurred regularly during the past three centuries, with major droughts indicated for 1768-1782, 1786-1796 and 1851-1864. In addition, analysis of the time series suggested that mean interannual variance (as standard deviation) was greater during the instrumental period (post-1884; SD = 0.080) than in the preceding period (1723-1883; SD = 0.063). This prominent increase in inter-annual variability may be related to a mid-19th century shift in the Pacific Decadal Oscillation (D’Arrigo et al., 2001). The consecutive years of summer aridity highlighted in Figure 8 are further evidence of prolonged aridity prior to Euro-Canadian settlement of the western Canadian plains and the instrumental observation of climate. The impacts of these droughts are described in Wolfe et al. (2001) and Sauchyn et al. (2003). Their relevance here is as possible analogues of anomalous precipitation expected with global warming in contrast to the climate of the 20th century, and particularly 1961-1990, the basis for our definition of the normal climate.
Paleoclimatic reconstructions based on the 2000-year salinity record from Chauvin Lake, Alberta, also suggest that past drought severity was greater than that recorded during the instrumental period. Validation using modern climatic records showed that detrended, diatom-inferred salinity was significantly and positively correlated (p < 0.05) with PET during 1910-1990 (Fig. 9). Correlations were also positive and significant at lags of one or two years. Using this observation, and the drought of 1988-1989 as a benchmark, analysis of past salinities revealed that climatic conditions were continuously arid during the 4th to 8th centuries and regularly dry during CE 1300-1750 (Fig. 10). Severe events were also recorded during 1860-1890, an interval of broadscale droughts in central North America. Salinities declined substantially during the 20th century, reaching values during the 1960s and 1970s that were the lowest of the past 750 years.
Other recent studies of the past hydroclimate of the northern Great Plains (Woodhouse and Overpeck, 1998; St. George and Nielsen, 2002) have similarly concluded that the droughts of the 20th century were of comparatively short duration. Much of the previous paleo drought research has focused on the United States (Cook et al., 1999). Our results fit with the timing of droughts (e.g., 1790s, 1860s, 1880s) in eastern Montana, immediately south of our study area, but generally do not correlate with tree-ring data from the eastern plains (St. George and Nielsen, 2002). The lower frequency variation in hydroclimate (i.e., aridity), on the other hand, tends to be consistent over a larger area. Our tree-ring and diatom data share decadal to century scale variability with dendroclimatic (Case and MacDonald, 1995; St. George and Nielsen, 2002) and lake sediment (Laird et al., 2003) records from across the northern Great Plains. Common periods of sustained aridity encompass much of the 16th and 18th centuries.
Conclusions and Implications
Analysis of climate normals (1961-1990) revealed that a large part of the southern Canadian Plains usually has dry subhumid conditions, and that P/PET was rarely < 0.5 (semiarid) when analyzed on a 50 x 50 km grid. In contrast, semiarid conditions were common, on average, one of every three years during the instrumental record, although the high degree of spatial variability suggested that any given site could experience a wide range of climate, from semiarid (P/PET < 0.5) to humid conditions (> 1.0). High temporal variability was evident in the tree-ring (Fig. 8) and lake sediment (Fig. 10) reconstructions, however, these records also demonstrated that more prolonged, intensely arid events were common in the past, and showed further that modern climate normals are not representative of the full range of potential climatic conditions, even in the absence of global warming.
The increasing aridity projected by the GCM scenarios represents a higher frequency of dry years over a larger area. As a result of an increase in mean annual temperature and, in turn, potential evapotranspiration, with even the coolest and wettest climate change scenario, there is a marginal increase in the area of subhumid climate in the Prairie Provinces by the 2080s. The Canadian GCM forecasts the least increase in precipitation and the largest increase in temperature and therefore an approximately 50 % increase in the area of subhumid climate and a significant area of semiarid climate by the 2050s. Past droughts may provide important benchmarks to evaluate potential impacts of increasing aridity forecast by GCMs. For example, comparison of climatic conditions observed in the pre-drought year 1987 shows a marked similarity to conditions forecast by the CGCM2 A21 scenario for 2050, in terms of the geographic extent of semiarid and dry subhumid climate (Fig. 11a). Dry conditions during 1987 are viewed as partly responsible for establishing more extensive drought during 1988-1989 (Fig. 4) and suggest that future increases in aridity may generate conditions that favour the regular development of large-scale droughts.
The implication of more frequent and widespread droughts is that they place a larger area at risk of desertification more often unless counteracting management strategies are enacted. The critical variable will be the duration of the semiarid conditions in any given area. The scenario of increased climate extremes includes unusually wet years during which social and biophysical systems may recover from drought. If aridity is prolonged, however, the resistance of systems could drop below critical thresholds, such that they will not recover during the next wet year. Based on evidence derived from climatic proxies, our analysis suggests that the climate of the 20th century was anomalous in terms of the absence of sustained drought. Further, comparison of instrumental records and climatic reconstructions suggests that the climate normal period of 1961-1990 may have been among the most benign of the past 750 years. Possibly, the climate of 1921-1950 (Fig. 11b) is more characteristic of the past and future climate of the Prairie Provinces than the “normal” climate of 1961-1990. Because both lake and tree-ring analyses recorded an abrupt amelioration of climatic conditions at near the start of the instrumental record, we suggest that the immediate impacts of future global warming may be to return the prairie environment to past conditions in which persistent aridity was recorded for intervals of decades or longer.
The results presented here are the product of research funded by the Natural Sciences and Engineering Research Council, the Prairie Adaptation Research Collaborative and the Canadian Foundation for Climate and Atmospheric Sciences. We thank Sam Kennedy and Brent Joss for compilation of the maps; Antoine Beriault, Jennifer Stroich and Julie Frischke for processing of the tree-ring data; Sybille Wunsam for diatom enumeration, Jim Rusak for diatom data management, and Brian Cumming for salinity reconstructions. We also thank Matthew Boyd and an anonymous referee for reviews of the manuscript.
- Ashmore, P. and Church, M., 2001. The Impact of Climate Change on Rivers and River Processes in Canada. Geological Survey of Canada, Ottawa, Bulletin 555, 58 p.
- Boer, G.J., Flato, G. and Ramsden, D., 2000. A transient climate change simulation with greenhouse gas and aerosol forcing: projected climate for the 21st century. Climate Dynamics, 16: 427-450.
- Campbell, I.A., 1998. Fluvial Processes, p. 39-48. In D.S. Lemmen, R.E. Vance, I.A. Campbell, P.P.David, D.J. Pennock, D.J. Sauchyn and S.A. Wolfe, eds., Geomorphic Systems of the Palliser Triangle: Description and Response to Changing Climate, Geological Survey of Canada, Ottawa, Bulletin 521, 72 p.
- Case, R.A. and MacDonald, G.M., 1995. A dendroclimatic reconstruction of annual precipitation on the western Canadian prairies since A.D. 1505 from Pinus flexilis James. Quaternary Research, 44: 267-275.
- Cook, E.R. and Holmes, R.L., 1999. Users Manual for Program ARSTAN. Laboratory of Tree-Ring Research, University of Arizona, Tucson, 16 p.
- Cook, E.R. and Kairiukstis, L.A. (eds.), 1990. Methods of Dendrochronology — Applications in the Environmental Sciences. Kluwer, Dordrecht, 394 p.
- Cook, E.R., Meko, D.M., Stahle, D.W. and Cleaveland, M.K., 1999. Drought reconstructions for the continental United States. Journal of Climate, 12: 1145-1162.
- D’Arrigo, R., Villalba, R. and Wiles, G., 2001. Tree-ring estimates of Pacific decadal climate variability. Climate Dynamics, 18: 219-224.
- Dunne, T. and Leopold, L.B., 1998. Water in Environmental Planning. W.H. Freeman, San Francisco, 818 p.
- Henderson, N., Hogg, E., Barrow, E. and Dolter, B., 2002. Climatic Change Impacts on the Island Forests of the Great Plains and the Implications for Nature Conservation Policy. Report to the Prairie Adaptation Research Collaborative (August, 2002), Regina, 116 p.
- Hengeveld, H., 2000. Projections for Canada’s Climate Future. Environment Canada, Climate Change Digest 00-01, Ottawa, 27 p.
- Herrington, R., Johnson, B. and Hunter, F., 1997. Responding to Global Climate Change in the Prairies - Volume III of the Canada Country Study: Climate Impacts and Adaptation. Environment Canada, Regina, 44 p.
- Hogg, E., 1994. Climate and the southern limit of the western Canadian boreal forest. Canadian Journal of Forest Research, 24: 1835-1845.
- ____ 1997. Temporal scaling of moisture and the forest-grassland boundary in western Canada. Agricultural and Forest Meteorology, 84: 115-122.
- Hostetler, S.W. and Giorgi, F., 1993. Use of output from high-resolution atmospheric models in landscape-scale hydrologic models: An assessment. Water Resources Research, 29: 1685-1695.
- Hulme, M., Barrow, E.M., Arnell, N.W., Harrison, P.A., Johns, T.C. and Downing, T.E., 1999. Relative impacts of human-induced climate change and natural climate variability. Nature, 397: 688-691.
- Hunter, F.G., Donald, D.B., Johnson, B.N., Hyde, W.D., Hanesiak, J.M., Kellerhals, M.O.B., Hopkinson, R.F. and Oegema, B.W., 2002. The Vanguard Torrential Storm (meteorology and hydrology). Canadian Water Resources Journal, 27: 213-228.
- IPCC-TGCIA, 1999. Guidelines on the Use of Scenario Data for Climate Impact and Adaptation Assessment. Version 1. Prepared by T.R. Carter, M. Hulme and M. Lal, Intergovernmental Panel on Climate Change, Task Group on Scenarios for Climate Impact Assessment, 69 p.
- Jones, D.K.C., 1993. Global warming and geomorphology. The Geographical Journal, 159: 124-130.
- Knox, J.C., 1984. Fluvial responses to small scale climate changes, p. 318-342. In J.E. Costa and P.J. Fleisher, eds., Developments and Applications of Geomorphology. Springer-Verlag, Berlin, 372 p.
- Laird, K.R., Cumming, B.F., Wunsam, S., Rusak, J.A., Oglesby, R.J., Fritz, S.C. and Leavitt, P.R., 2003. Lake sediments record large-scale shifts in moisture regimes across the northern prairies of North America during the past two millennia. Proceedings of the National Academy of Science, 100: 2483-2488.
- Laird, K.R., Fritz, S.C., Maasch, K.A. and Cumming, B.F., 1996. Greater drought intensity and frequency before AD 1200 in the northern Great Plains, USA. Nature, 384: 552-554.
- Last, W.M. and Vance, R.E., 2002. The Holocene history of Oro Lake, one of Western Canada’s longest continuous lacustrine records. Sedimentary Geology, 148: 161-184.
- Lemmen, D.S. and Vance, R.E., 1999. An overview of the Palliser Triangle Global Change Project, p. 7-22. In D.S. Lemmen and R.E. Vance, eds., Holocene Climate and Environmental Change in the Palliser Triangle: A Geoscientific Context for Evaluating the Impacts of Climate Change on the Southern Canadian Prairies. Geological Survey of Canada, Ottawa, Bulletin 534, 295 p.
- Middelton, N. and Thomas, D.S.G., 1992. World Atlas of Desertification. United Nations Environment Program, Edward Arnold, London, 69 p.
- Nakicenovic, N. and Swart, R., 2000. IPCC Special Report on Emission Scenarios. Cambridge University Press, Cambridge, 612 p.
- O’Hara, S L. and Campbell, I.A., 1993. Holocene geomorphology and stratigraphy of the Lower Falcon Valley, Dinosaur Provincial Park, Alberta, Canada. Canadian Journal of Earth Sciences, 30: 1846-1852.
- Rains, R.B., Burns, J.A. and Young, R.R., 1994. Postglacial alluvial terraces and an incorporated bison skeleton, Ghostpine Creek, southern Alberta. Canadian Journal of Earth Sciences, 31: 1501-1509.
- Sauchyn, D.J., 1998. Mass wasting processes, p. 48-54. In D.S. Lemmen, R.E. Vance, I.A. Campbell, P.P. David, D.J. Pennock, D.J. Sauchyn and S.A. Wolfe, eds., Geomorphic Systems of the Palliser Triangle: Description and Response to Changing Climate. Geological Survey of Canada, Ottawa, Bulletin 521, 72 p.
- Sauchyn, D.J. and Skinner, W.R., 2001. A proxy PDSI record for the southwestern Canadian plains. Canadian Water Resource Journal, 26: 253-272.
- Sauchyn, D.J., Stroich, J. and Beriault, A., 2003. A paleoclimatic context for the drought of 1999-2001 in the northern Great Plains of North America. The Geographical Journal, 169: 158-167.
- St. George, S. and Nielsen, E., 2002. Hydroclimatic change in southern Manitoba since A.D. 1409 inferred from tree rings. Quaternary Research, 58: 103-111.
- Thornes, J.B., 1985. The ecology of erosion. Geography, 70: 222-234.
- Thorpe, J., Wolfe, S., Campbell, J., Leblanc, J. and Molder, R., 2001. An Ecoregion Approach for Evaluating Land Use Management and Climate Change Adaptation Strategies on the Sand Dune Areas in the Prairie Provinces. Saskatchewan Research Council, Saskatoon, Publication 11368-1E01, 293 p.
- Wheaton, E.E., 1984. The Impacts of Climatic Change: The Potential Wind Erosion of Soil in the Canadian Prairies. The Case of Saskatchewan. Saskatchewan Research Council, Saskatoon, Publication E-906-29-D-84, 43 p.
- ____ 1990. Frequency and severity of drought and dust storms. Canadian Journal of Agricultural Economics, 38: 695-700.
- Wilson, S., Cumming, B. and Smol, J.P., 1996. Assessing the reliability of salinity inference models from diatom assemblages: An examination of a 219 lake dataset from western North America. Canadian Journal of Fisheries and Aquatic Science, 53: 1580-1594.
- Wolfe, S.A., 1997. Impact of increased aridity on sand dune activity in the Canadian Prairies. Journal of Arid Environments, 36: 421-432.
- Wolfe, S.A., Huntley, D.J., David, P.P., Ollerhead, J., Sauchyn, D.J. and Macdonald, G.M., 2001. Late 18th century drought-induced sand dune activity, Great Sand Hills, southwestern Saskatchewan. Canadian Journal of Earth Sciences, 38: 105-117.
- Wolfe, S.A., Huntley, D.J. and Ollerhead, J., 1995. Recent and late Holocene sand dune activity in southwest Saskatchewan. Geological Survey of Canada, Ottawa, Current Research, 1995-B: 131-140.
- Wolfe, S.A. and Nickling, W.G., 1997. Sensitivity of Eolian Processes to Climate Change in Canada. Geological Survey of Canada, Ottawa, Bulletin 421, 30 p.
- Woodhouse, C.A. and Overpeck, J.T., 1998. 2000 years of drought variability in the central United States. Bulletin of the American Meteorological Society, 79: 2693-2714.