Dynamics and controls of urban heat sink and island phenomena in a desert city: Development of a local climate zone scheme using remotely-sensed inputs

This study aims to determine the dynamics and controls of Surface Urban Heat Sinks (SUHS) and Surface Urban Heat Islands (SUHI) in desert cities, using Dubai as a case study. A Local Climate Zone (LCZ) schema was developed to subdivide the city into different zones based on similarities in land cover and urban geometry. Proximity to the Gulf Coast was also determined for each LCZ. The LCZs were then used to sample seasonal and daily imagery from the MODIS thermal sensor to determine Land Surface Temperature (LST) variations relative to desert sand. Canonical correlation techniques were then applied to determine which factors explained the variability between urban and desert LST.Our results indicate that the daytime SUHS effect is greatest during the summer months (typically ∼3.0 °C) with the strongest cooling effects in open high-rise zones of the city. In contrast, the night-time SUHI effect is greatest during the winter months (typically ∼3.5 °C) with the strongest warming effects in compact mid-rise zones of the city. Proximity to the Arabian Gulf had the largest influence on both SUHS and SUHI phenomena, promoting daytime cooling in the summer months and night-time warming in the winter months. However, other parameters associated with the urban environment such as building height had an influence on daytime cooling, with larger buildings promoting shade and variations in airflow. Likewise, other parameters such as sky view factor contributed to night-time warming, with higher temperatures associated with limited views of the sky.     

Remotely sensed surface temperature variation of an inland saline lake over the central Qinghai–Tibet Plateau

Research on surface water temperature (SWT) variations in large lakes over the Qinghai–Tibet Plateau (QTP) has been limited by lack of in situ measurements. By taking advantage of the increased availability of remotely sensed observations, this study investigated SWT variation of Siling Co in central QTP by processing complete MODIS Land surface temperature (LST) images over the lake covering from 2001 to 2013. The temporal (diurnal, intra-annul and inter-annul) variations of Siling Co SWT as well as the spatial patterns were analyzed. The results show that on average from late December to mid-April the lake is in a mixing state of water and ice and drastic diurnal temperature differences occur, especially along the shallow shoreline areas. The extent of spatial variations in monthly SWT ranges from 1.25 °C to 3.5 °C, and particularly large at nighttime and in winter months. The spatial patterns of annual average SWT were likely impacted by the cooling effect of river inflow from the west and east side of the lake. The annual cycle of spatial pattern of SWT is characterized by seasonal reversions between the shallow littoral regions and deep parts due to different heat capacity. Compared to the deep regions, the littoral shallow shoreline areas warms up quickly in spring and summer, and cool down drastically in autumn and winter, showing large diurnal and seasonal variation amplitudes of SWT. Two cold belt zones in the western and eastern side of the lake and warm patches along the southwestern and northeastern shorelines are shaped by the combined effects of the lakebed topography and river runoff. Overall, the lake-averaged SWT increased at a rate of 0.26 °C/decade during 2001–2013. Faster increase of temperature was found at nighttime (0.34 °C/decade) and in winter and spring, consistent with the asymmetric warming pattern over land areas reported in prior studies. The rate of temperature increase over Siling Co is remarkably lower than that over Bangoin station, which is probably attributable to the large heat capacity of water and partly reflects the sensitive of alpine saltwater lake to climate change.     

Analysis of urban built-up areas and surface urban heat island using downscaled MODIS derived land surface temperature data

Regional scale urban built-up areas and surface urban heat islands (SUHI) are important for urban planning and policy formation. Owing to coarse spatial resolution (1000 m), it is difficult to use Moderate Resolution Imaging Spectroradiometer (MODIS) Land surface temperature (LST) products for mapping urban areas and visualization, and SUHI-related studies. To overcome this problem, the present study downscaled MODIS (1000 m resolution)-derived LST to 250 m resolution to map and visualize the urban areas and identify the basic components of SUHI over 12 districts of Punjab, India. The results are compared through visual interpretation and statistical procedure based on similarity analysis. The increased entropy value in the downscaled LST signifies higher information content. The temperature variation within the built-up and its environs is due to difference in land use and is depicted better in the downscaled LST. The SUHI intensity analysis of four cities (Ludhiana, Patiala, Moga and Vatinda) indicates that mean temperature in urban built-up core is higher (38.87 °C) as compared to suburban (35.85 °C) and rural (32.41 °C) areas. The downscaling techniques demonstrated in this paper enhance the usage of open-source wide swath MODIS LST for continuous monitoring of SUHI and urban area mapping, visualisation and analysis at regional scale. Such initiatives are useful for the scientific community and the decision-makers.     

