Search

A large portion of Canada is covered by wetlands; mapping and monitoring them is of great importance for various applications. In this regard, Remote Sensing (RS) technology has been widely employed for wetland studies in Canada over the past 45 years. This study evaluates meta-data to investigate the status and trends of wetland studies in Canada using RS technology by reviewing the scientific papers published between 1976 and the end of 2020 (300 papers in total). Initially, a meta-analysis was conducted to analyze the status of RS-based wetland studies in terms of the wetland classification systems, methods, classes, RS data usage, publication details (e.g., authors, keywords, citations, and publications time), geographic information, and level of classification accuracies. The deep systematic review of 128 peer-reviewed articles illustrated the rising trend in using multi-source RS datasets along with advanced machine learning algorithms for wetland mapping in Canada. It was also observed that most of the studies were implemented over the province of Ontario. Pixel-based supervised classifiers were the most popular wetland classification algorithms. This review summarizes different RS systems and methodologies for wetland mapping in Canada to outline how RS has been utilized for the generation of wetland inventories. The results of this review paper provide the current state-of-the-art methods and datasets for wetland studies in Canada and will provide direction for future wetland mapping research. [ABSTRACT FROM AUTHOR]

Since individual neural networks have limited deep expressiveness and effectiveness, many learning frameworks face difficulties in the availability and balance of sample selection. As a result, in change detection, it is difficult to upgrade the hit rate of a high-performance model on both positive and negative pixels. Therefore, supposing that the sacrificed components coincide perfectly with the important evaluation objectives, such as positives, it would lose more than gain. To address this issue, in this paper, we propose a multi-visual collaborative deep network (MV-CDN) served by three collaborative network members that consists of three subdivision approaches, the CDN with one collaborator (CDN-C), CDN with two collaborators (CDN-2C), and CDN with three collaborators (CDN-3C). The purpose of the collaborator is to re-evaluate the feature elements in the network transmission, and thus to translate the group-thinking into a more robust field of vision. We use three sets of public double-temporal hyperspectral images taken by the AVIRIS and HYPERION sensors to show the feasibility of the proposed schema. The comparison results have confirmed that our proposed schema outperforms the existing state-of-the-art algorithms on the three tested datasets, which demonstrates the broad adaptability and progressiveness of the proposal. [ABSTRACT FROM AUTHOR]

This paper proposes an improved model-based forest height inversion method for airborne L-band dual-baseline repeat-pass polarimetric synthetic aperture radar interferometry (PolInSAR) collections. A two-layer physical model with various volumetric scattering attenuation and dynamic motion properties is first designed based on the traditional Random Motion over Ground (RMoG) model. Related PolInSAR coherence functions with both volumetric and temporal decorrelations incorporated are derived, where the impacts of homogenous and heterogeneous attenuation and dynamic motion properties on the performance of forest height inversion were investigated by the Linear Volume Attenuation (LVA), Quadratic Volume Attenuation (QVA), Linear Volume Motion (LVM), and Quadratic Volume Motion (QVM) depictions in the volume layer. Dual-baseline PolInSAR data were acquired to increase the degree of freedom (DOF) of the coherence observations and thereby provide extra constraints on the forest parameters to address the underdetermined problem. The experiments were carried out on a boreal forest in Canada and a tropical one in Gabon, where physical models with LVA + QVM (RMSE: 3.56 m) and QVA + LVM (RMSE: 6.83 m) exhibited better performances on the forest height inversion over the boreal and tropical forest sites, respectively. To leverage the advantages of LVA, QVA, LVM, and QVM, a pixel-wise optimization strategy was used to obtain the best forest height inversion performance for the range of attenuation and motion profiles considered. This pixel-wise optimization surpasses the best-performing single model and achieves forest height inversion results with an RMSE of 3.21 m in the boreal forest site and an RMSE of 6.48 m in the tropical forest site. [ABSTRACT FROM AUTHOR]

Crop identification and classification are of great significance to agricultural land use management. The physically constrained general model-based decomposition (PCGMD) has proven to be a promising method in comparison with the typical four-component decomposition methods in scattering mechanism interpretation and identifying vegetation types. However, the robustness of PCGMD requires further investigation from the perspective of final applications. This paper aims to validate the efficiency of the PCGMD method on crop classification for the first time. Seven C-band time-series RADARSAT-2 images were exploited, covering the entire growing season over an agricultural region near London, Ontario, Canada. Firstly, the response and temporal evolution of the four scattering components obtained by PCGMD were analyzed. Then, a forward selection approach was applied to achieve the highest classification accuracy by searching an optimum combination of multi-temporal SAR data with the random forest (RF) algorithm. For comparison, the general model-based decomposition method (GMD), the original and its three improved Yamaguchi four-component decomposition approaches (Y4O, Y4R, S4R, G4U), were used in all tests. The results reveal that the PCGMD method is highly sensitive to seasonal crop changes and matches well with the real physical characteristics of the crops. Among all test methods used, the PCGMD method using six images obtained the optimum classification performance, reaching an overall accuracy of 91.83%. [ABSTRACT FROM AUTHOR]

Canada's successful space-based earth-observation (EO) radar program has earned widespread and expanding user acceptance following the launch of RADARSAT-1 in 1995. RADARSAT-2, launched in 2007, while providing data continuity for its predecessor's imaging capabilities, added new polarimetric modes. Canada's follow-up program, the RADARSAT Constellation Mission (RCM), launched in 2019, while providing continuity for its two predecessors, includes an innovative suite of polarimetric modes. In an effort to make polarimetry accessible to a wide range of operational users, RCM uses a new method called hybrid compact polarization (HCP). There are two essential elements to this approach: (1) transmit only one polarization, circular; and (2) receive two orthogonal polarizations, for which RCM uses H and V. This configuration overcomes the conventional dual and full polarimetric system limitations, which are lacking enough polarimetric information and having a small swath width, respectively. Thus, HCP data can be considered as dual-pol data, while the resulting polarimetric classifications of features in an observed scene are of comparable accuracy as those derived from the traditional fully polarimetric (FP) approach. At the same time, RCM's HCP methodology is applicable to all imaging modes, including wide swath and ScanSAR, thus overcoming critical limitations of traditional imaging radar polarimetry for operational use. The primary image data products from an HCP radar are different from those of a traditional polarimetric radar. Because the HCP modes transmit circularly polarized signals, the data processing to extract polarimetric information requires different approaches than those used for conventional linearly polarized polarimetric data. Operational users, as well as researchers and students, are most likely to achieve disappointing results if they work with traditional polarimetric processing tools. New tools are required. Existing tutorials, older seminar notes, and reference papers are not sufficient, and if left unrevised, could succeed in discouraging further use of RCM polarimetric data. This paper is designed to provide an initial response to that need. A systematic review of studies that used HCP SAR data for environmental monitoring is also provided. Based on this review, HCP SAR data have been employed in oil spill monitoring, target detection, sea ice monitoring, agriculture, wetland classification, and other land cover applications. [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]

Wildfires play a crucial role in the transformation of forest ecosystems and exert a significant influence on the global climate over geological timescales. Recent shifts in climate patterns and intensified human–forest interactions have led to an increase in the incidence of wildfires. These fires are characterized by their extensive coverage, higher frequency, and prolonged duration, rendering them increasingly destructive. To mitigate the impact of wildfires on climate change, ecosystems, and biodiversity, it is imperative to conduct systematic monitoring of wildfire progression and evaluate their environmental repercussions on a global scale. Satellite remote sensing is a powerful tool, offering precise and timely data on terrestrial changes, and has been extensively utilized for wildfire identification, tracking, and impact assessment at both local and regional levels. The Canada Centre for Mapping and Earth Observation, in collaboration with the Canadian Forest Service, has developed a comprehensive National Burned Area Composite (NBAC). This composite serves as a benchmark for curating a bi-temporal multi-source satellite image dataset for change detection, compiled from the archives of Sentinel-2, Sentinel-1, and ALOS-2 PALSAR-2. To our knowledge, this dataset is the inaugural large-scale, multi-source, and multi-frequency satellite image dataset with 20 m spatial resolution for wildfire mapping, monitoring, and evaluation. It harbors significant potential for enhancing wildfire management strategies, building upon the profound advancements in deep learning that have contributed to the field of remote sensing. Based on our curated dataset, which encompasses major wildfire events in Canada, we conducted a systematic evaluation of the capability of multi-source satellite earth observation data in identifying wildfire-burned areas using statistical analysis and deep learning. Our analysis compares the difference between burned and unburned areas using post-event observation solely or bi-temporal (pre- and post-event) observations across diverse land cover types. We demonstrate that optical satellite data yield higher separability than C-Band and L-Band Synthetic Aperture Radar (SAR), which exhibit considerable overlap in burned and unburned sample distribution, as evidenced by SAR-based boxplots. With U-Net, we further explore how different input channels influence the detection accuracy. Our findings reveal that deep neural networks enhance SAR's performance in mapping burned areas. Notably, C-Band SAR shows a higher dependency on pre-event data than L-Band SAR for effective detection. A comparative analysis of U-Net and its variants indicates that U-Net works best with single-sensor data, while the late fusion architecture marginally surpasses others in the fusion of optical and SAR data. Accuracy across sensors is highest in closed forests, with sequentially lower performance in open forests, shrubs, and grasslands. Future work will extend the data from both spatial and temporal dimensions to encompass varied vegetation types and climate zones, furthering our understanding of multi-source and multi-frequency satellite remote sensing capabilities in wildfire detection and monitoring. [ABSTRACT FROM AUTHOR]

