Search

Arctic and Subarctic environments are among the most vulnerable regions to climate change. Increases in liquid precipitation and changes in snowmelt onset are cited as the main drivers of change in streamflow and water temperature patterns in some of the largest rivers of the Canadian Arctic. However, in spite of this evidence, there is still a lack of research on water temperature, particularly in the eastern Canadian Arctic. In this paper, we use the CEQUEAU hydrological‐water temperature model to derive consistent long‐term daily flow and stream temperature time series in Aux Mélèzes River, a non‐regulated basin (41 297 km2) in the eastern Canadian subarctic. The model was forced using reanalysis data from the fifth‐generation ECMWF atmospheric reanalyses (ERA5) from 1979 to 2020. We used water temperature derived from thermal infrared (TIR) images as reference data to calibrate CEQUEAU's water temperature model, with calibration performed using single‐site, multi‐site, and upscaling factors approaches. Our results indicate that the CEQUEAU model can simulate streamflow patterns in the river and shows excellent spatiotemporal performance with Kling‐Gupta Efficiency (KGE) metric >0.8. Using the best‐performing flow simulation as one of the inputs allowed us to produce synthetic daily water temperature time series throughout the basin, with the multi‐site calibration approach being the most accurate with root mean square errors (RMSE) <2.0°C. The validation of the water temperature simulations with a three‐year in situ data logger dataset yielded an RMSE = 1.38°C for the summer temperatures, highlighting the robustness of the calibrated parameters and the chosen calibration strategy. This research demonstrates the reliability of TIR imagery and ERA5 as sources of model calibration data in data‐sparse environments and underlines the CEQUEAU model as an assessment tool, opening the door to its use to assess climate change impact on the arctic regions of Canada. [ABSTRACT FROM AUTHOR]

This paper describes an approach for combining Landsat and Radarsat satellite images to generate national statistics for urban ecosystem accounting. These accounts will inform policy related to the development of mitigation measures for climatic and hydrologic events in Canada. Milton, Ontario was used as a test case for the development of an approach identifying urban ecosystem types and assessing change from 2001 to 2019. Methods included decomposition of Radarsat images into polarimetric parameters to test their usefulness in characterizing urban areas. Geographic object-based image analysis (GEOBIA) was used to identify urban ecosystem types following an existing classification of local climate zones. Three supervised classifiers: decision tree, random forest and support vector machine, were compared for their accuracy in mapping urban ecosystems. Ancillary geospatial datasets on roads, buildings, and Landsat-based vegetation were used to better characterize individual ecosystem assets. Change detection focused on the occurrence of changes that can impact ecosystem service supply – i.e., conversions from less to more built-up urban types. Results demonstrate that combining Radarsat polarimetric parameters with the Landsat images improved urban characterization using the GEOBIA random forest classifier. This approach for mapping urban ecosystem types provides a practical method for measuring and monitoring changes in urban areas. [ABSTRACT FROM AUTHOR]

Accurate and fine-scale forest data are essential to improve natural resource management, particularly in the face of climate change. Here, we present SCANFI, the Spatialized CAnadian National Forest Inventory, which provides coherent, 30 m resolution 2020 wall-to-wall maps of forest attributes (land cover type, canopy height, crown closure, aboveground tree biomass, and main species composition). These maps were developed using the NFI photo-plot dataset, a systematic regular sample grid of photo-interpreted high-resolution imagery covering all of Canada's non-arctic landmass. SCANFI was produced using temporally harmonized summer and winter Landsat spectral imagery along with hundreds of tile-level regional models based on a multiresponse k-nearest neighbours and random forest imputation method. This tile-level approach revealed the importance of radiometric variables in predicting vegetation attributes, namely winter radiometry, as the large-scale climate gradients were controlled at the tile-level. SCANFI was validated with rigorous cross-validation analyses, which revealed robust model performance for structural attributes (biomass R2 = 0.76; crown closure R2 = 0.82; height R2 = 0.78) and tree species cover (e.g., Douglas fir R2 = 0.60). SCANFI attributes were also validated with several independent external products, ranging from ground plot-based tree species cover (e.g., black spruce R2 = 0.53) to satellite LiDAR height data products (e.g., crown closure R2 = 0.71). SCANFI total aboveground biomass trends also followed those published by other studies. The methodology presented herein can be used to map time series of these attributes, identify the original training points used to make any given prediction, as well as map additional variables associated with the NFI photo-plots that are challenging to map using traditional remote sensing approaches. [ABSTRACT FROM AUTHOR]

Meltwater drainage systems beneath ice sheets are a poorly understood, yet fundamentally important environment for understanding glacier dynamics, which are strongly influenced by the nature and quantity of meltwater entering the subglacial system. Contemporary sub-ice sheet meltwater drainage systems are notoriously difficult to study, but we can utilise exposed beds of palaeo-ice sheets to further our understanding of subglacial drainage. In particular, eskers record deposition in glacial drainage channels and are widespread on the exposed beds of former ice sheets. This paper presents a 1:5,000,000 scale map of >20,000 large eskers (typically > 2 km long) deposited by the Laurentide Ice Sheet (LIS), mapped from Landsat imagery of Canada, in order to establish a dataset suitable for analysis of esker morphometry and drainage patterns at the ice sheet scale. Comparisons between eskers mapped from Landsat imagery and aerial photographs indicate that, in most areas, approximately 75% of eskers are detected using Landsat. The data presented in this map build on and extend previous work in providing a consistent map of an unprecedented sample of eskers for quantitative analysis. It offers an alternative perspective on the problems surrounding ice-sheet meltwater drainage and can be used for: (i) detailed investigations of esker morphometry and distribution from a large sample size; (ii), testing of numerical models of meltwater drainage routing that predict esker characteristics (e.g. channel spacing, sinuosity), (iii) assessment of the factors that control esker location and formation; and (iv), a refined understanding of ice margin configurations during retreat of the LIS. [ABSTRACT FROM AUTHOR]

