The role of prior model calibration on predictions with ensemble Kalman filter

This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.     

Rock salt mass in the Paleoproterozoic sequence of the Onega trough in Karelia (from the Onega parametric well data)

A thick (200 m) rock salt mass covering Late Archean granitoids was exposed for the first time in the Early Proterozoic volcanogenic-sedimentary formations in the Onega trough of the east Baltic Shield by the Onega parametric well. The mineral composition of salts, their geochemical features, and the isotopic composition of carbonate carbon and oxygen have been studied. After fluid inclusions present in salts, their metamorphism temperature and isotopic composition of helium and argon were determined. The obtained results give evidence of the fact that rock salts and magnesites associated with them were formed in an evaporate basin with participation of deep crustal processes. The age of the underlying granitoids (2.716 ± 9 Ma) is determined using the Pb—Pb method.     

Observation of Precipitating Systems over Complex Orography with Meteorological Doppler Radars: A Feasibility Study

Summary  This paper concerns the use of airborne or ground-based Doppler radars to observe precipitating systems over complex orography. As nearly all of the previous experiments involving Doppler radars were conducted over flat surfaces over the continents or the oceans, new techniques are needed firstly to separate ground clutter from meteorological signal and, in the case of airborne Doppler observations, to deduce navigational errors. Secondly, it is necessary to take the atmospheric circulation induced by orography into account in the three-dimensional wind field analysis. Variational techniques are presented to solve these problems. The proposed methods are tested with simulated ground-based and airborne Doppler radar observations for analytic flows over analytic terrains and for numerically simulated wind and reflectivity fields for the Brig event (22 September 1993) of heavy precipitation over the southern flank of the Alps (Cosma and Richard, 1998), and with actual airborne Doppler data relative to weak snow showers over the Rocky Mountains on 12 March 1995. Received March 22, 1999/Revised June 1, 1999     

Photoelectric device for monitoring the movement of a torsion balance beam

Manufactured and tested in his work is a photoelectric device that makes it possible to automatically determine the relative angular position of the beam of a miniature torsion balance, angular readings being discretized with various time step values and stored in a computer. Measurements were taken with this device that captured the occurrence of the partial solar eclipse on September 11, 2007, and confirmed the fact, established previously through visual observation, that the position of torsion balance beams changed at the instant an eclipse took place.     

Spectroscopic binaries near the North Galactic Pole paper 6: BD 33° 2206

Photoelectric radial-velocity measurements have confirmed O. C. Wilson’s finding that BD 33° 2206, the secondary star in the wide visual binary ADS 8470, is a spectroscopic binary. It has an eccentric orbit with a period of 100 days. Its γ-velocity is close to the constant radial velocity of the visual primary, confirming the physical association of the stars.     

Synchrotron spectra of short-period pulsars

A model with synchrotron radiation near the light cylinder is proposed to explain the observed spectra of short-period pulsars (P≤0.1 s). These spectra can be described if a power-law energy distribution of the emitting electrons with exponent γ=2–8 is assumed. For most pulsars, the peak frequency ν_m is below 10 MHz. The ν_m(γ) dependence is derived, and shows that the peak frequencies for pulsars with spectral indices α<1.5 may fall in the observable range. In particular, ν_m may be ν_m ～ 100 MHz for PSR J0751 + 1807 and PSR J1640 + 2224. The observed radio spectrum of Geminga (PSR J0633 + 1746) can be described by a synchrotron model with a monoenergetic or Maxwellian distribution of relativistic electrons and a small angle β between the spin axis and magnetic moment (β ～ 10°).     

Regularized derivatives of potential fields and their role in semi-automated interpretation methods

Evaluation of higher derivatives (gradients) of potential fields plays an important role in geophysical interpretation (qualitative and/or quantitative), as has been demonstrated in many approaches and methods. On the other hand, numerical evaluation of higher derivatives is an unstable process – it has the tendency to enlarge the noise content in the original data (to degrade the signal-to-noise ratio). One way to stabilize higher derivative evaluation is the utilization of the Tikhonov regularization. In the submitted contribution we present the derivation of the regularized derivative filter in the Fourier domain as a minimization task by means of using the classical calculus of variations. A very important part of the presented approach is the selection of the optimum regularization parameter – we are using the analysis of the C-norm function (constructed from the difference between two adjacent solutions, obtained for different values of regularization parameter). We show the influence of regularized derivatives on the properties of the classical 3D Euler deconvolution algorithm and apply it to high-sensitivity magnetometry data obtained from an unexploded ordnance detection survey. The solution obtained with regularized derivatives gives better focused depth-estimates, which are closer to the real position of sources (verified by excavation of unexploded projectiles).     

Assessment of regional air pollution distribution by point cumulative semivariogram method at Erzurum urban center, TURKEY

In the assessment of air quality, regional distribution and dispersion with distance are important, together with the variations of pollutants in time. On this occasion, the point cumulative semi-variogram (PCSV) method is used in order to find simply regional distribution of pollutants of Erzurum urban centre. This method is based simply on the summation of square differences in air pollutant concentrations between different sites. Monthly regional variation maps of Erzurum are constructed by finding radius of influence (for SO₂, from 1000 m to 3500 m and, for TSP, 1000–2000 m) and PCSV scattering diagram data at different levels by using monthly average sulphur dioxide (SO₂) and total suspended particulate (TSP) matter concentrations in 2001–2002 winter season. Consequently, the air pollution distribution of Erzurum is assessed.     

Thermal evolution of an extensional detachment as constrained by organic metamorphic data and thermal modeling: Graz Paleozoic Nappe Complex (Eastern Alps)

Following Early Cretaceous nappe stacking, the Eastern Alps were affected by late-orogenic extension during the Late Cretaceous. In the eastern segment of this range, a Late Cretaceous detachment separates a very low- to low-grade metamorphic cover (Graz Paleozoic Nappe Complex, GPNC) above a low- to high-grade metamorphic basement. Synchronously, the Kainach Gosau Basin (KGB) collapsed and subsided on top of the section.Metamorphism of organic material within this section has been investigated using vitrinite reflectance data and Raman spectra of extracted carbonaceous material. In the southern part of the GPNC, vitrinite reflectance indicates a decrease in organic maturity towards the stratigraphic youngest unit. The remaining part of the GPNC is characterized by an aureole of elevated vitrinite reflectance values and Raman R2 ratios that parallels the margins of the GPNC. Vitrinite reflectance in the KGB shows a steep coalification gradient and increases significantly towards the western basin margin. The observed stratigraphic trend in the southern GPNC is a result of deep Paleozoic to Early Cretaceous burial. This maturity pattern was overprinted along the margins by advective heat and convective fluids during Late Cretaceous to Paleogene exhumation of basement rocks.During shearing, the fault zone was heated up to ca. 500 °C. This overprint is explained by a two-dimensional thermal model with a ramp-flat fault geometry and a slip rate of 1 to 1.5 cm/year during 5 Ma fault movement. The collapse basin above the detachment subsided in a thermal regime which was characterized by relaxing isotherms.