Loss of forest cover and derived effects on forest ecosystems services has led to the establishment of land management policies and forest monitoring systems, and consequently to the demand for accurate and multitemporal data on forest extent and structure. In recent years, spaceborne Light Detection and Ranging (LiDAR) missions, such as the Global Ecosystem Dynamics Investigation (GEDI) instrument, have facilitated the repeated acquisition of data on the vertical structure of vegetation. In this study, we designed an approach incorporating GEDI and airborne LiDAR data, in addition to detailed forestry inventory data, for estimating tree-growth dynamics for the Laurentides wildlife reserve in Canada. We estimated an average tree-growth rate of 0.32 ± 0.23 (SD) m/year for the study site and evaluated our results against field data and a time series of NDVI from Landsat images. The results are in agreement with expected patterns in tree-growth rates related to tree species and forest stand age, and the produced dataset is able to track disturbance events resulting in the loss of canopy height. Our study demonstrates the benefits of using spaceborne-LiDAR data for extending the temporal coverage of forestry inventories and highlights the ability of GEDI data for detecting changes in forests' vertical structure. [ABSTRACT FROM AUTHOR]

Introducing ionospheric information into a precise point positioning (PPP) solution enables faster ambiguity resolution and significantly improves positioning accuracy. To compute such corrections over wide areas, sparse networks with potentially irregular station distributions are often used. This aspect brings a new level of complexity as ionospheric corrections should be weighted appropriately in the PPP filter. This paper presents a possible implementation of grid-based wide-area slant ionospheric delay corrections, with a focus on the reported uncertainties. A balance is obtained between obtaining corrections with formal errors small enough to enable fast convergence, while large enough to overbound most errors. Based on least-squares collocation, the method uses satellite-specific variograms based on the 99th percentile values in each distance bin. Tested in southern Canada over a 53-week period in 2020, ionospheric grids allowed dual-frequency receivers to obtain around 5 cm accuracy in each horizontal component within 5 min of static data collection. For single-frequency solutions using data from geodetic receivers, positioning errors were reduced by over 60% for both static and kinematic processing. [ABSTRACT FROM AUTHOR]

Predictive accuracy in wildland fire behavior is contingent on a thorough understanding of the 3D fuel distribution. However, this task is complicated by the complex nature of fuel forms and the associated constraints in sampling and quantification. In this study, twelve terrestrial laser scanning (TLS) plot scans were sampled within the mountain pine beetle-impacted forests of Jasper National Park, Canada. The TLS point clouds were delineated into eight classes, namely individual-tree stems, branches, foliage, downed woody logs, sapling stems, below-canopy branches, grass layer, and ground-surface points using a transformer-based deep learning classifier. The fine-scale 3D architecture of trees and branches was reconstructed using a quantitative structural model (QSM) based on the multi-class components from the previous step, with volume attributes extracted and analyzed at the branch, tree, and plot levels. The classification accuracy was evaluated by partially validating the results through field measurements of tree height, diameter-at-breast height (DBH), and live crown base height (LCBH). The extraction and reconstruction of 3D wood components enable advanced fuel characterization with high heterogeneity. The existence of ladder trees was found to increase the vertical overlap of volumes between tree branches and below-canopy branches from 8.4% to 10.8%. [ABSTRACT FROM AUTHOR]

"Snow gone" (SGN) stage is one of the critical variables that describe the start of the official forest fire season in the Canadian Province of Alberta. In this paper, our objective is to evaluate the potential of MODIS-based indices for determining the SGN stage. Those included: (i) enhanced vegetation index (EVI), (ii) normalized difference water index (NDWI) using the shortwave infrared (SWIR) spectral bands centered at 1.64 μm (NDWI1.64μm) and at 2.13 μm (NDWI2.13μm), and (iii) normalized difference snow index (NDSI). These were calculated using the 500 m 8-day gridded MODIS-based composites of surface reflectance data (i.e., MOD09A1 v.005) for the period 2006-08. We performed a qualitative evaluation of these indices over two forest fire prone natural subregions in Alberta (i.e., central mixedwood and lower boreal highlands). In the process, we generated and compared the natural subregion-specific lookout tower sites average: (i) temporal trends for each of the indices, and (ii) SGN stage using the ground-based observations available from Alberta Sustainable Resource Development. The EVI-values were found to have large uncertainty at the onset of the spring and unable to predict the SGN stages precisely. In terms of NDSI, it showed earlier prediction capabilities. On the contrary, both of the NDWI's showed distinct pattern (i.e., reached a minimum value before started to increase again during the spring) in relation to observed SGN stages. Thus further analysis was carried out to determine the best predictor by comparing the NDWI's predicted SGN stages with the ground-based observations at all of the individual lookout tower sites (approximately 120 in total) across the study area. It revealed that NDWI2.13μm demonstrated better prediction capabilities (i.e., on an average approximately 90% of the observations fell within ±2 periods or ±16 days of deviation) in comparison to NDWI1.64μm (i.e., on an average approximately 73% of the observations fell within ±2 periods or ±16 days of deviation). [ABSTRACT FROM AUTHOR]

This study examines the response of a cold-regions deltaic wetland ecosystem in northwestern Canada to two separate and differing seasonal wetting cycles. The goal of this paper was to examine the nature of reflected electromagnetic energy measured by earth observation (EO) satellites, and to assess whether seasonal wetland hydroperiod and episodic flooding events impact the information retrieved by the Sentinel-2 sensors. The year 2018 represents a year characterized by a large spring freshet and ice-jam flooding, while 2019 represents a year characterized more by summer open-water flooding. We applied the Modified Normalized Difference Wetness Index (MNDWI) to address the effects of the wetting cycles. The response of the vegetative cover was tracked using the fraction of the absorbed photosynthetically active radiation (fAPAR) and the Leaf Area Index (LAI). All three indices were viewed through the lens of cover classes as derived through a previously published study by the authors. The study provides a framework for designing longer-term studies where multiple intra- and inter-annual hydrological cycles can be accessed via EO data. Future studies will enable the examination of lag times inherent in the response to the various water sources applied to spectral response and incorporate this EO approach into a monitoring framework. [ABSTRACT FROM AUTHOR]

While many studies have been conducted regarding wind-driven Ka-band scattering on the ocean and sea surfaces, few have identified the impacts of Ka-band scattering on small inland water bodies, and fewer have identified the influence of wind on coherence over water. These previous studies have been limited in spatial scale, covering only large water bodies >25 km2. The recently launched Surface Water and Ocean Topography (SWOT) mission is the first Ka-band InSAR satellite designed for mapping water surface elevations and open water areas for rivers as narrow as 100 m and lakes as small as 0.0625 km2. Because measurements of these types are novel, there remains some uncertainty about expected backscatter amplitudes given wind-driven water surface roughness variability. A previous study using the airborne complement to SWOT, AirSWOT, found that low backscatter and low coherence values were indicative of higher errors in the water surface elevation products, recommending minimum thresholds for backscatter and coherence for filtering the data to increase the accuracy of averaged data for lakes and rivers. We determined that the global average wind speed over lakes is 4 m/s, and after comparing AirSWOT backscatter and coherence data with ERA-5 wind speeds, we found that the minimum required speed to retrieve high backscatter and coherence is 3 m/s. We examined 11,072 lakes across Canada and Alaska, with sizes ranging from 350 m2 to 156 km2, significantly smaller than what could be measured with previous Ka-band instruments in orbit. We found that small lakes (0.0625–0.25 km2) have significantly lower backscatter (3–5 dB) and 0.20–0.25 lower coherence than larger lakes (>1 km2). These results suggest that approximately 75% of SWOT observable lake areas around the globe will have consistently high-accuracy water surface elevations, though seasonal wind variability should remain an important consideration. Despite very small lakes presenting lower average backscatter and coherence, this study asserts that SWOT will be able to accurately resolve the water surface elevations and water surface extents for significantly smaller water bodies than have been previously recorded from satellite altimeters. This study additionally lays the foundation for future high-resolution inland water wind speed studies using SWOT data, when the data become available, as the relationships between wind speed and Ka-band backscatter reflect those of traditional scatterometers designed for oceanic studies. [ABSTRACT FROM AUTHOR]

