A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Contribution of storms to groundwater recharge in the semi‐arid region of Karnataka,India

This paper describes the application of environmental isotopes and injected tracer techniques in estimating the contribution of storms as well as annual precipitation to groundwater recharge and its circulation, in the semi‐arid region of Bagepalli, Kolar district, Karnataka. Environmental isotopes ²H, ¹⁸O and ³H were used to study the effect of storms on the hydrological system, and an isotope balance was used to compute the contribution of a storm component to the groundwater. Some of the groundwater samples collected during the post‐storm periods were highly depleted in stable isotope content with higher deuterium excess relative to groundwater from the pre‐storm periods. Significant variation in deuterium excess in groundwater from the same area, collected in two different periods, indicates the different origin of air masses. The estimated recharge component of a storm event of 600 mm to the groundwater was found to be in the range of 117–165 mm. There was no significant variation in environmental tritium content of post‐storm and pre‐storm groundwater, indicating the fast circulation of groundwater in the system. After completion of the environmental isotope work, an injected radiotracer ³H technique was applied to estimate the direct recharge of total precipitation to the groundwater. The estimated recharge to the groundwater is 33 mm of the 550 mm annual precipitation during 1992. Copyright © 2003 John Wiley & Sons, Ltd.     

Pb, Sr and Nd isotope study of two co-existing magmas in the Nízke Tatry Mountains, Western Carpathians (Slovakia)

Summary Two co-existing plutonic rocks (diorite and granodiorite) were studied from an intrusion of Variscan age in the Raztocna Valley – Nízke Tatry Mountains, Western Carpathians. Geochemical analyses of major and trace elements constrain a volcanic arc as emplacement environment and give the first hints of a mixture of two magmatic end-members: the so-called Prasivá granodiorite and the Raztocna diorite. The ⁸⁷Sr/⁸⁶Sr₍₀₎ ratios vary between 0.7075 and 0.7118, the ε Nd₍₀₎ values range from −1.4 to −5.0. Common Pb isotopes reveal a dominant crustal source with minor influences from a mantle and a lower crustal source. Modelling based on Sr and Nd isotope data and using three component mixing calculations indicates that mixing of 2/3 of upper mantle material with 1/3 upper crustal material can produce the isotopic composition of the Raztocna diorite. Very minor amounts of lower crust were incorporated in the diorite. For the Prasivá granodiorite, the mixing ratio of upper mantle and upper crust is similar, but a lower crustal reservoir contributed about 5–10% of the source material.     

Principal Components Analysis of Rainfall and Associated Synoptic Models of the Southwest Monsoon over India

Summary  The summer monsoon circulation shows various spatial and temporal oscillations and often a combination of systems produces an integrated effect. In this study phases of the southwest (SW) monsoon have been identified in an objective manner with the help of T-mode principal component analysis (PCA) of weekly rainfall anomalies. Mean composite charts have been prepared utilising all available upper air data (1977–1986) for each category of the SW monsoon epochs identified by the PCA. These sets of charts have been constructed for both the strong and weak phases associated with the first four significant principal components (PCs). A well defined east-west oriented trough system, extending from about 28° N Latitude/65° E Longitude to 20° N Latitude/90° E Longitude, in the lower levels, has been the main feature associated with the strong phase of the monsoon corresponding to PC I. The trough in the lower levels is more marked in the eastern half compared to the western half in both the sets of charts associated with strong phases of the monsoon related to the PC II and PC III. With PC II, the position of the troughs in the lower levels is further north of its location in PC III. The east-west trough system associated with the strong phase of PC IV has a large southward tilt with height. The charts corresponding to the weak phases of these PCs have synoptic features, such as the position of the trough close to the foothills of the Himalayas, and the shifting of middle and upper tropospheric anticyclones to the south. The study suggests an objective method of interpretation of principal components by utilising synoptic data. In addition, synoptic models and data sets corresponding to different phases of the monsoon can also be prepared in an objective manner by such PCA. Received July 18, 1997 Revised April 30, 1999     

Impact of Earth radiation pressure on GPS position estimates

GPS satellite orbits available from the International GNSS Service (IGS) show a consistent radial bias of up to several cm and a particular pattern in the Satellite Laser Ranging (SLR) residuals, which are suggested to be related to radiation pressure mismodeling. In addition, orbit-related frequencies were identified in geodetic time series such as apparent geocenter motion and station displacements derived from GPS tracking data. A potential solution to these discrepancies is the inclusion of Earth radiation pressure (visible and infrared) modeling in the orbit determination process. This is currently not yet considered by all analysis centers contributing to the IGS final orbits. The acceleration, accounting for Earth radiation and satellite models, is introduced in this paper in the computation of a global GPS network (around 200 IGS sites) adopting the analysis strategies from the Center for Orbit Determination in Europe (CODE). Two solutions covering 9 years (2000–2008) with and without Earth radiation pressure were computed and form the basis for this study. In previous studies, it has been shown that Earth radiation pressure has a non-negligible effect on the GPS orbits, mainly in the radial component. In this paper, the effect on the along-track and cross-track components is studied in more detail. Also in this paper, it is shown that Earth radiation pressure leads to a change in the estimates of GPS ground station positions, which is systematic over large regions of the Earth. This observed “deformation” of the Earth is towards North–South and with large scale patterns that repeat six times per GPS draconitic year (350 days), reaching a magnitude of up to 1 mm. The impact of Earth radiation pressure on the geocenter and length of day estimates was also investigated, but the effect is found to be less significant as compared to the orbits and position estimates.     

