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.     

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.     

InSAR to support sustainable urbanization over compacting aquifers: The case of Toluca Valley,Mexico

This paper illustrates how InSAR alone can be used to delineate potential ground fractures related to aquifer system compaction. An InSAR-derived ground fracturing map of the Toluca Valley, Mexico, is produced and validated through a field campaign. The results are of great interest to support sustainable urbanization and show that InSAR processing of open-access Synthetic Aperture Radar (SAR) data from the Sentinel-1 satellites can lead to reliable and cost-effective products directly usable by cities to help decision-making.The Toluca Valley Aquifer (TVA) sustains the water needs of two million inhabitants living within the valley, a growing industry, an intensively irrigated agricultural area, and 38% of the water needs of the megalopolis of Mexico City, located 40 km east of the valley. Ensuring water sustainability, infrastructure integrity, along with supporting the important economic and demographic growth of the region, is a major challenge for water managers and urban developers. This paper presents a long-term analysis of ground fracturing by interpreting 13 years of InSAR-derived ground displacement measurements. Small Baseline Subset (SBAS) and Persistent Scatterer Interferometry (PSI) techniques are applied over three SAR datasets totalling 93 acquisitions from Envisat, Radarsat-2, and Sentinel-1A satellites and covering the period from 2003 to 2016.From 2003 to 2016, groundwater level declines of up to 1.6 m/yr, land subsidence up to 77 mm/yr, and major infrastructure damages are observed. Groundwater level data show highly variable seasonal responses according to their connectivity to recharge areas. However, the trend of groundwater levels consistently range from −0.5 to −1.5 m/yr regardless of the well location and depth. By analysing the horizontal gradients of vertical land subsidence, we provide a potential ground fracture map to assist in future urban development planning in the Toluca Valley.     

20 years of active deformation on volcano caldera: Joint analysis of InSAR and AInSAR techniques

InSAR (interferometric synthetic aperture radar) techniques are applied to investigate last two decades of surface deformation of the Cerro Blanco/Robledo Caldera (CBRC). The objective is the identification of deforming patterns that alter the shape of these complex structures when they show low or null activity. The joint analysis between results by using different methods over a long time span, represents a unique opportunity to improve knowledge of volcanic structures located in remote area and, for this, poorly or not monitored.In this work we identify displacement patterns over the volcanic area, by using both classical differential InSAR analysis, and A-InSAR (advanced InSAR) analysis based on SAR data acquired by ERS-1/2 and ENVISAT sensors during the 1996–2010 time interval. The satellite-derived information allows us to characterize the deformation pattern that affected the CBRC and shows that the actively deforming CBRC is subsiding in the observed period. In order to figure out the deformation history of CBRC, we analyzed the four sub-periods 1992–1996, 1996–2000, and 2005–2010 by using standard differential InSAR technique, and the interval 2003–2007 by adopting an A-InSAR technique.Subsidence velocities of the CBRC caldera are about 2.6 cm/yr in the time interval 1992–1996 (measured with ERS descending data), 1.8 cm/yr in 1996–2000 (ERS descending data), 1.2 cm/yr in 2003–2007 (ENVISAT descending data), and finally, 0.87 cm/yr in 2005–2010 (ENVISAT ascending data). Moreover, outside the caldera and in particular in the NW area, we observe the presence of positive velocity values. Results show that: (a) a decreasing subsidence rate might be related to the reduction of volcanic activity in correspondence of the CBRC; (b) positive velocity signal, decreasing with time, might be interpreted as follows: – evidence of volcano structure lateral spreading, according to the velocity pattern distribution in this area and to the relative local flanks topographic convexity of the volcano structure; – uplift signal of this sector of mountain chain; – combination of the two mechanisms above.     

Modeling the spatial distribution of above-ground carbon in Mexican coniferous forests using remote sensing and a geostatistical approach