Rockfall constitutes a major threat to the safety and sustainability of transport corridors bordered by rocky cliffs. This research introduces a new approach to rockfall susceptibility modeling for the identification of potential rockfall source zones. This is achieved by developing a data-driven model to assess the local slope morphological attributes with respect to the rock slope evolution processes. The ability to address "where" a rockfall is more likely to occur via the analysis of historical event inventories with respect to terrain attributes and to define the probability of a given area producing a rockfall is a critical advance toward effective transport corridor management. The availability of high-quality digital volumetric change detection products permits new developments in rockfall assessment and prediction. We explore the potential of simulating the conceptualization of slope-scale rockfall susceptibility modeling using computer power and artificial intelligence (AI). We employ advanced 3D computer vision algorithms for analyzing point clouds to interpret high-resolution digital observations capturing the rock slope evolution via long-term, LiDAR-based 3D differencing. The approach has been developed and tested on data from three rock slopes: two in Canada and one in the UK. The results indicate clear potential for AI advances to develop local susceptibility indicators from local geometry and learning from recent rockfall activity. The resultant models produce slope-wide rockfall susceptibility maps in high resolution, producing up to 75% agreement with validated occurrences. [ABSTRACT FROM AUTHOR]

Wetlands are important ecosystems that are linked to climate change mitigation. As 25% of global wetlands are located in Canada, accurate and up-to-date wetland classification is of high importance, nationally and internationally. The advent of deep learning techniques has revolutionized the current use of machine learning algorithms to classify complex environments, specifically in remote sensing. In this paper, we explore the potential and possible limitations to be overcome regarding the use of ensemble deep learning techniques for complex wetland classification and discusses the potential and limitation of various solo convolutional neural networks (CNNs), including DenseNet, GoogLeNet, ShuffleNet, MobileNet, Xception, Inception-ResNet, ResNet18, and ResNet101 in three different study areas located in Newfoundland and Labrador, Canada (i.e., Avalon, Gros Morne, and Grand Falls). Moreover, to improve the classification accuracies of wetland classes of bog, fen, marsh, swamp, and shallow water, the results of the three best CNNs in each study area is fused using three supervised classifiers of random forest (RF), bagged tree (BTree), Bayesian optimized tree (BOT), and one unsupervised majority voting classifier. The results suggest that the ensemble models, in particular BTree, have a valuable role to play in the classification of wetland classes of bog, fen, marsh, swamp, and shallow water. The ensemble CNNs show an improvement of 9.63–19.04% in terms of mean producer's accuracy compared to the solo CNNs, to recognize wetland classes in three different study areas. This research indicates a promising potential for integrating ensemble-based learning and deep learning for operational large area land cover, particularly complex wetland type classification. [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]

The increasing frequency of flooding worldwide has driven research to improve near real-time flood mapping from remote-sensing data. Improved automation and processing speed to map both open water and vegetated area flooding have resulted from these research efforts. Despite these achievements, flood mapping in urban areas where a significant number of overall impacts are felt remains a challenge. Near real-time data availability, shadowing caused by manmade infrastructure, spatial resolution, and cloud cover inhibiting optical transmission, are all factors that complicate detailed urban flood mapping needed to inform response efforts. This paper uses numerous data sources collected during two major flood events that impacted the same region of Eastern Canada in 2017 and 2019 to test different urban flood mapping approaches presented as case studies in three separate urban boroughs. Cloud-free high-resolution 3 m PlanetLab optical data acquired near peak-flood in 2019 were used to generate a maximum flood extent product for that year. Approaches using new Lidar Digital Elevation Models (DEM)s and water height estimated from nineteen RADARSAT-2 flood maps, point-based flood perimeter observations from citizen geographic information, and simulated traffic camera or other urban sensor network data were tested and verified using independent data. Coherent change detection (CCD) using multi-temporal Interferometric Wide (IW) Sentinel-1 data was also tested. Results indicate that while clear-sky high-resolution optical imagery represents the current gold standard, its availability is not guaranteed due to timely coverage and cloud cover. Water height estimated from 8 to 12.5 m resolution RADARSAT-2 flood perimeters were not sufficiently accurate to flood adjacent urban areas using a Lidar DEM in near real-time, but all nineteen scenes combined captured boroughs that flooded at least once in both flood years. CCD identified flooded boroughs and roughly captured their flood extents, but lacked timeliness and sufficient detail to inform street-level decision-making in near real-time. Point-based flood perimeter observation, whether from in-situ sensors or high-resolution optical satellites combined with Lidar DEMs, can generate accurate full flood extents under certain conditions. Observed point-based flood perimeters on manmade features with low topographic variation produced the most accurate flood extents due to reliable water height estimation from these points. [ABSTRACT FROM AUTHOR]

Sustainable development of Canada's North requires an increased focus on renewable, zero-emission energy sources. Burwash Landing in Yukon is prospective for geothermal energy based on a high geothermal gradient, local occurrence of warm groundwater and proximity to the active, crustal-scale Denali fault. Uncertainties about the potential geothermal system include the nature and geometry of fluid pathways, and heat sources required to drive a hydrothermal system. In this study, we inverted three passive electromagnetic datasets—321 extremely low frequency electromagnetic, 33 audiomagnetotelluric and 51 magnetotelluric stations—to map the subsurface electrical structure to 8 km depth. Our new model reveals vertical conductive structures associated with the two main faults, Denali and Bock's Creek, which we interpret to represent fluid-deposited graphite and hydrothermal alteration, respectively. Our model supports an interpreted releasing bend on the main Denali fault strand. This is associated with the deepest conductivity anomaly along the fault and potential for deeper penetration of fluids. Enigmatic conductive bodies from 1 to > 6 km depth are associated with intermediate to mafic intrusions. Fluids released from these bodies may advect heat and provide a possible heat source to mobilize hot fluids and sustain a geothermal system in the region. [ABSTRACT FROM AUTHOR]

To monitor polar regions is of paramount importance for climatological studies. Climate change due to anthropogenic activities is inducing global warming that, for example, has resulted in glacier melting. This has had a significant impact on sea levels and ocean circulation. In this study, the temporal trend of the marine-terminated d'Iberville glacier (Ellesmere Island, Canada) is analysed using C-band synthetic aperture radar satellite imagery collected by the Radarsat-2 and Sentinel-1 missions. The data set consists of a time series of 10 synthetic aperture radar data collected from 2010 to 2022 in dual-polarimetric imaging mode, where a horizontally polarised electromagnetic wave was transmitted. An automatic approach based on a global threshold constant false alarm rate method is applied to the single- and dual-polarisation features, namely the HH-polarised normalised radar cross-section and a combination of the HH- and HV-polarised scattering amplitudes, with the aim of extracting the ice front of the glacier and, therefore, estimating its behaviour over time. Independent collocated satellite optical imagery from the Sentinel-2 multi-spectral instrument is also considered, where available, to support the experimental outcomes. The experimental results show that (1) the HH-polarised normalised radar cross-section achieved better performance with respect to the dual-polarised feature, especially under the most challenging case of a sea-ice infested sea surface; (2) when the HH-polarised normalised radar cross-section was considered, the ice front extraction methodology provided a satisfactory accuracy, i.e., a root mean square error spanning from about 1.1 pixels to 3.4 pixels, depending on the sea-surface conditions; and (3) the d'Iberville glacier exhibited, during the study period, a significant retreat whose average surface velocity was 160 m per year, resulting in a net ice area loss of 2.2 km2 (0.18 km2 per year). These outcomes demonstrate that the d'Iberville glacier is behaving as most of the marine-terminated glaciers in the study area while experiencing a larger ice loss. [ABSTRACT FROM AUTHOR]