Land cover information is necessary for a large range of environmental applications related to climate impacts and adaption, emergency response, wildlife habitat, etc. In Canada, a 2008 user survey indicated that the most practical land cover data is provided in a nationwide 30 m spatial resolution format, with an update frequency of five years. In response to this need, the Canada Centre for Remote Sensing (CCRS) has generated a 30 m land cover map of Canada for the base year 2010, as the first of a planned series of maps to be updated every five years, or more frequently. This land cover dataset is also the Canadian contribution to the 30 m spatial resolution 2010 Land Cover Map of North America, which is produced by Mexican, American and Canadian government institutions under a collaboration called the North American Land Change Monitoring System (NALCMS). This paper describes the mapping approach used for generating this land cover dataset for Canada from Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) Landsat sensor observations. The innovative part of the mapping approach is the local optimization of the land cover classifier, which has resulted in increased spatial consistency and accuracy. Training and classifying with locally confined reference samples over a large number of partially overlapping areas (i.e., moving windows) ensures the optimization of the classifier to a local land cover distribution, and decreases the negative effect of signature extension. A weighted combination of labels, which is determined by the classifier in overlapping windows, defines the final label for each pixel. Since the approach requires extensive computation, it has been developed and deployed using the Government of Canada's High-Performance Computing Center (HPC). An accuracy assessment based on 2811 randomly distributed samples shows that land cover data produced with this new approach has achieved 76.60% accuracy with no marked spatial disparities. [ABSTRACT FROM AUTHOR]

The monitoring of forest ecosystems is significantly affected by the lack of consistent historical data of low-severity (forest partially disturbed) or gradual disturbance (e.g. eastern spruce budworm epidemic). The goal of this paper is to explore the use of a subset of Landsat time series and deep learning models to identify both the type and the year of disturbances, including low-severity and gradual disturbances, in the boreal forest of eastern Canada at the pixel level. Remote sensing data such as the spectral information from Landsat time series are the best available option for large scale observations of disturbances that go back decades. Traditional modeling approaches, like LandTrendr, require substantial handcrafted pre-processing to remove noise and to extract temporal features from the image sequences before using them as input to a classical machine-learning model. Deep-learning models can autonomously discern which features are relevant within the coarse temporal and spectral information from the Landsat annual dense time series. We evaluated the performance of TempCNN and Transformer model in detecting and classifying the type and the year of the forest disturbance using Landsat time series subsequences. Our findings resulted in the generation of four disturbance maps outlining the forest history from 1986 to 2021 within the eastern Canadian boreal forest. Our experimental outcomes demonstrate several significant benefits of employing deep learning models. Firstly, using noisy Landsat time series they achieve comparable accuracy for classifying fire and total harvesting than existing publicly available disturbance maps. Secondly, the use of shorter time series subsequence with deep learning models enables to map adequately different overlapping disturbances occurring in the complete time series. Finally, they increase the number of distinguishable disturbance classes by adding partial harvesting, gradual disturbances, and forest recovery from older events, making them useful approaches for obtaining the first remote sensing-based map for areas affected by the eastern spruce budworm. • TempCNN, Transformer tested on Landsat time series for disturbance classification. • The type and the year of disturbance events were classified with good accuracy. • Deep learning effective for identifying partial and progressive disturbances. • Landsat subsequences extract overlapping disturbances. [ABSTRACT FROM AUTHOR]

Canada's forests play a critical role in the global carbon (C) cycle and are responding to unprecedented climate change as well as ongoing natural and anthropogenic disturbances. However, the representation of disturbance in boreal regions is limited in pre-existing land surface models (LSMs). Moreover, many LSMs do not explicitly represent subgrid-scale heterogeneity resulting from disturbance. To address these limitations, we implement harvest and wildfire forcings in the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) land surface model alongside dynamic tiling that represents subgrid-scale heterogeneity due to disturbance. The disturbances are captured using 30 m spatial resolution satellite data (Landsat) on an annual basis for 33 years. Using the pan-Canadian domain (i.e., all of Canada south of 76° N) as our study area for demonstration, we determine the model setup that optimally balances a detailed process representation and computational efficiency. We then demonstrate the impacts of subgrid-scale heterogeneity relative to standard average individual-based representations of disturbance and explore the resultant differences between the simulations. Our results indicate that the modeling approach implemented can balance model complexity and computational cost to represent the impacts of subgrid-scale heterogeneity resulting from disturbance. Subgrid-scale heterogeneity is shown to have impacts 1.5 to 4 times the impact of disturbance alone on gross primary productivity, autotrophic respiration, and surface energy balance processes in our simulations. These impacts are a result of subgrid-scale heterogeneity slowing vegetation re-growth and affecting surface energy balance in recently disturbed, sparsely vegetated, and often snow-covered fractions of the land surface. Representing subgrid-scale heterogeneity is key to more accurately representing timber harvest, which preferentially impacts larger trees on higher quality and more accessible sites. Our results show how different discretization schemes can impact model biases resulting from the representation of disturbance. These insights, along with our implementation of dynamic tiling, may apply to other tile-based LSMs. Ultimately, our results enhance our understanding of, and ability to represent, disturbance within Canada, facilitating a comprehensive process-based assessment of Canada's terrestrial C cycle. [ABSTRACT FROM AUTHOR]

