Printer Friendly

Climatology of the forest-tundra ecotone at a maritime subarctic-alpine site, Mealy Mountains, Labrador.

ABSTRACT. Climatological investigations were conducted from 2001 to 2009 in the central Mealy Mountains, Labrador (53.6[degrees] N), as part of interdisciplinary research on tree line ecology and climate change. The aim was to describe local climatic and edaphic variables along an altitudinal gradient from the closed forest edge to the alpine tundra and to relate recent changes in regional climate to potential changes in the forest tundra ecotone. Results show relatively warm, moist summers and cool winters, with abundant precipitation and moderate to gale-force winds. At ca. 600 m a.s.l., the tree line was characterized by 694 [+ or -] 85 growing degree-days and a mean July air temperature of 12.9 [+ or -] 0.8[degrees]C. Growing season soil temperatures of ca. 8[degrees]C were similar across the forest-tundra ecotone, but their seasonal regimes differed among subzones. Soil nutrient fluxes showed some variation but no consistent pattern that would suggest nutrients as a limiting factor. Snow depth in the forest-tundra subzone was more variable than in the forest, indicating that microtopography is an important factor for tree survival there in winter. Comparisons of the field data with long-term regional climatic and hydrographic records show that conditions have become warmer and drier in the most recent decade compared with the previous half-century or more. Concurrent vegetation studies indicate that changes expected with a warming climate are already occurring.

Key words: Labrador, subarctic-alpine climate, forest-tundra ecotone, growing degree-days, soil temperature, nutrient flux, water balance, climate warming

RESUME. De 2001 a 2009, des etudes climatologiques ont ete effectuees dans les monts centraux Mealy, au Labrador (53,6[degrees] N) dans le cadre d'une recherche interdisciplinaire sur l'ecologie de la limite forestiere et le changement climatique. Ces etudes avaient pour but de decrire les variables climatiques et edaphiques locales le long d'une pente altitudinale allant de la lisiere de la foret dense jusqu'a la toundra alpine. Elies avaient egalement pour but de faire le lien entre les recents changements enregistres dans le climat regional et les changements susceptibles de se produire dans l'ecotone de la toundra forestiere. Les resultats montrent des etes humides et relativement chauds ainsi que des hivers frais, avec des precipitations abondantes et des vents allant de moderes a coups de vent. A environ 600 m ASL, la limite forestiere etait caracterisee par 694 [+ or -] 85 degres-jours de croissance et une temperature de l'air moyenne de 12,9 [+ or -] 0,8 [degrees]C en juillet. En saison de croissance, les temperatures du sol d'environ 8 [degrees]C etaient semblables a la grandeur de Tecotone de la foret-toundra, bien que les regimes saisonniers differaient d'une sous-zone a l'autre. Les flux de nutriments du sol affichaient une certaine variation, mais il n'existait aucun modele uniforme laissant entrevoir que les nutriments pouvaient etre un facteur limitant. L'epaisseur de la neige dans la sous-zone foret-toundra variait plus que dans la foret, ce qui indique que la microtopographie est un facteur important pour la survie des arbres dans cet endroit-la l'hiver. La comparaison des donnees obtenues sur le terrain aux donnees climatiques et hydrographiques regionales a long terme laisse entrevoir que les conditions se sont rechauffees et assechees au cours de la toute derniere decennie comparativement au demi-siecle precedent ou plus. Les etudes concurrentes sur la vegetation indiquent que les changements qui etaient escomptes de meme que le rechauffement du climat se manifestent deja.

Mots cles: Labrador, climat alpin subarctique, ecotone foret-toundra, degres-jours de croissance, temperature du sol, flux des nutriments, equilibre hydrique, rechauffement du climat

Traduit pour la revue Arctic par Nicole Giguere.


Studies have shown that a simple model of poleward and upslope movement of the forest limit under a warming climate is inadequate to explain the complexity of boreal forest ecotone responses to a changing climate (Korner, 1998; Sveinbjornsson et al., 2002; Gamache and Payette, 2005; Holtmeier and Broil, 2005, 2010; Harsch et al., 2009). The measure of growing degree-days (GDD), calculated using the mean daily surface air temperature (TA) threshold of 5[degrees]C, is frequently used in bioclimatic studies, including those of high-altitude black spruce (Sirois et al., 1999; Sirois, 2000). Meunier et al. (2007) analyzed the thermal requirements for black spruce (Picea mariana), the dominant tree line species in the eastern Canadian boreal region, in relation to GDD. They found a threshold for seed viability of approximately 800 GDD and mapped this bioclimatic threshold for the Canadian boreal region in terms of recent and projected future climates. On the basis of the observations that root-zone temperature is strongly correlated to tree canopy temperatures and that both are important physiological criteria for trees, Korner and Paulsen (2004) analyzed 10 cm soil temperatures at 46 alpine tree line sites between latitude 68[degrees] N and 42[degrees] S. They found growing season temperatures averaged 6.7 [+ or -] 0.8[degrees]C overall and concluded that soil temperature was a better predictor of tree line on a global scale than either growing season length or GDD. Although their study included only three subarctic or boreal sites, it offers a global indicator of tree line climate for more intensive testing at the regional and local levels. Thawing degree-days (TDD), the cumulative number of daily heat units above 0[degrees]C, is a bioclimatic measure used in relation to vegetation zonation and climate change impacts on vegetation in Arctic tundra (Edlund and Alt, 1989; Walker, 2000; Hollister et al., 2005), the boreal forest-tundra transition zone (Rouse, 1991), and subarctic alpine regions (Callaghan et al., 2010). TDD can be a useful indicator of the thermal potential for forest species expansion into tundra areas (Sundqvist et al., 2008).

Along with growing season conditions, winter climate can also constrain tree line expansion through damage to shoots and seedlings by wind abrasion and desiccation (Kullman, 1998; Cairns, 2001), and the krummholz form of tree line is generally associated with such winter stresses (e.g., Harsch et al., 2009). No single climatic variable adequately captures winter severity, although mean winter TA is normally used as such a measure. Freezing degree-days (FDD), the cumulative degree-days below 0[degrees]C, provide a winter counterpart to GDD that has some explanatory value in relation to tree line environments (Gamache and Payette, 2004; Palmer et al., 2012). The fact that snow provides protection to seedlings and shoots in winter but may also cause damage by excessive weight and delay growth in spring (Lavoie and Payette, 1992; Gamache and Payette, 2004) suggests that snow cover should be considered in assessing tree line climates.

Hare (1951) first argued that snowfall in central Labrador was primarily the result of tropospheric flow from the south and west, with surface convergence in migrating low-pressure systems, even though surface winds into these lows at the Labrador Sea were easterly. Regionally, large precipitation amounts over the Labrador plateau were subsequently confirmed through analysis of runoff from major rivers in the region (Hare, 1980). A recent review of temperate and boreal rainforests in eastern Canada (Clayden et al., 2011) identified areas south of the Mealy Mountains as "per-humid" boreal forest, i.e., forest with annual precipitation greatly exceeding evapotranspiration, and noted that adjacent montane areas below the resolution of their analysis might have even higher precipitation amounts.

Soil development is a reflection of climate and may be expected to vary across a climatic gradient. Conifer species require soil macronutrients, particularly nitrogen (N), phosphorous (P), and potassium (K), in relatively large amounts, as well as small amounts of certain micronutrients, (Brady and Weil, 1996; Grossnickle, 2000), and differences across the ecotone may indicate a limiting nutrient factor in tree growth.

Research on tree line and tundra ecology has been conducted since 2001 in the Mealy Mountains, a subarctic-alpine location in south-central Labrador (Fig. 1). This work included obtaining baseline biophysical data to test hypotheses regarding tree line dynamics. The most recent regional climate records indicated that a warming trend was in progress (EC, 2013a); therefore, the research context was ecosystem change in a warming climate. Initial results demonstrated that transplanted conifer seedlings can survive in the Mealy Mountains above the altitude at which trees now occur, pointing to the possibility of future upward expansion of the forest limit (Munier et al., 2010). In 2007-09, this research was intensified and coordinated with that of other research groups in an International Polar Year project focusing on the circumpolar forest-tundra ecotone (Hofgaard and Harper, 2011). These studies have further elucidated the nature of the forest-tundra ecotone, both in the study area (Trindade et al., 2011; Wheeler et al., 2011; Trant et al., 2011; Jameson, 2012; Trant, 2013) and more generally in relation to local topography (and implicitly, microclimate) at this and other locations (Harper et al., 2011). Dendrochronological analysis (Trindade et al., 2011) showed that growth of conifer species is not limited by precipitation, supporting the case for abundant precipitation in the Mealy Mountains region.