Despite providing vital ecosystem services, wetlands are increasingly threatened across the globe by both anthropogenic activities and natural processes. Synthetic aperture radar (SAR) has emerged as a promising tool for rapid and accurate monitoring of wetland extent and type. By acquiring information on the roughness and moisture content of the surface, SAR offers unique potential for wetland monitoring. However, there are still challenges in applying SAR for mapping complex wetland environments. The backscattering similarity of different wetland classes is one of the challenges. Choosing the appropriate SAR specifications (incidence angle, frequency and polarization), based on the wetland type, is also a subject of debate and should be investigated more thoroughly. The geometric distortion of SAR imagery and loss of coherency are other remaining challenges in applying SAR and its processing techniques for wetland studies. Hence, this study provides a systematic meta-analysis based on compilation and analysis of indexed research studies that used SAR for wetland monitoring. This meta-analysis reviewed 172 papers and documented an upward trend in usage of SAR data, increasing usage of multi-sensor data, increasing integration of C- and L- bands over other configurations and higher classification accuracy with multi-frequency and multi-polarized SAR data. The highest number of wetland research studies using SAR data came from the USA, Canada and China. This meta-analysis highlighted the current challenges and solutions for wetland monitoring using SAR sensors. [ABSTRACT FROM AUTHOR]

Accurate visualization of air temperature distribution can be useful for various thermal analyses in fields such as human health and heat transfer of local area. This paper presents a novel approach to measuring air temperature from midwave hyperspectral Fourier transform infrared (FTIR) imaging in the carbon dioxide absorption band (between 4.25–4.35 μ m). In this study, the proposed visual air temperature (VisualAT) measurement is based on the observation that the carbon dioxide band shows zero transmissivity at short distances. Based on analysis of the radiative transfer equation in this band, only the path radiance by air temperature survives. Brightness temperature of the received radiance can provide the raw air temperature and spectral average, followed by a spatial median-mean filter that can produce final air temperature images. Experiment results tested on a database obtained by a midwave extended FTIR system (Telops, Quebec City, QC, Canada) from February to July 2018 show a mean absolute error of 1.25 ∘ K for temperature range of 2.6−26.4 ∘ C. [ABSTRACT FROM AUTHOR]

Land cover change (LCC) studies are increasingly using deep learning (DL) modeling techniques. Past studies have leveraged temporal or spatiotemporal sequences of historical LC data to forecast changes with DL models. However, these studies do not adequately assess the association between neighborhood size and DL model capability to forecast LCCs, where neighborhood size refers to the spatial extent captured by each data sample. The objectives of this research study were to: (1) evaluate the effect of neighborhood size on the capacity of DL models to forecast LCCs, specifically Temporal Convolutional Networks (TCN) and Convolutional Neural Networks (CNN-TCN), and (2) assess the effect of auxiliary spatial variables on model capacity to forecast LCCs. First, each model type and neighborhood setting configuration was assessed using data derived from multitemporal MODIS LC for the Regional District of Bulkley-Nechako, Canada, comparing subareas exhibiting different amounts of LCCs with trends obtained for the full region. Next, outcomes were compared with three other study regions. The modeling results were evaluated with three-map comparison measures, where the real-world LC for the next timestep, the real-world LC for the previous timestep, and the forecasted LC for the next year were used to calculate correctly transitioned areas. Across all regions explored, it was observed that increasing neighborhood sizes improved the DL model's capabilities to forecast short-term LCCs. CNN–TCN models forecasted the most correct LCCs for several regions while reducing error due to quantity when provided additional spatial variables. This study contributes to the systematic exploration of neighborhood sizes on selected spatiotemporal DL techniques for geographic applications. [ABSTRACT FROM AUTHOR]

We designed a workflow to investigate areas of potential neotectonic deformation, making use of well-developed techniques, but applied to a site characterized by low relief and low or moderate tectonic activity. In this pilot study, we targeted the Temiskaming Graben, in Eastern Canada, where recent and ongoing geophysical and sedimentological investigations have revealed recent activity along this ancient structure. The dataset compiled for this experimental study covers an area of nearly 147 square km across the provinces of Ontario and Quebec. For efficiency in terms of computational resources, we first performed cluster analysis on knickpoint location, identifying seven areas with a high density of disruptions along river profiles. We then performed more detailed morphometric analysis at 30 m resolution, identifying knickpoints along river profiles, calculating the hypsometric integral across the landscape with a moving window, and mapping and comparing lineaments with known structural features. The results of our workflow showed that these three techniques can be efficiently combined for neotectonic analysis, and the synergistic approach strengthens the reliability and accuracy of our results. Our research extends the application of morphometric analysis, commonly used for exploring areas with intense tectonism and high topography, to areas that are characterized by low relief and low or moderate tectonic activity. The new areas identified with the workflow proposed in this research require ground-truthing through mapping and shallow geophysical investigations. [ABSTRACT FROM AUTHOR]

This study examines the angular distribution of scattered solar radiation associated with wildfire smoke aerosols observed over boreal forests in Canada during the ARCTAS (Arctic Research of the Composition of the Troposphere from Aircraft and Satellites) campaign. First, it estimates smoke radiative parameters (550 nm optical depth of 3.9 and single scattering albedo of 0.90) using quasi-simultaneous multiangular and multispectral airborne measurements by the Cloud Absorption Radiometer (CAR). Next, the paper estimates the broadband top-of-atmosphere radiances that a satellite instrument such as the Clouds and the Earth's Radiant Energy System (CERES) could have observed, given the narrowband CAR measurements made from an aircraft circling about a kilometer above the smoke layer. This estimation includes both an atmospheric correction that accounts for the atmosphere above the aircraft and a narrowband-to-broadband conversion. The angular distribution of estimated radiances is found to be substantially different than the angular model used in the operational data processing of CERES observations over the same area. This is because the CERES model is a monthly average model that was constructed using observations taken under smoke-free conditions. Finally, a sensitivity analysis shows that the estimated angular distribution remains accurate for a fairly wide range of smoke and underlying surface parameters. Overall, results from this work suggest that airborne CAR measurements can bring some substantial improvements in the accuracy of satellite-based radiative flux 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]

Atmospheric correction of satellite optical imagery over inland waters is a key remaining challenge in aquatic remote sensing. This is due to numerous confounding factors such as the complexity of water optical properties, the surface glint, the heterogeneous nature of atmospheric aerosols, and the proximity of bright land surfaces. This combination of factors makes it difficult to retrieve accurate information about the system observed. Moreover, the impact of radiance coming from adjacent land (adjacency effects) in complex geometries further adds to this challenge, especially for small lakes. In this study, ten atmospheric correction algorithms were evaluated for high-resolution multispectral imagery of Landsat-8 Operational Land Imager and Sentinel-2 MultiSpectral Instrument using in situ optical measurements from ~300 lakes across Canada. The results of the validation show that the performance of the algorithms varies by spectral band and evaluation metrics. The dark spectrum fitting algorithm had the best performance in terms of similarity angle (spectral shape), while the neural network-based models showed the lowest errors and bias per band. However, none of the tested atmospheric correction algorithms meet a 30% retrieval accuracy target across all the visible bands, likely due to uncorrected adjacency effects. To quantify this process, three-dimensional radiative transfer simulations were performed and compared to satellite observations. These simulations show that up to 60% of the top of atmosphere reflectance in the near-infrared bands over the lake was from the adjacent lands covered with green vegetation. The significance of these adjacency effects on atmospheric correction has been analyzed qualitatively, and potential efforts to improve the atmospheric correction algorithms are discussed. [ABSTRACT FROM AUTHOR]

Machine learning (ML) algorithms have emerged as competent tools for identifying areas that are susceptible to flooding. The primary variables considered in most of these works include terrain models, lithology, river networks and land use. While several recent studies include average annual rainfall and/or temperature, other meteorological information such as snow accumulation and short-term intense rain events that may influence the hydrology of the area under investigation have not been considered. Notably, in Canada, most inland flooding occurs during the freshet, due to the melting of an accumulated snowpack coupled with heavy rainfall. Therefore, in this study the impact of several climate variables along with various hydro-geomorphological (HG) variables were tested to determine the impact of their inclusion. Three tests were run: only HG variables, the addition of annual average temperature and precipitation (HG-PT), and the inclusion of six other meteorological datasets (HG-8M) on five study areas across Canada. In HG-PT, both precipitation and temperature were selected as important in every study area, while in HG-8M a minimum of three meteorological datasets were considered important in each study area. Notably, as the meteorological variables were added, many of the initial HG variables were dropped from the selection set. The accuracy, F1, true skill and Area Under the Curve (AUC) were marginally improved when the meteorological data was added to the a parallel random forest algorithm (parRF). When the model is applied to new data, the estimated accuracy of the prediction is higher in HG-8M, indicating that inclusion of relevant, local meteorological datasets improves the result. [ABSTRACT FROM AUTHOR]