Determination of the crop row orientations from Formosat-2 multi-temporal and panchromatic images

This paper presents a technique developed for the retrieval of the orientation of crop rows, over anthropic lands dedicated to agriculture in order to further improve estimate of crop production and soil erosion management. Five crop types are considered: wheat, barley, rapeseed, sunflower, corn and hemp. The study is part of the multi-sensor crop-monitoring experiment, conducted in 2010 throughout the agricultural season (MCM’10) over an area located in southwestern France, near Toulouse. The proposed methodology is based on the use of satellite images acquired by Formosat-2, at high spatial resolution in panchromatic and multispectral modes (with spatial resolution of 2 and 8 m, respectively). Orientations are derived and evaluated for each image and for each plot, using directional spatial filters (45° and 135°) and mathematical morphology algorithms. “Single-date” and “multi-temporal” approaches are considered. The single-date analyses confirm the good performances of the proposed method, but emphasize the limitation of the approach for estimating the crop row orientation over the whole landscape with only one date. The multi-date analyses allow (1) determining the most suitable agricultural period for the detection of the row orientations, and (2) extending the estimation to the entire footprint of the study area. For the winter crops (wheat, barley and rapeseed), best results are obtained with images acquired just after harvest, when surfaces are covered by stubbles or during the period of deep tillage (0.27 > R² > 0.99 and 7.15° > RMSE > 43.02°). For the summer crops (sunflower, corn and hemp), results are strongly crop and date dependents (0 > R² > 0.96, 10.22° > RMSE > 80°), with a well-marked impact of flowering, irrigation equipment and/or maximum crop development. Last, the extent of the method to the whole studied zone allows mapping 90% of the crop row orientations (more than 45,000 ha) with an error inferior to 40°, associated to a confidence index ranging from 1 to 5 for each agricultural plot.     

Influence of neighbourhood information on ‘Local Climate Zone’ mapping in heterogeneous cities

Local climate zone (LCZ) mapping is an emerging field in urban climate research. LCZs potentially provide an objective framework to assess urban form and function worldwide. The scheme is currently being used to globally map LCZs as a part of the World Urban Database and Access Portal Tools (WUDAPT) initiative. So far, most of the LCZ maps lack proper quantitative assessment, challenging the generic character of the WUDAPT workflow. Using the standard method introduced by the WUDAPT community difficulties arose concerning the built zones due to high levels of heterogeneity. To overcome this problem a contextual classifier is adopted in the mapping process. This paper quantitatively assesses the influence of neighbourhood information on the LCZ mapping result of three cities in Belgium: Antwerp, Brussels and Ghent. Overall accuracies for the maps were respectively 85.7 ± 0.5, 79.6 ± 0.9, 90.2 ± 0.4%. The approach presented here results in overall accuracies of 93.6 ± 0.2, 92.6 ± 0.3 and 95.6 ± 0.3% for Antwerp, Brussels and Ghent. The results thus indicate a positive influence of neighbourhood information for all study areas with an increase in overall accuracies of 7.9, 13.0 and 5.4%. This paper reaches two main conclusions. Firstly, evidence was introduced on the relevance of a quantitative accuracy assessment in LCZ mapping, showing that the accuracies reported in previous papers are not easily achieved. Secondly, the method presented in this paper proves to be highly effective in Belgian cities, and given its open character shows promise for application in other heterogeneous cities worldwide.     

An integrated remote sensing and GIS approach for monitoring areas affected by selective logging: A case study in northern Mato Grosso,Brazilian Amazon

