On fine orbit selection for particular geodetic and oceanographic missions involving passage through resonances

 The identification of mean semi-major axes (suitably defined) for satellite orbits to satisfy a variety of requirements for geodesy, geophysics and oceanography, in terms of repeat orbits (with orbital resonances), is investigated. Various options for the definition of semi-major axis, from the viewpoint of satellite dynamics, are described. Simple simulations of the expected resonant changes in inclination are presented, and tools for the analysis of orbit resonances to extract certain lumped harmonic coefficients of the geopotential (e.g. from the very precise CHAMP orbit) are resurrected. Finally, a preliminary example of the 46th-order resonance analysis possible for CHAMP, based on the mean orbital elements produced by GFZ (GeoForschungs Zentrum) for ephemeris prediction, is presented. Received: 10 July 2001 / Accepted: 17 July 2002 Correspondence to: J. Klokočník at Ondřejov Observatory Acknowledgements. We thank Prof. Dr. Ch. Reigber, Dr. P. Schwintzer, Dr. T. Gruber and Dr. R. K?nig from GFZ Potsdam for various consultations and discussions, and for the CHAMP two-line mean elements. This investigation was performed under the aegis of CEDR (Center for Earth's Dynamics Research, Prague-Ondřejov); it has been supported by project LN00A005 (provided by the Ministry of Education of the Czech Republic) and by grant A 3004 of the Grant Agency of the Academy of Sciences of the Czech Republic.     

Large Scale Rockfall Reach Susceptibility Maps in La Cabrera Sierra (Madrid) performed with GIS and Dynamic Analysis at 1:5,000

For the testing of the effect on the tsunami prevention facilities, a simplified methodfor tsunami risk assessment was suggested without wave run-up analysis. This methodis proposed using calculated offshore tsunami waveform and field reconnaissance suchas the seawall height, time necessary for residents' evacuation and tsunami warninginsurance. Then, two normalized values are evaluated; one is the ratio of calculatedmaximum tsunami height to seawall height, the other is the ratio of time betweentsunami over-topping and evacuation completion to total time required for evacuation.These two values are used to qualitatively estimate the safety of residents and the effectof tsunami prevention facilities, eliminating the necessity to compute complicatedtsunami run-up onshore.     

The transition from an active to a passive margin (SW end of the South Shetland Trench, Antarctic Peninsula)

The lateral ending of the South Shetland Trench is analysed on the basis of swath bathymetry and multichannel seismic profiles in order to establish the tectonic and stratigraphic features of the transition from an northeastward active to a southwestward passive margin style. This trench is associated with a lithospheric-scale thrust accommodating the internal deformation in the Antarctic Plate and its lateral end represents the tip-line of this thrust. The evolutionary model deduced from the structures and the stratigraphic record includes a first stage with a compressional deformation, predating the end of the subduction in the southwestern part of the study area that produced reverse faults in the oceanic crust during the Tortonian. The second stage occurred during the Messinian and includes distributed compressional deformation around the tip-line of the basal detachment, originating a high at the base of the slope and the collapse of the now inactive accretionary prism of the passive margin. The initial subduction of the high at the base of the slope induced the deformation of the accretionary prism and the formation of another high in the shelf—the Shelf Transition High. The third stage, from the Early Pliocene to the present-day, includes the active compressional deformation of the shelf and the base-of-slope around the tip-line of the basal detachment, while extensional deformations are active in the outer swell of the trench.     

P and S Mach Waves Generated by the Detonation of a Cylindrical Explosive Charge - Experiments and Simulations

During blasting operations in a competent rock mass without open joints or faults, the structure of the wave fields generated by the detonation of a cylindrical charge depends mainly on the relation between the velocity of detonation ( c d ) and the velocity of the P and S waves, c P and c S , respectively, in the medium. If c d is larger than c P , detonation is supersonic and Mach P and S waves appear. If c d is between c P and c S , detonation is transonic, and only an S Mach wave is generated. If c d is less than c S , subsonic detonation occurs. In this paper, we describe an experiment conducted at the Andina mine (Chile) in order to demonstrate the differences between the wave fields generated by supersonic and transonic detonation. We also present a simulation of the detonation of a cylindrical charge using the finite differences Itasca software FLAC3.4. In particular, we show how important it is to secure a good adjustment of the rheological behaviour of the medium, in our case, a strain softening elasto-plastic model.     

Numerical modeling of the transient hydrogeological response produced by tunnel construction in fractured bedrocks

Groundwater inflows into tunnels constructed in fractured bedrocks not only constitute an important factor controlling the rate of advancement in driving the tunnel but may pose potential hazards. Drawdowns caused by tunnel construction may also induce geotechnical and environmental impacts. Here we present a numerical methodology for the dynamic simulation of the hydrogeological transient conditions induced by the tunnel front advance. The methodology is based on the use of a Cauchy boundary condition at the points lying along the tunnel according to which water discharge, Q, is computed as the product of a leakage coefficient, , and the head difference, (H−h), where H is the prescribed head at the tunnel wall and h is the hydraulic head in the fractured rock in the close vicinity of the tunnel. At a given position of the tunnel, is zero until the tunnel reaches such position when it is assigned a positive value. The use of step-wise time functions for allows an efficient and accurate simulation of the transient hydrogeological conditions at and around the tunnel during the excavation process. The methodology has been implemented in TRANMEF-3, a finite element computer code for groundwater flow in 3D fractured media developed at the University of A Coruña, Spain, and has been used to simulate the impact of a tunnel on the groundwater system at the Äspö island (Sweden). This tunnel was constructed to access an underground laboratory for research on radioactive waste disposal. The large amount of available data at this site provides a unique opportunity to test the performance of the numerical model and the proposed methodology for tunnel advance. With just minor calibration, the numerical model is able to reproduce accurately the measurements of inflows into the tunnel at several reaches and hydraulic heads at surface-drilled boreholes. These results obtained at the Äspö site lead us to conclude that accurate predictions of the transient hydrogeological responses induced by tunneling works in fractured bedrocks, can be achieved provided that a sound hydrogeological characterization of large-scale fracture zones is available.     

On Modelling Discrete Geological Structures as Markov Random Fields

The purpose of this paper is to extend the locally based prediction methodology of BayMar to a global one by modelling discrete spatial structures as Markov random fields. BayMar uses one-dimensional Markov-properties for estimating spatial correlation and Bayesian updating for locally integrating prior and additional information. The methodology of this paper introduces a new estimator of the field parameters based on the maximum likelihood technique for one-dimensional Markov chains. This makes the estimator straightforward to calculate also when there is a large amount of missing observations, which often is the case in geological applications. We make simulations (both unconditional and conditional on the observed data) and maximum a posteriori predictions (restorations) of the non-observed data using Markov chain Monte Carlo methods, in the restoration case by employing simulated annealing. The described method gives satisfactory predictions, while more work is needed in order to simulate, since it appears to have a tendency to overestimate strong spatial dependence. It provides an important development compared to the BayMar-methodology by facilitating global predictions and improved use of sparse data.     

最大效果测度值法研究矿井突水水源

为了快速准确地查清矿井突水水源，排除水患对生产的影响，利用最大效果测度值燕子山矿14-2＃煤层顶板砂岩裂隙水，其导水通道为工作面所在盘区发育的F23和F38正断层。最大效果测度值为矿井防治水工作提供了一个新的思路。     

Late Neoproterozoic crustal growth in the European Variscides: Nd isotope and geochemical evidence from the Sierra de Córdoba Andesites (Ossa-Morena Zone, Southern Spain)

Because the Hercynian overprint was extremely weak, the Sierra de Córdoba (southeastern Ossa-Morena Zone, OMZ) provides an excellent opportunity to study the tectonic evolution of sequences deposited close to the Late Neoproterozoic–Early Palaeozoic boundary. In order to put constraints on the sources and geodynamic significance of the Late Proterozoic magmatism, a representative set of 18 igneous rocks, and 3 interbedded sedimentary rocks from the San Jerónimo Formation have been studied for major and trace element geochemistry and for the Sm–Nd isotopic systematics. The igneous rocks are generally porphyritic to microporphyritic andesites, with abundant plagioclase (±amphibole) phenocrysts. With the exception of two intrusive rocks, possibly not related to the Late Proterozoic episode, all the samples display positive _{Nd550 Ma} values, ranging from +2.9 to +7.6. Most of them, with +4<_{Nd550 Ma}<+6, exhibit LREE enrichment, high La/Nb ratios, and elevated Zr/Nb ratios ranging from 21 to 32. There is no obvious correlation between the shape of REE patterns, La/Nb ratios and _{Nd550 Ma} values, precluding simple models of late-stage interaction with typical crustal components having low _Nd and high LREE/HREE and La/Nb ratios. Based on their major element composition and enriched, continental crust-like trace element characteristics, combined with distinctly positive _Nd initial values, the Córdoba andesites document an episode of crustal growth through the addition of calc-alkaline magmas, extracted from a mantle reservoir which was strongly depleted in LREE on a time-integrated basis. The occurrence of interlayered sediments of continental provenance (negative _Nd values) does not favour a purely ensimatic arc setting, remote from continental land masses, for this subduction-related magmatism, but the geochemical data suggest an active margin environment located on relatively juvenile crust. In any case, the Córdoba andesites document the addition of materials chemically similar to the bulk continental crust which were extracted from mantle sources with strong time-integrated LREE depletion. Therefore, they provide evidence for crustal growth related to Cadomian orogenic events during Late Proterozoic times.