Electrical structure of the upper mantle beneath Central Europe: Results of the CEMES project

In the years 2001–2003, we accomplished the experimental phase of the project CEMES by collecting long-period magnetotelluric data at positions of eleven permanent geomagnetic observatories situated within few hundreds kilometers along the south-west margin of the East European Craton. Five teams were engaged in estimating independently the magnetotelluric responses by using different data processing procedures. The conductance distributions at the depths of the upper mantle have been derived individually beneath each observatory. By averaging the individual cross-sections, we have designed the final model of the geoelectrical structure of the upper mantle beneath the CEMES region. The results indicate systematic trends in the deep electrical structure of the two European tectonic plates and give evidence that the electrical structure of the upper mantle differs between the East European Craton and the Phanerozoic plate of west Europe, with a separating transition zone that generally coincides with the Trans-European Suture Zone.     

Estimation of phosphorus loads with sparse data for agricultural watersheds in the Czech Republic

Abstract

Agricultural watersheds in the Czech Republic are one of the primary sources of non-point-source phosphorus (P) loads in receiving waters. Since such non-point sources are generally located in headwater catchments, streamflow and P concentration data are sparse. We show how very short daily streamflow and P concentration records can be combined with nearby longer existing daily streamflow records to result in reliable estimates of daily and annual P concentrations and loads. Maintenance of variance streamflow record extension methods (MOVE) can be employed to extend short streamflow records. Constituent load regressions are used to predict daily P constituent loads from streamflow and other time varying characteristics. Annual P loads are then estimated for individual watersheds. Resulting annual P load estimates ranged from 0.21 to 95.4 kg year^-1 with a mean value of 11.77 kg year^-1. Similarly annual P yield estimates ranged from 0.01 to 0.3 kg ha^-1 year^-1 with an average yield of 0.07 kg ha^-1 year^-1. We document how short records of daily streamflow and P concentrations can be combined with a national network of daily streamflow records in the Czech Republic to arrive at meaningful and reliable estimates of annual P loads for small agricultural watersheds.

Citation Beránková, T., Vogel, R. M., Fiala, D. & Rosendorf, P. (2010) Estimation of phosphorus loads with sparse data for agricultural watersheds in the Czech Republic. Hydrol. Sci. J. 55(8), 1417–1426.     

Gravimetric monitoring of the first field‐wide steam injection in a fractured carbonate field in Oman – a feasibility study

Gas‐Oil Gravity Drainage is to be enhanced by steam injection in a highly fractured, low permeability carbonate field in Oman. Following a successful pilot, field‐wide steam injection is being implemented, first of this type in the world. A dedicated monitoring program has been designed to track changes in the reservoir. Various observations are to be acquired, including, surface deformation, temperature measurements, microseismic, well logs, pressure and saturation measurements to monitor the reservoir. Because significant changes in the reservoir density are expected, time‐lapse gravimetry is also being considered. In this paper we investigate the feasibility of gravimetric monitoring of the thermally enhanced gravity drainage process at the carbonate field in Oman. For this purpose, forward gravity modelling is performed. Based on field groundwater measurements, the estimates of the hydrological signal are considered and it is investigated under what conditions the groundwater influences can be minimized. Using regularized inversion of synthetic gravity data, we analyse the achievable accuracy of heat‐front position estimates. In case of large groundwater variations at the field, the gravity observations can be significantly affected and, consequently, the accuracy of heat‐front monitoring can be deteriorated. We show that, by applying gravity corrections based on local observations of groundwater, the hydrological influences can to a large extent be reduced and the accuracy of estimates can be improved. We conclude that gravimetric monitoring of the heat‐front evolution has a great potential.     

Induction Sounding of the Earth’s Mantle at a New Russian Geophysical Observatory

Deep magnetotelluric (MT) sounding data were collected and processed in the western part of the East European Craton (EEC). The MT sounding results correspond well with impedances obtained by magnetovariation (MV) sounding on the new geophysical observatory situated not far from the western border of Russia. Inversion based on combined data of both induction soundings let us evaluate geoelectrical structure of the Earth’s crust and upper and mid-mantle at depths up to 2000 km, taking into account the harmonics of 11-year variations. Results obtained by different authors and methods are compared with similar investigations on the EEC such as international projects CEMES in Central Europe and BEAR in Fennoscandia.     