Forest conservation is considered an option for mitigating the effect of greenhouse gases on global climate, hence monitoring forest carbon pools at global and local levels is important. The present study explores the capability of remote-sensing variables (vegetation indices and textures derived from SPOT-5; backscattering coefficient and interferometric coherence of ALOS PALSAR images) for modeling the spatial distribution of above-ground biomass in the Environmental Conservation Zone of Mexico City. Correlation and spatial autocorrelation coefficients were used to select significant explanatory variables in fir and pine forests. The correlation for interferometric coherence in HV polarization was negative, with correlations coefficients r = −0.83 for the fir and r = −0.75 for the pine forests. Regression-kriging showed the least root mean square error among the spatial interpolation methods used, with 37.75 tC/ha for fir forests and 29.15 tC/ha for pine forests. The results showed that a hybrid geospatial method, based on interferometric coherence data and a regression-kriging interpolator, has good potential for estimating above-ground biomass carbon.     

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.     

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.     

Identifying natural and anthropogenically-induced geohazards from satellite ground motion and geospatial data: Stoke-on-Trent,UK

Determining the location and nature of hazardous ground motion resulting from natural and anthropogenic processes such as landslides, tectonic movement and mining is essential for hazard mitigation and sustainable resource use. Ground motion estimates from satellite ERS-1/2 persistent scatterer interferometry (PSI) were combined with geospatial data to identify areas of observed geohazards in Stoke-on-Trent, UK. This investigation was performed within the framework of the EC FP7-SPACE PanGeo project which aimed to provide free and open access to geohazard information for 52 urban areas across Europe. Geohazards identified within the city of Stoke-on-Trent and neighbouring rural areas are presented here alongside an examination of the PanGeo methodology.A total of 14 areas experiencing ground instability caused by natural and anthropogenic processes have been defined, covering 122.35 km². These are attributed to a range of geohazards, including landslides, ground dissolution, made ground and mining activities. The dominant geohazard (by area) is ground movement caused by post-mining groundwater recharge and mining-related subsidence (93.19% of total geohazard area), followed by landsliding (5.81%). Observed ground motions along the satellite line-of-sight reach maxima of +35.23 mm/yr and −22.57 mm/yr. A combination of uplift, subsidence and downslope movement is displayed.‘Construction sites’ and ‘continuous urban fabric’ (European Urban Atlas land use types) form the land uses most affected (by area) by ground motion and ‘discontinuous very low density urban fabric’ the least. Areas of ‘continuous urban fabric’ also show the highest average velocity towards the satellite (5.08 mm/yr) and the highest PS densities (1262.92 points/km²) along with one of the lowest standard deviations. Rural land uses tend to result in lower PS densities and higher standard deviations, a consequence of fewer suitable reflectors in these regions. PSI is also limited in its ability to identify especially rapid ground motion. As a consequence the supporting geospatial data proved especially useful for the identification of landslides and some areas of ground dissolution. The mapped areas of instability are also compared with modelled potential geohazards (the BGS GeoSure dataset).     

Estimating temporal changes in soil carbon stocks at ecoregional scale in Madagascar using remote-sensing

Soil organic carbon (SOC) plays an important role in climate change regulation notably through release of CO2 following land use change such a deforestation, but data on stock change levels are lacking. This study aims to empirically assess SOC stocks change between 1991 and 2011 at the landscape scale using easy-to-access spatially-explicit environmental factors. The study area was located in southeast Madagascar, in a region that exhibits very high rate of deforestation and which is characterized by both humid and dry climates. We estimated SOC stock on 0.1 ha plots for 95 different locations in a 43,000 ha reference area covering both dry and humid conditions and representing different land cover including natural forest, cropland, pasture and fallows. We used the Random Forest algorithm to find out the environmental factors explaining the spatial distribution of SOC. We then predicted SOC stocks for two soil layers at 30 cm and 100 cm over a wider area of 395,000 ha. By changing the soil and vegetation indices derived from remote sensing images we were able to produce SOC maps for 1991 and 2011. Those estimates and their related uncertainties where combined in a post-processing step to map estimates of significant SOC variations and we finally compared the SOC change map with published deforestation maps. Results show that the geologic variables, precipitation, temperature, and soil-vegetation status were strong predictors of SOC distribution at regional scale. We estimated an average net loss of 10.7% and 5.2% for the 30 cm and the 100 cm layers respectively for deforested areas in the humid area. Our results also suggest that these losses occur within the first five years following deforestation. No significant variations were observed for the dry region. This study provides new solutions and knowledge for a better integration of soil threats and opportunities in land management policies.     

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.     