The fate of the carbon stored in frozen peatlands in Canada's Hudson Bay Lowlands (HBL) depends strongly on water, with wetter conditions favoring long-term storage. Dynamic surface water products generated from historical Landsat data exist to inform surface water trends in the HBL based on binary classifications of land vs water at 30 m spatial resolution. However, the HBL contains many water features smaller than 30 m, including streams and patterned fens that require a sub-pixel mapping approach to fully capture their dynamics. In this paper, we leverage an existing binary dynamic surface water product to generate a spatially and temporally comprehensive Landsat sub-pixel surface water time-series over the HBL. We first generate Landsat Red, NIR and SWIR composites from 1985 to 2021 and evaluate data-driven machine learning and physical models to derive sub-30 m water fractions across the entire mapping domain. Data-driven models were calibrated on 30 m water fractions obtained from 16 of 17 WorldView scene water fraction maps and assessed on the 17th held-out scene in a jackknife validation. Data-driven models were also calibrated by relating scaled Landsat reflectance to water fractions obtained from binary water masks at corresponding scales from 90 to 270 m. Physical models involved simple linear unmixing using a new approach that samples land endmembers in the vicinity of each pixel to represent land reflectance, and a global water endmember sampled beneath permanent water in each image. The assumption of linear scaling was verified at sub-30 m to 270 m resolution for NIR, SWIR and NDVI, which were all highly correlated with water fraction and had no significant differences among water fraction – spectral feature relationships across the full range of scales. Each of 11 different methods and input spectral features were evaluated against WorldView water fractions and ranked based on nine criteria. We found that machine learning using all spectral features including local NIR land endmembers performed best among data-driven approaches, and linear unmixing of the NIR band using local land endmembers performed best among physical. Both best-performing methods were applied to the entire time-series and surface water area trends were calculated and compared with trends obtained from binary water masks. All three methods agree on the direction of surface water change towards wetting, with sub-pixel inter-annual variation agreeing with stream discharge noted in other studies. • Physical and data-driven Landsat water fraction methods tested in the Hudson Bay Lowlands. • Physical methods perform linear unmixing between local land and global water endmembers. • Data-driven methods consist of machine learning calibrated using high-resolution water. • Best methods applied to 1985–2021 Landsat time-series reveal overall wetting trends. • Spatial trends generally suggest drying in bogs and wetting in fens. [ABSTRACT FROM AUTHOR]

Information on biomass distribution is needed to estimate GHG emissions and removals from land use changes in Canada's north for UNFCCC reporting. This paper reports aboveground biomass measurements along the Dempster Highway transect in 2004, and around Yellowknife and the Lupin Gold Mine in 2005. The measured aboveground biomass ranges are 10-100 t ha-1 for woodlands, 1-100 t ha-1 for shrub sites, and 0.5-10 t ha-1 for grass/herbs sites. The root mean squared error (RMSE) of measurements is 21%, and the median absolute percentage error (MedAPE) is 14%. The combination of JERS backscatter and Landsat TM4/TM5 gives the best biomass equation for the Dempster Highway transect, with r 2 = 0.72 when using a one-step approach (i.e. using all points) and 0.78 when using a two-step approach (i.e. stratifying data into three classes: grass, shrub, and woodlands). The two-step approach reduces the MedAPE from 53% to 33%. The validation against Yellowknife & Lupin data indicates that the equations have good transferability. The improvement of two-step approach over the one-step approach, however, is not significant for the validation dataset, suggesting that the one-step approach is as good as the two-step approach when applied over areas outside where the equations are developed. The relationships and error analysis of this study, as well as the final estimate of GHG emission/removal over Canada's north have been incorporated into Canada's 2006 UNFCCC report. [ABSTRACT FROM AUTHOR]

This paper describes single-date and multi-date land-cover classification accuracy results using segment-based, gap-filled Landsat 7 Enhanced Thematic Mapper data compared with Landsat 5 Thematic Mapper data captured one day apart. Maximum likelihood and Decision tree classification algorithms were evaluated. The same training and verification sets of ground data were used for each classification evaluation. For the comparison with the single-date classification, an average decrease of 2.8% in the classification accuracy was obtained with the use of the gap-filled Landsat data. Area estimates for the mid-summer images differed, on average, from 0.6% to 1.9% for a four-class and eight-class classification, respectively. A multi-date land-cover classification was also completed with the addition of a late spring Landsat 5 image, resulting in an average decrease in classification accuracy of 1.8%. [ABSTRACT FROM AUTHOR]

