Using Imaging Spectroscopy to study soil properties

Imaging Spectroscopy (IS) is a promising tool for studying soil properties in large spatial domains. Going from point to image spectrometry is not only a journey from micro to macro scales, but also a long stage where problems such as dealing with data having a low signal-to-noise level, contamination of the atmosphere, large data sets, the BRDF effect and more are often encountered. In this paper we provide an up-to-date overview of some of the case studies that have used IS technology for soil science applications. Besides a brief discussion on the advantages and disadvantages of IS for studying soils, the following cases are comprehensively discussed: soil degradation (salinity, erosion, and deposition), soil mapping and classification, soil genesis and formation, soil contamination, soil water content, and soil swelling. We review these case studies and suggest that the 15 data be provided to the end-users as real reflectance and not as raw data and with better signal-to-noise ratios than presently exist. This is because converting the raw data into reflectance is a complicated stage that requires experience, knowledge, and specific infrastructures not available to many users, whereas quantitative spectral models require good quality data. These limitations serve as a barrier that impedes potential end-users, inhibiting researchers from trying this technique for their needs. The paper ends with a general call to the soil science audience to extend the utilization of the IS technique, and it provides some ideas on how to propel this technology forward to enable its widespread adoption in order to achieve a breakthrough in the field of soil science and remote sensing. (C) 2009 Elsevier Inc. All rights reserved.

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and interannually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R(2) = 0.80 and R(2) = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages. (C) 2010 Elsevier Inc. All rights reserved.

MODIS NDVI time-series allow the monitoring of Eucalyptus plantation biomass

The use of remote sensing is necessary for monitoring forest carbon stocks at large scales. Optical remote sensing, although not the most suitable technique for the direct estimation of stand biomass, offers the advantage of providing large temporal and spatial datasets. In particular, information on canopy structure is encompassed in stand reflectance time series. This study focused on the example of Eucalyptus forest plantations, which have recently attracted much attention as a result of their high expansion rate in many tropical countries. Stand scale time-series of Normalized Difference Vegetation Index (NDVI) were obtained from MODIS satellite data after a procedure involving un-mixing and interpolation, on about 15,000 ha of plantations in southern Brazil. The comparison of the planting date of the current rotation (and therefore the age of the stands) estimated from these time series with real values provided by the company showed that the root mean square error was 35.5 days. Age alone explained more than 82% of stand wood volume variability and 87% of stand dominant height variability. Age variables were combined with other variables derived from the NDVI time series and simple bioclimatic data by means of linear (Stepwise) or nonlinear (Random Forest) regressions. The nonlinear regressions gave r-square values of 0.90 for volume and 0.92 for dominant height, and an accuracy of about 25 m(3)/ha for volume (15% of the volume average value) and about 1.6 m for dominant height (8% of the height average value). The improvement including NDVI and bioclimatic data comes from the fact that the cumulative NDVI since planting date integrates the interannual variability of leaf area index (LAI), light interception by the foliage and growth due for example to variations of seasonal water stress. The accuracy of biomass and height predictions was strongly improved by using the NDVI integrated over the two first years after planting, which are critical for stand establishment. These results open perspectives for cost-effective monitoring of biomass at large scales in intensively-managed plantation forests. (C) 2011 Elsevier Inc. All rights reserved.

Relating MODIS vegetation index time-series with structure, light absorption and stem production of fast-growing Eucalyptus plantations

By allowing the estimation of forest structural and biophysical characteristics at different temporal and spatial scales, remote sensing may contribute to our understanding and monitoring of planted forests. Here, we studied 9-year time-series of the Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) on a network of 16 stands in fast-growing Eucalyptus plantations in Sao Paulo State, Brazil. We aimed to examine the relationships between NDVI time-series spanning entire rotations and stand structural characteristics (volume, dominant height, mean annual increment) in these simple forest ecosystems. Our second objective was to examine spatial and temporal variations of light use efficiency for wood production, by comparing time-series of Absorbed Photosynthetically Active Radiation (APAR) with inventory data. Relationships were calibrated between the NDVI and the fractions of intercepted diffuse and direct radiation, using hemispherical photographs taken on the studied stands at two seasons. APAR was calculated from the NDVI time-series using these relationships. Stem volume and dominant height were strongly correlated with summed NDVI values between planting date and inventory date. Stand productivity was correlated with mean NDVI values. APAR during the first 2 years of growth was variable between stands and was well correlated with stem wood production (r(2) = 0.78). In contrast, APAR during the following years was less variable and not significantly correlated with stem biomass increments. Production of wood per unit of absorbed light varied with stand age and with site index. In our study, a better site index was accompanied both by increased APAR during the first 2 years of growth and by higher light use efficiency for stem wood production during the whole rotation. Implications for simple process-based modelling are discussed. (C) 2009 Elsevier B.V. All rights reserved.

Mapping changes in coastline geomorphic features using Landsat TM and ETM+ imagery: examples in southeastern Brazil