The lowest elevation of spring snow ("snowline") is an important factor influencing recruitment and survival of wildlife in alpine areas. In this study, we assessed the spatial and temporal variability of alpine spring snowline across major Dall sheep mountain areas in Alaska and northwestern Canada. We used a daily MODIS snow fraction product to estimate the last day of 2000-2016 spring snow for each 500-m pixel within 28 mountain areas. We then developed annual (2000-2016) regression models predicting the elevation of alpine snowline during mid-May for each mountain area. MODIS-based regression estimates were compared with estimates derived using a Normalized Difference Snow Index from Landsat-8 Operational Land Imager (OLI) surface reflectance data. We also used 2000-2009 decadal climate grids to estimate total winter precipitation and mean May temperature for each of the 28 mountain areas. Based on our MODIS regression models, the 2000-2016 mean 15 May snowline elevation ranged from 339 m in the cold arctic class to 1145 m in the interior mountain class. Spring snowline estimates from MODIS and Landsat OLI were similar, with a mean absolute error of 106 m. Spring snowline elevation was significantly related to mean May temperature and total winter precipitation. The late spring of 2013 may have impacted some sheep populations, especially in the cold arctic mountain areas which were snow-covered in mid-May, while some interior mountain areas had mid-May snowlines exceeding 1000 m elevation. We found this regional (>500,000 km²) remote sensing application useful for determining the inter-annual and regional variability of spring alpine snowline among 28 mountain areas. [ABSTRACT FROM AUTHOR]

The Terrain Observation by Progressive Scans (TOPS) acquisition mode of the Sentinel-1 mission provides a wide coverage per acquisition with resolutions of 5 m in range and 20 m in azimuth, which makes this acquisition mode attractive for glacier velocity monitoring. Here, we retrieve surface velocities from the southern Ellesmere Island ice caps (Canadian Arctic) using both offset tracking and Differential Interferometric Synthetic Aperture Radar (D-InSAR) techniques and combining ascending and descending passes. We optimise the offset tracking technique by omitting the azimuth offsets. By doing so, we are able to improve the final resolution of the velocity product, as Sentinel-1 shows a lower resolution in the azimuth direction. Simultaneously, we avoid the undesired ionospheric effect manifested in the data as azimuth streaks. The D-InSAR technique shows its merits when applied to slow-moving areas, while offset tracking is more suitable for fast-moving areas. This research shows that the methods used here are complementary and the use of both to determine glacier velocities is better than only using one or the other. We observe glacier surface velocities of up to 1200 m year-1 for the fastest tidewater glaciers. The land-terminating glaciers show typical velocities between 12 and 33 m year-1, though with peaks up to 150 m year-1 in narrowing zones of the confining valleys. [ABSTRACT FROM AUTHOR]

Wildfires are burning increasingly extensive areas of forest in Canada, reducing their capacity as carbon sinks. Here we compare the answers that four independent land cover datasets, produced from different satellite images (SPOT, Landsat, and MODIS), provide for the question of what burned in Canada in recent years. We harmonized the different datasets into a common, simpler legend consisting of three classes of forest (needle-leaf, broadleaf, and mixed) plus non-forest, and resampled them to a common pixel size (250 m). Then we used annual maps of burned area to count, for each map and year from 2011 to 2014, the number of burned pixels of each class, and we tabulated them by terrestrial ecozone and Canada-wide. While all four maps agree that needle-leaf forest is the most frequently burned class in Canada, there is great disparity in the results from each map regarding the proportion of burned area that each class represents. Proportions reported by one map can be more than double those reported by another map, and more than four times at the ecozone level. We discuss the various factors that can explain the observed discrepancies and conclude that none of the maps provides a sufficiently accurate answer for applications such as carbon accounting. There is a need for better information in areas lacking forest inventory, especially in the vast unmanaged forest of Canada. [ABSTRACT FROM AUTHOR]

The winter regime of river-ice covers in high northern latitude regions is often a determining factor in the management of water resources, conservation of aquatic ecosystems and preservation of traditional and cultural lifestyles of local peoples. As ground-based monitoring of river-ice regimes in high northern latitudes is expensive and restricted to a few locations due to limited accessibility to most places along rivers from shorelines, remote sensing techniques are a suitable approach for monitoring. This study developed a RADARSAT-2 based method to monitor the spatio-temporal variation of ice covers, as well as ice types during the freeze-up period, along the main channel of the Slave River Delta in the Northwest Territories of Canada. The spatio-temporal variation of ice covers along the river was analyzed using the backscatter-based coefficient of variation (CV) in the 2013-2014 and 2014-2015 winters. As a consequence of weather and flow conditions, the ice cover in the 2013-2014 winter had the higher variation than the 2014-2015 winter, particularly in the potential areas of flooded/cracked ice covers. The river sections near active channels (e.g., Middle Channel and Nagle Channel), Big Eddy, and Great Slave Lake also yielded higher intra-annual variation of ice cover characteristics during the winters. With the inclusion of backscatter and texture analysis from RADARSAT-2 data, four water and ice cover classes consisting of open water, thermal ice, juxtaposed ice, and consolidated ice, were discriminated in the images acquired between November and March in both the studied winters. In addition to river geomorphology and climatic conditions such as river width, sinuosity or air temperature, the fluctuation of water flows during the winter has a significant impact on the variation of ice cover as well as the formation of different ice types in the Slave River. The RADARSAT-2 based monitoring algorithm can also be applied to other river systems in high latitude ecosystems to annually monitor their river-ice variation and formation during the freeze-up and ice cover progression period. [ABSTRACT FROM AUTHOR]

Satellite forest inventories are the only feasible way to map Canada's vast, remote forest regions, such as those in the Northwest Territories (NWT). A method used to create such inventories is the k-nearest neighbour (k-NN) algorithm, which spatially extends information from forest inventory (FI) plots to the entire forest land base using wall-to-wall features typically derived from Landsat data. However, the benefits of integrating L-band synthetic aperture radar (SAR) data, strongly correlated to forest biomass, have not been assessed for Canadian northern boreal forests. Here we describe an optimized multivariate k-NN implementation of a 151,700 km2 area in southern NWT that included ca. 2007 Landsat and dual-polarized Phased Array type L-band SAR (PALSAR) data on board the Advanced Land Observing Satellite (ALOS). Five forest attributes were mapped at 30 m cells: stand height, crown closure, stand/total volume and aboveground biomass (AGB). We assessed accuracy gains compared to Landsat-based maps. To circumvent the scarcity of FI plots, we used 3600 footprints from the Geoscience Laser Altimeter System (GLAS) as surrogate FI plots, where forest attributes were estimated using Light Detection and Ranging (LiDAR) metrics as predictors. After optimization, k-NN predicted forest attribute values for each pixel as the average of the 4 nearest (k = 4) surrogate FI plots within the Euclidian space of 9 best features (selected among 6 PALSAR, 10 Landsat, and 6 environmental features). Accuracy comparisons were based on 31 National Forest Inventory ground plots and over 1 million airborne LiDAR plots. Maps that included PALSAR HV backscatter resulted in forest attribute predictions with higher goodness of fit (adj. R2), lower percent mean error (ME%), and percent root mean square error (RMSE%), and lower underestimation for larger attribute values. Predictions were most accurate for conifer stand height (RMSE% = 32.1%, adj. R2 = 0.58) and AGB (RMSE% = 47.8%, adj. R2 = 0.74), which is much more abundant in the area than mixedwood or broadleaf. Our study demonstrates that optimizing k-NN parameters and feature space, including PALSAR, Landsat, and environmental variables, is a viable approach for inventory mapping of the northern boreal forest regions of Canada. [ABSTRACT FROM AUTHOR]