Fire is the dominant disturbance agent in Alaskan and Canadian boreal ecosystems and releases large amounts of carbon into the atmosphere. Burned area and carbon emissions have been increasing with climate change, which have the potential to alter the carbon balance and shift the region from a historic sink to a source. It is therefore critically important to track the spatiotemporal changes in burned area and fire carbon emissions over time. Here we developed a new burned-area detection algorithm between 2001–2019 across Alaska and Canada at 500 m (meters) resolution that utilizes finer-scale 30 m Landsat imagery to account for land cover unsuitable for burning. This method strictly balances omission and commission errors at 500 m to derive accurate landscape- and regional-scale burned-area estimates. Using this new burned-area product, we developed statistical models to predict burn depth and carbon combustion for the same period within the NASA Arctic–Boreal Vulnerability Experiment (ABoVE) core and extended domain. Statistical models were constrained using a database of field observations across the domain and were related to a variety of response variables including remotely sensed indicators of fire severity, fire weather indices, local climate, soils, and topographic indicators. The burn depth and aboveground combustion models performed best, with poorer performance for belowground combustion. We estimate 2.37×106 ha (2.37 Mha) burned annually between 2001–2019 over the ABoVE domain (2.87 Mha across all of Alaska and Canada), emitting 79.3 ± 27.96 Tg (±1 standard deviation) of carbon (C) per year, with a mean combustion rate of 3.13 ± 1.17 kg C m -2. Mean combustion and burn depth displayed a general gradient of higher severity in the northwestern portion of the domain to lower severity in the south and east. We also found larger-fire years and later-season burning were generally associated with greater mean combustion. Our estimates are generally consistent with previous efforts to quantify burned area, fire carbon emissions, and their drivers in regions within boreal North America; however, we generally estimate higher burned area and carbon emissions due to our use of Landsat imagery, greater availability of field observations, and improvements in modeling. The burned area and combustion datasets described here (the ABoVE Fire Emissions Database, or ABoVE-FED) can be used for local- to continental-scale applications of boreal fire science. [ABSTRACT FROM AUTHOR]

Forest managers across the Canadian boreal require detailed fire pattern information to support disturbance-based management. However, there are no consistent classifications of post-fire patterns, and those that exist rely on field-data that is both expensive and lacking in spatial representation. As a result, across the managed boreal forest there is limited fire pattern information, no standardized protocols to derive the quantitative fire pattern information nor agreement on the most appropriate metrics needed to implement management planning based on emulating natural disturbances. In this paper we first derived three classes of tree mortality based on Landsat spectral data and aerially photo-interpreted (API) polygons. Next we used these derived classes to compute seven fire pattern metrics over 14 fires that occurred from 1984 to 2006 in the Canadian provinces of Alberta and Saskatchewan, Canada. We then compared the metrics derived from the Landsat data to those computed from the API coverage. Lastly, we assessed the influence of a suite of environmental and fire pattern variables on the Landsat-derived indices. Overall we found very close correspondence for the two event-scale indices: the total area affected by fire and the complexity of the perimeter. The more detailed within-fire event metrics presented more varied results and were somewhat less precise. For example, the total amount of residuals or the largest disturbed patch were accurately captured. Other indices such as the number of disturbed patches or total amount of island remnants presented moderate systematic biases, but still might be serviceable given that the biases were predictable in direction. Fire size and seasonal drought variables were the most correlated with the Landsat derived pattern indices and their inclusion as model variables is likely to increase overall prediction accuracy. A pan-boreal geospatial database of detailed fire pattern metrics would be an invaluable tool to support the implementation of disturbance-based management approaches. Here we demonstrate that a cost-effective Landsat methods produces comparable fire pattern information to conventional approaches based on manual interpretation of aerial photographs. Building upon the standardized methods proposed here hundreds of fires could be mapped, to potentially, create a comprehensive national fire pattern database. [ABSTRACT FROM AUTHOR]

Avulsions change river courses and transport water and sediment to new channels impacting infrastructure, floodplain evolution, and ecosystems. Abrupt avulsion events (occurring over days to weeks) are potentially catastrophic to society and thus receive more attention than slow avulsions, which develop over decades to centuries and can be challenging to identify. Here, we examine gradual channel changes of the Peace‐Athabasca River Delta (PAD), Canada using in situ measurements and 37 years of Landsat satellite imagery. A developing avulsion of the Athabasca River is apparent along the Embarras River–Mamawi Creek (EM) distributary. Its opening and gradual enlargement since 1982 are evident from multiple lines of observation: Between 1984 and 2021 the discharge ratio between the EM and the Athabasca River more than doubled, increasing from 9% to 21%. The EM has widened by +53% since 1984, whereas the Athabasca River channel width has remained stable. The downstream Mamawi Creek delta is growing at a discharge‐normalized rate roughly twice that of the Athabasca River delta in surface area. Longitudinal global navigation satellite systems field surveys of water surface elevation reveal the EM possesses a ∼2X slope advantage (8 × 10−5 vs. 4 × 10−5) over the Athabasca River, and unit stream power and bed shear stress suggest enhanced sediment transport and erosional capacity through the evolving flow path. Our findings: (a) indicate that a slow avulsion of the Athabasca River is underway with potentially long‐term implications for inundation patterns, ecosystems, and human use of the PAD; and (b) demonstrate an observational approach for identifying other slow avulsions at river bifurcations globally. Plain Language Summary: Avulsions shift river courses and move water and sediment to new channels, which affect infrastructure, floodplains, and ecosystems. Slow avulsions take decades to develop and are more difficult to identify. Using on‐the‐ground measurements and 37 years of Landsat satellite imagery, we analyze gradual channel changes in the Peace‐Athabasca River Delta (PAD), Canada. The Athabasca River is changing course such that more of its water enters its westernmost outlet, the Embarras River–Mamawi Creek (EM) channel. Multiple lines of evidence demonstrate that the EM channel has been gradually opening since 1982. Between 1984 and 2021, the water entering the EM channel increased from 9% to 21% of the river's total flow. Since 1984, the EM channel has widened by 53%, while the Athabasca River channel has remained stable. The delta forming at the EM mouth (i.e., Mamawi Creek delta) has grown twice as fast as the Athabasca River delta. Field measurements of water surface elevation show the slope of the EM channel is twice as steep as the slope of the lower Athabasca River (8 × 10−5 vs. 4 × 10−5). Because water tends to flow down the steepest slope, we expect more water to flow down the EM channel in the future. Our findings indicate a slow capture of Athabasca River water into its EM channel, with potential long‐term implications for the delta's inundation pattern, ecosystems, and traditional Indigenous activities. Key Points: We assess a potential avulsion of the Athabasca River in the Peace‐Athabasca Delta, Canada using field measurements and remote sensingAnalysis of hydrological and morphological observations affirm that a slow avulsion is currently underwayThe avulsion may accelerate in the future and cause transformative effects on the delta's vegetation, habitat, and ecosystems [ABSTRACT FROM AUTHOR]