Here we apply results from investigations of the climate of the Mealy Mountains (Chan, 2010) toward understanding the constraints on forest-tundra ecotone change. The forest-tundra ecotone is defined here as "the transition zone between forest and tundra at high elevation or latitude" (Harper et al., 2011:480). We examine relevant abiotic factors, including air and soil temperature and moisture, precipitation, and soil nutrients (Timoney, 1995; Grace et al., 2002), and compare the results with biogeographic generalizations regarding tree line climate (e.g., Korner and Paulsen, 2004; Meunier et al., 2007). The overall objectives are 1) to provide a quantitative description of the local climate in the context of the long-term regional climate and climate change, 2) to assess the spatial variability of specific climatic and edaphic variables along an altitudinal gradient from the closed forest edge to the alpine tundra as factors determining tree line character, and 3) to relate observed changes in climate to reported and potential changes in the forest-tundra ecotone.


Central Labrador is characterized by a broad glaciated plateau that rises westward from the Labrador Sea coast to elevations of more than 500 m above sea level (a.s.l.). The coast is incised by many narrow inlets, the most prominent being Lake Melville, a saltwater body that extends some 170 km inland from the Labrador Sea. The Mealy Mountains rise abruptly from a narrow coastal plain south of Lake Melville (Fig. 1). The local bedrock geology consists mainly of granite and anorthosite (Wardle et al., 1997). With summit elevations approaching 1100 m a.s.l., these highly glaciated mountains have a different climatic regime compared with the Labrador Sea coast and the interior lowlands adjacent to Lake Melville.

The bioclimatic zones of Labrador can be represented in terms of tundra, forest-tundra, open woodland, and closed forest (Hare and Ritchie, 1972). The northernmost limits of scattered, coniferous trees, predominantly black spruce and white spruce (Picea glauca), define the boundary between tundra and forest-tundra on the regional scale. In northern Quebec and Labrador, this limit dips southward from ca. 58.5[degrees] N at Ungava Bay to about 55[degrees] N over the northern Labrador Plateau before returning northward to ca. 58[degrees] N at the Labrador Sea (Hare and Ritchie, 1972; Elliott and Short, 1979; Payette, 2007).

According to Meades (2007), the central Mealy Mountains are the most southerly outlier of the High Subarctic Tundra ecoregion, fringed by elements of the High Boreal Lorest ecoregion and abutted to the south by extensive peatlands of the Eagle River Plateau. Altitude results in a local bioclimatic gradient that encompasses a range of vegetation zones, from isolated patches of closed forest to alpine tundra. At latitude 53.6[degrees] N, the Mealy Mountains are more than 100 km south of the present latitudinal tree limit, but 200 km north of the 800 GDD isoline mapped by Meunier et al. (2007) as the northern limit of black spruce seedling viability.

The study area (53.6[degrees] N, 58.8[degrees] W) is centered in an eastward-trending valley below an unnamed 1057 m a.s.l. summit, and it extends from alpine tundra at 1000 m a.s.l. to forest east of a chain of headwater lakes at ca. 500 m a.s.l. (Pig. 1). A complex glacial history is evident in the surface features and geomorphology of the valley (Fig. 2), where extensive glacial and glaciofluvial deposits, glacial erratics, and ground moraine overlay bedrock. On the basis of regional temperatures, this area is considered to be in the zone of sporadic discontinuous permafrost (Smith and Riseborough, 2002).

The forest in this area is mainly open lichen woodland made up of four boreal conifer species: black spruce, white spruce, balsam fir (Abies balsamea), and eastern larch (Larix laricina). Surveys in the study area (Trant, 2013) found black spruce to be dominant in the closed forest and krummholz, while in the transition between those zones eastern larch is dominant, followed by balsam fir. Patches of shrub birch (Betula spp.) and mountain alder (Alnus crispa) occur in open areas, becoming dense in riparian floodplains. Locally, the forest-tundra ecotone extends from the highest open woodland patches at ca. 550 m a.s.l. to the limit of tree species in krummholz form near 750 m a.s.l. (Trant et al., 2011). Above that level, where soils and substrate permit, the vegetation consists of typical Arctic-alpine low shrubs, herbaceous plants, moss and lichen (Munier et ah, 2010). Summits are predominately bare bedrock and boulder fields. Soils are observed to be generally shallow and include cryogenic mineral soils, poorly developed podzols, and some hygric organic soils. Except in well-drained tills and gravels, soils are assumed to be near saturation for most of the year, as is the case over much of northern and central Labrador (Hare, 1980).

