Preliminary study of crust-upper mantle structure of the Tibetan Plateau by using broadband teleseismic body waveforms

As part of a joint Sino-U.S. research project to study the deep structure of the Tibetan Plateau, 11 broadband digital seismic recorders were deployed on the Plateau for one year of passive seismic recording. In this report we use teleseimic P waveforms to study the seismic velocity structure of crust and upper mantle under three stations by receiver function inversion. The receiver function is obtained by first rotating two horizontal components of seismic records into radial and tangential components and then deconvolving the vertical component from them. The receiver function depends only on the structure near the station because the source and path effects have been removed by the deconvolution. To suppress noise, receiver functions calculated from events clustered in a small range of back-azimuths and epicentral distances are stacked. Using a matrix formalism describing the propagation of elastic waves in laterally homogeneous stratified medium, a synthetic receiver function and differential receiver functions for the parameters in each layer can be calculated to establish a linearized inversion for one-dimensional velocity structure. Preliminary results of three stations, Wen-quan, Golmud and Xigatze (Coded as WNDO, TUNL and XIGA), located in central, northern and southern Plateau are given in this paper. The receiver functions of all three stations show clear P-S converted phases. The time delays of these converted phases relative to direct P arrivals are: WNDO 7.9s (for NE direction) and 8.3s (for SE direction), TUNL 8.2s, XIGA 9.0s. Such long time delays indicate the great thickness of crust under the Plateau. The differences between receiver function of these three station shows the tectonic difference between southern and north-central Plateau. The waveforms of the receiver functions for WNDO and TUNL are very simple, while the receiver function of XIGA has an additional midcrustal converted phase. The S wave velocity structures at these three stations are estimated from inversions of the receiver function. The crustal shear wave velocities at WNDO and TUNL are vertically homogeneous, with value between 3.5–3.6 km/s down to Moho. This value in the lower crust is lower than the normal value for the lower crust of continents, which is consistent with the observed strong Sn attenuation in this region. The velocity structure at XIGA shows a velocity discontinuity at depth of 20 km and high velocity value of 4.0 km/s in the midcrust between 20–30 km depth. Similar results are obtained from a DSS profile in southern Tibet. The velocity under XIGA decreases below a depth of 30 km, reaching the lowest value of 3.2 km/s between 50–55 km. depth. This may imply that the Indian crust underthrusts the low part of Tibetan crust in the southern Plateau, forming a “double crust”. The crustal thickness at each of these sites is: WNDO, 68 km; TUNL, 70 km; XI-GA, 80 km. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 581–592, 1992.     

Guest editorial

The Jason-1 satellite was launched on 7 December 2001 with the primary objective of continuing the high accuracy time series of altimeter measurements that began with the TOPEX/Poseidon mission in 1992. To achieve this goal, it is necessary to validate the performance of the Jason-1 measurement system, and to verify that its error budget is at least at the same level as that of the TOPEX/Poseidon mission. The article reviews the main components of the Jason-1 altimetric error budget from instrument characterization to the geophysical use of the data. Using the Interim Geophysical Data Records (16DR) that were distributed to the Jason-1 Science Working Team during the verification phase of the mission, it is shown that the Jason-1 mission is performing well enough to continue studies of the large-scale features of the ocean, and especially to continue time series of mean sea-level variations with an accuracy comparable to TOPEX/Poseidon.     

The Poseidon Operational Tool for the Prediction of Floating Pollutant Transport

In this work the development and the application of an operational management tool for the Greek Seas is described. This tool consists of a three-dimensional floating pollutant prediction model coupled with a weather, a hydrodynamic and a wave model in order to track the movements and the spreading of the pollutants and indicate those coastal areas which might be affected. The tool is part of the Poseidon system which has been designed to provide real time data and forecasts for marine environmental conditions in the Greek Seas. In this paper, we present four case studies based on realistic scenarios that show the value of the application for long-term strategic planning and short-term decision making in oil spill accidents.     

Landscape transform and spatial metrics for mapping spatiotemporal land cover dynamics using Earth Observation data-sets

Analysis of Earth observation (EO) data, often combined with geographical information systems (GIS), allows monitoring of land cover dynamics over different ecosystems, including protected or conservation sites. The aim of this study is to use contemporary technologies such as EO and GIS in synergy with fragmentation analysis, to quantify the changes in the landscape of the Rajaji National Park (RNP) during the period of 19 years (1990–2009). Several statistics such as principal component analysis (PCA) and spatial metrics are used to understand the results. PCA analysis has produced two principal components (PC) and explained 84.1% of the total variance, first component (PC1) accounted for the 57.8% of the total variance while the second component (PC2) has accounted for the 26.3% of the total variance calculated from the core area metrics, distance metrics and shape metrics. Our results suggested that notable changes happened in the RNP landscape, evidencing the requirement of taking appropriate measures to conserve this natural ecosystem.     