Evidence for short term changes in the intensity and spectrum of hard X-rays from Cyg X-1

The paper presents experimental evidence for the existence of fast intensity fluctuations with time scales of the order of a minute in the X-ray emission from Cyg X-1 at energies greater than 29 keV. Spectral variations over time intervals of 20–25 min are also observed in the same energy range. Whereas, similar intensity and spectral fluctuations have been reported earlier at lower energies the observations presented here is the first evidence for the existence of similar fluctuations at high energies.     

Airborne gravimetry used in precise geoid computations by ring integration

Two detailed geoids have been computed in the region of North Jutland. The first computation used marine data in the offshore areas. For the second computation the marine data set was replaced by the sparser airborne gravity data resulting from the AGMASCO campaign of September 1996. The results of comparisons of the geoid heights at on-shore geometric control showed that the geoid heights computed from the airborne gravity data matched in precision those computed using the marine data, supporting the view that airborne techniques have enormous potential for mapping those unsurveyed areas between the land-based data and the off-shore marine or altimetrically derived data. Received: 7 July 1997 / Accepted: 22 April 1998     

Cascaded image analysis for dynamic crack detection in material testing

Concrete probes in civil engineering material testing often show fissures or hairline-cracks. These cracks develop dynamically. Starting at a width of a few microns, they usually cannot be detected visually or in an image of a camera imaging the whole probe. Conventional image analysis techniques will detect fissures only if they show a width in the order of one pixel. To be able to detect and measure fissures with a width of a fraction of a pixel at an early stage of their development, a cascaded image analysis approach has been developed, implemented and tested. The basic idea of the approach is to detect discontinuities in dense surface deformation vector fields. These deformation vector fields between consecutive stereo image pairs, which are generated by cross correlation or least squares matching, show a precision in the order of 1/50 pixel. Hairline-cracks can be detected and measured by applying edge detection techniques such as a Sobel operator to the results of the image matching process. Cracks will show up as linear discontinuities in the deformation vector field and can be vectorized by edge chaining. In practical tests of the method, cracks with a width of 1/20 pixel could be detected, and their width could be determined at a precision of 1/50 pixel.     

Reducing the draconitic errors in GNSS geodetic products

Systematic errors at harmonics of the GPS draconitic year have been found in diverse GPS-derived geodetic products like the geocenter $Z$ -component, station coordinates, $Y$ -pole rate and orbits (i.e. orbit overlaps). The GPS draconitic year is the repeat period of the GPS constellation w.r.t. the Sun which is about 351 days. Different error sources have been proposed which could generate these spurious signals at the draconitic harmonics. In this study, we focus on one of these error sources, namely the radiation pressure orbit modeling deficiencies. For this purpose, three GPS+GLONASS solutions of 8 years (2004–2011) were computed which differ only in the solar radiation pressure (SRP) and satellite attitude models. The models employed in the solutions are: (1) the CODE (5-parameter) radiation pressure model widely used within the International GNSS Service community, (2) the adjustable box-wing model for SRP impacting GPS (and GLONASS) satellites, and (3) the adjustable box-wing model upgraded to use non-nominal yaw attitude, specially for satellites in eclipse seasons. When comparing the first solution with the third one we achieved the following in the GNSS geodetic products. Orbits: the draconitic errors in the orbit overlaps are reduced for the GPS satellites in all the harmonics on average 46, 38 and 57 % for the radial, along-track and cross-track components, while for GLONASS satellites they are mainly reduced in the cross-track component by 39 %. Geocenter $Z$ -component: all the odd draconitic harmonics found when the CODE model is used show a very important reduction (almost disappearing with a 92 % average reduction) with the new radiation pressure models. Earth orientation parameters: the draconitic errors are reduced for the $X$ -pole rate and especially for the $Y$ -pole rate by 24 and 50 % respectively. Station coordinates: all the draconitic harmonics (except the 2nd harmonic in the North component) are reduced in the North, East and Height components, with average reductions of 41, 39 and 35 % respectively. This shows, that part of the draconitic errors currently found in GNSS geodetic products are definitely induced by the CODE radiation pressure orbit modeling deficiencies.