On the accuracy of the bathymetry-generated gravitational field quantities for a depth-dependent seawater density distribution

In geophysical studies investigating the lithosphere structure, the gravitational field generated by the ocean density contrast (i.e., bathymetry-generated gravitational field) represents a significant amount of the signal to be modelled and subsequently removed from the Earth’s gravity field. The ocean density contrast is typically calculated as the difference between the mean density values of the Earth’s crust and seawater. The approximation of the actual seawater density distribution by its mean value yields relative errors up to about 2% in computed quantities of the gravitational field. To reduce these errors, a more realistic model of the seawater density distribution is utilized based on the analysis of existing oceanographic data of salinity, temperature, and pressure (depth). We study the accuracy of the bathymetry-generated gravitational field quantities formulated for a depth-dependent model of the seawater density distribution. This density distribution approximates the seawater density variations due to an increasing pressure with depth, whereas smaller lateral density variations caused by salinity, temperature, and other oceanographic factors are not taken into consideration. The error analysis reveals that the approximation of the seawater density by the depth-dependent density model reduces the maximum errors to less than 0.6%. The corresponding depth-averaged errors are below 0.1%. The depth-dependent seawater density model is further facilitated in expressions for computing the bathymetry-generated gravitational field quantities by means of the spherical bathymetric (ocean bottom depth) functions. The numerical realization reveals large differences in the results obtained with and without consideration of the depth-dependent seawater density distribution. The maxima of absolute differences reach 201 m²/s² and 16.5 mGal in computed values of the potential and attraction, respectively. The application of the depth-dependent seawater density model thus significantly improves the accuracy in the forward modelling of the bathymetric gravitational field quantities.     

Quantitative analysis of travel-time curves of seismic events in the Ostrava-Karviná coal mining district

Summary The paper is intended as a contribution to the quantitative analysis of travel-time curves of seismic events recorded in the Ostrava-Karviná District (OKD). The input data represent a set of 2621 seismic events, recorded by the local seismological network of 26 mine stations DSLA and a regional diagnostic polygon consisting of five surface Lennartz stations. All the events were processed automatically in the Operational Seismological Centre of the Czechoslovak Army Mine in Karviná and stored in the seismological data base. The results are presented in the form of graphs of arrival times versus distance for the whole OKD, for two mines and one tectonic block.Travel-time curves of direct P and S waves, as well as of reflected and refracted waves are given. The direct P and S waves propagate well practically throughout the whole region studied, but their apparent velocities of propagation are affected by the properties of the rock medium.As a result of the complicated geological conditions, the recorded wave image is quite complicated. Methods of mathematical modelling, using kinematic and dynamic parameters of seismic waves, will have to be applied to identify the separate wave groups uniquely.     

Optimal Model for Geoid Determination from Airborne Gravity

Two different approaches for transformation of airborne gravity disturbances, derived from gravity observations at low-elevation flying platforms, into geoidal undulations are formulated, tested and discussed in this contribution. Their mathematical models are based on Green's integral equations. They are in these two approaches defined at two different levels and also applied in a mutually reversed order. While one of these approaches corresponds to the classical method commonly applied in processing of ground gravity data, the other approach represents a new method for processing of gravity data in geoid determination that is unique to airborne gravimetry. Although theoretically equivalent in the continuous sense, both approaches are tested numerically for possible numerical advantages, especially due to the inverse of discretized Fredholm's integral equation of the first kind applied on different data. High-frequency synthetic gravity data burdened by the 2-mGal random noise, that are expected from current airborne gravity systems, are used for numerical testing. The results show that both approaches can deliver for the given data a comparable cm-level accuracy of the geoidal undulations. The new approach has, however, significantly higher computational efficiency. It would be thus recommended for real life geoid computations. Additional errors related to regularization of gravity data and the geoid, and to accuracy of the reference field, that would further deteriorate the quality of estimated geoidal undulations, are not considered in this study.