Estimate of heavy metals in soil and streams using combined geochemistry and field spectroscopy in Wan-sheng mining area,Chongqing, China

Heavy metals contaminated soils and water will become a major environmental issue in the mining areas. This paper intends to use field hyper-spectra to estimate the heavy metals in the soil and water in Wan-sheng mining area in Chongqing. With analyzing the spectra of soil and water, the spectral features deriving from the spectral of the soils and water can be found to build the models between these features and the contents of Al, Cu and Cr in the soil and water by using the Stepwise Multiple Linear Regression (SMLR). The spectral features of Al are: 480 nm, 500 nm, 565 nm, 610 nm, 680 nm, 750 nm, 1000 nm, 1430 nm, 1755 nm, 1887 nm, 1920 nm, 1950 nm, 2210 nm, 2260 nm; The spectral features of Cu are: 480 nm, 500 nm, 610 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1920 nm, 2150 nm, 2260 nm; And the spectral features of Cr are: 480 nm, 500 nm, 610 nm, 715 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1755 nm, 1920 nm, 1950 nm. With these features, the best models to estimate the heavy metals in the study area were built according to the maximal R². The R² of the models of estimating Al, Cu and Cr in the soil and water are 0.813, 0.638, 0.604 and 0.742, 0.584, 0.513 respectively. And the gradient maps of these three types of heavy metals’ concentrations can be created by using the Inverse distance weighted (IDW).The gradient maps indicate that the heavy metals in the soil have similar patterns, but in the North-west of the streams in the study area, the contents are of great differences. These results show that it is feasible to predict contaminated heavy metals in the soils and streams due to mining activities by using the rapid and cost-effective field spectroscopy.     

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.     

Using remote sensing to monitor the influence of river discharge on watershed outlets and adjacent coral Reefs: Magdalena River and Rosario Islands,Colombia

Worldwide, coral reef ecosystems are being increasingly threatened by sediments loads from river discharges, which in turn are influenced by changing rainfall patterns due to climate change and by growing human activity in their watersheds. In this case study, we explored the applicability of using remote sensing (RS) technology to estimate and monitor the relationship between water quality at the coral reefs around the Rosario Islands, in the Caribbean Sea, and the rainfall patterns in the Magdalena River watershed. From the Moderate Resolution Imaging Spectroradiometer (MODIS), this study used the water surface reflectance product (MOD09GQ) to estimate water surface reflectance as a proxy for sediment concentration and the land cover product (MCD12Q1 V51) to characterize land cover of the watershed. Rainfall was estimated by using the 3B43 V7 product from the Tropical Rainforest Measuring Mission (TRMM). For the first trimester of each year, we investigated the inter-annual temporal variation in water surface reflectance at the Rosario Islands and at the three main mouths of the Magdalena River watershed. No increasing or decreasing trends of water surface reflectance were detected for any of the sites for the study period 2001–2014 (p > 0.05) but significant correlations were detected among the trends of each site at the watershed mouths (r = 0.57–0.90, p < 0.05) and between them and the inter-annual variation in rainfall on the watershed (r = 0.63–0.67, p < 0.05). Those trimesters with above-normal water surface reflectance at the mouths and above-normal rainfall at the watershed coincided with La Niña conditions while the opposite was the case during El Niño conditions. Although, a preliminary analysis of inter-annual land cover trends found only cropland cover in the watershed to be significantly correlated with water surface reflectance at two of the watershed mouths (r = 0.58 and 0.63, p < 0.05), the validation analysis draw only a 40.7% of accuracy in this land cover classification. This requires further analysis to confirm the impact of the cropland on the water quality at the watershed outlets. Spatial analysis with MOD09GQ imagery detected the overpass of river plumes from Barbacoas Bay over the Rosario Islands waters.     

