Geology and metamorphism of the Ladakh Terrane and Shyok Suture Zone in the Chogo Lungma – Turmik area (northern Pakistan)

Abstract

Metamorphism of the Askore Amphibolite, metabasaltic and metasedimentary medium-grade hornblendebearing schists at the northernmost portion of the Ladakh Terrane and of the Shyok Suture Zone, mainly a low-grade volcano-sedimentary series, has been studied in the area between the Chogo Lungma glacier and the Indus river halfway between Skardu and Rondu.

In the Askore Amphibolite the peak assemblage in the amphibolite facies defines the regional metamorphic foliation, and is overprinted by a later static recrystallization at comparable P-T conditions. In spite of similar peak temperatures (630–650°C), geobarometry based on amphibole composition reveals a marked difference between garnet – epidote – andesine amphibolites exposed just above the Main Mantle Thrust at the head of Turmik valley, which equilibrated at high pressures (about 10 kbar) in late Miocene (Tortonian), and biotite – epidote – oligoclase amphibolites outcropping at the mouth of Turmik valley, which equilibrated at pressures of c. 6 kbar before late Eocene (Priabonian).

The Dasu Ultramafite and other smaller lens-shaped bodies of low- to medium-grade metaperidotite separate the Ladakh Terrane from the Shyok Suture Zone. They are antigorite serpentinites, often with talc and magnesite, in which relict cumulitic structures are locally recognisable. The ultramafites may represent remnants of oceanic lithosphere separating the Ladakh-Kohistan island arc from the Asian plate, or they may be deep crustal rocks stripped from the basement of the arc.

The mostly greenschist-facies Shyok Suture Zone shows the lithology of a calc-alkaline volcano-sedimentary series. It is supposed to be a remnant of a back arc basin of early Cretaceous age, separating the arc from the southern margin of Asia. Chloritoid, kyanite and biotite have been found in individual thrust sheets occurring at different structural levels and totally subordinate in volume to very low- and low-grade rocks. Such sharp differences in mineral paragenesis, together with field evidence of local shear, suggest a complex internal structure for the Shyok Suture Zone. From the head of Chogo Lungma glacier to the Basha valley, close to the contact with the Karakorum Metamorphic Complex, the rocks of the Shyok Suture Zone record a late Miocene metamorphic event at medium pressures and temperatures. Thermobarometric and geochronological evidence suggests that this event can be related to the exhumation and thrusting of the Karakorum metamorphic core over the Shyok Suture Zone.     

Evaluating the Poroelastic Effect on Anisotropic, Organic-Rich, Mudstone Systems

Understanding the poroelastic effect on anisotropic organic-rich mudstones is of high interest and value for evaluating coupled effects of rock deformation and pore pressure, during drilling, completion and production operations in the oilfield. These applications include modeling and prevention of time-dependent wellbore failure, improved predictions of fracture initiation during hydraulic fracturing operations (Suarez-Rivera et al. Presented at the Canadian Unconventional Resources Conference held in Calgary, Alberta, Canada, 15–17 November 2011. CSUG/SPE 146998 2011), improved understanding of the evolution of pore pressure during basin development, including subsidence and uplift, and the equilibrated effective in situ stress (Charlez, Rock mechanics, vol 2 1997; Katahara and Corrigan, Pressure regimes in sedimentary basins and their prediction: AAPG Memoir, vol 76, pp 73–78 2002; Fjær et al. Petroleum related rock mechanics. 2nd edn 2008). In isotropic rocks, the coupled poro-elastic deformations of the solid framework and the pore fluids are controlled by the Biot and Skempton coefficients. These are the two fundamental properties that relate the rock framework and fluid compressibility and define the magnitude of the poroelastic effect. In transversely isotropic rocks, one desires to understand the variability of these coefficients along the directions parallel and longitudinal to the principal directions of material symmetry (usually the direction of bedding). These types of measurements are complex and uncommon in low-porosity rocks, and particularly problematic and scarce in tight shales. In this paper, we discuss a methodology for evaluating the Biot’s coefficient, its variability along the directions parallel and perpendicular to bedding as a function of stress, and the homogenized Skempton coefficient, also as a function of stress. We also predict the pore pressure change that results during undrained compression. Most importantly, we provide values of transverse and longitudinal Biot’s coefficients and the homogenized Skempton coefficient for two important North American, gas-producing, organic-rich mudstones. These results could be used for petroleum-related applications.     