COMPUTATIONAL GEOMETRY AND GEOGRAPHY *

New developments in computational geometry are applicable to geographic data processing. Computationally efficient procedures incorporating logarithmic searches over appropriate data structures, and other shortcuts, can yield important savings in computation time and data storage space. There are “fast”methods for two-dimensional problems such as point inclusion, nearest-neighbor and two-nearest-neighbors determination, and grid overlay. Worst-case analysis of algorithms provides useful performance bounds for the developers of geographic production software.     

Continental rift systems and anorogenic magmatism

Precambrian Laurentia and Mesozoic Gondwana both rifted along geometric patterns that closely approximate truncated-icosahedral tessellations of the lithosphere. These large-scale, quasi-hexagonal rift patterns manifest a least-work configuration. For both Laurentia and Gondwana, continental rifting coincided with drift stagnation, and may have been driven by lithospheric extension above an insulated and thermally expanded mantle. Anorogenic magmatism, including flood basalts, dike swarms, anorthosite massifs and granite-rhyolite provinces, originated along the Laurentian and Gondwanan rift tessellations. Long-lived volcanic regions of the Atlantic and Indian Oceans, sometimes called hotspots, originated near triple junctions of the Gondwanan tessellation as the supercontinent broke apart. We suggest that some anorogenic magmatism results from decompression melting of asthenosphere beneath opening fractures, rather than from random impingement of hypothetical deep-mantle plumes.     

Travel time controls the magnitude of nitrate discharge in groundwater bypassing the riparian zone to a stream on Virginia's coastal plain

Groundwater that bypasses the riparian zone by travelling along deep flow paths may deliver high concentrations of fertilizer‐derived NO₃^? to streams, or it may be impacted by the NO₃^? removal process of denitrification in streambed sediments. In a study of a small agricultural catchment on the Atlantic coastal plain of Virginia's eastern shore, we used seepage meters deployed in the streambed to measure specific discharge of groundwater and its solute concentrations for various locations and dates. We used values of Cl^? concentration to discriminate between bypass water recharged distal to the stream and that contained high NO₃^? but low Cl^? concentrations and riparian‐influenced water recharged proximal to the stream that contained low NO₃^? and high Cl^? concentrations. The travel time required for bypass water to transit the 30‐cm‐thick, microbially active denitrifying zone in the streambed determined the extent of NO₃^? removal, and hydraulic conductivity determined travel time through the streambed sediments. At all travel times greater than 2 days, NO₃^? removal was virtually complete. Comparison of the timescales for reaction and transport through the streambed sediments in this system confirmed that the predominant control on nitrate flux was travel time rather than denitrification rate coefficients. We conclude that extensive denitrification can occur in groundwater that bypasses the riparian zone, but a residence time in biologically active streambed sediments sufficient to remove a large fraction of the NO₃^? is only achieved in relatively low‐conductivity porous media. Instead of viewing them as separate, the streambed and riparian zone should be considered an integrated NO₃^? removal unit. Copyright © 2011 John Wiley & Sons, Ltd.     

Dynamic stochastic estimation of physical variables

A fundamental problem facing the physical sciences today is analysis of natural variations and mapping of spatiotemporal processes. Detailed maps describing the space/time distribution of groundwater contaminants, atmospheric pollutant deposition processes, rainfall intensity variables, external intermittency functions, etc. are tools whose importance in practical applications cannot be overestimated. Such maps are valuable inputs for numerous applications including, for example, solute transport, storm modeling, turbulent-nonturbulent flow characterization, weather prediction, and human exposure to hazardous substances. The approach considered here uses the spatiotemporal random field theory to study natural space/time variations and derive dynamic stochastic estimates of physical variables. The random field model is constructed in a space/time continuum that explicitly involves both spatial and temporal aspects and provides a rigorous representation of spatiotemporal variabilities and uncertainties. This has considerable advantages as regards analytical investigations of natural processes. The model is used to study natural space/time variations of springwater calcium ion data from the Dyle River catchment area, Belgium. This dataset is characterized by a spatially nonhomogeneous and temporally nonstationary variability that is quantified by random field parameters, such as orders of space/time continuity and random field increments. A rich class of covariance models is determined from the properties of the random field increments. The analysis leads to maps of continuity orders and covariances reflecting space/time calcium ion correlations and trends. Calcium ion estimates and the associated statistical errors are calculated at unmeasured locations/instants over the Dyle region using a space/time kriging algorithm. In practice, the interpretation of the results of the dynamic stochastic analysis should take into consideration the scale effects.