Forest cover disturbances due to processes such as logging and forest fires are a widespread issue especially in the tropics, and have heavily affected forest biomass and functioning in the Brazilian Amazon in the past decades. Satellite remote sensing has played a key role for assessing logging activities in this region; however, there are still remaining challenges regarding the quantification and monitoring of these processes affecting forested lands. In this study, we propose a new method for monitoring areas affected by selective logging in one of the hotspots of Mato Grosso state in the Brazilian Amazon, based on a combination of object-based and pixel-based classification approaches applied on remote sensing data. Logging intensity and changes over time are assessed within grid cells of 300 m × 300 m spatial resolution. Our method encompassed three main steps: (1) mapping forest/non-forest areas through an object-based classification approach applied to a temporal series of Landsat images during the period 2000–2015, (2) mapping yearly logging activities from soil fraction images on the same Landsat data series, and (3) integrating information from previous steps within a regular grid-cell of 300 m × 300 m in order to monitor disturbance intensities over this 15-years period. The overall accuracy of the baseline forest/non-forest mask (year 2000) and of the undisturbed vs disturbed forest (for selected years) were 93% and 84% respectively. Our results indicate that annual forest disturbance rates, mainly due to logging activities, were higher than annual deforestation rates during the whole period of study. The deforested areas correspond to circa 25% of the areas affected by forest disturbances. Deforestation rates were highest from 2001 to 2005 and then decreased considerably after 2006. In contrast, the annual forest disturbance rates show high temporal variability with a slow decrease over the 15-year period, resulting in a significant increase of the ratio between disturbed and deforested areas. Although the majority of the areas, which have been affected by selective logging during the period 2000–2014, were not deforested by 2015, more than 70% of the deforested areas in 2015 had been at least once identified as disturbed forest during that period.     

SAR interferometry and optical remote sensing for analysis of co-seismic deformation,source characteristics and mass wasting pattern of Lushan (China,April 2013) earthquake

Co-seismic deformation associated with the Lushan (China) earthquake that occurred along the south-western segment of the Longmenshan Fault Zone (LFZ) on the 20th April 2013 has been estimated by differential interferometric SAR (DInSAR) technique using Radarsat-2 data. The Lushan earthquake resulted in the deformation of the Sichuan basin and the Longmenshan ranges in proximity to the LFZ. The line of sight (LOS) displacement values obtained from DInSAR technique mainly range between −4.0 cm to +3.0 cm. The western Sichuan basin shows oblique westward movement with predominant downward component in areas farther from LFZ and predominant westward component over the downward movement in areas closer to the source fault. Inversion modelling has been used to derive the seismic source characteristics from DInSAR derived deformation values using elastic dislocation source type. The linear inversion model converged at a double-fault source solution consisting of a deeper, steep, NW dipping fault plane-1 of 60 km × 16 km dimension and a shallower, gentle, NW dipping fault plane-2 of 60 km × 15 km dimension, with distributed slip values varying between 0 to 2.26 m. These fault planes (fault planes-1 and -2) coincide with the Dachuan-Shuangshi fault and the buried Range Front Fault, respectively. The inversion model gives a moment magnitude of 6.81 and the geodetic moment of 2.07 × 10¹⁹ Nm, comparable to those given in literature, derived using teleseismic body wave data. Thus DInSAR technique helped to quantify the co-seismic deformation and to retrieve the source characteristics from the estimated deformation values. The study also evaluated the distribution pattern of earthquake induced landslides (EIL) triggered fresh or re-activated during the Lushan earthquake and found that they show spatial association with the seismic source zone and also with various pre-conditioning factors of slope instability.     

Mapping seagrass coverage and spatial patterns with high spatial resolution IKONOS imagery

Seagrass habitats in subtidal coastal waters provide a variety of ecosystem functions and services and there is an increasing need to acquire information on spatial and temporal dynamics of this resource. Here, we explored the capability of IKONOS (IKO) data of high resolution (4 m) for mapping seagrass cover [submerged aquatic vegetation (%SAV) cover] along the mid-western coast of Florida, USA. We also compared seagrass maps produced with IKO data with that obtained using the Landsat TM sensor with lower resolution (30 m). Both IKO and TM data, collected in October 2009, were preprocessed to calculate water depth invariant bands to normalize the effect of varying depth on bottom spectra recorded by the two satellite sensors and further the textural information was extracted from IKO data. Our results demonstrate that the high resolution IKO sensor produced a higher accuracy than the TM sensor in a three-class % SAV cover classification. Of note is that the OA of %SAV cover mapping at our study area created with IKO data was 5–20% higher than that from other studies published. We also examined the spatial distribution of seagrass over a spatial range of 4–240 m using the Ripley’s K function [L(d)] and IKO data that represented four different grain sizes [4 m (one IKO pixel), 8 m (2 × 2 IKO pixels), 12 m (3 × 3 IKO pixels), and 16 m (4 × 4 IKO pixels)] from moderate-dense seagrass cover along a set of six transects. The Ripley’s K metric repeatedly indicated that seagrass cover representing 4 m × 4 m pixels displayed a dispersed (or slightly dispersed) pattern over distances of <4–8 m, and a random or slightly clustered pattern of cover over 9–240 m. The spatial pattern of seagrass cover created with the three additional grain sizes (i.e., 2 × 24 m IKO pixels, 3 × 34 m IKO pixels, and 4 × 4 m IKO pixels) show a dispersed (or slightly dispersed) pattern across 4–32 m and a random or slightly clustered pattern across 33–240 m. Given the first report on using satellite observations to quantify seagrass spatial patterns at a spatial scale from 4 m to 240 m, our novel analyses of moderate-dense SAV cover utilizing Ripley’s K function illustrate how data obtained from the IKO sensor revealed seagrass spatial information that would be undetected by the TM sensor with a 30 m pixel size. Use of the seagrass classification scheme here, along with data from the IKO sensor with enhanced resolution, offers an opportunity to synoptically record seagrass cover dynamics at both small and large spatial scales.     

