The Global Impact of Clouds on the Production of MODIS Bidirectional Reflectance Model-Based Composites for Terrestrial Monitoring

A global data set of cloud occurrence probability derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua gridded daily data is analyzed to investigate the probability of obtaining at least a minimum number of cloud-free observations within various compositing periods. The probabilities derived from Terra and Aqua, with morning and afternoon overpass times, respectively, are similar and increase with compositing period. Compositing both Terra and Aqua observations results in considerably higher probabilities of obtaining a sufficient number of observations for bidirectional reflectance model-based compositing. Given that the only alternative to obtaining sufficient samples is to extend the observation period, which can cause significant problems when the surface state changes, it is concluded that using data from the two MODIS sensors provides the most effective way of generating composited products. Findings with respect to the availability of cloud-free composites when n-day composites are generated on a temporally overlapping daily rolling basis, i.e., every day, rather than every n-days, are also discussed for regional and global applications     

Supervised SAR Image MPM Segmentation Based on Region-Based Hierarchical Model

This letter presents a novel method of supervised multiresolution segmentation for synthetic aperture radar images. The method uses a region-based half-tree hierarchical Markov random field model for multiresolution segmentation. To form the region-based multilayer model, the watershed algorithm is employed at each resolution level independently. The nodes of a quadtree in the proposed model are defined as regions instead of pixels. The relationship over scale is studied, and the region-based upward and downward maximization of posterior marginal estimations are deduced. The experimental results for the segmentation of homogeneous areas prove the region-based model much better in terms of robustness to speckle and preservation of edges compared to the pixel-based hierarchical model and the Gibbs sampler with the single-resolution model     

Geomasking sensitive health data and privacy protection: an evaluation using an E911 database

Geomasking is used to provide privacy protection for individual address information while maintaining spatial resolution for mapping purposes. Donut geomasking and other random perturbation geomasking algorithms rely on the assumption of a homogeneously distributed population to calculate displacement distances, leading to possible under-protection of individuals when this condition is not met. Using household data from 2007, we evaluated the performance of donut geomasking in Orange County, North Carolina. We calculated the estimated k-anonymity for every household based on the assumption of uniform household distribution. We then determined the actual k-anonymity by revealing household locations contained in the county E911 database. Census block groups in mixed-use areas with high population distribution heterogeneity were the most likely to have privacy protection below selected criteria. For heterogeneous populations, we suggest tripling the minimum displacement area in the donut to protect privacy with a less than 1% error rate.     

Soil erosion modelled with USLE and PESERA using QuickBird derived vegetation parameters in an alpine catchment

The focus of soil erosion research in the Alps has been in two categories: (i) on-site measurements, which are rather small scale point measurements on selected plots often constrained to irrigation experiments or (ii) off-site quantification of sediment delivery at the outlet of the catchment. Results of both categories pointed towards the importance of an intact vegetation cover to prevent soil loss. With the recent availability of high-resolution satellites such as IKONOS and QuickBird options for detecting and monitoring vegetation parameters in heterogeneous terrain have increased. The aim of this study is to evaluate the usefulness of QuickBird derived vegetation parameters in soil erosion models for alpine sites by comparison to Cesium-137 (Cs-137) derived soil erosion estimates. The study site (67 km²) is located in the Central Swiss Alps (Urseren Valley) and is characterised by scarce forest cover and strong anthropogenic influences due to grassland farming for centuries. A fractional vegetation cover (FVC) map for grassland and detailed land-cover maps are available from linear spectral unmixing and supervised classification of QuickBird imagery. The maps were introduced to the Pan-European Soil Erosion Risk Assessment (PESERA) model as well as to the Universal Soil Loss Equation (USLE). Regarding the latter model, the FVC was indirectly incorporated by adapting the C factor. Both models show an increase in absolute soil erosion values when FVC is considered. In contrast to USLE and the Cs-137 soil erosion rates, PESERA estimates are low. For the USLE model also the spatial patterns improved and showed “hotspots” of high erosion of up to 16 t ha⁻¹ a⁻¹. In conclusion field measurements of Cs-137 confirmed the improvement of soil erosion estimates using the satellite-derived vegetation data.     

Combining EGM2008 and SRTM/DTM2006.0 residual terrain model data to improve quasigeoid computations in mountainous areas devoid of gravity data

A global geopotential model, like EGM2008, is not capable of representing the high-frequency components of Earth’s gravity field. This is known as the omission error. In mountainous terrain, omission errors in EGM2008, even when expanded to degree 2,190, may reach amplitudes of 10 cm and more for height anomalies. The present paper proposes the utilisation of high-resolution residual terrain model (RTM) data for computing estimates of the omission error in rugged terrain. RTM elevations may be constructed as the difference between the SRTM (Shuttle Radar Topography Mission) elevation model and the DTM2006.0 spherical harmonic topographic expansion. Numerical tests, carried out in the German Alps with a precise gravimetric quasigeoid model (GCG05) and GPS/levelling data as references, demonstrate that RTM-based omission error estimates improve EGM2008 height anomaly differences by 10 cm in many cases. The comparisons of EGM2008-only height anomalies and the GCG05 model showed 3.7 cm standard deviation after a bias-fit. Applying RTM omission error estimates to EGM2008 reduces the standard deviation to 1.9 cm which equates to a significant improvement rate of 47%. Using GPS/levelling data strongly corroborates these findings with an improvement rate of 49%. The proposed RTM approach may be of practical value to improve quasigeoid determination in mountainous areas without sufficient regional gravity data coverage, e.g., in parts of Asia, South America or Africa. As a further application, RTM omission error estimates will allow refined validation of global gravity field models like EGM2008 from GPS/levelling data.     

Progressive Transformation of a Flux Rope to an ICME

The solar wind conditions at one astronomical unit (AU) can be strongly disturbed by interplanetary coronal mass ejections (ICMEs). A subset, called magnetic clouds (MCs), is formed by twisted flux ropes that transport an important amount of magnetic flux and helicity, which is released in CMEs. At 1 AU from the Sun, the magnetic structure of MCs is generally modeled by neglecting their expansion during the spacecraft crossing. However, in some cases, MCs present a significant expansion. We present here an analysis of the huge and significantly expanding MC observed by the Wind spacecraft during 9 – 10 November 2004. This MC was embedded in an ICME. After determining an approximate orientation for the flux rope using the minimum variance method, we obtain a precise orientation of the cloud axis by relating its front and rear magnetic discontinuities using a direct method. This method takes into account the conservation of the azimuthal magnetic flux between the inbound and outbound branches and is valid for a finite impact parameter (i.e., not necessarily a small distance between the spacecraft trajectory and the cloud axis). The MC is also studied using dynamic models with isotropic expansion. We have found (6.2±1.5)×10²⁰ Mx for the axial flux and (78±18)×10²⁰ Mx for the azimuthal flux. Moreover, using the direct method, we find that the ICME is formed by a flux rope (MC) followed by an extended coherent magnetic region. These observations are interpreted by considering the existence of a previously larger flux rope, which partially reconnected with its environment in the front. We estimate that the reconnection process started close to the Sun. These findings imply that the ejected flux rope is progressively peeled by reconnection and transformed to the observed ICME (with a remnant flux rope in the front part).     

Geochronology of anthropogenic pollutants in riparian wetland sediments of the Lippe River (Germany)

Anthropogenic pollutants were determined in a dated sediment core, collected from a riparian wetland of the Lippe River (Ruhr district, Germany). The historical trend in heavy metals, polycyclic aromatic hydrocarbons, organochlorines, polychlorinated biphenyls and linear alkylbenzenes as well as more recent contaminants such as industrial additives, organotins, synthetic musks, methyltriclosan and some other compounds were determined for the time period between 1930 and 1986. Emission sources, information on technical production and usage, as well as on the individual pollution pathways, with appropriate environmental stability, were considered in the interpretation of the sediment contamination over the past 50 years.Contaminants were analysed and interpreted according to two different criteria: (a) the origin of the main contaminants as related to mining and industrial activities as well as municipal sewage and agricultural effluents and (b) the pollution history. Due to a significant appearance of formerly missing contaminants in sediments deposited since 1970, we suggest classifying contaminants as either common (predating 1970) or modern (postdating 1970).In summary, the study provided a comprehensive reconstruction of the pollution history of the Lippe River system.     

Spatial and temporal patterns of dry spell lengths in the Iberian Peninsula for the second half of the twentieth century

Summary Daily pluviometric records of 43 meteorological stations across the Iberian Peninsula have permitted a detailed analysis of dry spell patterns for the period 1951–2000 by distinguishing daily amount thresholds of 0.1, 1.0, 5.0 and 10.0 mm/day. The analyses are based on three annual series, namely the number of dry spells, N, the average dry spell length, L, and the extreme dry spell length, L _max. First, the statistical significance of local trends for the annual series of N, L and L _max has been investigated by means of the Mann-Kendall test and significant field trends have been established by means of Monte Carlo simulations. Clear signs of negative field trends are detected for N (1.0 and 10.0 mm/day) and L (0.1 mm/day). Second, the Weibull model fits well the empirical distributions of dry spell lengths for all the rain gauges, whatever the daily amount threshold, with a well ranged spatial distribution of their parameters u and k. On the basis of the Weibull distribution, return period maps for 2, 5, 10, 25 and 50 years have been obtained for dry spell lengths with respect to the four daily threshold levels. While for 0.1 and 1.0 mm/day the longest dry spells are expected at the south of the Iberian Peninsula, for 5.0 and 10.0 mm/day they are mostly detected at the southeast. Finally, the elapsed time between consecutive dry spells has been analysed by considering the same rain amount thresholds and different dry spell lengths at increasing intervals of 10 days. This analysis makes evident a significant negative field trend of the elapsed time between consecutive dry spells of lengths ranging from 10 to 20 days for daily amount thresholds of 1.0, 5.0 and 10.0 mm/day. Authors’ addresses: X. Lana, C. Serra, Departament de Física i Enginyeria Nuclear, ETSEIB, Universitat Politècnica de Catalunya, Av. Diagonal 647 planta 11, 08028 Barcelona, Spain; M. D. Marínez, Departament de Física Aplicada, Universitat Politècnica de Catalunya, 08028 Barcelona, Spain; A. Burgue?o, Departament de Meteorologia i Astronomia, Universitat de Barcelona, 08028 Barcelona, Spain; J. Martín-Vide, L. Gómez, Grup de Climatologia, Universitat de Barcelona, 08028 Barcelona, Spain.     

Comparison of Vegetation Indices from Two Ground Based Sensors

The best and commonly used ground-based sensor to monitor crop growth, ASD FieldSpecPro Spectroradiometer (Analytical Spectral Devices, Boulder, CO, USA) is a passive sensor, which can be used under adequate light condition. However, now-a-days active sensors such as GreenSeeker? (GS) handheld crop response (Trimble Agriculture division, USA) are used for monitoring crop growth and are flexible in terms of timeliness and illumination conditions besides being cheaper than the ASD. Before its wide use, the suitability and accuracy of GS should be assessed by comparing the NDVI measured by this instrument with that by ASD, under diverse wheat growing conditions of India. Keeping this in view, the present experiment was undertaken with the following objectives: (1) to find out the temporal variation of NDVI measured both by ASD and GS treatments, (2) to find out relationship between the NDVI measured through ASD and GS and, (3) to evaluate the suitability of GS for NDVI measurements. It was observed that the numerical value of NDVI as measured by GS was always significantly (P < 0.05) lower than that measured by ASD for all the experiments under study. The NDVI-ASD and NDVI-GS were significantly positively correlated (P < 0.01) with the correlation coefficients being +0.94, +0.88 and +0.87 for irrigation and nitrogen experiment, irrigation and cultivars experiment, and tillage, residue and nitrogen experiments, respectively. Further, the regression equation developed between the NDVI-ASD and NDVI-GS: [NDVI-GS = 1.070 × (NDVI-ASD ? 0.292] can be successfully used to compute the NDVI of ASD from that computed by GS.