Multitemporal Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) imagery was used to assess coastline morphological changes in southeastern Brazil. A spectral linear mixing approach (SLMA) was used to estimate fraction imagery representing amounts of vegetation, clean water (a proxy for shade) and soil. Fraction abundances were related to erosive and depositional features. Shoreline, sandy banks (including emerged and submerged banks) and sand spits were highlighted mainly by clean water and soil fraction imagery. To evaluate changes in the coastline geomorphic features, the fraction imagery generated for each data set was classified in a contextual approach using a segmentation technique and ISOSEG, an unsupervised classification. Evaluation of the classifications was performed visually and by an error matrix relating ground-truth data to classification results. Comparison of the classification results revealed an intense transformation in the coastline, and that erosive and depositional features are extremely dynamic and subject to change in short periods of time.

Geographical distribution of American cutaneous leishmaniasis and its phlebotomine vectors (Diptera: Psychodidae) in the state of Sao Paulo, Brazil

Background: American cutaneous leishmaniasis (ACL) is a re-emerging disease in the state of Sao Paulo, Brazil. It is important to understand both the vector and disease distribution to help design control strategies. As an initial step in applying geographic information systems (GIS) and remote sensing (RS) tools to map disease-risk, the objectives of the present work were to: (i) produce a single database of species distributions of the sand fly vectors in the state of Sao Paulo, (ii) create combined distributional maps of both the incidence of ACL and its sand fly vectors, and (iii) thereby provide individual municipalities with a source of reference material for work carried out in their area. Results: A database containing 910 individual records of sand fly occurrence in the state of Sao Paulo, from 37 different sources, was compiled. These records date from between 1943 to 2009, and describe the presence of at least one of the six incriminated or suspected sand fly vector species in 183/645 (28.4%) municipalities. For the remaining 462 (71.6%) municipalities, we were unable to locate records of any of the six incriminated or suspected sand fly vector species (Nyssomyia intermedia, N. neivai, N. whitmani, Pintomyia fischeri, P. pessoai and Migonemyia migonei). The distribution of each of the six incriminated or suspected vector species of ACL in the state of Sao Paulo were individually mapped and overlaid on the incidence of ACL for the period 1993 to 1995 and 1998 to 2007. Overall, the maps reveal that the six sand fly vector species analyzed have unique and heterogeneous, although often overlapping, distributions. Several sand fly species - Nyssomyia intermedia and N. neivai - are highly localized, while the other sand fly species - N. whitmani, M. migonei, P. fischeri and P. pessoai - are much more broadly distributed. ACL has been reported in 160/183 (87.4%) of the municipalities with records for at least one of the six incriminated or suspected sand fly vector species, while there are no records of any of these sand fly species in 318/478 (66.5%) municipalities with ACL. Conclusions: The maps produced in this work provide basic data on the distribution of the six incriminated or suspected sand fly vectors of ACL in the state of Sao Paulo, and highlight the complex and geographically heterogeneous pattern of ACL transmission in the region. Further studies are required to clarify the role of each of the six suspected sand fly vector species in different regions of the state of Sao Paulo, especially in the majority of municipalities where ACL is present but sand fly vectors have not yet been identified.

The land-atmosphere water flux in the tropics

Tropical vegetation is a major source of global land surface evapotranspiration, and can thus play a major role in global hydrological cycles and global atmospheric circulation. Accurate prediction of tropical evapotranspiration is critical to our understanding of these processes under changing climate. We examined the controls on evapotranspiration in tropical vegetation at 21 pan-tropical eddy covariance sites, conducted a comprehensive and systematic evaluation of 13 evapotranspiration models at these sites, and assessed the ability to scale up model estimates of evapotranspiration for the test region of Amazonia. Net radiation was the strongest determinant of evapotranspiration (mean evaporative fraction was 0.72) and explained 87% of the variance in monthly evapotranspiration across the sites. Vapor pressure deficit was the strongest residual predictor (14%), followed by normalized difference vegetation index (9%), precipitation (6%) and wind speed (4%). The radiation-based evapotranspiration models performed best overall for three reasons: (1) the vegetation was largely decoupled from atmospheric turbulent transfer (calculated from X decoupling factor), especially at the wetter sites; (2) the resistance-based models were hindered by difficulty in consistently characterizing canopy (and stomatal) resistance in the highly diverse vegetation; (3) the temperature-based models inadequately captured the variability in tropical evapotranspiration. We evaluated the potential to predict regional evapotranspiration for one test region: Amazonia. We estimated an Amazonia-wide evapotranspiration of 1370 mm yr(-1), but this value is dependent on assumptions about energy balance closure for the tropical eddy covariance sites; a lower value (1096 mm yr(-1)) is considered in discussion on the use of flux data to validate and interpolate models.

MODIS aerosol optical depth Retrievals with high spatial resolution over an urban area using the critical reflectance