Spectral vegetation indices have been the primary resources for characterizing grassland vegetation based on remotely sensed data. However, the use of spectral indices for vegetation characterization in grasslands has been challenged by the confounding effects from external factors, such as soil properties, dead materials, and shadowing of vegetation canopies. Dead materials refer to the dead component of vegetation, including fallen litter and standing dead grasses accumulated from previous years. The abundant dead materials have been presenting challenges to accurately estimate green vegetation using spectral vegetation indices (VIs) derived from remote sensing data in mixed grasslands. Therefore, a close investigation of the relationship between VIs and dead materials is needed. The identified relationships could provide better insight into not only using remote sensing data for quantitative estimation of dead materials, but also the improvement of green vegetation estimation in the mixed grassland that has a high proportion of dead materials. In this article, the spectral reflectance of dead materials and green vegetation mixtures and dead material cover were measured in mixed grasslands located in Grassland National Park (GNP), Saskatchewan, Canada. Nine VIs were derived from the measured spectral reflectance. The relationship between dead material cover and VIs was quantified using the regression model and sensitivity analysis. Results indicated that the relationship between dead material cover and VIs is a function of the amount of dead material cover. Weak positive relationship was found between VIs and dead materials where the cover was less than 50%, and a significant high negative relationship was evident when cover was greater than 50%. When the combined exponential and linear model was applied to fit the negative relationships, more than 90% variation in dead material cover could be explained by VIs. Sensitivity analysis was further applied to the developed models, indicating that sensitivities of all VIs were significant over the entire range of dead material cover except for the triangular vegetation index (TVI), which has insignificant sensitivity when dead material cover was greater than 94%. Among all VIs, the weighted difference vegetation index (WDVI) had the highest sensitivity to changes in dead material cover higher than 50%. The results from this study indicated that vegetation indices based on combination of reflectance in red and NIR bands can be used to estimate dead material cover that is greater than 50%. [ABSTRACT FROM AUTHOR]

Simple Summary: This study aims to quantify the effects of smoke originating from boreal biomass-burning fires on solar radiation propagating through the atmosphere. The wildfires that took place in summer 2017 along Greenland's west coast and northern Canada produced a considerable amount of particles that were transported north and northeast, respectively, and were detected at the Thule High Arctic Atmospheric Observatory (THAAO; 76.53 ∘ N, 68.74 ∘ W). Solar radiation measurements carried out at THAAO, satellite data, and modeled radiations allowed the estimation of surface cooling and the warming aloft at two atmospheric layers (at altitudes of approximately 5 and 11 km asl) due to the presence of the biomass-burning aerosol particles. Boreal fires have increased during the last years and are projected to become more intense and frequent as a consequence of climate change. Wildfires produce a wide range of effects on the Arctic climate and ecosystem, and understanding these effects is crucial for predicting the future evolution of the Arctic region. This study focuses on the impact of the long-range transport of biomass-burning aerosol into the atmosphere and the corresponding radiative perturbation in the shortwave frequency range. As a case study, we investigate an intense biomass-burning (BB) event which took place in summer 2017 in Canada and subsequent northeastward transport of gases and particles in the plume leading to exceptionally high values (0.86) of Aerosol Optical Depth (AOD) at 500 nm measured in northwestern Greenland on 21 August 2017. This work characterizes the BB plume measured at the Thule High Arctic Atmospheric Observatory (THAAO; 76.53 ∘ N, 68.74 ∘ W) in August 2017 by assessing the associated shortwave aerosol direct radiative impact over the THAAO and extending this evaluation over the broader region ( 60 ∘ N– 80 ∘ N, 110 ∘ W– 0 ∘ E). The radiative transfer simulations with MODTRAN6.0 estimated an aerosol heating rate of up to 0.5 K/day in the upper aerosol layer (8–12 km). The direct aerosol radiative effect (ARE) vertical profile shows a maximum negative value of −45.4 Wm − 2 for a 78 ∘ solar zenith angle above THAAO at 3 km altitude. A cumulative surface ARE of −127.5 TW is estimated to have occurred on 21 August 2017 over a portion (∼ 3.1 × 10 6 km 2 ) of the considered domain ( 60 ∘ N– 80 ∘ N, 110 ∘ W– 0 ∘ E). ARE regional mean daily values over the same portion of the domain vary between −65 and −25 Wm − 2 . Although this is a limited temporal event, this effect can have significant influence on the Arctic radiative budget, especially in the anticipated scenario of increasing wildfires. [ABSTRACT FROM AUTHOR]

This research aims at modeling the microwave backscatter of corn fields by coupling an incoherent, interaction-based scattering model with a semi-empirical bulk vegetation dielectric model. The scattering model is fitted to co-polarized phase difference measurements over several corn fields imaged with fully polarimetric synthetic aperture radar (SAR) images with incidence angles ranging from 20° to 60°. The dataset comprised two field campaigns, one over Canada with the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR, 1.258 GHz) and the other one over Argentina with Advanced Land Observing Satellite 2 (ALOS-2) Phased Array type L-band Synthetic Aperture Radar (PALSAR-2) (ALOS-2/PALSAR-2, 1.236 GHz), totaling 60 data measurements over 28 grown corn fields at peak biomass with stalk gravimetric moisture larger than 0.8 g/g. Co-polarized phase differences were computed using a maximum likelihood estimation technique from each field's measured speckled sample histograms. After minimizing the difference between the model and data measurements for varying incidence angles by a nonlinear least-squares fitting, well agreement was found with a root mean squared error of 24.3° for co-polarized phase difference measurements in the range of −170.3° to −19.13°. Model parameterization by stalk gravimetric moisture instead of its complex dielectric constant is also addressed. Further validation was undertaken for the UAVSAR dataset on earlier corn stages, where overall sensitivity to stalk height, stalk gravimetric moisture, and stalk area density agreed with ground data, with the sensitivity to stalk diameter being the weakest. This study provides a new perspective on the use of co-polarized phase differences in retrieving corn stalk features through inverse modeling techniques from space. [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]

The population of loggerhead shrike (Lanius ludovicianus excubutirudes) in Grassland National Park of Canada (GNPC) has undergone a severe decline due to habitat loss and limitation. Shrike habitat availability is highly impacted by the biophysical characteristics of grassland landscapes. This study was conducted in the west block of GNPC. The overall purpose was to extract important biophysical and topographical variables from both SPOT satellite imagery and in situ measurements. Statistical analysis including Analysis of Variance (ANOVA), measuring Coefficient Variation (CV), and regression analysis were applied to these variables obtained from both imagery and in situ measurement. Vegetation spatial variation and heterogeneity among active, inactive and control nesting sites at 20 m × 20 m, 60 m × 60 m and 100 m × 100 m scales were investigated. Results indicated that shrikes prefer to nest in open areas with scattered shrubs, particularly thick or thorny species of smaller size, to discourage mammalian predators. The most important topographical characteristic is that active sites are located far away from roads at higher elevation. Vegetation index was identified as a good indicator of vegetation characteristics for shrike habitats due to its significant relation to most relevant biophysical factors. Spatial variation analysis showed that at all spatial scales, active sites have the lowest vegetation abundance and the highest heterogeneity among the three types of nesting sites. For all shrike habitat types, vegetation abundance decreases with increasing spatial scales while habitat heterogeneity increases with increasing spatial scales. This research also indicated that suitable shrike habitat for GNPC can be mapped using a logistical model with ATSAVI and dead material in shrub canopy as the independent variables. [ABSTRACT FROM AUTHOR]

In order to better constrain the utility of multispectral datasets in the characterization of surface materials, Landsat Thematic Mapper (TM) data were evaluated in the discrimination of geological classes in the Cape Smith Belt of Quebec, a greenstone belt that hosts Early Proterozoic units including those of the Purtuniq ophiolite. Ground-based measurements collected for the study area highlight the importance of chemical alteration in controlling the reflectance properties of key geological classes. The spatial distribution of exposed lithologies in the study area was determined through (1) image classification using a feedforward backpropagation neural network classifier; and (2) generation of fraction images for spectral end members using a linear unmixing algorithm and ground reflectance data. Despite some shortcomings, the database of surface cover generated by the neural network classifier is a useful representation of the spatial distribution of exposed geological materials in the study area, with an overall agreement with ground truth of 87.7%. In contrast, the fraction images generated through unmixing are poor representations of ground truth for several key lithological classes. These results underscore both the considerable utility and marked limitations of Landsat TM data in the mapping of igneous and metamorphic lithologies. [ABSTRACT FROM AUTHOR]

