Wavelet analysis of MODIS time series to detect expansion and intensification of row-crop agriculture in Brazil

Since 2000, the southwestern Brazilian Amazon has undergone a rapid transformation from natural vegetation and pastures to row-crop agricultural with the potential to affect regional biogeochemistry. The goals of this research are to assess wavelet algorithms applied to MODIS time series to determine expansion of row-crops and intensification of the number of crops grown. MODIS provides data from February 2000 to present, a period of agricultural expansion and intensification in the southwestern Brazilian Amazon. We have selected a study area near Comodoro, Mato Grosso because of the rapid growth of row-crop agriculture and availability of ground truth data of agricultural land-use history. We used a 90% power wavelet transform to create a wavelet-smoothed time series for five years of MODIS EVI data. From this wavelet-smoothed time series we determine characteristic phenology of single and double crops. We estimate that over 3200 km(2) were converted from native vegetation and pasture to row-crop agriculture from 2000 to 2005 in our study area encompassing 40,000 km(2). We observe an increase of 2000 km(2) of agricultural intensification, where areas of single crops were converted to double crops during the study period. (C) 2007 Elsevier Inc. All rights reserved.

Facies, delta(13)C, delta(15)N and C/N analyses in a late Quaternary compound estuarine fill, northern Brazil and relation to sea level

The present work integrates sedimentary facies, (14)C dating, delta(13)C, delta(15)N, and C/N with geologic and geomorphologic data available from literature. The aim was to characterize the depositional settings of a late Quaternary estuary in northeastern Marajo Island and analyze its evolution within the context of relative sea level fluctuations. The data derive from four continuous cores along a proximal-to-distal transect of a paleoestuary, previously recognized using remote sensing information. Fifteen sediment samples recorded ages ranging from 42,580 +/- 1430 to 3184 +/- 37 (14)C yr B.P. Fades analysis indicated fine- to coarse-grained sands with parallel lamination or cross stratification, massive or laminated muds and heterolithic deposits. delta(13)C (-28.1 parts per thousand to -19.7 parts per thousand, mean = -23.0 parts per thousand), delta(15)N (+ 14.8 parts per thousand to + 4.7 parts per thousand, mean = + 9.2 parts per thousand) and C/N (14.5 to 1.5, mean = 7.9) indicate mostly marine and freshwater phytoplankton sources for the organic matter. The results confirm a large late Quaternary paleoestuary in northeastern Marajo Island. The distribution of delta(13)C, delta(15)N, and C/N, together with fades associations, led to identify depositional settings related to fluvial channel, floodplain, tidal channel/tidal flat, central basin, tidal delta, and tidal inlet/sand barrier. These deposits are consistent with a wave-dominated estuary. Variations in stratigraphy and geochemistry are controlled by changes in relative sea level, revealing a main transgression from an undetermined time around 42,000 (14)C yr B.P. and 29,340 (+/- 200) (14)C yr B.P., which is synchronous to the overall drop in sea level after the last interglacial. Following this period, and probably until 9110 +/- 37 (14)C yr B.P., i.e., during a time interval encompassing two glacial episodes including the Last Glacial and the Younger Dryas, there was a pronounced drop in sea level, recorded by the development of a major erosional discontinuity due to valley re-incision. Sea level rose again until 5464 +/- 40 (14)C yr B.P, just before the main worldwide mid-Holocene transgressive peak. Mid to late Holocene coastal progradation ended the Marajo paleoestuarine history, and promoted the establishment of continental conditions throughout the island. The divergence comparing the Marajo sea level behavior with the eustatic curve allows hypothesizing that post-rifting tectonics along the Brazilian Equatorial margin influenced the sedimentary evolution of the studied paleoestuary. Considering that sedimentary facies in estuarine settings are highly variable both laterally and vertically, the present integration of facies with isotope and elemental analyses was crucial to provide a more precise interpretation of the Late Pleistocene and Holocene Marajo paleoestuary, and analyze its sea level history within the eustatic and tectonic context. (C) 2010 Elsevier B.V. All rights reserved.

Soil density evaluated by spectral reflectance as an evidence of compaction effects

Soil compaction, reflected by high bulk density, is an environmental degradation process and new technologies are being developed for its detection. Despite the proven efficiency of remote sensing, it has not been widely used for soil density. Our objective was to evaluate the density of two soils: a Typic Quartzpisament (TQ) and a Rhodic Paleudalf (RP), using spectral reflectance obtained by a laboratory spectroradiometer between 450 and 2500 nm. Undisturbed samples were taken at two depths (0-20 and 60-80 cm), and were artificially compacted. Spectral data, obtained before and after compaction, were compared for both wet and dried compacted samples. Results demonstrated that soil density was greater in RP than in TQ at both depths due to its clayey texture. Spectral data detected high density (compacted) from low density (non-compacted) clayey soils under both wet and dry conditions. The detection of density in sandy soils by spectral reflectance was not possible. The intensity of spectral reflectance of high soil bulk density (compacted) samples was higher than for low density (non-compacted) soils due to changes in soil structure and porosity. Dry samples with high bulk density showed differences in the spectral intensity, but not in the absorption features. Wet samples in equal condition had statistically higher reflectance intensity than that of the low soil bulk density (non-compacted), and absorption differences at 1920 nm, which was due to the altered position of the water molecules. Soil line and spectral reflectance used together could detect soil bulk density variations for the clay soil. This technique could assist in the detection of high soil density in the laboratory by providing new soil information.

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.

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.