Devastating floods occur regularly around the world. Recently, machine learning models have been used for flood susceptibility mapping. However, even when these algorithms are provided with adequate ground truth training samples, they can fail to predict flood extends reliably. On the other hand, the height above nearest drainage (HAND) model can produce flood prediction maps with limited accuracy. The objective of this research is to produce an accurate and dynamic flood modeling technique to produce flood maps as a function of water level by combining the HAND model and machine learning. In this paper, the HAND model was utilized to generate a preliminary flood map; then, the predictions of the HAND model were used to produce pseudo training samples for a R.F. model. To improve the R.F. training stage, five of the most effective flood mapping conditioning factors are used, namely, Altitude, Slope, Aspect, Distance from River and Land use/cover map. In this approach, the R.F. model is trained to dynamically estimate the flood extent with the pseudo training points acquired from the HAND model. However, due to the limited accuracy of the HAND model, a random sample consensus (RANSAC) method was used to detect outliers. The accuracy of the proposed model for flood extent prediction, was tested on different flood events in the city of Fredericton, NB, Canada in 2014, 2016, 2018, 2019. Furthermore, to ensure that the proposed model can produce accurate flood maps in other areas as well, it was also tested on the 2019 flood in Gatineau, QC, Canada. Accuracy assessment metrics, such as overall accuracy, Cohen's kappa coefficient, Matthews correlation coefficient, true positive rate (TPR), true negative rate (TNR), false positive rate (FPR) and false negative rate (FNR), were used to compare the predicted flood extent of the study areas, to the extent estimated by the HAND model and the extent imaged by Sentinel-2 and Landsat satellites. The results confirm that the proposed model can improve the flood extent prediction of the HAND model without using any ground truth training data. [ABSTRACT FROM AUTHOR]

A variety of evidence suggests that the boreal forests of Canada are responding to climate change. Specifically, several studies have inferred that widespread browning trends detected in time series of the Normalized Difference Vegetation Index (NDVI) from the Advanced Very High Resolution Radiometer (AVHRR) reflect the response of boreal forests to longer growing seasons, increased summer drought stress, and higher frequency of fires. Data from the Thematic Mapper (TM5) and Enhanced Thematic Mapper Plus (ETM+) sensors onboard Landsat 5 and 7, respectively, span essentially the same time period as the AVHRR record, but provide data with substantially higher radiometric and spatial fidelity, and by extension, a much improved basis for evaluating decadal-scale trends in spectral vegetation indices such as the NDVI. However, detection of trends, which are often subtle, requires careful attention to ensure that artifacts associated with the quality and stability of inter- and intra-sensor calibration do not lead to spurious conclusions in results from time series analyses. In this paper, we use time series of TM5 and ETM+ images for fifteen sites distributed across the Canadian boreal forest zone to explore if and how sensor geometry and inter- and intra-sensor calibration affect trends in spectral vegetation indices derived from multi-decadal Landsat time series. To do this, we created annual cloud-free composites for each Landsat spectral band based on peak summer NDVI at each site from 1984 to 2011 using all available TM5 and ETM+ data. To distinguish trends arising from long term climate change from those related to disturbance, we isolated areas within each site that were undisturbed during the Landsat record, and used these locations to analyze sources of variance in time series of red reflectance, near-infrared (NIR) reflectance, the NDVI, and the Enhanced Vegetation Index (EVI). Our results highlight the challenges involved in distinguishing trends in surface properties from data artifacts caused by undetected atmospheric effects, changes in sensor view angles, and subtle radiometric differences between the TM5 and ETM+ sensors. In particular, differences in sensor view geometry across adjacent overlapping Landsat scenes cause vegetated pixels in the eastern portion of Landsat scenes to have higher reflectances in the red and NIR bands (by 5 and 6 percent, respectively) than pixels in the western portion of scenes. While this effect does not significantly change NDVI values, it does affect EVI values. We also found modest, but potentially significant, differences between the red band reflectance of each sensor, with TM5 data having 14 percent higher red reflectance on average for vegetated pixels, which can introduce spurious trends in time series that combine TM5 and ETM+ data. More generally, the results from this work demonstrate that while the 30 + year Landsat archive provides unprecedented opportunities for studying changes to the Earth's terrestrial biosphere over the last three decades, care must be taken when inferring trends in these data without considering how sources of variance unrelated to surface processes affect the integrity of Landsat time series. [ABSTRACT FROM AUTHOR]

New 'big data' streams such as street-level imagery are offering unprecedented possibilities for developing health-relevant data on the urban environment. Urban environmental features derived from street-level imagery have been used to assess pedestrian-friendly neighbourhood design and to predict active commuting, but few such studies have been conducted in Canada. Using 1.15 million Google Street View (GSV) images in seven Canadian cities, we applied image segmentation and object detection computer vision methods to extract data on persons, bicycles, buildings, sidewalks, open sky (without trees or buildings), and vegetation at postal codes. The associations between urban features and walk-to-work rates obtained from the Canadian Census were assessed. We also assessed how GSV-derived urban features perform in predicting walk-to-work rates relative to more widely used walkability measures. Results showed that features derived from street-level images are better able to predict the percent of people walking to work as their primary mode of transportation compared to data derived from traditional walkability metrics. Given the increasing coverage of street-level imagery around the world, there is considerable potential for machine learning and computer vision to help researchers study patterns of active transportation and other health-related behaviours and exposures. [ABSTRACT FROM AUTHOR]