Since 2002, global snow water equivalent (SWE) estimates have been generated using Advanced Microwave Scanning Radiometer (AMSR-E)/Aqua data. Accurate estimates of SWE are important to improve monitoring and managing of water resources in specific regions. SWE and snow map product accuracy are functions of topography and of land cover type because landscape characteristics have a strong influence on redistribution and physical properties of snow cover, and influence the microwave properties of the surface. Here we evaluate the AMSR-E SWE and derived snow map products in the Mackenzie River Basin (MRB), Canada, which is characterized by complex topography and varying land cover types from tundra to boreal forest. We compare in situ snow depth observations and Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover maps from January 2003 to December 2007 with passive microwave remotely sensed SWE from AMSR-E and derived snow cover maps. In the MRB the mean absolute error ranges from 12 mm in the early winter season to 50 mm in the late winter season and overestimations of snow cover maps based on a 1 mm threshold of AMSR-E SWE varies from 4% to 8%. The optimal threshold for AMSR-E SWE to classify the pixels as snow ranges from 6 mm to 9 mm. The overall accuracy of new snow cover maps from AMSR-E varies from 91% to 94% in different sub-basins in the MRB. [ABSTRACT FROM AUTHOR]

The Sporadic-E (Es) layer is an often-observed phenomenon at high latitudes; however, our understanding of the polar cap Es layer is severely limited due to the scarce number of measurements. Here, the first comprehensive study of the polar cap Es layer associated with Global Positioning System (GPS) Total Electron Content (TEC) variations and scintillations is presented with multiple measurements at Resolute, Canada (Canadian Advanced Digital Ionosonde (CADI), Northward-looking face of Resolute Incoherent-Scatter Radar (RISR-N), and GPS receiver). According to the joint observations, the polar cap Es layer is a thin patch structure with variously high electron density, which gradually develops into the lower E region (~100 km) and horizontally extends >200 km. Moreover, the TEC variations produced by the polar cap Es layer are pulse-like enhancements with a general amplitude of ~0.5 TECu and are followed by smaller but rapid TEC perturbations. Furthermore, the possible scintillation effects likely associated with the polar cap Es layer are also discussed. As a consequence, the results widely expand our understanding on the polar cap Es layer, in particular on TEC variations. [ABSTRACT FROM AUTHOR]

Four burned area tools were implemented in Google Earth Engine (GEE), to obtain regular processes related to burned area (BA) mapping, using medium spatial resolution sensors (Landsat and Sentinel-2). The four tools are (i) the BA Cartography tool for supervised burned area over the user-selected extent and period, (ii) two tools implementing a BA stratified random sampling to select the scenes and dates for validation, and (iii) the BA Reference Perimeter tool to obtain highly accurate BA maps that focus on validating coarser BA products. Burned Area Mapping Tools (BAMTs) go beyond the previously implemented Burned Area Mapping Software (BAMS) because of GEE parallel processing capabilities and preloaded geospatial datasets. BAMT also allows temporal image composites to be exploited in order to obtain BA maps over a larger extent and longer temporal periods. The tools consist of four scripts executable from the GEE Code Editor. The tools' performance was discussed in two case studies: in the 2019/2020 fire season in Southeast Australia, where the BA cartography detected more than 50,000 km2, using Landsat data with commission and omission errors below 12% when compared to Sentinel-2 imagery; and in the 2018 summer wildfires in Canada, where it was found that around 16,000 km2 had burned. [ABSTRACT FROM AUTHOR]

Satellite altimetry is an important contributor for measuring the water level of continental water bodies. The technique has been applied for almost three decades. In this period the data quality has increased and the applications have evolved from the study of a few large lakes and rivers, to near global applications at various scales. Products from current satellite altimetry missions should be validated to continuously improve the measurements. Sentinel-3A has been operating since 2016 and is the first mission operating in synthetic aperture radar (SAR) mode globally. Here we evaluate its performance in capturing lake level variations based on a physical and an empirical retracker provided in the official level 2 product. The validation is performed for more than 100 lakes in the United States and Canada where the altimetry based water levels are compared with in situ data. As validation measures we consider the root mean squared error, the Pearson correlation, and the percentage of outliers. For the US sites the median of the RMSE value is 25 cm and 19 cm and the median of the Pearson correlations are 0.86 and 0.93 for the physical and empirical retracker, respectively. The percentage of outliers (median) is 11% for both retrackers. The validations measures are slightly poorer for the Canadian sites; the median RMSE is approximately 5 cm larger, the Pearson correlation 0.1 lower, and the percentage of outliers 5% larger. The poorer performance for the Canadian sites is mainly related to the presence of lake ice in the winter period where the surface elevations are not able to map the surface correctly. The validation measures improve considerably when evaluated for summer data only. For both areas we show that the reconstruction of the water level variations based on the empirical retracker is significantly better compared to that of the physical retracker in terms of the RMSE and the Pearson correlation. [ABSTRACT FROM AUTHOR]

Wildfires are a concerning issue in Canada due to their immediate impact on people's lives, local economy, climate, and environment. Studies have shown that the number of wildfires and affected areas in Canada has increased during recent decades and is a result of a warming and drying climate. Therefore, identifying potential wildfire risk areas is increasingly an important aspect of wildfire management. The purpose of this study is to investigate if remotely sensed soil moisture products from the Soil Moisture and Ocean Salinity (SMOS) satellite can be used to identify potential wildfire risk areas for better wildfire management. We used the National Fire Database (NFDB) fire points and polygons to group the wildfires according to ecozone classifications, as well as to analyze the SMOS soil moisture data over the wildfire areas, between 2010–2017, across fourteen ecozones in Canada. Timeseries of 3-day, 5-day, and 7-day soil moisture anomalies prior to the onset of each wildfire occurrence were examined over the ecozones individually. Overall, the results suggest, despite the coarse-resolution, SMOS soil moisture products are potentially useful in identifying soil moisture anomalies where wildfire hot-spots may occur. [ABSTRACT FROM AUTHOR]

When severe flooding occurs in Canada, the Emergency Geomatics Service (EGS) is tasked with creating and disseminating maps that depict flood extents in near real time. EGS flood mapping methods were created with efficiency and robustness in mind, to allow maps to be published quickly, and therefore have the potential to generate high-repeat water products that can enhance frequent wetland monitoring. The predominant imagery currently used is synthetic aperture radar (SAR) from RADARSAT-2 (R2). With the commissioning phase of the RADARSAT Constellation Mission (RCM) complete, the EGS is adapting its methods for use with this new source of SAR data. The introduction of RCM's circular-transmit linear-receive (CTLR) beam mode provides the option to exploit compact polarimetric (CP) information not previously available with R2. The aim of this study was to determine the most effective CP parameters for use in mapping open water and flooded vegetation, using current EGS methodologies, and compare these products to those created by using R2 data. Nineteen quad-polarization R2 scenes selected from three regions containing wetlands prone to springtime flooding were used to create reference flood maps, using existing EGS tools. These scenes were then used to simulate 22 RCM CP parameters at different noise floors and spatial resolutions representative of the three RCM beam modes. Using multiple criteria, CP parameters were ranked in order of importance and entered into a stepwise classification procedure, for evaluation against reference R2 products. The top four CP parameters —m-chi-volume or m-delta-volume, RR intensity, Shannon Entropy intensity (SEi), and RV intensity—achieved a maximum agreement with baseline R2 products of upward of 98% across all 19 scenes and three beam modes. Separability analyses between flooded vegetation and other land-cover classes identified four candidate CP parameters—RH intensity, RR intensity, SEi, and the first Stokes parameter (SV0)—suitable for flooded-vegetation-region growing. Flooded-vegetation-region-growing CP thresholds were found to be dependent on incidence angle for each of these four parameters. After region growing using each of the four candidate CP parameters, RH intensity was deemed best to map flooded vegetation, based on our evaluations. The results of the study suggest a set of suitable CP parameters to generate flood maps from RCM data, using current EGS methodologies that must be validated further as real RCM data become available. [ABSTRACT FROM AUTHOR]

