A Lost Realm in the Internal Domains of the Betic-Rif Orogen (Spain and Morocco): Evidence from Conglomerates and Consequences for Alpine Geodynamic Evolution

The Malaguide-Ghomaride Complex is capped by Upper Oligocene-Aquitanian clastic deposits postdating early Alpine orogenesis but predating the main tectonic-metamorphic evolution, end of nappe emplacement, unroofing, and exhumation of the metamorphic units of the Betic-Rif Orogen. Two conglomerate intervals within these deposits are characterized by clasts of sedimentary, epimetamorphic, and mafic volcanic rocks derived from Malaguide-Ghomaride units and by clasts of acidic magmatic and orthogneissic rocks of unknown provenance, here studied. Magmatic rocks originated from late-Variscan two-mica cordierite-bearing granitoids and, subordinately, from aplitic dikes. Orthogneisses derive from similar plutonic rocks but are affected by an Alpine metamorphic overprint evolving from greenschist (T=510&j0;-530 degrees C and P=5-6 kbar) to low-temperature amphibolite facies (T>550&j0;C and P<3 kbar). Such a plutonic rock suite is unknown in any Betic-Rif unit or in the basement of the Alboran Sea, and the metamorphic evolution in the orthogneisses is different from (and older than) that of Alpujarride-Sebtide rocks to which they were formerly ascribed. Magmatic and metamorphic rocks very similar to those studied characterize the basements of some Kabylia and Calabria-Peloritani units. Therefore, the source area is a currently lost continental-crust realm of Calabria-Peloritani-Kabylia type, located to the ESE of the Malaguide-Ghomaride Domain and affected by a pre-latest Oligocene Alpine metamorphism. Increasingly active tectonics transformed this realm into rising areas from which erosion fed small subsiding synorogenic basins formed on the Malaguide-Ghomaride Complex. This provenance analysis demonstrates that all these domains constituted a single continental-crust block until Aquitanian-Burdigalian times, before its dispersal around nascent western Mediterranean basins.     

Glacier development and topographic context

This paper analyses the topographic context of the remaining glaciated areas in the Maladeta Massif (Central Spanish Pyrenees). These ice‐covered surfaces have been incorporated into a geographic information system (GIS) in an attempt at correlating the presence of ice with a range of topographic variables obtained from a digital elevation model. The use of generalized additive models and binary regression tree models enabled us (i) to quantify the spatial variability in the distribution of glaciers attributable to characteristics of the local terrain, (ii) to investigate the interaction between the variables that account for the ice cover distribution and (iii) to map the probability of glacier development. Our results show that although the development of glaciers depends on regional climate conditions, the topographic context is of paramount importance in determining the location, extent, shape and recent evolution of each glacial body. Thus, the joint effect of altitude, exposure to incoming solar radiation, slope and mean curvature is able to explain more than 70 per cent of the observed variance. Copyright © 2006 John Wiley & Sons, Ltd.     

Enforcing smoothness and assessing uncertainty in non-linear one-dimensional prestack seismic inversion

Estimation of elastic properties of rock formations from surface seismic amplitude measurements remains a subject of interest for the exploration and development of hydrocarbon reservoirs. This paper develops a global inversion technique to estimate and appraise 1D distributions of compressional‐wave velocity, shear‐wave velocity and bulk density, from normal‐moveout‐corrected PP prestack surface seismic amplitude measurements. Specific objectives are: (a) to evaluate the efficiency of the minimization algorithm (b) to appraise the impact of various data misfit functions, and (c) to assess the effect of the degree and type of smoothness criterion enforced by the inversion. Numerical experiments show that very fast simulated annealing is the most efficient minimization technique among alternative approaches considered for global inversion. It is also found that an adequate choice of data misfit function is necessary for a reliable and efficient match of noisy and sparse seismic amplitude measurements. Several procedures are considered to enforce smoothness of the estimated 1D distributions of elastic parameters, including predefined quadratic measures of length, flatness and roughness. Based on the general analysis of global inversion techniques, we introduce a new stochastic inversion algorithm that initializes the search for the minimum with constrained random distributions of elastic parameters and enforces predefined autocorrelation functions (semivariograms). This strategy readily lends itself to the assessment of model uncertainty. The new global inversion algorithm is successfully tested on noisy synthetic amplitude data. Moreover, we present a feasibility analysis of the resolution and uncertainty of prestack seismic amplitude data to infer 1D distributions of elastic parameters measured with wireline logs in the deepwater Gulf of Mexico. The new global inversion algorithm is computationally more efficient than the alternative global inversion procedures considered here.     