Many of the methods used for estimating population size from ecological surveys have limitations on precision, cost, and/or applicability. The CountEm method was proposed recently for estimating the number of individuals in large groups from single images. It is simple and efficient, and can be applied to any species. Here we present a case study by applying CountEm to a real ecological survey with 278 images of Greater Snow Geese (Anser caerulescens atlanticus) and Common Eiders (Somateria mollissima) flocks taken from fixed‐wing aircraft in Eastern Canada. First, we evaluated the precision and counting time of CountEm on single images. Second, we developed and tested a new multi‐image version of the CountEm software. We show that flock sizes of N > 35,000 can be estimated on single images in ∼5 min, from counting a sample of ∼200 birds, yielding relative SEs in the 5%−10% range. Processing times increased to 10–20 min when simultaneously processing large numbers of images that contained over half a million birds with only modest increases in relative SE (range: 10%−15%). Our results suggest that CountEm may be used to save time and resources if incorporated into monitoring programs that utilize imagery in the abundance estimates. [ABSTRACT FROM AUTHOR]

Wildfires are a major natural disturbance in Canada that are postulated to increase under a warming climate. To derive accurate trends in burned area and to quantify the effects of fire frequency, duration, and extent, a sufficiently long time series of reliable burned area maps is required. With that in mind, we extended Canada's National Burned Area Composite (NBAC) dataset from its previous start year (2004) back to 1986. NBAC consists of annual maps in polygon format where the area burned in each fire event is represented by the best available delineation among various mapping methods and data sources of varying quality. Ordered from more to less reliability, in the new 35-year time series (1986 to 2020), 10% of the total burned area was derived from airborne and high-resolution (<5 m) satellite imagery, 81% from change detection methods using 30 m Landsat satellite imagery, and the remaining 9% was largely from aerial surveys. Total (Canada-wide) annual burned area estimates ranged from 215,797 ha in 2020 to 6.7 million ha (Mha) in 1989. We computed 95% confidence intervals for the estimate of each year from 1986 to 2020 based on the accuracy and relative contribution in that year of the different data sources, for both the new NBAC time series and the polygon version of the Canadian National Fire Database (CNFDB), a commonly used source of spatially explicit data on burned area in Canada. NBAC confidence intervals were on average ±9.7% of the annual figure, about one-third the width of the confidence intervals derived for CNFDB. The NBAC time series also included nearly 5000 fire events (totalling 4 Mha, with the largest event being 120,661 ha in size) that are missing in the CNFDB. In a regional analysis for the Northwest Territories, retroactive fire mapping from Landsat imagery reduced historical estimates by 3 Mha (16%), which would result in a 1.6 Mha increase in the reported undisturbed critical habitat for threatened woodland caribou. The NBAC dataset is freely downloadable from the Canadian Wildland Fire Information System. [ABSTRACT FROM AUTHOR]

Understanding the development of landscape patterns over broad spatial and temporal scales is a major contribution to ecological sciences and is a critical area of research for forested land management. Boreal forests represent an excellent case study for such research because these forests have undergone significant changes over recent decades. We analyzed the temporal trends of four widely-used landscape pattern indices for boreal forests of Canada: forest cover, largest forest patch index, forest edge density, and core (interior) forest cover. The indices were computed over landscape extents ranging from 5,000 ha (n = 18,185) to 50,000 ha (n = 1,662) and across nine major ecozones of Canada. We used 26 years of Landsat satellite imagery to derive annualized trends of the landscape pattern indices. The largest declines in forest cover, largest forest patch index, and core forest cover were observed in the Boreal Shield, Boreal Plain, and Boreal Cordillera ecozones. Forest edge density increased at all landscape extents for all ecozones. Rapidly changing landscapes, defined as the 90th percentile of forest cover change, were among the most forested initially and were characterized by four times greater decrease in largest forest patch index, three times greater increase in forest edge density, and four times greater decrease in core forest cover compared with all 50,000 ha landscapes. Moreover, approximately 18% of all 50,000 ha landscapes did not change due to a lack of disturbance. The pattern database results provide important context for forest management agencies committed to implementing ecosystem-based management strategies. [ABSTRACT FROM AUTHOR]

Boreal woodland caribou (Rangifer tarandus caribou) are designated as threatened under Canada's Species at Risk Act. A Recovery Strategy for the boreal population of caribou identified critical habitat for all but 1 of 51 caribou ranges - Saskatchewan's Boreal Shield (SK1). The strategy identified 65% undisturbed habitat as the threshold below which a local population was not likely to be self-sustaining. Disturbance was measured as the combined effects of fires, 40 years and anthropogenic land use. The fire component of the total disturbance model used fire polygons that were delineated using traditional mapping methods. Our study maps fire from 1988-2013 using the differenced Normalized Burn Ratio analysis of Landsat Thematic mapper and Operational Land Imager. Annual burned areas based on fire perimeters were similar between traditionally and Landsat-derived inventory approaches, but the traditional methods overestimated within-burn areas by 31.8%, as a result of including post-fire residuals and water bodies as burned. The federal recovery model assumes that all lands within provincial fire polygons (,40 years) are inadequate as caribou habitat, and ignores the potential value of post-fire residuals and water bodies as habitat. For some Boreal Shield ranges including SKI, where fire comprises the majority of the total disturbance and residual patches are abundant, total disturbance calculations, critical habitat designation and range planning decisions should take into account residuals, including water bodies. [ABSTRACT FROM AUTHOR]