Lake ice is a dominant component of Canada's landscape and can act as an indicator for how freshwater aquatic ecosystems are changing with warming climates. While lake ice monitoring through government networks has decreased in the last three decades, the increased availability of remote sensing images can help to provide consistent spatial and temporal coverage for areas with annual ice cover. Synthetic aperture radar (SAR) data are commonly used for lake ice monitoring, due to the acquisition of images in any condition (time of day or weather). Using Sentinel-1 A/B images, a high-density time series of SAR images was developed for Lake Hazen in Nunavut, Canada, from 2015–2018. These images were used to test two different methods of monitoring lake ice phenology: one method using the first difference between SAR images and another that applies the Otsu segmentation method. Ice phenology dates determined from the two methods were compared with visual interpretation of the Sentinel-1 images. Mean errors for the pixel comparison of the first difference method ranged 3–10 days for ice-on and ice-off, while average error values for the Otsu method ranged 2–10 days. Mean errors for comparisons of different sections of the lake ranged 0–15 days for the first difference method and 2–17 days for the Otsu method. This research demonstrates the value of temporally consistent image acquisition for improving the accuracy of lake ice monitoring. [ABSTRACT FROM AUTHOR]

Understanding the ecological effects of human activities on an ecosystem is integral to the implementation of conservation management plans. The plasticity of plant functional traits presents an opportunity to examine the capacity for intraspecific functional trait variations to be indicators of anthropogenic landscape modifications. The presence of intraspecific trait variation would indicate that plants of a single species could to be used to evaluate and map functional diversity, a common metric used to measure biodiversity. This study uses leaf spectroscopy, light detection and ranging (LiDAR) and partial least squares regression (PLSR) to examine the intraspecific variation of functional traits in a population of 40 Quercus garryana experiencing varying levels of anthropogenic influence at the site level (<0.3 km2) in Duncan, B.C., Canada. These individuals vary in their spatial relationship to roads, agricultural land use change and an encroaching Coastal Douglas-fir forest. A total of 14 functional traits were estimated using pre-determined PLSR coefficients from a multi-species dataset. LiDAR data for each tree and were organized into functional categories based on their influence of plant lifeform, leaf growth or leaf structure. Principal components analysis was performed on each functional category to determine the relative influence of each trait. Results show that leaf growth and lifeform functional trait categories express significant variation in relation to three anthropogenic landscape modifications, while traits associated to leaf structure only varied between land use types (p = 0.05). Diameter at breast height (DBH), mass-based chlorophyll and leaf mass per area (LMA) showed the strongest variation across treatments. These findings support the hypothesis that trait variation exists in small populations of the same species and illustrate that spectroscopy can be used to indirectly sense land use via the leaf functional traits of a single tree species. [ABSTRACT FROM AUTHOR]

Land cover change (LCC) is typically characterized by infrequent changes over space and time. Data-driven methods such as deep learning (DL) approaches have proven effective in many domains for predictive and classification tasks. When applied to geospatial data, sequential DL methods such as long short-term memory (LSTM) have yielded promising results in remote sensing and GIScience studies. However, the characteristics of geospatial datasets selected for use with these methods have demonstrated important implications on method performance. The number of data layers available, the rate of LCC, and inherent errors resulting from classification procedures are expected to influence model performance. Yet, it is unknown how these can affect compatibility with the LSTM method. As such, the main objective of this study is to explore the capacity of LSTM to forecast patterns that have emerged from LCC dynamics given varying temporal resolutions, persistent land cover classes, and auxiliary data layers pertaining to classification confidence. Stacked LSTM modeling approaches are applied to 17-year MODIS land cover datasets focused on the province of British Columbia, Canada. This geospatial data is reclassified to four major land cover (LC) classes during pre-processing procedures. The evaluation considers the dataset at variable temporal resolutions to demonstrate the significance of geospatial data characteristics on LSTM method performance in several scenarios. Results indicate that LSTM can be utilized for forecasting LCC patterns when there are few limitations on temporal intervals of the datasets provided. Likewise, this study demonstrates improved performance measures when there are classes that do not change. Furthermore, providing classification confidence data as ancillary input also demonstrated improved results when the number of timesteps or temporal resolution is limited. This study contributes to future applications of DL and LSTM methods for forecasting LCC. [ABSTRACT FROM AUTHOR]

This research aimed at exploiting the joint use of machine learning and polarimetry for improving the retrieval of surface soil moisture content (SMC) from synthetic aperture radar (SAR) acquisitions at C-band. The study was conducted on two agricultural areas in Canada, for which a series of RADARSAT-2 (RS2) images were available along with direct measurements of SMC from in situ stations. The analysis confirmed the sensitivity of RS2 backscattering (σ°) to SMC. The comparison of SMC with the compact polarimetry (CP) parameters, computed from the RS2 acquisitions by the CP data simulator, pointed out that some CP parameters had a sensitivity to SMC equal or better than σ°, with correlation coefficients up to R ≃ 0.4. Based on these results, the potential of machine learning (ML) for SMC retrieval was exploited by implementing and testing on the available data an artificial neural network (ANN) algorithm. The algorithm was implemented using several combinations of σ° and CP parameters. Validation results of the algorithm with in situ observations confirmed the promising capabilities of the ML techniques for SMC monitoring. Furthermore, results pointed out the potential of CP in improving the SMC retrieval accuracy, especially when used in combination with linearly polarized σ°. Depending on the considered input combination, the ANN algorithm was able to estimate SMC with Root Mean Square Error (RMSE) between 3% and 7% of SMC and R between 0.7 and 0.9. [ABSTRACT FROM AUTHOR]

Observations of vegetation phenology at regional-to-global scales provide important information regarding seasonal variation in the fluxes of energy, carbon, and water between the biosphere and the atmosphere. Numerous algorithms have been developed to estimate phenological transition dates using time series of remotely sensed spectral vegetation indices. A key challenge, however, is that different algorithms provide inconsistent results. This study provides a comprehensive comparison of start of season (SOS) and end of season (EOS) phenological transition dates estimated from 500 m MODIS data based on two widely used sources of such data: the TIMESAT program and the MODIS Global Land Cover Dynamics (MLCD) product. Specifically, we evaluate the impact of land cover class, criteria used to identify SOS and EOS, and fitting algorithm (local versus global) on the transition dates estimated from time series of MODIS enhanced vegetation index (EVI). Satellite-derived transition dates from each source are compared against each other and against SOS and EOS dates estimated from PhenoCams distributed across the Northeastern United States and Canada. Our results show that TIMESAT and MLCD SOS transition dates are generally highly correlated (r = 0.51-0.97), except in Central Canada where correlation coefficients are as low as 0.25. Relative to SOS, EOS comparison shows lower agreement and higher magnitude of deviations. SOS and EOS dates are impacted by noise arising from snow and cloud contamination, and there is low agreement among results from TIMESAT, the MLCD product, and PhenoCams in vegetation types with low seasonal EVI amplitude or with irregular EVI time series. In deciduous forests, SOS dates from the MLCD product and TIMESAT agree closely with SOS dates from PhenoCams, with correlations as high as 0.76. Overall, our results suggest that TIMESAT is well-suited for local-to-regional scale studies because of its ability to tune algorithm parameters, which makes it more flexible than the MLCD product. At large spatial scales, where local tuning is not feasible, the MLCD product provides a readily available data set based on a globally consistent approach that provides SOS and EOS dates that are comparable to results from TIMESAT. [ABSTRACT FROM AUTHOR]

Information on crop seeding date is required in many applications such as crop management and yield forecasting. This study presents a novel method to estimate crop seeding date at the field level from time-series 250-m Moderate Resolution Imaging Spectroradiometer (MODIS) data and growing degree days (GDD; base 5 °C; °C-days). The start of growing season (SOS) was first derived from time-series EVI2 (two-band Enhanced Vegetation Index) calculated from a MODIS 8-day composite surface reflectance product (MOD09Q1; Collection 6). Based on GDD calculated from the Daymet gridded estimates of daily weather parameters, a simple model was developed to establish a linkage between the observed seeding date and the SOS. Calibration and validation of the model was conducted on three major crops, spring wheat, canola and oats in the Province of Manitoba, Canada. The estimated SOS had a strong linear correlation with the observed seeding date; with a deviation of a few days depending on the year. The seeding date of the three crops can be calculated from the SOS by adjusting the number of days needed to accumulate GDD (AGDD) for emergence. The overall root-mean-square-difference (RMSD) of the estimated seeding date was less than 10 days. Validation showed that the accuracy of the estimated seeding date was crop-type independent. The developed method is useful for estimating the historical crop seeding date from remote sensing data in Canada to support studies of the interactions among seeding date, crop management and crop yield under climate change. It is anticipated that this method can be adapted to other crops in other locations using the same or different satellite data. [ABSTRACT FROM AUTHOR]