Although there are a few earlier records, the regional instrumental record of climate in central Labrador begins ca. 1940 with meteorological stations established at Cartwright in 1936 and Goose Bay in 1941 (Fig. 1). Warming over much of the Arctic in the mid-to-late 20th century (ACIA, 2004) did not initially include the eastern Canadian Arctic and Labrador. From the mid-1970s through the early 1990s, a pronounced winter cooling occurred over northeastern Canada, including Labrador (Banfield and Jacobs, 1998; Feldstein, 2002; Chouinard et ah, 2007); however, by 2000 a warming trend in all seasons was underway in this region (Vose et ah, 2005). As is generally the case around the North Atlantic, Labrador temperature records show decadal-scale variability that is strongly associated with the North Atlantic Oscillation and Arctic Oscillation (NAO/ AO) (D'Arrigo et ah, 1993, 2003). Over the 20th century, a significant increase in annual precipitation in Labrador was accompanied by an increase in the rain-to-snow ratio, particularly in the spring (Zhang et ah, 2000).

Climate variability is also reflected in the regional hydrology. The Eagle River has a drainage area of 10900 [km.sup.2] and flows south and eastward from the Mealy Mountains to the Labrador Sea (Fig. 1). Gauges near its mouth have produced a mostly uninterrupted record from 1969 to the present of the combined runoff from the Mealy Mountains and the adjacent Eagle River plateau. In Whitfield and Cannon's (2000) analysis of hydrologic records across Canada, the Eagle River fell into the class of stations characteristic of the boreal forest zone, with low winter flow, a snowmelt-driven spring flood, and a relatively high sustained summer flow diminishing into fall. Comparison of differences in flow between the decades 1976-85 and 1986-95 showed the Eagle River to be in a group that experienced lower winter temperatures, slight decreases in precipitation, and decreased flow in the latter period (Whitfield and Cannon, 2000). Our investigations provide an opportunity to extend the analysis of the Eagle River record, while providing new information on the hydroclimatology of the Mealy Mountains.


Field studies were carried out from July 2001 through September 2009, and in 2008 and 2009, surveys and measurements were standardized according to IPY-PPSA protocols (Hofgaard and Rees, 2008; Harper et al., 2011). For purposes of this analysis, the ecotone is divided into three subzones: F, the continuous forest below 550 m a.s.l.; FT, a subzone with isolated erect trees (stem height > 2 m) and dense krummholz transitioning upslope increasingly into krummholz and shrub; and T, the alpine tundra above the uppermost limit of tree species, observed to be ca. 750 m a.s.l. (Trant et al., 2011). Examples of vegetation in these subzones are shown in Figure 2.

Automatic climate stations (Campbell Scientific, 2014) were installed in 2001 near the lower margin of the forest-tundra ecotone ("Lower") at 570 m a.s.l. and at a tundra site near the summit ("Upper") at 995 m a.s.l. These stations recorded the following variables at hourly intervals: air temperature (TA) and relative humidity (RH) at 1.5 m above the ground, solar irradiance (KJ,), and ground temperature (T0) at the mineral soil-bedrock interface (1.0 m at Upper and 0.7 m at Lower). At both sites, a storage type of precipitation collector, charged with antifreeze and mineral oil to inhibit freezing and evaporation, was used to measure annual total precipitation (P). In 2005, a third station ("Base") was set up at the base campsite at 600 m a.s.l. to record wind speed and direction, TA, RH, and vapor pressure deficit (VPD), as well as soil temperature at 10, 30, and 60 cm depths. During summer periods when the camp was occupied, rainfall (R) was measured daily at Base using a manual gauge.

Short-term and seasonal data, including air temperature and humidity, soil temperature, and soil moisture, were collected along the altitudinal transect. Portable TA and RH sensors (Onset Hobo ProV2[TM], Onset Computer Corp., 2012) mounted in solar radiation shields on posts at 1.5 m height served for these temporary installations. Soil temperature at 10 cm depth ([T.sub.10]) was recorded with buried miniature temperature loggers (Onset Tidbit[TM] and Onset Pendant[TM]), supplemented by manual measurements using a digital thermometer. Soil moisture was sampled initially using gypsum block sensors, which proved to be relatively insensitive in this near-saturated environment, and subsequently with a Delta-T[TM] radio frequency soil moisture probe (Delta-T, 2005).

Soils were sampled for chemistry initially in 2005 (Yurich, 2007) and more intensively along the transect in 2007 and 2008. Soil pits were dug for bulk samples and to determine the depth of the organic layer (defined as the layer between the base of loose litter and the humus-mineral interface). To provide a standard measure of availability of the principal and minor nutrients (Hofgaard and Rees, 2008), fluxes of major ions were sampled during the growing seasons of 2007 and 2008 using ion exchange membranes buried to a depth of 10 cm. The membranes were of a commercial type (PRS[TM]) obtained from and analyzed by Western Ag Innovations Inc. (2006).

Snow surveys were conducted in late March in 2008, 2009, and 2010 (Leblanc et al., 2010). On each occasion, pits were dug at seven sites from forest to tundra to determine depth, density, and structure of the snowpack. Snow depth was further sampled at 1 m intervals over a radius of 10 m around each snow pit. Digital photography was used to document the late-winter landscape.

Annual thawing degree-days (TDD) were calculated using cumulative daily mean [T.sub.A] above 0[degrees]C; growing degree-days (GDD), using [T.sub.A] above 5[degrees]C; and freezing degree-days (FDD), using [T.sub.A] below 0[degrees]C. Soil temperature at 10 cm depth ([T.sub.10]) was analyzed with respect to two soil temperature thresholds: 0[degrees]C as indicating root-zone freeze and thaw and 3.2[degrees]C to define growing season length, following Korner and Paulsen (2004).

Climatological measurements provide the data for estimating evaporation and transpiration (evapotranspiration, ET) when combined with hydrologic data. Discharge from the Eagle River (Fig. 1) is a response to hydroclimatic conditions within its drainage basin, which includes the study area. Basin runoff, Ro, is given by the water balance equation, Ro = P - ET [+ or -] [DELTA]S, where P is precipitation, ET is evapotranspiration, and [DELTA]S is the change in storage. For annual Ro, and given a sufficiently long record, it can be assumed to a first approximation that [DELTA]S = 0 (Hare, 1980). Discharge data for the Eagle River were obtained from the Water Survey of Canada (EC, 2013b).

To understand the results in a regional context, we compared climate data from Goose Bay and Cartwright with data from the field site, which is located midway between them. Homogenized and adjusted surface air temperature and precipitation datasets of the Canadian Climate Program (Mekis and Vincent, 2011; Vincent et al., 2012; EC, 2013c) were used in this analysis. To estimate long-term climate series for the Mealy Mountains, midpoint series were interpolated using arithmetic means of Goose Bay and Cartwright data (monthly and seasonal [T.sub.A] and annual P, GDD, TDD, and FDD). Further adjustments to series for [T.sub.A], TDD, GDD, and FDD were made using linear regression between data from the study area (dependent) and the interpolated midpoint values for the eight years of observations. ANOVA was applied to test the significance of the results, with the Durbin-Watson statistic used to check for autocorrelation effects (Bowerman and O'Connell, 1993). Regional temperature anomalies were analyzed using gridded NCAR/NCEP reanalysis data (Kalnay et al., 1996) obtained from the interactive mapping site of the National Oceanic and Atmospheric Administration (NOAA, 2012a). Statistical and graphical analyses were done using Microsoft Excel and IBM SPSS Statistics Version 19 (IBM, 2012). Wind frequency plots were constructed using Rose Works (UAI Environmental, 2003).


Tables 1 and 2 summarize climatic and soil data, respectively, obtained at each ecotone subzone in our study area in the Mealy Mountains during the period 2001 to 2009, as well as data from Cartwright and Goose Bay from the same period for comparison. Table 3 shows the statistics of the regression analysis used to estimate long-term climate series for our tree line site, and Table 4 shows the results.

Solar Irradiance

K[down arrow] was consistently higher (11% annually) at Lower station than at Upper (Table 1). This difference was greater than the maximum [+ or -] 5% difference in hourly fluxes found in a post-field comparison of the instruments and the [+ or -] 3% maximum instrumental error specified by the manufacturer. Lower K[down arrow] near the summit can be attributed to more frequent cloud cover there, as was often observed during the field season, but possibly also to rime in colder periods. During the June-September growing season, K[down arrow] at Lower station was only 4% higher than at the summit, i.e., within the range of instrumental uncertainty.


The prevailing winds were westerly, with wind speed most frequently between 5 and 10 m [s.sup.-1]. Seasonal plots (Fig. 3) indicate a wind shift to northwest from winter to spring, with a slight increase in the frequency of northerly through northeasterly winds. Wind speeds are

consistently lower in the summer months (by about 5 m [s.sup.-1]) and increase in the fall. Peak daily winds averaged more than 12 m [s.sup.-1] from December through March, with an hourly maximum of 28 m [s.sup.-1] recorded in January 2008.


Monthly and seasonal means of [T.sub.A], TDD, GDD, and FDD were calculated from eight years of daily values (sample size from 224 to 248 days for each month). Growing season (June-September) [T.sub.A] averaged 11[degrees]C at the forest edge (500 m), 10[degrees]C within the transition, and 8[degrees]C at the upper tundra site (995 m), with corresponding GDD of 764, 694, and 455, respectively (Table 1). Mean annual [T.sub.A] was below freezing at all sites, as was [T.sub.G] at 1 m depth at the Tundra site, indicating the presence of permafrost at higher elevations in the area. The altitudinal gradient of [T.sub.A] between the Lower and Upper stations averaged -0.6 [+ or -] 0.2[degrees]C/100 m annually, -0.4 [+ or -] 0.1[degrees]C/100 m in winter (December-February), and -0.8 [+ or -] 0.1[degrees]C/100 m in the June-September growing period. Nighttime and winter hourly data for these gradients were positively skewed, indicating more frequent temperature inversions at those times.


Annual precipitation at Lower station (3265 mm) was nearly twice that at Upper (1682 mm). Both gauges were of the same design with wind shields to reduce blow-by and were on locally elevated sites where little snow accumulated on the ground. However, the summit ridge location of Upper station was more exposed to prevailing westerly winds than sites lower in the valley, and it is likely that these winds further reduced snow and rain catch in that gauge relative to the catch at the valley site. Snow blowing from the ridge tops and summit may have added to the catch in the gauge lower in the valley. Given that these sites were unmonitored for most of the year, the measurements reported here can only be taken as approximate. A conservative estimate of the annual precipitation is between 2000 and 3000 mm, compared with about 1000 mm at the two permanent stations (Table 1), indicating a strong local orographic enhancement of precipitation in the central Mealy Mountains. This enhancement was also seen in rainfall observations during the summer field season (total of 229 field-days from 2001 to 2009) at Base Camp, where July-August totals averaged 329 mm, or 1.5 times the amounts recorded in the same summer periods at Goose Bay and Cartwright. March snow survey data for 2008-10 show large variation in depth within subzones and between years (Fig. 4). Snowpack density for the study area ranged from 400 to 450 kg [m.sup.-3], resulting in a late-March snow water equivalent (SWE) of 0.65 m.


RH was recorded hourly at all climate stations, while Base station also recorded vapor pressure (VP) and vapor pressure deficit (VPD) derived from RH and TA. Average RH increased with elevation (from 75% at Lower to 78% at Base and 88% at Upper station), as would be expected with decreasing temperature. At Upper station, RH exceeded 80% during 90% of the hourly observations. Hourly VPD was 0.52 [+ or -] 0.45 kPa at Base station (600 m a.s.l.) during the summer months (June-August) between 0800 and 2200 local time.

Analysis across the Forest-Tundra Ecotone

Climate and soil data were analyzed according to ecotone subzone (Tables 1 and 2). Upper station was taken to represent conditions in the T subzone, while the data for Base and Lower stations were averaged to provide a single estimate for the FT subzone. There was no climate station within the continuous forest, but soil temperature was measured there. Air temperatures and degree-day variables were estimated by extrapolation to 500 m in the F subzone using seasonal altitudinal gradients observed between Lower and Upper stations.

To provide further detail on the spatial variability of [T.sub.A], a temporary station in the upper FT subzone at 700 m was operated for 10 days in the summer of 2009. The full array of stations provided a sample of 240 hourly measurements at 570, 600, 700, and 995 m a.s.l. The comparison showed the altitudinal gradient of [T.sub.A] during typical midsummer days to be non-linear, falling initially at about 1[degrees]C/100 m between Lower and Base stations, then less steeply at about 0.3[degrees]C/100 m to the upper part of the FT zone at 700 m. Between 700 m and the tundra site at 995 m, the gradient was about 0.9[degrees]C/100 m. Post-field comparisons found differences between sensors to be less than 0.3[degrees]C. This analysis shows that there is a non-linear relationship between [T.sub.A] and altitude within the FT subzone, with similar growing season temperatures over a relatively broad ecotone area.

Late-March snow depths were taken over three consecutive years (Fig. 4). Aggregated over the three years, mean depths (m) by subzone were 1.98 [+ or -] 0.30 at F, 1.42 [+ or -] 0.58 at FT, and 1.28 [+ or -] 0.51 at T. Differences in depth between F and FT and between F and T were significant (Independent Samples Median Test with Multiple Comparisons, N = 748, p < 0.01), but those between FT and T were not (p = 0.08). Spatial variability was least in the F subzone.

Analysis of two years of measurements found soil temperature [T.sub.10] to be significantly different between subzones in winter, but not in summer and fall (Table 2). The annual cycle of [T.sub.10] differed between subzones (Fig. 5): no deep freezing occurred at the F sites because of their more uniformly deep snow cover (Fig. 4) and because soil moisture was greater in the forest. The period with [T.sub.10] above 0[degrees]C was nearly a month longer at F sites than at the FT and T sites (Table 2). Growing season length, based on [T.sub.10] above 3.2[degrees]C (sensu Korner and Paulsen, 2004), was similar, and growing season [T.sub.10] means of 7.9[degrees] to 8.2[degrees]C did not differ significantly between subzones across an altitudinal range of nearly 500 m.

Soil moisture was greater at F sites than at FT and T sites, while the latter two did not differ significantly (Table 2). Humic layer thickness was significantly less at F sites than at FT or T. The macronutrients total N, N[H.sub.4]-N, P, and S were not significantly different (p > 0.01), but significant differences were found for N-N[O.sub.3], Ca, Mg, and K and the micronutrients Mn and Zn (p < 0.01). Overall, there was no consistent pattern to indicate a more favorable nutrient environment at the forest edge than in the other subzones.

Regional Water Balance

An estimate of basin runoff (Ro) was obtained from long-term (1967-2010) discharge records for the Eagle River. The total annual discharge near the river mouth was divided by the drainage basin area to get the annual Ro. Annual P for the basin was estimated by the arithmetic mean of the annual total P reported at Goose Bay and Cartwright. Over the 44-year record, Eagle River Ro correlated significantly (p < 0.001) with P at Goose Bay (r = 0.62) and Cartwright (r = 0.52), and more strongly (r = 0.66) with P interpolated at their midpoint. The estimated annual basin P (1967-2010) was 1187 [+ or -] 133 mm, and the observed Ro was 734 [+ or -] 124 mm, giving an estimated basin ET of 452 [+ or -] 106 mm. The actual value for ET is likely larger, considering that locally enhanced precipitation effects were not considered in this estimate.

Representativeness of Study Period Conditions in Relation to the Long-Term Climate

Monthly averages of [T.sub.A], TDD, GDD, FDD, and P were highly correlated between Goose Bay and Cartwright over the 1942-2010 period of record. No significant trend was found in those variables over that period (N = 68, p > 0.05). For the period 2001 to 2009, monthly averages of [T.sub.A], GDD, TDD, and FDD in the study area were highly correlated with those at Goose Bay and Cartwright ([R.sup.2] from 0.57 to 0.90, p < 0.01). Differences between stations are attributed to altitude and distance from the Labrador Sea. Eight years of observations is a relatively small sample for climatological estimates. However, comparisons of mean monthly and seasonal values among stations are strengthened by the fact that the data derive from daily values.

As noted in the preceding section, annual P interpolated between the two permanent stations was taken as an estimate for the Mealy Mountains region, without further adjustment for orographic effects. For [T.sub.A], TDD, GDD, and FDD, linear regression was applied between the midpoint series and the study area data for 2001-09, and the resulting equations were used to provide long-term estimates for the study area at tree line elevation. All regression models showed a high degree of explanation of the estimated variables (Table 3). Effects of autocorrelation were ruled out at the [alpha] = 0.01 significance level using the Durbin-Watson statistic for small samples (Bowerman and O'Connell, 1993).

The estimated long-term climatic series for the Mealy Mountains study area are shown in Figure 6. There was no significant trend in growing season [T.sub.A] or GDD until the late 1990s, when the current regional warming started. The water balance terms show a positive trend in P at an average rate of 3.3 mm [y.sup.-1] ([R.sup.2] = 0.16, F = 13.1, and p = 0.001), a slight decrease in Ro (not significant), and a corresponding increase in ET. Table 4 shows comparative summary statistics from those series for the periods 1942-2000, 1971-2000, and 2001-10. The most recent decade was the warmest over the 68 years of record, with significantly higher temperatures in all seasons and higher TDD and GDD. In contrast, differences in P, Ro, and ET were not statistically significant. Winter temperatures and FDD, indicators of winter severity, were less negative in the most recent period than during the colder 1971-2000 period. Gridded NCEP/NCAR reanalysis data for the same periods confirmed an increase in growing season [T.sub.A] of approximately 1[degrees]C in the Mealy Mountains area relative to the preceding tri-decade (Fig. 7), while annual P over the region increased by about 6%. These interpolated gridded data have a relatively coarse resolution of 2.5[degrees] latitude by 2.5[degrees] longitude (Kalnay et al., 1996; NOAA, 2012a) and therefore do not show anomalies on the scale of the Mealy Mountains. Concurrently, during most of 2001-10, sea surface temperatures in the adjacent Labrador Sea were 1[degrees] to 3[degrees]C above their 1971-2000 averages (NOAA, 2012b).


This study has examined climatic and soil factors relevant to vegetation growth across an altitudinal range of ca. 500 m, where outliers of closed-canopy boreal forest transition through open woodland and krummholz forms into alpine tundra. Several factors can be ruled out as constraining forest advance. Observations of solar irradiance agreed with climatological estimates for this region (Gullett, 1987) and did not differ significantly across the ecotone during the growing season. Annual precipitation, summer rainfall, and soil moisture were sufficient for tree growth at all levels, and the availability of most soil nutrients did not differ significantly from forest to tundra. Growing season climate and winter severity are very likely the controlling factors. Additionally, some vegetation changes might also be associated with the regional warming that has been shown to be underway.

Mean [T.sub.A] in the FT subzone for the warmest month (July) was 13[degrees]C, which is well above 10[degrees]C, the value often associated with the northern latitudinal limit of erect trees (Sveinbjornsson, 1992), but within the range for altitudinal tree lines reviewed by Korner (1998). An estimated 764 GDD for the F subzone (Table 1) was near the thermal limit for black spruce seed viability of 800 to 940 GDD found by Meunier et al. (2007). The value within the FT zone was only 694 GDD; however, eastern larch seedlings were found there (Trant, 2013), indicating that viable seeds were being produced for that species. For winter conditions in the FT subzone, the 2001-09 average of-2059 FDD is well below the value of -1816 FDD reported for the same period at Goose Bay (Table 1). However, continuous boreal forest at Goose Bay (Meades, 2007) is a response to the longer-term climate and is better represented by the 1971-2000 normals value of -2025 FDD (EC, 2013a). Therefore, winter severity based on FDD could explain the absence of continuous coniferous forest above 600 m a.s.l. in the Mealy Mountains, only if long-term FDD were much lower than present, as was the case for most of the 20th century (Table 4) and probably longer.

From the forest edge through the ecotone to the tundra site, 10 cm soil temperatures for the growing season were about 2[degrees]C above the average value reported by Korner and Paulsen (2004) in their global survey. Our observations took place in an anomalously warm decade, but the tree line we observed was largely a response to the climate of the previous several decades, when air and soil temperatures were 1[degrees] to 2[degrees]C cooler than at present (Fig. 6), but still within the 6[degrees] to 7[degrees]C range. Therefore, long-term growing season soil temperature does not contraindicate the presence of trees in the FT subzone, a conclusion supported by the presence there of isolated erect conifers in some sheltered areas.

VPD, coupled with temperature, is a predictor for leaf stomatal water vapor conductance and net photosynthesis. Stomatal conductance and C[O.sub.2] assimilation decrease when VPD increases (Oke, 1987; Aber and Federer, 1992), i.e., with higher VPD, a tree is more water-stressed and photosynthesis is slowed. Judging by the relatively low summer daytime VPD observed in the FT subzone, drought or moisture stress do not seem to be significant factors constraining tree growth, although more detailed observations would be needed to confirm this.

Our observations have confirmed an orographic enhancement of P and particularly snowfall in the Mealy Mountains, as shown by Flare (1951, 1980) at a larger scale for central Labrador. Snowfall makes up more than half of annual P, while rainfall is substantial in summer and is likely significant in late spring and early fall. The snow pit stratigraphy (Leblanc et al., 2010) showed evidence of midwinter freezing rain events. Such icing conditions are a major factor inhibiting the survival of arboreal vegetation (e.g., Carlson et al., 2011), and they are likely to become more frequent in a warming climate.

There was large spatial variability in late-winter snow depth in the ecotone compared with sites inside the forest edge (Fig. 4). Prevailing westerly winds, persistent and sometimes gale-force in winter (Fig. 3), resulted in redistribution of snow in exposed areas, giving protection to vegetation in the lee of obstacles, but exposing unsheltered vegetation to abrasion and desiccation. This kind of severe winter environment cannot be captured in the FDD measure, but is well recognized as limiting growth and survival of seedlings and saplings and inhibiting the formation of erect stems in krummholz (Harsch et al., 2009; Holtmeier and Broil, 2010). In the Mealy Mountains, a glaciated landscape with a variety of surficial deposits, spatial differences in snow accumulation, and shelter from or exposure to winter winds, as well as localized summer heating facilitated by microtopography, all contribute to an extended and heterogeneous forest-tundra ecotone.

In the absence of other limiting factors, our analysis indicates that the regional warming of the most recent decade, which is projected to continue (Meehl et al, 2007; Meunier et al., 2007), is sufficient to cause an upslope expansion of forest in the Mealy Mountains. But will it be the same perhumid forest, or will it become drier? The 1942-2010 series for P interpolated from Goose Bay and Cartwright to the study area (Fig. 6) had a positive trend of 3.3 mm [y.sup.-1], in agreement with a regional analysis by Zhang et al. (2000). However, our analysis for 2001-10 shows that while P, [T.sub.A], and ET increased, annual Ro in the Eagle River basin decreased by 8% relative to the 1971-2000 average (Fig. 6), continuing a decrease observed in the mid-1990s (Whitfield and Cannon, 2000). The decrease in Ro was most pronounced in the spring, which is likely a reflection of less winter snowfall. For the same period, snowfall as a proportion of total P at both Goose Bay and Cartwright decreased from 48% to 40% (EC, 2013a). Decreased snowfall means less soil moisture stored into the growing season. Our conclusion is that although precipitation in the area is increasing, the growing season climate is becoming drier as a consequence of increased evapotranspiration driven by climate warming. If sustained, this trend will result in significant changes to forest growth patterns and feedbacks to climate (Meunier et al., 2007).

Warmer, drier summers increase the likelihood of disturbance due to insect herbivory and wildfires (Allen et al., 2010; Peng et al., 2011). Cone collections in the study area during the warm summer of 2008 revealed a high proportion of insect damage to seeds (Jameson, 2012). Foster (1983) investigated wildfire patterns in southeastern Labrador using dendrochronology to date the fire history of areas burned between ca. 1870 and 1979. The central plateau area bordered by the south slopes of the Mealy Mountains was found to have very infrequent small fires, with a fire interval estimated to be over 500 years. The perhumid climate was identified as the dominant factor in this regime. Ignition was attributed almost exclusively to lightning, while the tendency for fires to burn upslope was inhibited by the patchy nature of the glaciated terrain, with its many natural firebreaks. Although wildfires were associated with summer drought, fire did not occur in many years of below-average summer precipitation, possibly for lack of ignition (Foster, 1983). During the 2001-09 field seasons, occasional thunderstorms with lightning were observed in the Mealy Mountains, although no near strikes were observed. Our observations of summer rainfall and soil moisture do not support the conclusion that a shift to a more frequent fire regime has yet occurred in the Mealy Mountains.

It is expected that the warming climate will produce a response in the vegetation. In northern Quebec, where warming is also occurring, Gamache and Payette (2004) found a positive correlation at a number of sites between stem elongation in black spruce krummholz and regional GDD. In the Mealy Mountains FT subzone, conifer species appeared mostly in krummholz form, and reproduction was by layering. However, erect stems were observed emerging from black spruce krummholz, and true seedlings of eastern larch were found (Trant et al., 2011; Trant, 2013). Mealy Mountains GDD for 2001-10 were 20% higher than in the climatological baseline or "normal" period 1971-2000 and 22% higher than the average for the longer-term record (Table 4). This means that during the previous half-century, and almost certainly longer, the climate of this tree line zone was less favorable than at present for seedling germination and survival and for tree growth generally. Winter severity, as measured by lower [T.sub.A] and FDD, was significantly greater in the preceding periods, particularly during the 1971-2000 tri-decade. Long-term studies elsewhere (e.g., Kullman, 2007; Payette, 2007; Harsch et ah, 2009; Harper et ah, 2011; Tremblay and Boudreau, 2011) indicate that boreal forest ecosystems may take several decades to adjust to a change in climate. Our results show that local climate warming is well underway in the tree line ecotone in the Mealy Mountains, and according to concurrent vegetation surveys (Trant et ah, 2011; Trant, 2013), that change is first becoming evident in the krummholz.


We acknowledge the assistance of many members of the Labrador Highlands Research Group over a decade of interdisciplinary research in the Mealy Mountains. Thanks in particular to Trevor Bell and Philippe Leblanc for the snow surveys, Andrew Trant and Ryan Jameson for sharing results of their vegetation studies, and Luise Hermanutz for initiating and maintaining the momentum of this project. Thanks also to Charles Conway for assistance with the graphics and to two anonymous reviewers for many helpful suggestions. We are grateful to the Labrador Innu Nation for encouraging and supporting our work in their traditional territory. Funding has come from Environment Canada (Northern Ecosystem Initiative), the Meteorological Service of Canada--Atlantic Region, the Northern Scientific Training Program, and the Government of Canada Program for the International Polar Year, with in-kind support from the Government of Newfoundland and Labrador (Wildlife Division).


Aber, J.D., and Federer, C.A. 1992. A generalized lumped-parameter model of photosynthesis, evapotranspiration and net primary production in temperate and boreal forest ecosystems. Oecologica 92(4):463-474.

ACIA. 2004. Arctic climate impact assessment: Impacts of a warming Arctic. Cambridge: Cambridge University Press.

Allen, C.D., Macalady, A.K., Chenchouni, H., Bachelet, D., McDowell, N., Vennetier, M., Kitzberger, T., et al. 2010. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Forest Ecology and Management 259(4):660-684. http://dx.doi.Org/10.1016/j.foreco.2009.09.001

Banfield, C.E., and Jacobs, J.D. 1998. Regional patterns of temperature and precipitation for Newfoundland and Labrador during the past century. The Canadian Geographer 42(4):354-364. http://dx.doi.Org/10.llll/j.1541-0064.1998.tb01351.x

Bowerman, B., and O'Connell, R.T. 1993. Forecasting and time series: An applied approach, 3rd ed. Belmont, California: Duxbury Press.

Brady, N.C., and Weil, R.R. 1996. The nature and properties of soils, 11th ed. Upper Saddle River, New Jersey: Prentice Hall.

Cairns, D.M. 2001. Patterns of winter desiccation in krummholz forms of Abies lasiocarpa at treeline sites in Glacier National Park, Montana, USA. Geografiska Annaler, Series A 83(3): 157-168. http://dx.doi.Org/10.llll/j.0435-3676.2001.00151.x

Callaghan, T.V., Bergholm, F., Christensen, T.R., Jonasson, C., Kokfelt, U., and Johansson, M. 2010. A new climate era in the sub-Arctic: Accelerating climate changes and multiple impacts. Geophysical Research Letters 37, L14705.

Campbell Scientific. 2014. Weather& climate: Campbell automated weather stations (AWS) and meteorologic instruments.

Carlson, B.Z., Monroe, J.S., and Hegman, B. 2011. Distribution of alpine tundra in the Adirondack Mountains of New York, U.S.A. Arctic, Antarctic, and Alpine Research 43(3):331-342. http://dx.doi.Org/10.1657/1938-4246-43.3.331

Chan, S. 2010. Regional and local climatology of a subarctic alpine treeline, Mealy Mountains, Labrador. MSc thesis, Department of Geography, Memorial University, St. John's, Newfoundland and Labrador.

Chouinard, C., Fortier, R., and Mareschal, J.-C. 2007. Recent climate variations in the subarctic inferred from three borehole temperature profiles in northern Quebec, Canada. Earth and Planetary Science Letters 263(3-4):355-369. http://dx.doi.Org/10.l016/j.epsl.2007.09.017

Clayden, S.R., Cameron, R.P., and McCarthy, J.W. 2011. Perhumid boreal and hemiboreal forests of eastern Canada. In: DellaSala, D.A., ed. Temperate and boreal rainforests of the world: Ecology and conservation. Washington, D.C.: Island Press. 111-131.

D'Arrigo, R.D., Cook E.R., Jacoby, G.C., and Briffa, K.R. 1993. NAO and sea surface temperature signatures in tree-ring records from the North Atlantic sector. Quaternary Science Reviews 12(6):431-440. )80007-1

D'Arrigo, R., Buckley, B., Kaplan, S., and Woollett, J. 2003. Interannual to multidecadal modes of Labrador climate variability inferred from tree rings. Climate Dynamics 20(23):219-228.

Delta-T Devices Ltd. 2005. User manual for the moisture meter Type HH2. Cambridge: Delta-T Devices.

Edlund, S.A., and Alt, B.T. 1989. Regional congruence of vegetation and summer climate patterns in the Queen Elizabeth Islands, Northwest Territories, Canada. Arctic 42(1):3-23.

Elliott, D.L., and Short, S.K. 1979. The northern limit of trees in Labrador: A discussion. Arctic 32(3):201-206.

EC (Environment Canada). 2013a. Historical climate data.

--. 2013b. Water Survey of Canada.

--. 2013c. Adjusted and homogenized Canadian climate data (AHCCD). A-l

Feldstein, S.B. 2002. The recent trend and variance increase of the annular mode. Journal of Climate 15:88-94.<0088:TRTAVI> 2.0.CO;2

Foster, D.R. 1983. The history and pattern of fire in the boreal forest of southeastern Labrador. Canadian Journal of Botany 61(9):2459-2471.

Gamache, I., and Payette, S. 2004. Height growth response of tree line black spruce to recent climate warming across the forest-tundra of eastern Canada. Journal of Ecology 92:835-845. http://dx.doi.Org/10.llll/j.0022-0477.2004.00913.x.

--. 2005. Latitudinal response of subarctic tree lines to recent climate change in eastern Canada. Journal of Biogeography 32(5):849-862. http://dx.doi.Org/10.llll/j.1365-2699.2004.0H82.x

Grossnickle, S.C. 2000. Ecophysiology of northern spruce species: The performance of planted seedlings. Ottawa: NRC Research Press.

Grace, J., Berninger, F., and Nagy, L. 2002. Impacts of climate change on the tree line. Annals of Botany 90(4):537-544.

Gullett, D. 1987. Climatic Atlas Canada: A series of maps portraying Canada's climate. Map Series 4: Bright sunshine and solar radiation. Downsview, Ontario: Atmospheric Environment Service. 102 maps, 32 p.

Hare, F.K. 1951. The present-day snowfall of Labrador-Ungava. American Journal of Science 249:654-670. http://dx.doi.Org/10.2475/ajs.249.9.654

--. 1980. Long-term annual surface heat and water balances over Canada and the United States south of 60[degrees]N: Reconciliation of precipitation, run-off and temperature fields. Atmosphere-Ocean 18(2):127-153.

Hare, F.K., and Ritchie, J.C. 1972. The boreal bioclimates. Geographical Review 62(3):333-365.

Harper, K.A., Danby, R.K., De Fields, D.L., Lewis, K.P., Trant, A.J., Starzomski, B.M., Savidge, R., and Hermanutz, L. 2011. Tree spatial pattern within the forest-tundra ecotone: A comparison of sites across Canada. Canadian Journal of Forest Research 41(3):479-489.

Harsch, M.A., Hulme, P.E., McGlone, M.S., and Duncan, R.P. 2009. Are treelines advancing? A global meta-analysis of treeline response to climate warming. Ecology Letters 12(10):1040-1049. http://dx.doi.Org/10.l 11 l/j,1461-0248.2009.01355.x

Hofgaard, A., and Harper, K.A. 2011. Tree recruitment, growth, and distribution at the circumpolar forest-tundra transition: Introduction. Canadian Journal of Forest Research 41(3):435-436.

Hofgaard, A., and Rees, G., eds. 2008. PPS Arctic manual: Common protocols for field measurements and handling of collected material.

Hollister, R.D., Webber, P.J., and Bay, C. 2005. Plant response to temperature in northern Alaska: Implications for predicting vegetation change. Ecology 86(6): 1562-1570.

Holtmeier, F.-K., and Broil, G. 2005. Sensitivity and response of Northern Hemisphere altitudinal and polar treelines to environmental change at landscape and local scales. Global Ecology and Biogeography 14(5):395-410.! 11 l/j,1466-822X.2005.00168.x

--. 2010. Wind as an ecological agent at treelines in North America, the Alps, and the European Subarctic. Physical Geography 31(3):203-233.

IBM. 2012. SPSS statistics. http://www-01 statistics/

Jameson, R.G. 2012. Conifer seed production, seed viability and relative potentials for upslope advance at a multispecies treeline, central Labrador, Canada. MSc thesis, Memorial University, St. John's, Newfoundland and Labrador.

Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., et al. 1996. The NCEP/NCAR 40-Year Reanalysis Project. Bulletin of the American Meteorological Society 77(3):437-472. http://dx.doi.0rg/lO.l 175/1520-0477(1996)077<0437:TNYRP> 2.0.CO;2

Korner, C. 1998. A re-assessment of high elevation treeline positions and their explanation. Oecologia 115:445-459.

Korner, C., and Paulsen, J. 2004. A world-wide study of high altitude tree-line temperatures. Journal of Biogeography 31:713-732. 111/j. 1365-2699.2003.01043.x

Kullman, L. 1998. Tree-limits and montane forests in the Swedish Scandes: Sensitive biomonitors of climate change and variability. Ambio 27(4):312-321.

--. 2007. Tree line population monitoring of Pinus sylvestris in the Swedish Scandes, 1973-2005: Implications for tree line theory and climate change ecology. Journal of Ecology 95(1):41 -52. http://dx.doi.Org/10.l 111/j. 1365-2745.2006.01190.x

Lavoie, C., and Payette, S. 1992. Black spruce growth forms as a record of a changing winter environment at treeline, Quebec, Canada. Arctic and Alpine Research 24(l):40-49.

Leblanc, P., Bell, T., Trant, A., and Jacobs, J. 2010. Snow surveys Mealy Mountains. Unpubl. ms. available at the Department of Geography, Memorial University, St. John's, Newfoundland and Labrador, A1B 3X9, Canada.

Meades, S.J. 2007. Ecoregions of Newfoundland and Labrador. Deer Lake, NL, Parks and Natural Areas Division, Department of Environment and Conservation, Government of Newfoundland and Labrador.

Meehl, G.A., Stocker, T.F., Collins, W.D., Friedlingstein, P., Gaye, A.T., Gregory, J.M., Kitoh, A., et al. 2007. Global climate projections. In: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K.B., Tignor M., and Miller, H.L., eds. Climate change 2007: The physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. New York: Cambridge University Press. 747-845.

Mekis, E., and Vincent, L.A. 2011. An overview of the second generation adjusted daily precipitation dataset for trend analysis in Canada. Atmosphere-Ocean 49(2):163-177.

Meunier, C., Sirois, L., and Begin, Y. 2007. Climate and Picea mariana seed maturation relationships: A multi-scale perspective. Ecological Monographs 77(3):361-376. http://dx.d0i.0rg/l 0.1890/06-1543.1

Munier, A., Hermanutz, L., Jacobs, J.D., and Lewis, K. 2010. The interacting effects of temperature, ground disturbance, and herbivory on seedling establishment: Implications for treeline advance with climate warming. Plant Ecology 210(1):19-30.

NOAA (National Oceanic and Atmospheric Administration). 2012a. Monthly/seasonal climate composites. Earth System Research Laboratory, Physical Sciences Division.

--. 2012b. State of the climate: Global analysis. National Climate Data Center.

Oke, T.R. 1987. Boundary layer climates, 2nd ed. New York: Methuen.

Onset Computer Corp. 2012. Specifications for the Onset Hobo and Tidbit temperature data loggers. Onset Computer Corp., Bourne, Massachusetts 02532.

Palmer, M.J., Burn, C.R., and Kokelj, S.V. 2012. Factors influencing permafrost temperatures across tree line in the uplands east of the Mackenzie Delta, 2004-2010. Canadian Journal of Earth Sciences 49(8):877-894.

Payette, S. 2007. Contrasted dynamics of northern Labrador tree lines caused by climate change and migrational lag. Ecology 88(3):770-780.

Peng, C., Ma, Z., Lei, X., Zhu, Q., Chen, H., Wang, W., Liu, S., Li, W., Fang, X., and Zhou, X. 2011. A drought-induced pervasive increase in tree mortality across Canada's boreal forests. Nature Climate Change 1:467-471. http://dx.doi.Org/10:1038/nclimatel293

Rouse, W.R. 1991. Impacts of Hudson Bay on the terrestrial climate of the Hudson Bay Lowlands. Arctic and Alpine Research 23(1):24-30.

Sirois, L. 2000. Spatiotemporal variation in black spruce cone and seed crops along a boreal forest--tree line transect. Canadian Journal of Forest Research 30(6):900-909. http://dx.doi.Org/10.l 139/x00-015

Sirois, L., Begin, Y., and Parent, J. 1999. Female gametophyte and embryo development of black spruce along a shore-hinterland climatic gradient of a recently created reservoir, northern Quebec. Canadian Journal of Botany 77(l):61-69.

Smith, M.W., and Riseborough, D.W. 2002. Climate and the limits of permafrost: A zonal analysis. Permafrost and Periglacial Processes 13(1):1-15.

Sundqvist, M.K., Bjork, R.G., and Molau, U. 2008. Establishment of boreal forest species in alpine dwarf-shrub heath in subarctic Sweden. Plant Ecology & Diversity l(l):67-75.

Sveinbjornsson, B. 1992. Arctic tree line in a changing climate. In: Chapin, F.S., III, Jefferies, R.L., Reynolds, J.F., Shaver, G.R., and Svoboda, J., eds. Arctic ecosystems in a changing climate. San Diego: Academic Press. 239-256.

Sveinbjornsson, B., Hofgaard, A., and Lloyd, A. 2002. Natural causes of the tundra-taiga boundary. Ambio Special Report 12:23-29.

Timoney, K. 1995. Tree and tundra cover anomalies in the subarctic forest-tundra of northwestern Canada. Arctic 48(1):13-21.

Trant, A. 2013. Changing treelines: How variability in scale and approach improve our understanding. PhD Thesis, Memorial University. St. John's, Newfoundland and Labrador.

Trant, A.J., Jameson, R.G., and Hermanutz, L. 2011. Persistence at the tree line: Old trees as opportunists. Arctic 64(3):367-370.

Tremblay, G.D., and Boudreau, S. 2011. Black spruce regeneration at the treeline ecotone: Synergistic impacts of climate change and caribou activity. Canadian Journal of Forest Research 41:460-468.

Trindade, M., Bell, T., Laroque, C.R, Jacobs, J.D., and Flermanutz, L. 2011. Dendroclimatic response of a coastal alpine treeline ecotone: A multispecies perspective from Labrador. Canadian Journal of Forest Research 41:469-478.

UAI Environmental. 2003. Rose Works, Version 1.0. Lancaster, Pennsylvania: UAI Environmental Inc.

Vincent, L.A., Wang, X.L., Milewska, E.J., Wan, H., Yang, F., and Swail, V. 2012. A second generation of homogenized Canadian monthly surface air temperature for climate trend analysis. Journal of Geophysical Research 117, D18110.

Vose, R., Easterling, D., and Gleason, B. 2005. Maximum and minimum temperature trends for the globe: An update through 2004. Geophysical Research Letters 32, L23822.

Walker, D.A. 2000. Hierarchical subdivision of Arctic tundra based on vegetation response to climate, parent material and topography. Global Change Biology 6(Suppl. 1):19-34. http://dx.doi.Org/10.1046/j.1365-2486.2000.06010.x

Wardle, R., Gower, C., Ryan, B., James, D., and Kerr, A. 1997. Geological map of Labrador; 1:1 million scale. Geological Survey Map 97-07. St. John's: Department of Mines and Energy, Government of Newfoundland and Labrador.

Western Ag Innovations Inc. 2006. PRS[TM] Probe operations manual. Saskatoon, Saskatchewan: Western Ag Innovations Inc.

Wheeler, J.A., Hermanutz, L., and Marino, P.M. 2011. Feathermoss seedbeds facilitate black spruce seedling recruitment in the forest-tundra ecotone (Labrador, Canada). Oikos 120(8): 1263-1271. http://dx.doi.Org/10.llll/j.1600-0706.2010.18966.x

Whitfield, P.H., and Cannon, A.J. 2000. Recent variations in climate and hydrology in Canada. Canadian Water Resources Journal 25(1):19-65.

Yurich, C. 2006. The effects of altitude and site on soil chemistry, seed bank, and soil temperature in the Mealy Mountains. Unpubl. research report submitted in partial fulfillment of the requirements for the degree of Master of Environmental Science, School of Graduate Studies, Memorial University, St. John's, Newfoundland and Labrador.

Zhang, X., Vincent, L.A., Hogg, W.D., and Niitsoo, A. 2000. Temperature and precipitation trends in Canada during the 20th century. Atmosphere-Ocean 38(3):395-429.


(Received 17 April 2012; accepted in revised form 22 April 2013)

(1) Department of Geography, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador A1B 3X9, Canada;

(2) Present address: Yukon Environmental and Socio-Economic Assessment Board, PO Box 2126, Haines Junction, Yukon Y0B 1L0, Canada

(3) Present address: 101-2050 Scotia St., Vancouver, British Columbia V5T 4T1, Canada

TABLE 1. Summary of observational annual and seasonal
climatic data for the Mealy Mountains study site (2001-09),
by ecotone subzone. For comparison, Cartwright (WCA) and
Goose Bay (YYR) data for the same period are also shown.
[T.sub.G] = ground temperature ([degrees]C), [T.sub.A] =
surface air temperature ([degrees]C), GDD = growing degree-
days, TDD = thawing degree-days, FDD = freezing degree-
days, P = precipitation (mm), R = rainfall (mm), and K[down
arrow] = solar irradiance (MJ [m.sup.-2] [d.sup.-1]). "n.d."
means no data.

Ecotone subzone                   F                      FT

Elevation (m)                    500                     600
[T.sub.G] (Annual)               n.d.              -0.1 [+ or -] 0.7
[T.sub.A] (Annual)       (-1.2 [+ or -] 0.8)(1)    -1.9 [+ or -] 0.8
[T.sub.A]                (11.0 [+ or -] 0.6)       10.1 [+ or -] 0.6
[T.sub.A]               (-14.8 [+ or -] 1.5)      -15.2 [+ or -] 1.6
[T.sub.A] (July)         (13.6 [+ or -] 0.9)       12.9 [+ or -] 0.8
Days [T.sub.A] > 0        (119 [+ or -] 13)         113 [+ or -] 12
Days [T.sub.A] > 5         (71 [+ or -] 9)           65 [+ or -] 7
GDD (Annual)              (764 [+ or -] 93)         694 [+ or -] 85
TDD (Annual)             (1502 [+ or -] 157)       1399 [+ or -] 149
FDD (October-May)       (-1968 [+ or -] 227)      -2059 [+ or -] 231
P (Annual)                       n.d               3265 [+ or -] 853
R (July-August)                  n.d.               329 [+ or -] 161
K[down arrow]                    n.d.              11.8 [+ or -] 6.9
K[down arrow]                    n.d.              15.8 [+ or -] 4.3

Ecotone subzone                 T                   WCA

Elevation (m)                  995                   14
[T.sub.G] (Annual)       -1.9 [+ or -] 0.9           n.d.
[T.sub.A] (Annual)       -4.3 [+ or -] 0.8     1.1 [+ or -] 0.8
[T.sub.A]                 7.7 [+ or -] 0.7    11.8 [+ or -] 0.7
[T.sub.A]               -16.7 [+ or -] 1.5   -10.7 [+ or -] 1.8
[T.sub.A] (July)         10.5 [+ or -] 0.8    13.8 [+ or -] 1.7
Days [T.sub.A] > 0         90 [+ or -] 10      179 [+ or -] 5
Days [T.sub.A] > 5         42 [+ or -] 6        83 [+ or -] 12
GDD (Annual)              455 [+ or -] 62      922 [+ or -] 108
TDD (Annual)             1029 [+ or -] 119    1806 [+ or -] 176
FDD (October-May)       -2393 [+ or -] 234   -1386 [+ or -] 229
P (Annual)               1682 [+ or -] 521    1055 [+ or -] 117
R (July-August)                n.d.            186 [+ or -] 48
K[down arrow]            10.6 [+ or -] 6.9           n.d.
K[down arrow]            15.2 [+ or -] 4.3           n.d.

Ecotone subzone                YYR

Elevation (m)                   48
[T.sub.G] (Annual)             n.d.
[T.sub.A] (Annual)        0.9 [+ or -] 0.8
[T.sub.A]                14.0 [+ or -] 0.7
[T.sub.A]               -14.2 [+ or -] 1.5
[T.sub.A] (July)         16.7 [+ or -] 1.6
Days [T.sub.A] > 0        180 [+ or -] 5
Days [T.sub.A] > 5        104 [+ or -] 10
GDD (Annual)             1227 [+ or -] 121
TDD (Annual)             2133 [+ or -] 175
FDD (October-May)       -1816 [+ or -] 202
P (Annual)                917 [+ or -] 67
R (July-August)           226 [+ or -] 39
K[down arrow]                  n.d.
K[down arrow]                  n.d.

(1) Values in parentheses were estimated by extrapolation.

TABLE 2. Soil variables by ecotone subzone. Values are
spatial averages and standard deviations from observations
at multiple (n) sample sites. Comparisons used one-way ANOVA
with post hoc tests. Differences between subzones are
considered significant when p [less than or equal to] 0.01.

Ecotone subzone            Forest                  Transition
Soil variable                 F                        FT

Soil temperature            n = 8                    n = 10
  at 10 cm
Annual (1)            3.0 [+ or -] 0.2 *         2.1 [+ or -] 0.7 *
Jun-Nov               6.2 [+ or -] 0.4           6.2 [+ or -] 0.9
Dec-May              -0.2 [+ or -] 0.4 *        -1.9 [+ or -] 1.4 *
Growing               8.2 [+ or -] 0.6           7.9 [+ or -] 1.5
  season (2)
Days                         111                       115
  [T.sub.10] >
Days                         220                       193
  [T.sub.10] >
Soil moisture              n = 26                    n = 131
  % volumetric       77.2 [+ or -] 17.0 *       42.2 [+ or -] 21.2
  capacity (3)
Soil humic                 n = 40                    n = 64
  layer               1.5 [+ or -] 1.6 *         6.9 [+ or -] 4.7
Soil nutrient              n = 66                    n = 113
  flux (4)
Total N               19.1 [+ or -] 8.4         17.8 [+ or -] 9.2
N-N[O.sub.3]          2.40 [+ or -] 2.88        3.24 [+ or -] 2.79
N[H.sub.4]-N          16.7 [+ or -] 8.4         14.6 [+ or -] 9.5
Ca                     379 [+ or -] 329 *        455 [+ or -] 292 *
Mg                    90.1 [+ or -] 75.6        96.9 [+ or -] 61.0
K                    137.9 [+ or -] 116.9 *     76.7 [+ or -] 47.7 *
P                     1.55 [+ or -] 1.62        1.67 [+ or -] 2.78
S                     17.4 [+ or -] 14.8        19.2 [+ or -] 15.6
Mn                    5.00 [+ or -] 5.82        3.00 [+ or -] 6.16
Zn                    1.89 [+ or -] 1.67        1.72 [+ or -] 1.24
B                     0.97 [+ or -] 0.45        1.00 [+ or -] 0.53
Fe                    4.29 [+ or -] 5.37        34.0 [+ or -] 188
Al                    34.0 [+ or -] 20.4        39.6 [+ or -] 23.5

Ecotone subzone            Tundra             ANOVA
Soil variable                 T              Sig. (p)

Soil temperature           n = 6
  at 10 cm
Annual (1)             1.4 [+ or -] 5.9 *      <0.01
Jun-Nov                5.9 [+ or -] 0.6         0.66
Dec-May               -3.2 [+ or -] 1.3 *      <0.01
Growing                8.0 [+ or -] 1.2         0.90
  season (2)
Days                        106                 n.d.
  [T.sub.10] >
Days                        179                 n.d.
  [T.sub.10] >
Soil moisture              n = 33
  % volumetric        55.2 [+ or -] 24.5       <0.01
  capacity (3)
Soil humic                 n = 40
  layer                5.6 [+ or -] 4.5        <0.01
Soil nutrient              n = 42
  flux (4)
Total N               18.0 [+ or -] 12.0        0.68
N-N[O.sub.3]          4.69 [+ or -] 4.52 *     <0.01
N[H.sub.4]-N          13.3 [+ or -] 10.5        0.16
Ca                     177 [+ or -] 58 *       <0.01
Mg                    35.6 [+ or -] 18.7 *     <0.01
K                     49.2 [+ or -] 22.0 *     <0.01
P                     1.03 [+ or -] 0.55        0.27
S                     23.6 [+ or -] 16.0        0.13
Mn                    0.94 [+ or -] 1.04 *     <0.01
Zn                    0.71 [+ or -] 0.44 *     <0.01
B                     0.86 [+ or -] 0.38        0.25
Fe                    17.8 [+ or -] 87.2        0.38
Al                    43.4 [+ or -] 12.8        0.06

* Post-hoc test shows value differs significantly from those
in other subzones at p < 0.01.

(1) Soil temperatures based on daily means over two years
at multiple sites.

(2) Based on a threshold of 3.2[degrees]C (Korner and Paulsen, 2004).

(3) Non-parametric test used on soil moisture.

(4) Units are [micro]g/10 [cm.sup.2]--four weeks.

TABLE 3. Statistics of regression models for predicting [T.sub.A],
GDD, TDD, and FDD in the Mealy Mountains study area from midpoint
averages of data from Goose Bay and Cartwright. N = 8 years.
[R.sup.2] is the coefficient of determination, F is the ANOVA
F-statistic for the regression, p is the significance level of F,
and d is the Durbin-Watson statistic, which for all the variables
significantly exceeded the critical value for positive
autocorrelation (p < 0.01).

Variable              [R.sup.2]   SE est.     F       p       d

[T.sub.A] (Annual)      0.92       0.28      53.5    0.001   2.3
[T.sub.A] (January)     0.97       0.62     190.5   <0.001   1.5
[T.sub.A] (Dec-Feb)     0.98       0.22     352.7   <0.001   2.7
[T.sub.A] (July)        0.94       0.24      87.3   <0.001   2.6
[T.sub.A] (Jun-Sep)     0.77       0.36      19.7    0.004   3.2
GDD                     0.88      33.6       43.9    0.001   1.9
TDD                     0.83      70.2       18.8    0.004   1.7
FDD                     0.91      75.6       36.1    0.001   1.3

TABLE 4. Climate and water balance variables during three
different periods for the Mealy Mountains tree line at 600
m a.s.l, calculated from regression on Goose Bay and
Cartwright data. ET is estimated from midpoint P and Eagle
annual Ro, all in mm. All temperatures are [degrees]C. Values
marked by an asterisk (*) were significantly greater in the
most recent decade than in the previous intervals (Student's
t, p < 0.01), and "n.d." indicates no data.

Variable\Period           1942-2000            1971-2000

[T.sub.A] (Annual)     -2.6 [+ or -] 0.9    -2.8 [+ or -] 1.0
[T.sub.A] (January)   -17.2 [+ or -] 3.5   -17.9 [+ or -] 3.0
[T.sub.A] (Dec-Feb)   -16.2 [+ or -] 2.5   -16.9 [+ or -] 1.9
[T.sub.A] (July)       12.4 [+ or -] 0.8    12.3 [+ or -] 1.0
[T.sub.A] (Jun-Sep)     9.0 [+ or -] 0.8     9.2 [+ or -] 0.8
GDD                     591 [+ or -] 67      600 [+ or -] 68
TDD                    1110 [+ or -] 137    1132 [+ or -] 140
FDD                   -2285 [+ or -] 298   -2360 [+ or -] 333
P                      1128 [+ or -] 172    1201 [+ or -] 139
Ro                            n.d.           752 [+ or -] 136
ET                            n.d.           450 [+ or -] 113

Variable\Period             2001-10

[T.sub.A] (Annual)     -1.3 [+ or -] 1.0 *
[T.sub.A] (January)   -16.1 [+ or -] 3.5 *
[T.sub.A] (Dec-Feb)   -14.8 [+ or -] 2.3 *
[T.sub.A] (July)       13.2 [+ or -] 1.1 *
[T.sub.A] (Jun-Sep)    10.5 [+ or -] 0.6 *
GDD                     724 [+ or -] 76 *
TDD                    1350 [+ or -] 137 *
FDD                   -2053 [+ or -] 330
P                      1181 [+ or -] 103
Ro                      694 [+ or -] 82
ET                      487 [+ or -] 58
COPYRIGHT 2014 Arctic Institute of North America of the University of Calgary
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Jacobs, John D.; Chan, Sarah; Sutton, Elizabeth
Article Type:Report
Geographic Code:1CNEW
Date:Mar 1, 2014
Previous Article:Invasive species management programs in Alaska: a survey of statewide expenditures, 2007--11.
Next Article:The Boulder Patch (North Alaska, Beaufort Sea) and its benthic algal flora.

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters