Review of resource analysis and related studies: Working group,Mathematical Geology,freie Universität Berlin

The working group Mathematical Geology of the Freie Universität Berlin was formed in 1971. We review quantitative methods used by the working group since 1983 to treat mineral exploration problems. The methods applied vary from elementary statistical analysis of multivariate exploration data to optimal strategies for selecting favorable targets, and from multiobjective decision-making for additional drill hole locations to expert systems in exploration.The methods applied are directly related to the level of information at each stage of the exploration process. Special emphasis was placed on the relationship between and evaluation of subjective and objective data. Case studies illustrating the various methods are presented for different kinds of mineral deposits and exploration environments.     

Applications of AnnAGNPS model for soil loss estimation and nutrient loading for Malaysian conditions

The study evaluated the performance and suitability of AnnAGNPS model in assessing runoff, sediment loading and nutrient loading under Malaysian conditions. The watershed of River Kuala Tasik in Malaysia, a combination of two sub-watersheds, was selected as the area of study. The data for the year 2004 was used to calibrate the model and the data for the year 2005 was used for validation purposes. Several input parameters were computed using methods suggested by other researchers and studies carried out in Malaysia. The study shows that runoff was predicted well with an overall R² value of 0.90 and E value of 0.70. Sediment loading was able to produce a moderate result of R² = 0.66 and E = 0.49, nitrogen loading predictions were slightly better with R² = 0.68 and E = 0.53, and phosphorus loading performance was slightly poor with an R² = 0.63 and E = 0.33. The erosion map developed was in agreement with the erosion risk map produced by the Department of Agriculture, Malaysia. Rubber estates and urban areas were found to be the main contributors to soil erosion. The simulation results showed that AnnAGNPS has the potential to be used as a valuable tool for planning and management of watersheds under Malaysian conditions.     

Using spectral information from the NIR water absorption features for the retrieval of canopy water content

Canopy water content (CWC) is important for mapping and monitoring the condition of the terrestrial ecosystem. Spectral information related to the water absorption features at 970 nm and 1200 nm offers possibilities for deriving information on CWC. In this study, we compare the use of derivative spectra, spectral indices and continuum removal techniques for these regions. Hyperspectral reflectance data representing a range of canopies were simulated using the combined PROSPECT + SAILH model. Best results in estimating CWC were obtained by using spectral derivatives at the slopes of the 970 nm and 1200 nm water absorption features. Real data from two different test sites were analysed. Spectral information at both test sites was obtained with an ASD FieldSpec spectrometer, whereas at the second site HyMap airborne imaging spectrometer data were also acquired. Best results were obtained for the derivative spectra. In order to avoid the potential influence of atmospheric water vapour absorption bands the derivative of the reflectance on the right slope of the canopy water absorption feature at 970 nm can best be used for estimating CWC.     

The use of MODIS data to derive acreage estimations for larger fields: A case study in the south-western Rostov region of Russia

Recent developments in remote sensing technology, in particular improved spatial and temporal resolution, open new possibilities for estimating crop acreage over larger areas. Remotely sensed data allow in some cases the estimation of crop acreage statistics independently of sub-national survey statistics, which are sometimes biased and incomplete. This work focuses on the use of MODIS data acquired in 2001/2002 over the Rostov Oblast in Russia, by the Azov Sea. The region is characterised by large agricultural fields of around 75 ha on average. This paper presents a methodology to estimate crop acreage using the MODIS 16-day composite NDVI product. Particular emphasis is placed on a good quality crop mask and a good quality validation dataset. In order to have a second dataset which can be used for cross-checking the MODIS classification a Landsat ETM time series for four different dates in the season of 2002 was acquired and classified. We attempted to distinguish five different crop types and achieved satisfactory and good results for winter crops. Three hundred and sixty fields were identified to be suitable for the training and validation of the MODIS classification using a maximum likelihood classification. A novel method based on a pure pixel field sampling is introduced. This novel method is compared with the traditional hard classification of mixed pixels and was found to be superior.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

GPS observations of the ionospheric F2-layer behavior during the 20th November 2003 geomagnetic storm over South Korea

The ionospheric F2-layer peak density (NmF2) and its height (hmF2) are of great influence on the shape of the ionospheric electron density profile Ne (h) and may be indicative of other physical processes within the ionosphere, especially those due to geomagnetic storms. Such parameters are often estimated using models such as the semiempirical international reference ionosphere (IRI) models or are measured using moderately priced to expensive instrumentation, such as ionosondes or incoherent scatter radars. Global positioning system (GPS) observations have become a powerful tool for mapping high-resolution ionospheric structures, which can be used to study the ionospheric response to geomagnetic storms. In this paper, we describe how 3-D ionospheric electron density profiles were produced from data of the dense permanent Korean GPS network using the tomography reconstruction technique. These profiles are verified by independent ionosonde data. The responses of GPS-derived parameters at the ionospheric F2-layer to the 20th November 2003 geomagnetic storm over South Korea are investigated. A fairly large increase in the electron density at the F2-layer peak (the NmF2) (positive storm) has been observed during this storm, which is accompanied by a significant uplift in the height of the F2 layer peak (the hmF2). This is confirmed by independent ionosonde observations. We suggest that the F2-layer peak height uplift and NmF2 increase are mainly associated with a strong eastward electric field, and are not associated with the increase of the O/N₂ ratio obtained from the GUVI instruments aboard the TIMED satellite. It is also inferred that the increase in NmF2 is not caused by the changes in neutral composition, but is related to other nonchemical effects, such as dynamical changes of vertical ion motions induced by winds and E × B drifts, tides and waves in the mesosphere/lower thermosphere region, which can be dynamically coupled upward to generate ionospheric perturbations and oscillations.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.