Using Java for distributed computing in the Gaia satellite data processing

In recent years Java has matured to a stable easy-to-use language with the flexibility of an interpreter (for reflection etc.) but the performance and type checking of a compiled language. When we started using Java for astronomical applications around 1999 they were the first of their kind in astronomy. Now a great deal of astronomy software is written in Java as are many business applications. We discuss the current environment and trends concerning the language and present an actual example of scientific use of Java for high-performance distributed computing: ESA’s mission Gaia. The Gaia scanning satellite will perform a galactic census of about 1,000 million objects in our galaxy. The Gaia community has chosen to write its processing software in Java. We explore the manifold reasons for choosing Java for this large science collaboration. Gaia processing is numerically complex but highly distributable, some parts being embarrassingly parallel. We describe the Gaia processing architecture and its realisation in Java. We delve into the astrometric solution which is the most advanced and most complex part of the processing. The Gaia simulator is also written in Java and is the most mature code in the system. This has been successfully running since about 2005 on the supercomputer “Marenostrum” in Barcelona. We relate experiences of using Java on a large shared machine. Finally we discuss Java, including some of its problems, for scientific computing.     

Spatial and temporal variations of atmospheric aerosol optical thickness in northwestern Mexico

The purpose of this paper was to study aerosol particles in the Northwestern region of Mexico (NWM) through Aerosol Optical Thickness (AOT) parameter in the atmosphere. This parameter represents one of the extinction coefficients of solar radiation and the rate of suspended particles in the atmosphere. For determination of AOT, we considered the use of remote sensors outside of the atmosphere. In particular, Moderate Resolution Imaging Spectroradiometer (MODIS) which can measure the atmospheric AOT thickness. Data from the MODIS sensor must be validated before they are considered reliable. For this task, we required surface measurements to obtain a correlation with the data acquired with the remote radiometer. The paper describes the validation process performed for data obtained with MODIS through measurements provided by an AErosol RObotic NETwork (AERONET) photometer located in the city of Hermosillo, Sonora, NWM. Additionally, we carried out a temporal analysis based on the behavior of the AOT graphics and spatial analysis supported in maps with sufficient information.     

Factors contributing to heavy metal accumulation in sediments and in the intertidal mussel Perna perna in the Gulf of Annaba (Algeria)

This paper presents the results of a seasonal survey of heavy metals accumulated in sediments and in the soft parts of the body of the mussel Perna perna at four stations in the Gulf of Annaba (Algeria). Pooled soft tissues from 10 mussels representing the entire range of sizes were digested in nitric acid. Statistical analysis reveals a significant seasonal effect on all the measured metals, the highest values being recorded in winter. With the exception of Cr, the levels for all metals were significantly higher in the east, at the outlet of the Seybouse River, than at all other monitoring stations. The study also shows that north-western waters are subject to a significantly lower degree of heavy metal pollution than elsewhere in the gulf. Levels were nevertheless within the limits of public health standards.     

Seabird aggregative patterns: A new tool for offshore wind energy risk assessment

The emerging development of offshore wind energy has raised public concern over its impact on seabird communities. There is a need for an adequate methodology to determine its potential impacts on seabirds. Environmental Impact Assessments (EIAs) are mostly relying on a succession of plain density maps without integrated interpretation of seabird spatio-temporal variability. Using Taylor’s power law coupled with mixed effect models, the spatio-temporal variability of species’ distributions can be synthesized in a measure of the aggregation levels of individuals over time and space. Applying the method to a seabird aerial survey in the Ebro Delta, NW Mediterranean Sea, we were able to make an explicit distinction between transitional and feeding areas to define and map the potential impacts of an offshore wind farm project. We use the Ebro Delta study case to discuss the advantages of potential impacts maps over density maps, as well as to illustrate how these potential impact maps can be applied to inform on concern levels, optimal EIA design and monitoring in the assessment of local offshore wind energy projects.     