[1] The retrieval of aerosol optical depth (Ta) over land by satellite remote sensing is still a challenge when a high spatial resolution is required. This study presents a tool that uses satellite measurements to dynamically identify the aerosol optical model that best represents the optical properties of the aerosol present in the atmosphere. We use aerosol critical reflectance to identify the single scattering albedo of the aerosol layer. Two case studies show that the Sao Paulo region can have different aerosol properties and demonstrates how the dynamic methodology works to identify those differences to obtain a better T a retrieval. The methodology assigned the high single scattering albedo aerosol model (pi o( lambda = 0.55) = 0.90) to the case where the aerosol source was dominated by biomass burning and the lower pi(o) model (pi(o) (lambda = 0.55) = 0.85) to the case where the local urban aerosol had the dominant influence on the region, as expected. The dynamic methodology was applied using cloud-free data from 2002 to 2005 in order to retrieve Ta with Moderate Resolution Imaging Spectroradiometer ( MODIS). These results were compared with collocated data measured by AERONET in Sao Paulo. The comparison shows better results when the dynamic methodology using two aerosol optical models is applied (slope 1.06 +/- 0.08 offset 0.01 +/- 0.02 r(2) 0.6) than when a single and fixed aerosol model is used (slope 1.48 +/- 0.11 and offset - 0.03 +/- 0.03 r(2) 0.6). In conclusion the dynamical methodology is shown to work well with two aerosol models. Further studies are necessary to evaluate the methodology in other regions and under different conditions.

Multiscale Analysis of Topographic Surface Roughness in the Midland Valley, Scotland

Surface roughness is an important geomorphological variable which has been used in the Earth and planetary sciences to infer material properties, current/past processes, and the time elapsed since formation. No single definition exists; however, within the context of geomorphometry, we use surface roughness as an expression of the variability of a topographic surface at a given scale, where the scale of analysis is determined by the size of the landforms or geomorphic features of interest. Six techniques for the calculation of surface roughness were selected for an assessment of the parameter`s behavior at different spatial scales and data-set resolutions. Area ratio operated independently of scale, providing consistent results across spatial resolutions. Vector dispersion produced results with increasing roughness and homogenization of terrain at coarser resolutions and larger window sizes. Standard deviation of residual topography highlighted local features and did not detect regional relief. Standard deviation of elevation correctly identified breaks of slope and was good at detecting regional relief. Standard deviation of slope (SD(slope)) also correctly identified smooth sloping areas and breaks of slope, providing the best results for geomorphological analysis. Standard deviation of profile curvature identified the breaks of slope, although not as strongly as SD(slope), and it is sensitive to noise and spurious data. In general, SD(slope) offered good performance at a variety of scales, while the simplicity of calculation is perhaps its single greatest benefit.

Palaeodrainage on Marajo Island, northern Brazil, in relation to Holocene relative sea-level dynamics

A large area in northeastern Marajo Island, northern Brazil, has been characterized geomorphologically, applying information acquired from Landsat imagery. This study was combined with detailed sedimentologic analysis of continuous cores, which provided a record of depositional settings developed in this area through the Holocene. The results revealed well-preserved, meandering to anastomosed drainage networks of wide palaeochannels that were superimposed by a narrower palaeochannel system. In both cases, the sedimentary record consists of sands, heterolithic deposits and muds, locally rich in plant debris. The strata are organized into fining upward successions that reach approximately 18 m thick in the wide channels and 4 m thick in the narrow channels. Sedimentary features suggestive of a coastal location for the wider palaeochannels and reworking of sediments by tidal currents include the prevalence of well to moderately sorted, rounded to sub-rounded, fine- to medium-grained sands displaying foreset packages separated by mud couplets, suggestive of tidal cycles. The data presented herein point to a rise in relative sea level reaching the Lake Arari area during the early to late/mid Holocene. This event was followed by a relative sea level drop. Tectonics seem to have contributed to an overall lowering in relative sea level in the study area since the mid-Holocene, which does not follow the same pattern recorded in other areas along the northern Brazilian coast.

Quaternary tectonics in a passive margin: Marajo Island, northern Brazil

Marajo Island is located in a passive continental margin that evolved from rifting associated with the opening of the Equatorial South Atlantic Ocean in the Late Jurassic/Early Cretaceous period. This study, based on remote sensing integrated with sedimentology, as well as subsurface and seismographic data available from the literature, allows discussion of the significance of tectonics during the Quaternary history of marginal basins. Results show that eastern Marajo Island contains channels with evidence of tectonic control. Mapping of straight channels defined four main groups of lineaments (i.e. NNE-SSW, NE-SW, NW-SE and E-W) that parallel main normal and strike-slip fault zones recorded for the Amazon region. Additionally, sedimentological studies of late Quaternary and Holocene deposits indicate numerous ductile and brittle structures within stratigraphic horizons bounded by undeformed strata, related to seismogenic deformation during or shortly after sediment deposition. This conclusion is consistent with subsurface Bouguer mapping suggestive of eastern Marajo Island being still part of the Marajo graben system, where important fault reactivation is recorded up to the Quaternary. Together with the recognition of several phases of fault reactivation, these data suggest that faults developed in association with rift basins might remain active in passive margins, imposing important control on development of depositional systems. Copyright (C) 2007 John Wiley & Sons, Ltd.