Trends in soil temperature are important, but rarely reported, indicators of climate change. On the basis of the soil temperature data from 30 climate stations across Canada during 1958–2008, trends ...in soil temperatures at 5, 10, 20, 50, 100, and 150 cm depths were analyzed, together with atmospheric variables, such as air temperature, precipitation, and depth of snow on the ground, observed at the same locations. There was a significant positive trend with soil temperatures in spring and summer means, but not for the winter and annual means. A positive trend with time in soil temperature was detected at about two‐thirds of the stations at all depths below 5 cm. A warming trend of 0.26–0.30°C/decade was consistently detected in spring (March–April–May) at all depths between 1958 and 2008. The warming trend in soil temperatures was associated with trends in air temperatures and snow cover depth over the same period. A significant decreasing trend in snow cover depth in winter and spring was associated with increasing air temperatures. The combined effects of the higher air temperature and reduced snow depth probably resulted in an enhanced increasing trend in spring soil temperatures, but no significant trends in winter soil temperatures. The thermal insulation by snow cover appeared to play an important role in the response of soil temperatures to climate change and must be accounted for in projecting future soil‐related impacts of climate change.
Accurate measurement of leaf chlorophyll concentration (LChl) in the field using a portable chlorophyll meter (PCM) is crucial to support methodology development for mapping the spatiotemporal ...variability of crop nitrogen status using remote sensing. Several PCMs have been developed to measure LChl instantaneously and non-destructively in the field, however, their readings are relative quantities that need to be converted into actual LChl values using conversion functions. The aim of this study was to investigate the relationship between actual LChl and PCM readings obtained by three PCMs: SPAD-502, CCM-200, and Dualex-4. Field experiments were conducted in 2016 on four crops: corn (Zea mays L.), soybean (Glycine max L. Merr.), spring wheat (Triticum aestivum L.), and canola (Brassica napus L.), at the Central Experimental Farm of Agriculture and Agri-Food Canada in Ottawa, Ontario, Canada. To evaluate the impact of other factors (leaf internal structure, leaf pigments other than chlorophyll, and the heterogeneity of LChl distribution) on the conversion function, a global sensitivity analysis was conducted using the PROSPECT-D model to simulate PCM readings under different conditions. Results showed that Dualex-4 had a better performance for actual LChl measurement than SPAD-502 and CCM-200, using a general conversion function for all four crops tested. For SPAD-502 and CCM-200, the error in the readings increases with increasing LChl. The sensitivity analysis reveals that deviations from the calibration functions are more induced by non-uniform LChl distribution than leaf architectures. The readings of Dualex-4 can have a better ability to restrict these influences than those of the other two PCMs.
Climate warming has significant impacts on permafrost, infrastructure and soil organic carbon at the northern high latitudes. These impacts are mainly driven by changes in soil temperature (TS). Snow ...insulation can cause significant differences between TS and air temperature (TA), and our understanding about this effect through space and time is currently limited. In this study, we compiled soil and air temperature observations (measured at about 0.2 m depth and 2 m height, respectively) at 588 sites from climate stations and boreholes across the northern high latitudes. Analysis of this circumpolar dataset demonstrates the large offset between mean TS and TA in the low arctic and northern boreal regions. The offset decreases both northward and southward due to changes in snow conditions. Correlation analysis shows that the coupling between annual TS and TA is weaker, and the response of annual TS to changes in TA is smaller in boreal regions than in the arctic and the northern temperate regions. Consequently, the inter-annual variation and the increasing trends of annual TS are smaller than that of TA in boreal regions. The systematic and significant differences in the relationship between TS and TA across the circumpolar north is important for understanding and assessing the impacts of climate change and for reconstruction of historical climate based on ground temperature profiles for the northern high latitudes.
Synthetic aperture radar (SAR) is more sensitive to the dielectric properties and structure of the targets and less affected by weather conditions than optical sensors, making it more capable of ...detecting changes induced by management practices in agricultural fields. In this study, the capability of C-band SAR data for detecting crop seeding and harvest events was explored. The study was conducted for the 2019 growing season in Temiskaming Shores, an agricultural area in Northern Ontario, Canada. Time-series SAR data acquired by Sentinel-1 constellation with the interferometric wide (IW) mode with dual polarizations in VV (vertical transmit and vertical receive) and VH (vertical transmit and horizontal receive) were obtained. interferometric SAR (InSAR) processing was conducted to derive coherence between each pair of SAR images acquired consecutively in time throughout the year. Crop seeding and harvest dates were determined by analyzing the time-series InSAR coherence and SAR backscattering. Variation of SAR backscattering coefficients, particularly the VH polarization, revealed seasonal crop growth patterns. The change in InSAR coherence can be linked to change of surface structure induced by seeding or harvest operations. Using a set of physically based rules, a simple algorithm was developed to determine crop seeding and harvest dates, with an accuracy of 85% (n = 67) for seeding-date identification and 56% (n = 77) for harvest-date identification. The extra challenge in harvest detection could be attributed to the impacts of weather conditions, such as rain and its effects on soil moisture and crop dielectric properties during the harvest season. Other factors such as post-harvest residue removal and field ploughing could also complicate the identification of harvest event. Overall, given its mechanism to acquire images with InSAR capability at 12-day revisiting cycle with a single satellite for most part of the Earth, the Sentinel-1 constellation provides a great data source for detecting crop field management activities through coherent or incoherent change detection techniques. It is anticipated that this method could perform even better at a shorter six-day revisiting cycle with both satellites for Sentinel-1. With the successful launch (2019) of the Canadian RADARSAT Constellation Mission (RCM) with its tri-satellite system and four polarizations, we are likely to see improved system reliability and monitoring efficiency.
To prevent potentially unsuitable activities during vegetation restoration, it is important to examine the impact of historical restoration activities on the target ecological system to inform future ...restoration policies. Taking the Loess Plateau of China as an example, a regionalization method and corresponding scheme were proposed to select suitable vegetation types (forested lands, woody grasslands/bushlands, grasslands, or xerophytic shrublands and semi-shrublands) for a given location using remote sensing technology in order to analyze the vegetation growth status before and after the largest ecological conservation project in the country: The Grain for Green Program (GTGP). To design the scheme, remote sensing data covering the periods before and after the implementation of the GTGP (the 1980s and 2001-2013) were collected, along with soil, meteorological, and topographic data. The net primary production (NPP) values for 2001-2013 were calculated using the Carnegie-Ames-Stanford Approach (CASA) model. Locations representing the native vegetation and the restored vegetation were first recognized using maps of vegetation cover. Then, for the restored vegetation area, the places suitable for planting the covered vegetation type were selected by comparing the NPP value of the corresponding vegetation type in the native vegetation area to the NPP value in the site under consideration. Third, half of these sites were uniformly selected based on their NPP value, and these areas and the native vegetation area were used as training regions. Based on weather, soil, and topographic data, a new regionalization scheme was designed using standardized Euclidean distances. Finally, data from the remainder of the Loess Plateau were used to validate the new regionalization scheme, which was also compared to an existing Chinese eco-geographical regionalization scheme. The results showed that the new regionalization scheme performed well, with an average potential classification accuracy of 81.81%. Compared with the eco-geographical regionalization scheme, the new scheme exhibited improved the consistency of vegetation dynamics, reflecting the potential to better guide vegetation restoration activities on the Loess Plateau.
Green leaf area index (LAI) is an important variable related to crop growth. Accurate and timely information on LAI is essential for developing suitable field management strategies to mitigate risk ...and boost yield. Several remote sensing (RS) based methods have been recently developed to estimate LAI at the regional scale. However, the performance of these methods tends to be affected by the quality of RS data, especially when time-series LAI are required. For crop LAI estimation, supplementary growth information from crop model is helpful to address this issue. In this study, we focus on the regional-scale LAI estimations of spring maize for the entire growth season. Using time-series multispectral RS data acquired by an unmanned aerial vehicle (UAV) and the World Food Studies (WOFOST) crop model, three methods were applied at different crop growth stages: empirical method using vegetation index (VI), data assimilation method and hybrid method. The VI-based method and assimilation method were used to generate time-series LAI estimations for the whole crop growth season. Then, a hybrid method specially for the late-stage LAI retrieval was developed by integrating WOFOST model and data assimilation. Using field-collected LAI data in Hongxing Farm in 2014, the performances of these three methods were evaluated. At the early stage, the VI-based method (R
= 0.63, RMSE = 0.16,
= 36) achieved higher accuracy than the assimilation method (R
= 0.54, RMSE = 0.52,
= 36), whereas at the mid stage, the assimilation method (R
= 0.63, RMSE = 0.46,
= 28) showed higher accuracy than the VI-based method (R
= 0.41, RMSE = 0.51,
= 28). At the late stage, the hybrid method yielded the highest accuracy (R
= 0.63, RMSE = 0.46,
= 29), compared with the VI-based method (R
= 0.19, RMSE = 0.43,
= 28) and the assimilation method (R
= 0.20, RMSE = 0.44,
= 29). Based on the results above, we considered a combination of the three methods, i.e., the VI-based method for the early stage, the assimilation method for the mid stage, and the hybrid method for the late stage, as an ideal strategy for spring-maize LAI estimation for the entire growth season of 2014 in Hongxing Farm, and the accuracy of the combined method over the whole growth season is higher than that of any single method.
This study investigated the estimation of grain yields of three major annual crops in Ontario (corn, soybean, and winter wheat) using MODIS reflectance data extracted with a general cropland mask and ...crop-specific masks. Time-series two-band enhanced vegetation index (EVI2) was derived from the 8 day composite 250 m MODIS reflectance data from 2003 to 2016. Using a general cropland mask, the strongest positive linear correlation between crop yields and EVI2 was observed at the end of July to early August, whereas a negative correlation was observed in spring. Using crop-specific masks, the time of the strongest positive linear correlation for winter wheat was found between mid-May and early June, corresponding to peak growth stages of the crop. EVI2 derived at peak growth stages of a crop provided good predictive capability for grain yield estimation, with considerable inter-annual variation. A multiple linear regression model was established for county-level yield estimation using EVI2 at peak growth stages and the year as independent variables. The model accounted for the spatiotemporal variability of grain yields of about 30% and 47% for winter wheat, 63% and 65% for corn, and 59% and 64% for soybean using the general cropland mask and crop-specific masks, respectively. A negative correlation during the spring indicated that vegetation index extracted using a general cropland mask should be used with caution in regions with mixed crops, as factors other than the growth conditions of the targeted crops may also be captured by remote sensing data.
We evaluated the utility of Terra/MODIS-derived crop metrics for yield estimation across the Canadian Prairies. This study was undertaken at the Census Agriculture Region (CAR) and the Rural ...Municipality (RM) of the province of Saskatchewan, in three prairie agro-climate zones. We compared MODIS-derived vegetation indices, gross primary productivity (GPP), and net primary productivity (NPP) to the known yields for barley, canola, and spring wheat. Multiple linear regressions were used to assess the relationships between the metrics and yield at the CAR and RM levels for the years 2000 to 2016. Models were evaluated using a leave-one-out cross validation (LOOCV) approach. Results showed that vegetation indices at crop peak growing stages were better predictors of yield than GPP or NPP, and EVI2 was better than NDVI. Using seasonal maximum EVI2, CAR-level crop yields can be estimated with a relative root-mean-square-error (RRMSE) of 14-20% and a Nash-Sutcliffe model efficiency coefficient (NSE) of 0.53-0.70, though the exact relationship varies by crop type and agro-climate zone. LOOCV showed the stability of the models across different years, although interannual fluctuations of estimation accuracy were observed. Assessments using RM-level yields showed slightly reduced accuracy, with NSE of 0.37-0.66, and RRMSE of 18-28%. The best performing models were used to map annual crop yields at the Soil Landscapes of Canada (SLC) polygon level. The results indicated that the models could perform well at both spatial scales, and thus, could be used to disaggregate coarse resolution crop yields to finer spatial resolutions using MODIS data.
Representative subsets of global climate models (GCMs) are often used in climate change impact studies to account for uncertainty in ensemble climate projections. However, the effectiveness of such ...subsets has seldom been assessed for the estimations of either the mean or the spread of the full ensembles. We assessed two different approaches that were employed to select 5 GCMs from a 20-member ensemble of GCMs from the CMIP5 ensemble for projecting canola and spring wheat yields across Canada under RCP 4.5 and 8.5 emission scenarios in the periods 2040-2069 and 2070-2099, based on crop simulation models. Averages and spreads of the simulated crop yields using the 5-GCM subsets selected by T&P and KKZ approaches were compared with the full 20-GCM ensemble. Our results showed that the 5-GCM subsets selected by the two approaches could produce full-ensemble means with a relative absolute error of 2.9-4.7% for canola and 1.5-2.2% for spring wheat, and covers 61.8-91.1% and 66.1-80.8% of the full-ensemble spread for canola and spring wheat, respectively. Our results also demonstrated that both approaches were very likely to outperform a subset of randomly selected 5 GCMs in terms of a smaller error and a larger range.
ABSTRACT
In addition to the uncertainty associated with crop models, climate scenarios are still a major source of uncertainty in projecting crop yield changes under climate change. Regional climate ...models (RCMs) are used as tools for dynamic downscaling of climate scenarios from global climate models (GCMs) to regional scales for climate change impact studies. It is known that running an RCM is more expensive and time consuming at a higher resolution than at a lower resolution. Therefore, it is interesting to investigate how resolutions of an RCM and statistical processing of RCM data may result in differences in projected crop yield changes. We used the decision support system for agrotechnology transfer (DSSAT)–CERES‐Wheat model to simulate yield changes of spring wheat at 13 locations across the Canadian Prairies, with climate scenarios from a Canadian RCM (CanRCM4) driven by a Canadian earth system model (CanESM2) with the forcing scenarios RCP4.5 and RCP8.5 at 25 and 50 km resolutions. Bias correction and a stochastic weather generator referred to as AAFC‐WG were used as statistical processing tools to develop future climate scenarios as input to the crop model. The results showed that when changes were averaged across the locations, whether 25‐ or 50‐km resolution CanRCM4 data were used, the projected yield changes were fairly consistent with those based on its driving GCM–CanESM2, especially if AAFC‐WG was used to develop future climate scenarios. The spatial patterns of the projected yield changes were also similar for the two resolutions of CanRCM4, although the magnitude of the projected changes had a relatively large range among the scenarios at some locations. These results indicated that using future climate scenarios based on climate change simulations by GCMs might be sufficient for projecting regional crop yield changes on the Canadian Prairies.