Vegetation chlorophyll estimates in the Amazon from multi-angle MODIS observations and canopy reflectance model

As a preparatory study for future hyperspectral missions that can measure canopy chemistry, we introduce a novel approach to investigate whether multi-angle Moderate Resolution Imaging Spectroradiometer (MODIS) data can be used to generate a preliminary database with long-term estimates of chlorophyll. MODIS monthly chlorophyll estimates between 2000 and 2015, derived from a fully coupled canopy reflectance model (ProSAIL), were inspected for consistency with eddy covariance fluxes, tower-based hyperspectral images and chlorophyll measurements. MODIS chlorophyll estimates from the inverse model showed strong seasonal variations across two flux-tower sites in central and eastern Amazon. Marked increases in chlorophyll concentrations were observed during the early dry season. Remotely sensed chlorophyll concentrations were correlated to field measurements (r² = 0.73 and r² = 0.98) but the data deviated from the 1:1 line with root mean square errors (RMSE) ranging from 0.355 μg cm⁻² (Tapajós tower) to 0.470 μg cm⁻² (Manaus tower). The chlorophyll estimates were consistent with flux tower measurements of photosynthetically active radiation (PAR) and net ecosystem productivity (NEP). We also applied ProSAIL to mono-angle hyperspectral observations from a camera installed on a tower to scale modeled chlorophyll pigments to MODIS observations (r² = 0.73). Chlorophyll pigment concentrations (Chl_A+B) were correlated to changes in the amount of young and mature leaf area per month (0.59 ≤  r² ≤  0.64). Increases in MODIS observed Chl_A+B were preceded by increased PAR during the dry season (0.61 ≤ r² ≤  0.62) and followed by changes in net carbon uptake. We conclude that, at these two sites, changes in LAI, coupled with changes in leaf chlorophyll, are comparable with seasonality of plant productivity. Our results allowed the preliminary development of a 15-year time series of chlorophyll estimates over the Amazon to support canopy chemistry studies using future hyperspectral sensors.     

Upscaling plot-scale soil respiration in winter wheat and summer maize rotation croplands in Julu County,North China

Soil respiration (R_s) data from 45 plots were used to estimate the spatial patterns of R_s during the peak growing seasons of winter wheat and summer maize in Julu County, North China, by combining satellite remote sensing data, field-measured data, and a support vector regression (SVR) model. The observed R_s values were well reproduced by the model at the plot scale, with a root-mean-square error (RMSE) of 0.31 μmol CO₂ m⁻² s⁻¹ and a coefficient of determination (R²) of 0.73. No significant difference was detected between the prediction accuracy of the SVR model for winter wheat and summer maize. With forcing from satellite remote sensing data and gridded soil property data, we used the SVR model to predict the spatial distributions of R_s during the peak growing seasons of winter wheat and summer maize rotation croplands in Julu County. The SVR model captured the spatial variations of R_s at the county scale. The satellite-derived enhanced vegetation index was found to be the most important input used to predict R_s. Removal of this variable caused an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.42 μmol CO₂ m⁻² s⁻¹. Soil properties such as soil organic carbon (SOC) content and soil bulk density (SBD) were the second most important factors. Their removal led to an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.37 μmol CO₂ m⁻² s⁻¹. The SVR model performed better than multiple regression in predicting spatial variations of R_s in winter wheat and summer maize rotation croplands, as shown by the comparison of the R² and RMSE values of the two algorithms. The spatial patterns of R_s are better captured using the SVR model than performing multiple regression, particularly for the relatively high and relatively low R_s values at the center and northeast study areas. Therefore, SVR shows promise for predicting spatial variations of R_s values on the basis of remotely sensed data and gridded soil property data at the county scale.     

Optimising three-band spectral indices to assess aerial N concentration,N uptake and aboveground biomass of winter wheat remotely in China and Germany