Past emission of metal-bearing particulate matter, sulphur dioxide (SO2), and sulphuric acid by base metal smelters in the Sudbury region led to widespread loss of vegetation, contamination of soils, and formation of black coatings on rock surfaces. These black coatings formed through the incorporation of smelter-borne particulate matter into the partly dissolved uppermost layers of siliceous minerals on exposed rock, and are characterized by high heavy-metal content. This study involved assessment of the reflectance properties of black coatings in the Sudbury region, and determination of the geographic distribution of coatings through supervised classification of reflectance data derived from a Landsat Enhanced Thematic Mapper Plus (ETM+) image. Classifications involved the use of the Spectral Angle Mapper (SAM), Maximum Likelihood, and Feedforward Backpropagation Neural Network algorithms. The reflectance spectra of black coatings in the Sudbury region are relatively flat and featureless, and are characterized by reflectance values less than ~13% across the visible, near-infrared, and short-wave infrared. Spectral properties are similar to those of magnetite, a spinel-group mineral known to be present in Sudbury coatings. The presence of carbon-rich soot particles may be an important influence on the reflectance properties of coatings. SAM classification results are characterized by the widespread mislabelling of uncoated urban and open-pit sites as mantled by black coatings, and neural network results problematically mislabel some uncoated wetland sites as coated. Results generated by the Maximum Likelihood algorithm most usefully depict the distribution of exposed black coatings in the Sudbury region. The mapping of black coatings using remote-sensing methods can provide useful information on the spatial character of environmental degradation in the vicinity of smelters, and should be helpful in the monitoring of environmental recovery where emissions have been reduced or eliminated. [ABSTRACT FROM AUTHOR]

The Nagagamisis Central Plateau (located in Northern Ontario, Canada) is an area of distinct natural and cultural significance. The importance of this land was officially recognized in 1957 through the establishment of the Nagagamisis Provincial Park Reserve. The park has experienced significant expansion since its inception and is currently under development as one of Ontario Parks 'Signature Sites'. Since the 1980s, timber harvest activity has led to widespread forest disturbance just outside of the park boundaries. This research is focused on the detection of stand level forest disturbances associated with timber harvest occurring near Nagagamisis Provincial Park. The image time-series data selected for this project were Landsat TM and ETM+; spanning a twenty-five year period from 1984 to 2009. The Tasselled Cap Transformation and Normalized Difference Moisture Index were derived for use in unsupervised image classification to determine the land cover for each image in the time-series. Image band differencing and raster arithmetic were performed to create maps illustrating the size and spatial distribution of stand level forest disturbances between image dates. A total area of 1649 km² or 26.1% of the study area experienced stand level disturbance during the analysis period. [ABSTRACT FROM AUTHOR]

We recently completed the accuracy assessment of a Landsat-derived landcover polygon layer covering the entire province of Alberta (660,000 km2), Canada, for which we gathered reference information for nearly 5000 randomly selected polygons ranging from two hectares to thousands of hectares in size. This gave us the unique opportunity to quantify, for the first time, how the probability of correctly classifying a landcover object varies with its size. Irrespective of whether they are represented as polygons or as sets of connected pixels with the same label, the classification accuracy of landcover objects decreases as their size decreases, steadily for large and medium sizes, and more dramatically when they are within two orders of magnitude of the pixel size of the input image. We show that this size-dependency is bound to occur whenever the size distribution of landcover objects follows an inverse power law. Our results are consistent with previous studies on related issues, confirm the need to account for size when assessing the accuracy of object-based landcover maps, and cast doubts on the validity of (1) recently proposed object-based accuracy estimators, and (2) landscape pattern analyses where the minimum patch size is close to the pixel size. [ABSTRACT FROM AUTHOR]

The analysis of a series of five normalized difference vegetation index (NDVI) images produced information about a Labrador (Canada) portion of the tundra-taiga interface. The twenty-five year observation period ranges from 1983 to 2008. The series composed of Landsat, SPOT and ASTER images, provided insight into regional scale characteristics of the tundra-taiga interface that is usually monitored from coarse resolution images. The image set was analyzed by considering an ordinal classification of the NDVI to account for the cumulative effect of differences of near-infrared spectral resolutions, the temperature anomalies, and atmospheric conditions. An increasing trend of the median values in the low, intermediate and high NDVI classes is clearly marked while accounting for variations attributed to cross-sensor radiometry, phenology and atmospheric disturbances. An encroachment of the forest on the tundra for the whole study area was estimated at 0 to 60 m, depending on the period of observation, as calculated by the difference between the median retreat and advance of an estimated location of the tree line. In small sections, advances and retreats of up to 320 m are reported for the most recent four- and seven-year periods of observations. [ABSTRACT FROM AUTHOR]