DEVEX-disdrometer evaluation experiment: Basic results and implications for hydrologic studies

Increasing our understanding of the small scale variability of drop size distributions (DSD), and therefore of several bulk characteristics of rainfall processes, has major implications for our interpretation of the remote sensing based estimates of precipitation and its uncertainty. During the spring and summer of 2002 the authors conducted the DEVEX experiment (disdrometer evaluation experiment) to compare measurements of natural rain made with three different types of disdrometers collocated at the Iowa City Municipal Airport in Iowa City, Iowa in the Midwestern United States. This paper focuses on the evaluation of the instruments rather than analysis of the hydrometeorological aspects of the observed events. The comparison demonstrates discrepancies between instruments. The authors discuss the systematic and random effects in terms of rainfall quantities, drop size distribution properties, and the observed drop size vs. velocity relationships. Since the instruments were collocated, the effects of the natural variability of rain are reduced some with time integration, isolating the instrumental differences. The authors discuss the status of DSD measurement technologies and the implications for a range of hydrologic applications from remote sensing of rainfall to atmospheric deposition to soil erosion and sediment transport in the environment. The data set collected during the DEVEX experiment is made available to the research community.     

Uncertainty in GPS Networks due to Remaining Systematic Errors: The Interval Approach

A realistic assessment of the total uncertainty budget of Global Positioning System (GPS) observations and its adequate mathematical treatment is a basic requirement for all analysis and interpretation of GPS-derived point positions, in particular GPS heights, and their respective changes. This implies not only the random variability but also the remaining systematic errors. At present in geodesy, the main focus is on stochastic approaches in which errors are modeled by means of random variables. Here, an alternative approach based on interval mathematics is presented. It allows us to model and to quantify the impact of remaining systematic errors in GPS carrier-phase observations on the final results using deterministic error bands. In this paper, emphasis is given to the derivation of the observation intervals based on influence parameters and to the study of the complex linear transfer of this type of uncertainty to estimated point positions yielding zonotopes. From the presented simulation studies of GPS baselines, it turns out that the uncertainty due to remaining systematic effects dominates the total uncertainty budget for baselines longer than 200 km.     

The relation between rigorous and Helmert’s definitions of orthometric heights

Following our earlier definition of the rigorous orthometric height [J Geod 79(1-3):82–92 (2005)] we present the derivation and calculation of the differences between this and the Helmert orthometric height, which is embedded in the vertical datums used in numerous countries. By way of comparison, we also consider Mader and Niethammer’s refinements to the Helmert orthometric height. For a profile across the Canadian Rocky Mountains (maximum height of ~2,800 m), the rigorous correction to Helmert’s height reaches ~13 cm, whereas the Mader and Niethammer corrections only reach ~3 cm. The discrepancy is due mostly to the rigorous correction’s consideration of the geoid-generated gravity disturbance. We also point out that several of the terms derived here are the same as those used in regional gravimetric geoid models, thus simplifying their implementation. This will enable those who currently use Helmert orthometric heights to upgrade them to a more rigorous height system based on the Earth’s gravity field and one that is more compatible with a regional geoid model.     

A comparative study of the design spectra defined by Eurocode 8, UBC, IBC and Turkish Earthquake Code on R/C sample buildings