Methodology for estimating availability of cloud-free image composites: A case study for southern Canada

Image composites are often used for earth surface phenomena studies at regional or national level. The compromise between residual clouds and the length of compositing period is a necessary corollary to the choice of satellite optical data for monitoring earth surface phenomena dynamics. This paper introduced a methodology for estimating availability of cloud-free image composites for optical sensors with various revisiting intervals, using MODIS MOD06 L2 cloud fraction product in the period of 2000–2008. The methodology starts with downscaling of the cloud fraction product to 1 km × 1 km cloud cover binary images. The binary images are then used for the exploration of spatial and temporal characteristics of cloud dynamics, and subsequently for the simulation of cloud-free composite availability with various revisiting intervals of optical sensors. Using Canada's southern provinces as an application case, the study explored several factors important for the design of environmental monitoring system using optical sensors of earth observation, in particular, cloud dynamics and its inter-annual variability, sensors’ revisiting intervals, and cloud-free threshold for targeting composites. While the cloud images used in the analysis are at 1 km × 1 km resolution, our analysis suggests that the simulated availabilities of cloud-free image composites may also provide reasonable estimates for optical sensors with higher than 1 km × 1 km resolution, though the closer to 1 km × 1 km resolution the optical sensor, the more pertinent the application. Also, the methodology can be parameterised to different temporal period and different spatial region, depending on applications.     

Detecting forest damage after a low-severity fire using remote sensing at multiple scales