Chronology and paleoenvironmental interpretation of talus flatiron sequences in a sub‐humid mountainous area: Tremp Depression,Spanish Pyrenees

Talus flatirons are debris‐covered relict slopes, disconnected from the source area, which are relatively common in arid and semi‐arid areas. Talus flatiron sequences record the alternation of accumulation and incision phases. These chronosequences may be used for infer temporal changes in the morphogenetic processes acting on the slopes as well as information on the local paleoclimatic history. Talus flatiron sequences developed in the Tremp Depression, eastern Spanish Pyrenees, are analysed from the geomorphological, chronological and paleoenvironmental perspective. The two groups of relict slopes differentiated by means of detailed geomorphological mapping have been dated by optically stimulated luminescence (OSL) and radiocarbon dating at 25–20 kyr (S3) and 5·4–1·7 cal kyr (S2). The talus flatiron group S3 is correlated with a fluvial terrace of the Noguera Pallaresa River (c. 23 kyr bp ). The comparison of the ages obtained in the Tremp Depression with chronologies published for talus flatiron sequences in semi‐arid areas and other paleoclimatic proxies suggests that the aggradation phases in the slopes occurred during periods with higher humidity and vegetation cover. The chronological differences observed between semi‐arid Spain and the Tremp Depression may be partly related to the more humid climate of the latter mountain area. Copyright © 2013 John Wiley & Sons, Ltd.     

Finding appropriate interpolation techniques for topographic surface generation for mudslide risk zonation

Mudslides and debris flows are now more common problems in tropical regions than landslides. This article addresses the issue of modelling mudslide-susceptible locations using an information value approach. Landsat 7 (ETM+) was employed to create a land-use map with limited field checks. Other parameters considered were lineament, road, soil, stream network and lithology. Topographic parameters such as slope and aspect play a dominant role in slope stability studies. Maps for slope and aspects were developed from a digital elevation model using statistical surface interpolation techniques. This article offers insight into the importance of the selection of suitable surface interpolation techniques. The three surface interpolation techniques evaluated in the study were: inverse distance weighted, Kriging and Spline. These were found to have varied accuracies of interpolation surfaces for all parameters, including elevation, slope and aspect. Surfaces offering the best accuracy were adopted for the information value approach for mudslide susceptibility zonation.     

Current Behaviour and Dynamics of the Lowermost Italian Glacier (Montasio Occidentale,Julian Alps)

Smaller glaciers (<0.5 km²) react quickly to environmental changes and typically show a large scatter in their individual response. Accounting for these ice bodies is essential for assessing regional glacier change, given their high number and contribution to the total loss of glacier area in mountain regions. However, studying small glaciers using traditional techniques may be difficult or not feasible, and assessing their current activity and dynamics may be problematic. In this paper, we present an integrated approach for characterizing the current behaviour of a small, avalanche‐fed glacier at low altitude in the Italian Alps, combining geomorphological, geophysical and high‐resolution geodetic surveying with a terrestrial laser scanner. The glacier is still active and shows a detectable mass transfer from the accumulation area to the lower ablation area, which is covered by a thick debris mantle. The glacier owes its existence to the local topo‐climatic conditions, ensured by high rock walls which enhance accumulation by delivering avalanche snow and reduce ablation by providing topographic shading and regulating the debris budget of the glacier catchment. In the last several years the glacier has displayed peculiar behaviour compared with most glaciers of the European Alps, being close to equilibrium conditions in spite of warm ablation seasons. Proportionally small relative changes have also occurred since the Little Ice Age maximum. Compared with the majority of other Alpine glaciers, we infer for this glacier a lower sensitivity to air temperature and a higher sensitivity to precipitation, associated with important feedback from increasing debris cover during unfavourable periods.