Earthquake codes have been revised and updated depending on the improvements in the representation of ground motions, soils and structures. These revisions have been more frequently seen in recent years. One of the key changes in earthquake codes has been performed on the design spectra. In this paper, the design spectra recommended by Turkish Earthquake Code and three other well known codes (Uniform Building Code, Eurocode 8, and International Building Code) are considered for comparison. The main purpose of this study is to investigate the differences caused by the use of different codes in the dynamic analysis and seismic verification of given types of buildings located at code defined different sites. The differences in expressions and some important points for elastic and inelastic spectra defined by the codes are briefly illustrated in tables and figures. Periods, base shears, lateral displacements and interstory drifts for the analyzed buildings located at code defined ground type are comparatively presented.     

Effect of common point selection on coordinate transformation parameter determination

The use of satellite positioning techniques commonly requires a transformation from a Conventional Terrestrial coordinate system to a Geodetic coordinate system, or vice versa. For such a transformation, the main problem is the determination of transformation parameters between these coordinate systems. The transformation parameters are estimated by a least-squares process using “common” points, i.e., those points whose coordinates are known in both systems. Therefore, the precision of so estimated transformation parameters is closely related to certain characteristics of the common points. In this contribution, we have formulated some theoretical relations between the transformation parameters and the number and the distribution of common points, and corroborated the theoretical results numerically, using a simulated geodetic network.     

Seismic Moment Tensor Resolution on a Local Scale: Simulated Rockburst and Mine-induced Seismic Events in the Kopanang Gold Mine, South Africa

Seismic records contain information about the effect of the source as well as the effect of wave propagation through the rock mass. The effect of wave propagation is usually not well known as only simplified models of geological structures are available. Therefore, the information about the source retrieved by inverting seismograms may include errors due to incomplete knowledge of the rock mass along the propagation path, which in turn cause a distortion in the calculated moment tensor (MT). The distortion of the MT on a local scale was observed by inverting records of a simulated rockburst conducted at the Kopanang gold mine in South Africa. A dominant isotropic component of the explosive characteristics was found from the inversion. The deviatoric components retrieved from the blast are spurious. A test of their stability indicated that they are not significant, assuming an uncertainty above 5% for velocities and 10% for attenuation within the homogeneous model available for the mine. Thus, the retrieval of the MT from records of local networks in mines using a homogeneous model of the rock mass seems to be feasible. However, the homogeneous model of the rock mass can only be applied to close stations, within a few kilometers of the source. The seismic records from distant stations were too complex to be modelled by a homogeneous rock mass. Records of six mine-induced seismic events recorded at the Kopanang gold mine were also inverted. A vertical linear dipole along the pressure (P) axis was found for three of the events, suggesting a pillar burst. The mechanism of two events contains an isotropic implosion together with a nearly vertical dip-slip, and seems to indicate a combination of a cavity collapse with a down dip-slip along a nearly vertical fault. One event corresponds to a dipole along the tensional (T) axis. However, it is vertical, thus its association with tensile faulting of the hangingwall is uncertain.     

Explicit formula for the geoid-quasigeoid separation

The explicit formula for the geoid-to-quasigeoid correction is derived in this paper. On comparing the geoidal height and height anomaly, this correction is found to be a function of the mean value of gravity disturbance along the plumbline within the topography. To evaluate the mean gravity disturbance, the gravity field of the Earth is decomposed into components generated by masses within the geoid, topography and atmosphere. Newton’s integration is then used for the computation of topography-and atmosphere-generated components of the mean gravity, while the combined solution for the downward continuation of gravity anomalies and Stokes’ boundary-value problem is utilized in computing the component of mean gravity disturbance generated by mass irregularities within the geoid. On application of this explicit formulism a theoretical accuracy of a few millimetres can be achieved in evaluation of the geoid-to-quasigeoid correction. However, the real accuracy could be lower due to deficiencies within the numerical methods and to errors within the input data (digital terrain and density models and gravity observations).