Remote sensing technologies are an ideal platform to examine the extent and impact of fire on the landscape. In this study we assess that capacity of the RapidEye constellation and Landsat (Thematic Mapper and Operational Land Imager to map fine-scale burn attributes for a small, low severity prescribed fire in a dry Western Canadian forest. Estimates of burn severity from field data were collated into a simple burn index and correlated with a selected suite of common spectral vegetation indices. Burn severity classes were then derived to map fire impacts and estimate consumed woody surface fuels (diameter ≥2.6 cm). All correlations between the simple burn index and vegetation indices produced significant results (p < 0.01), but varied substantially in their overall accuracy. Although the Landsat Soil Adjusted Vegetation Index provided the best regression fit (R² = 0.56), results suggested that RapidEye provided much more spatially detailed estimates of tree damage (Soil Adjusted Vegetation Index, R² = 0.51). Consumption estimates of woody surface fuels ranged from 3.38 ± 1.03 Mg ha⁻¹ to 11.73 ± 1.84 Mg ha⁻¹, across four derived severity classes with uncertainties likely a result of changing foliage moisture between the before and after fire images. While not containing spectral information in the short wave infrared, the spatial variability provided by the RapidEye imagery has potential for mapping and monitoring fine scale forest attributes, as well as the potential to resolve fire damage at the individual tree level.     

Intercomparison and validation of MODIS and GLASS leaf area index (LAI) products over mountain areas: A case study in southwestern China

The validation study of leaf area index (LAI) products over rugged surfaces not only gives additional insights into data quality of LAI products, but deepens understanding of uncertainties regarding land surface process models depended on LAI data over complex terrain. This study evaluated the performance of MODIS and GLASS LAI products using the intercomparison and direct validation methods over southwestern China. The spatio-temporal consistencies, such as the spatial distributions of LAI products and their statistical relationship as a function of topographic indices, time, and vegetation types, respectively, were investigated through intercomparison between MODIS and GLASS products during the period 2011–2013. The accuracies and change ranges of these two products were evaluated against available LAI reference maps over 10 sampling regions which standed for typical vegetation types and topographic gradients in southwestern China.The results show that GLASS LAI exhibits higher percentage of good quality data (i.e. successful retrievals) and smoother temporal profiles than MODIS LAI. The percentage of successful retrievals for MODIS and GLASS is vulnerable to topographic indices, especially to relief amplitude. Besides, the two products do not capture seasonal dynamics of crop, especially in spring over heterogeneously hilly regions. The yearly mean LAI differences between MODIS and GLASS are within ±0.5 for 64.70% of the total retrieval pixels over southwestern China. The spatial distribution of mean differences and temporal profiles of these two products are inclined to be dominated by vegetation types other than topographic indices. The spatial and temporal consistency of these two products is good over most area of grasses/cereal crops; however, it is poor for evergreen broadleaf forest. MODIS presents more reliable change range of LAI than GLASS through comparison with fine resolution reference maps over most of sampling regions. The accuracies of direct validation are obtained for GLASS LAI (r = 0.35, RMSE = 1.72, mean bias = −0.71) and MODIS LAI (r = 0.49, RMSE = 1.75, mean bias = −0.67). GLASS performs similarly to MODIS, but may be marginally inferior to MODIS based on our direct validation results. The validation experience demonstrates the necessity and importance of topographic consideration for LAI estimation over mountain areas. Considerable attention will be paid to the improvements of surface reflectance, retrieval algorithm and land cover types so as to enhance the quality of LAI products in topographically complex terrain.     

Modeling biophysical properties of broad-leaved stands in the hyrcanian forests of Iran using fused airborne laser scanner data and ultraCam-D images

Inventories of mixed broad-leaved forests of Iran mainly rely on terrestrial measurements. Due to rapid changes and disturbances and great complexity of the silvicultural systems of these multilayer forests, frequent repetition of conventional ground-based plot surveys is often cost prohibitive. Airborne laser scanning (ALS) and multispectral data offer an alternative or supplement to conventional inventories in the Hyrcanian forests of Iran. In this study, the capability of a combination of ALS and UltraCam-D data to model stand volume, tree density, and basal area using random forest (RF) algorithm was evaluated. Systematic sampling was applied to collect field plot data on a 150 m × 200 m sampling grid within a 1100 ha study area located at 36°38′- 36°42′N and 54°24′–54°25′E. A total of 308 circular plots (0.1 ha) were measured for calculation of stand volume, tree density, and basal area per hectare. For each plot, a set of variables was extracted from both ALS and multispectral data. The RF algorithm was used for modeling of the biophysical properties using ALS and UltraCam-D data separately and combined. The results showed that combining the ALS data and UltraCam-D images provided a slight increase in prediction accuracy compared to separate modeling. The RMSE as percentage of the mean, the mean difference between observed and predicted values, and standard deviation of the differences using a combination of ALS data and UltraCam-D images in an independent validation at 0.1-ha plot level were 31.7%, 1.1%, and 84 m³ ha⁻¹ for stand volume; 27.2%, 0.86%, and 6.5 m² ha⁻¹ for basal area, and 35.8%, −4.6%, and 77.9 n ha⁻¹ for tree density, respectively. Based on the results, we conclude that fusion of ALS and UltraCam-D data may be useful for modeling of stand volume, basal area, and tree density and thus gain insights into structural characteristics in the complex Hyrcanian forests.     

Geologic remote sensing for geothermal exploration: A review

This paper is a comprehensive review of the potential for remote sensing in exploring for geothermal resources. Temperature gradients in the earth crust are typically 25–30 °C per kilometer depth, however in active volcanic areas situated in subduction or rift zones gradients of up to 150 °C per kilometer depth can be reached. In such volcanic areas, meteoric water in permeable and porous rocks is heated and hot water is trapped to form a geothermal reservoir. At the Earth's surface hot springs and fumaroles are evidence of hot geothermal water. In low enthalpy systems the heat can be used for heating/cooling and drying while in high enthalpy systems energy is generated using hot water or steam. In this paper we review the potential of remote sensing in the exploration for geothermal resources. We embark from the traditional suite of geophysical and geochemical prospecting techniques to arrive at parameters at the Earth surface that can be measured by earth observing satellites. Next, we summarize direct and indirect detection of geothermal potential using alteration mineralogy, temperature anomalies and heat fluxes, geobotanical anomalies and Earth surface deformation. A section of this paper is dedicated to published remote sensing studies illustrating the principles of mapping: surface deformation, gaseous emissions, mineral mapping, heat flux measurements, temperature mapping and geobotany. In a case study from the La Pacana caldera (Chili) geothermal field we illustrate the cross cutting relationships between various surface manifestations of geothermal reservoirs and how remotely sensed indicators can contribute to exploration. We conclude that although remote sensing of geothermal systems has not reached full maturity, there is great potential for integrating these surface measurements in a exploration framework. A number of recommendations for future research result from our analysis of geothermal systems and the present contributions of remote sensing to studying these systems. These are grouped along a number of question lines: ‘how reproducible are remote sensing products’, ‘can long term monitoring of geothermal systems be achieved’ and ‘do surface manifestations link to subsurface features’?     

3-D voxel-based solid modeling of a broad-leaved tree for accurate volume estimation using portable scanning lidar

We developed a method to produce a 3-D voxel-based solid model of a tree based on portable scanning lidar data for accurate estimation of the volume of the woody material. First, we obtained lidar measurements with a high laser pulse density from several measurement positions around the target, a Japanese zelkova tree. Next, we converted lidar-derived point-cloud data for the target into voxels. The voxel size was 0.5 cm × 0.5 cm × 0.5 cm. Then, we used differences in the spatial distribution of voxels to separate the stem and large branches (diameter > 1 cm) from small branches (diameter ⩽ 1 cm). We classified the voxels into sets corresponding to the stem and to each large branch and then interpolated voxels to fill out their surfaces and their interiors. We then merged the stem and large branches with the small branches. The resultant solid model of the entire tree was composed of consecutive voxels that filled the outer surface and the interior of the stem and large branches, and a cloud of voxels equivalent to small branches that were discretely scattered in mainly the upper part of the target. Using this model, we estimated the woody material volume by counting the number of voxels in each part and multiplying the number of voxels by the unit voxel volume (0.13 cm³). The percentage error of the volume of the stem and part of a large branch was 0.5%. The estimation error of a certain part of the small branches was 34.0%.     

Timber production assessment of a plantation forest: An integrated framework with field-based inventory,multi-source remote sensing data and forest management history

Timber production is the purpose for managing plantation forests, and its spatial and quantitative information is critical for advising management strategies. Previous studies have focused on growing stock volume (GSV), which represents the current potential of timber production, yet few studies have investigated historical process-harvested timber. This resulted in a gap in a synthetical ecosystem service assessment of timber production. In this paper, we established a Management Process–based Timber production (MPT) framework to integrate the current GSV and the harvested timber derived from historical logging regimes, trying to synthetically assess timber production for a historical period. In the MPT framework, age-class and current GSV determine the times of historical thinning and the corresponding harvested timber, by using a “space-for-time” substitution. The total timber production can be estimated by the historical harvested timber in each thinning and the current GSV. To test this MPT framework, an empirical study on a larch plantation (LP) with area of 43,946 ha was conducted in North China for a period from 1962 to 2010. Field-based inventory data was integrated with ALOS PALSAR (Advanced Land-Observing Satellite Phased Array L-band Synthetic Aperture Radar) and Landsat-8 OLI (Operational Land Imager) data for estimating the age-class and current GSV of LP. The random forest model with PALSAR backscatter intensity channels and OLI bands as input predictive variables yielded an accuracy of 67.9% with a Kappa coefficient of 0.59 for age-class classification. The regression model using PALSAR data produced a root mean square error (RMSE) of 36.5 m³ ha⁻¹. The total timber production of LP was estimated to be 7.27 × 10⁶ m³, with 4.87 × 10⁶ m³ in current GSV and 2.40 × 10⁶ m³ in harvested timber through historical thinning. The historical process-harvested timber accounts to 33.0% of the total timber production, which component has been neglected in the assessments for current status of plantation forests. Synthetically considering the RMSE for predictive GSV and misclassification of age-class, the error in timber production were supposed to range from −55.2 to 56.3 m³ ha⁻¹. The MPT framework can be used to assess timber production of other tree species at a larger spatial scale, providing crucial information for a better understanding of forest ecosystem service.     

Understanding long-term (1982–2013) patterns and trends in winter wheat spring green-up date over the North China Plain

Monitoring the spring green-up date (GUD) has grown in importance for crop management and food security. However, most satellite-based GUD models are associated with a high degree of uncertainty when applied to croplands. In this study, we introduced an improved GUD algorithm to extract GUD data for 32 years (1982–2013) for the winter wheat croplands on the North China Plain (NCP), using the third-generation normalized difference vegetation index form Global Inventory Modeling and Mapping Studies (GIMMS3g NDVI). The spatial and temporal variations in GUD with the effects of the pre-season climate and soil moisture conditions on GUD were comprehensively investigated. Our results showed that a higher correlation coefficient (r = 0.44, p < 0.01) and lower root mean square error (22 days) and bias (16 days) were observed in GUD from the improved algorithm relative to GUD from the MCD12Q2 phenology product. In spatial terms, GUD increased from the southwest (less than day of year (DOY) 60) to the northeast (more than DOY 90) of the NCP, which corresponded to spatial reductions in temperature and precipitation. GUD advanced in most (78%) of the winter wheat area on the NCP, with significant advances in 37.8% of the area (p < 0.05). GUD occurred later at high altitudes and in coastal areas than in inland areas. At the interannual scale, the average GUD advanced from DOY 76.9 in the 1980s (average 1982–1989) to DOY 73.2 in the 1990s (average 1991–1999), and to DOY 70.3 after 2000 (average 2000–2013), indicating an average advance of 1.8 days/decade (r = 0.35, p < 0.05). Although GUD is mainly controlled by the pre-season temperature, our findings underline that the effect of the pre-season soil moisture on GUD should also be considered. The improved GUD algorithm and satellite-based long-term GUD data are helpful for improving the representation of GUD in terrestrial ecosystem models and enhancing crop management efficiency.     

Agriculture pest and disease risk maps considering MSG satellite data and land surface temperature

Pest risk maps for agricultural use are usually constructed from data obtained from in-situ meteorological weather stations, which are relatively sparsely distributed and are often quite expensive to install and difficult to maintain. This leads to the creation of maps with relatively low spatial resolution, which are very much dependent on interpolation methodologies. Considering that agricultural applications typically require a more detailed scale analysis than has traditionally been available, remote sensing technology can offer better monitoring at increasing spatial and temporal resolutions, thereby, improving pest management results and reducing costs. This article uses ground temperature, or land surface temperature (LST), data distributed by EUMETSAT/LSASAF (with a spatial resolution of 3 × 3 km (nadir resolution) and a revisiting time of 15 min) to generate one of the most commonly used parameters in pest modeling and monitoring: “thermal integral over air temperature (accumulated degree-days)”. The results show a clear association between the accumulated LST values over a threshold and the accumulated values computed from meteorological stations over the same threshold (specific to a particular tomato pest). The results are very promising and enable the production of risk maps for agricultural pests with a degree of spatial and temporal detail that is difficult to achieve using in-situ meteorological stations.     

Land subsidence in major cities of Central Mexico: Interpreting InSAR-derived land subsidence mapping with hydrogeological data

Significant structural damages to urban infrastructures caused by compaction of over-exploited aquifers are an important problem in Central Mexico. While the case of Mexico City has been well-documented, insight into land subsidence problems in other cities of Central Mexico is still limited. Among the cities concerned, we present and discuss the cases of five of them, located within the Trans-Mexican Volcanic Belt (TMVB): Toluca, Celaya, Aguascalientes, Morelia, and Queretaro. Applying the SBAS-InSAR method to C-Band RADARSAT-2 data, five high resolution ground motion time-series were produced to monitor the spatio-temporal variations of displacements and fracturing from 2012 to 2014. The study presents recent changes of land subsidence rates along with concordant geological and water data. It aims to provide suggestions to mitigate future damages to infrastructure and to assist in groundwater resources management.Aguascalientes, Celaya, Morelia and Queretaro (respectively in order of decreasing subsidence rates) are typical cases of fault-limited land subsidence of Central Mexico. It occurs as a result of groundwater over-exploitation in lacustrine and alluvial deposits covering highly variable bedrock topography, typical of horst-graben geological settings. Aguascalientes and Toluca show high rates of land subsidence (up to 10 cm/yr), while Celaya and Morelia show lower rates (from 2 to 5 cm/yr). Comparing these results with previous studies, it is inferred that the spatial patterns of land subsidence have changed in the city of Toluca. This change appears to be mainly controlled by the spatial heterogeneity of compressible sediments since no noticeable change occurred in groundwater extraction and related drawdown rates. While land subsidence of up to 8 cm/yr has been reported in the Queretaro Valley before 2011, rates inferior to 1 cm/yr are measured in 2013–2014. The subsidence has been almost entirely mitigated by major changes in the water management practices of the city, i.e., a 122 km long pipeline bringing surface water from an adjacent state allowed to cease pumping in half of the wells.     

Spatio-temporal variation of vegetation coverage and its response to climate change in North China plain in the last 33 years

Global climate change has led to significant vegetation changes in the past half century. North China Plain, the most important grain production base of china, is undergoing a process of prominent warming and drying. The vegetation coverage, which is used to monitor vegetation change, can respond to climate change (temperature and precipitation). In this study, GIMMS (Global Inventory Modelling and Mapping Studies)-NDVI (Normalized Difference Vegetation Index) data, MODIS (Moderate-resolution Imaging Spectroradiometer) – NDVI data and climate data, during 1981–2013, were used to investigate the spatial distribution and changes of vegetation. The relationship between climate and vegetation on different spatial (agriculture, forest and grassland) and temporal (yearly, decadal and monthly) scales were also analyzed in North China Plain. (1) It was found that temperature exhibiting a slight increase trend (0.20 °C/10a, P < 0.01). This may be due to the disappearance of 0 °C isotherm, the rise of spring temperature. At the same time, precipitation showed a significant reduction trend (−1.75 mm/10a, P > 0.05). The climate mutation period was during 1991–1994. (2) Vegetation coverage slight increase was observed in the 55% of total study area, with a change rate of 0.00039/10a. Human activities may not only accelerate the changes of the vegetation coverage, but also c effect to the rate of these changes. (3) Overall, the correlation between the vegetation coverage and climatic factor is higher in monthly scale than yearly scale. The correlation analysis between vegetation coverage and climate changes showed that annual vegetation coverage was better correlatend with precipitation in grassland biome; but it showed a better correlated with temperature i the agriculture biome and forest biome. In addition, the vegetation coverage had sensitive time-effect respond to precipitation. (4) The vegetation coverage showed the same increasing trend before and after the climatic variations, but the rate of increase slowed down. From the vegetation coverage point of view, the grassland ecological zone had an obvious response to the climatic variations, but the agricultural ecological zones showed a significant response from the vegetation coverage change rate point of view. The effect of human activity in degradation region was higher than that in improvement area. But after the climate abruptly changing, the effect of human activity in improvement area was higher than that in degradation region, and the influence of human activity will continue in the future.     

Augmenting forest inventory attributes with geometric optical modelling in support of regional susceptibility assessments to bark beetle infestations

Assessment of the susceptibility of forests to mountain pine beetle (Dendroctonus ponderosae Hopkins) infestation is based upon an understanding of the characteristics that predispose the stands to attack. These assessments are typically derived from conventional forest inventory data; however, this information often represents only managed forest areas. It does not cover areas such as forest parks or conservation regions and is often not regularly updated resulting in an inability to assess forest susceptibility. To address these shortcomings, we demonstrate how a geometric optical model (GOM) can be applied to Landsat-5 Thematic Mapper (TM) imagery (30 m spatial resolution) to estimate stand-level susceptibility to mountain pine beetle attack. Spectral mixture analysis was used to determine the proportion of sunlit canopy and background, and shadow of each Landsat pixel enabling per pixel estimates of attributes required for model inversion. Stand structural attributes were then derived from inversion of the geometric optical model and used as basis for susceptibility mapping. Mean stand density estimated by the geometric optical model was 2753 (standard deviation ± 308) stems per hectare and mean horizontal crown radius was 2.09 (standard deviation ± 0.11) metres. When compared to equivalent forest inventory attributes, model predictions of stems per hectare and crown radius were shown to be reasonably estimated using a Kruskal–Wallis ANOVA (p < 0.001). These predictions were then used to create a large area map that provided an assessment of the forest area susceptible to mountain pine beetle damage.