Remotely and accurately quantifying the canopy nitrogen status in crops is essential for regional studies of N budgets and N balances. In this study, we optimised three-band spectral algorithms to estimate the N status of winter wheat. This study extends previous work to optimise the band combinations further and identifies the optimised central bands and suitable bandwidths of the three-band nitrogen planar domain index (NPDI) for estimating the aerial N uptake, N concentration and aboveground biomass. Analysis of the influence of bandwidth change on the accuracy of estimating the canopy N status and aboveground biomass indicated that the suitable bandwidths for optimised central bands were 37 nm at 846 nm, 13 nm at 738 nm and 57 nm at 560 nm for assessing the aerial N uptake and were 37 nm at 958 nm, 21 nm at 696 nm and 73 nm at 578 nm for the assessment of the aerial N concentration and were 49 nm at 806 nm, 17 nm at 738 nm and 57 nm at 560 nm for the estimation of aboveground biomass. The optimised three-band NPDI could consistently and stably estimate the aerial N uptake and aboveground biomass of winter wheat in the vegetative stage and the aerial N concentration in the reproductive stage compared to the fixed band combinations. With suitable bandwidths, the broadband NPDI demonstrated excellent performance in estimating the aerial N concentration, N uptake and biomass. We conclude that the band-optimised algorithm represents a promising tool to measure the improved performance of the NPDI in estimating the aerial N uptake and biomass in the vegetative stage and the aerial N concentration in the reproductive stage, which will be useful for designing improved nitrogen diagnosis systems and for enhancing the applications of ground- and satellite-based sensors.     

The July 11, 1995 Myanmar–China earthquake: A representative event in the bookshelf faulting system of southeastern Asia observed from JERS-1 SAR images

On July 11, 1995, an Mw 6.8 earthquake struck eastern Myanmar near the Chinese border; hereafter referred to as the 1995 Myanmar–China earthquake. Coseismic surface displacements associated with this event are identified from JERS-1 (Japanese Earth Resources Satellite-1) SAR (Synthetic Aperture Radar) images. The largest relative displacement reached 60 cm in the line-of-sight direction. We speculate that a previously unrecognized dextral strike-slip subvertical fault striking NW–SE was responsible for this event. The coseismic slip distribution on the fault planes is inverted based on the InSAR-derived deformation. The results indicate that the fault slip was confined to two lobes. The maximum slip reached approximately 2.5 m at a depth of 5 km in the northwestern part of the focal region. The inverted geodetic moment was approximately Mw = 6.69, which is consistent with seismological results. The 1995 Myanmar–China earthquake is one of the largest recorded earthquakes that has occurred around the “bookshelf faulting” system between the Sagaing fault in Myanmar and the Red River fault in southwestern China.     

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.     

Global rain-fed,irrigated, and paddy croplands: A new high resolution map derived from remote sensing,crop inventories and climate data

Irrigation accounts for 70% of global water use by humans and 33–40% of global food production comes from irrigated croplands. Accurate and timely information related to global irrigation is therefore needed to manage increasingly scarce water resources and to improve food security in the face of yield gaps, climate change and extreme events such as droughts, floods, and heat waves. Unfortunately, this information is not available for many regions of the world. This study aims to improve characterization of global rain-fed, irrigated and paddy croplands by integrating information from national and sub-national surveys, remote sensing, and gridded climate data sets. To achieve this goal, we used supervised classification of remote sensing, climate, and agricultural inventory data to generate a global map of irrigated, rain-fed, and paddy croplands. We estimate that 314 million hectares (Mha) worldwide were irrigated circa 2005. This includes 66 Mha of irrigated paddy cropland and 249 Mha of irrigated non-paddy cropland. Additionally, we estimate that 1047 Mha of cropland are managed under rain-fed conditions, including 63 Mha of rain-fed paddy cropland and 985 Mha of rain-fed non-paddy cropland. More generally, our results show that global mapping of irrigated, rain-fed, and paddy croplands is possible by combining information from multiple data sources. However, regions with rapidly changing irrigation or complex mixtures of irrigated and non-irrigated crops present significant challenges and require more and better data to support high quality mapping of irrigation.