Many countries undertake a national forest inventory to enable statistically valid monitoring in support of national and international reporting of forest conditions and change. Canada's National Forest Inventory (NFI) program is designed to operate on a 10-year remeasurement cycle, with an interim report produced at the 5-year mid-point. The NFI is a sample-based inventory, with approximately 18,850 2 ×2-km photo plots across the country, distributed on a 20×20-km grid of sample points; these photo plots are the primary data source for the NFI. Capacity to provide annual monitoring information is required to keep policy and decision makers apprised of current forest conditions. In this study, we implemented a multistage monitoring framework and used a Moderate Resolution Imaging Spectroradiometer (MODIS) change product to successfully identify 78% of the changes in forest cover area that were captured with a Landsat change detection approach. Of the NFI photo plots that were identified by both the Landsat and MODIS approaches as having changes in forest cover, the proportion of change area within the plots was similar ( R = 0.78). Approximately 70% of the Landsat-derived change events occupied less than 40% of a single MODIS pixel, and more than 90% of the change events of this size were successfully detected with the MODIS product. Finally, MODIS estimates of the proportion of forest cover change at the NFI photo plot level were comparable to change estimates for the ecoregions as a whole ( R = 0.95). High-temporal, low-spatial resolution imagery such as MODIS, in combination with other remotely sensed data sources, can provide information on disturbance events within a national forest inventory remeasurement cycle, thereby satisfying the interim information needs of policy and decision makers as well as the requirements of national and international reporting commitments. [ABSTRACT FROM AUTHOR]

The present extent of glacier ice on Bylot Island, Arctic Canada, is mapped using high-resolution Landsat 7 ETM+ satellite imagery. The island is 43% ice covered, with 4783 km² of ice. Most ice is centered on the northwest-southeast--trending Byam Martin Mountains, flowing outward as radial valley glaciers and piedmont lobes. The largest glacier is 49 km long and 6.5 km wide. The majority of glaciers terminate on land, but many have margins ending in lakes and two calve into the sea. The late summer snowline, mapped from satellite imagery, is highest along the southern and central parts of the island at about 1050 m, with lower values along the east-northeastern margin of the ice down to about 700 m. These snowline-elevation differences suggest a predominant moisture source from the northeast. Several valley glaciers and piedmont lobes have deformed medial moraines and ice-surface foliation suggesting past surge activity. Ten glaciers are interpreted to be of possible surgetype. The modern extent of glaciers is compared with that of two earlier time intervals. First, we have mapped glacier margins in several areas of Bylot Island from aerial photographs acquired in 1958 and 1961. Secondly, former positions of ice fronts are mapped from moraine systems deposited during the Neoglacial maximum and identified on satellite data. Glaciers have retreated from 0.9 to 1.8 km since the Neoglacial maximum about 120 years ago, with most retreat occurring between 1958/1961 and 2001. Approximately 253 km² or 5% of the 1958/1961 ice-covered area has been lost. Overall, marked glacier retreat has occurred, although a few glaciers, possibly of surge-type, show small readvances. This retreat is consistent with observed climate warming in the Canadian Arctic, especially since the 1960s. [ABSTRACT FROM AUTHOR]

Posterior probabilities of occurrence for Zn-Pb Mississippi Valley Type (MVT) mineralization were calculated based on evidence maps derived from regional geology, Landsat-TM, RADARSAT-1, a digital elevation model and aeromagnetic data sets in the Borden Basin of northern Baffin Island, Canada. The vector representation of geological contacts and fault traces were refined according to their characteristics identified in Landsat-TM, RADARSAT-1, DEM, slope, aspect, and shaded relief data layers. Within the study area, there is an association between the occurrence of MVT mineralization and proximity to the contact of platformal carbonates and shale units of the adjacent geological formation. A spatial association also tends to exist between mineralization and proximity to E-W and NW-SE trending faults. The relationships of known MVT occurrences with the geological features were investigated by spatial statistical techniques to generate evidence maps. Supervised classification and filtering were applied to Landsat-TM data to divide the Society Cliffs Formation into major stratigraphic subunits. Because iron oxides have been observed at some of the MVT occurrences within the Borden Basin, Landsat-TM data band ratio (3/1) was calculated to highlight the potential presence of iron-oxides as another evidence map. Processed Landsat-TM data and other derived geological evidence maps provided useful indicators for identifying areas of potential MVT mineralization. Weights of evidence and logistic regression were used independently to integrate and generate posterior probability maps showing areas of potential mineralization based on all derived evidence maps. Results indicate that in spite of the lack of important data sets such as stream or lake sediment geochemistry, Landsat-TM data and regional geological data can be useful for MVT mineral-potential mapping. [ABSTRACT FROM AUTHOR]

AimTo detect possible expansion of boreal forest stands in response to recent warming. Previous modelling studies have concluded that major shifts in vegetation patterns, including changes in boreal forest extent, could arise during the next two centuries under global warming scenarios. However, field investigations of tree stands at ecotones have so far revealed little indication of stand response to warming during the last 100 years. This study uses a c. 25-year record of Landsat satellite observations to quantify changes in forest stand extent in two areas of northern Canada. LocationTwo regions of northern Canada, near Richmond Gulf, Quebec, and Great Slave Lake, north-west Territories. MethodsNormalized-difference vegetation index (NDVI) plots across forest-tundra boundaries were obtained from radiometrically corrected Landsat imagery acquired during the 1970s and 1990s. These curves were evaluated to look for changes over the c. 25-year period related to forest stand expansion. ResultsAlthough forest-tundra boundaries could be clearly mapped from the satellite data, no obvious change in forest boundaries was apparent during the duration of the image time series, constraining recent geographical expansion rates to <200–300 m per century. Also, no evidence for local expansion of forest stands (e.g. within sheltered valleys) was found. Main conclusionsThe results are consistent with field observations, and suggest that, at the moment, boreal forest extents remain basically stable. This may reflect inherent lags between forest response and climate change, or competitive pressures between tree stands and surrounding tundra and herbaceous vegetation. [ABSTRACT FROM AUTHOR]