Microearthquake seismicity and active tectonics of northwestern Greece

We carried out a microearthquake survey lasting for six weeks in northwest Greece using 18 portable seismograph stations to examine a region in which normal and thrust faulting have been reported in close proximity to one another. With this array we located 148 events and determined fault plane solutions for eight events using only rays radiated upwards. The seismicity of the region is diffuse with events extending to depths of nearly 30 km, and there is a minimum in activity near a depth of 15 km. The fault plane solutions exhibit a wide spectrum of fault types and orientations and are not consistent with simple zones of shortening or extension. Neither tractions applied to the edges or bottom of the region nor deviatoric stresses that compensate for lateral variations in crustal thickness can account for the variety of fault plane solutions. We think that the complicated behavior is a manifestation of inhomogeneous deformation due, at least in part, to a pre-existing complicated juxtaposition of structures and formations.     

Numerical simulations of katabatic jumps in coats land,Antartica

A non-hydrostatic numerical model, the Regional Atmospheric Modeling System (RAMS), has been used to investigate the development of katabatic jumps in Coats Land, Antarctica. In the control run with a 5 m s^-1downslope directed initial wind, a katabatic jump develops near the foot of the idealized slope. The jump is manifested as a rapid deceleration of the downslope flow and a change from supercritical to subcritical flow, in a hydraulic sense, i.e., the Froude number (Fr) of the flow changes from Fr > 1 to Fr> 1. Results from sensitivity experiments show that an increase in the upstream flow rate strengthens the jump, while an increase in the downstream inversion-layer depth results in a retreat of the jump. Hydraulic theory and Bernoulli's theorem have been used to explain the surface pressure change across the jump. It is found that hydraulic theory always underestimates the surface pressure change, while Bernoulli's theorem provides a satisfactory estimation. An analysis of the downs balance for the katabatic jump indicates that the important forces are those related to the pressure gradient, advection and, to a lesser extent, the turbulent momentum divergence. The development of katabatic jumps can be divided into two phases. In phase I, the t gradient force is nearly balanced by advection, while in phase II, the pressure gradient force is counterbalanced by turbulent momentum divergence. The upslope pressure gradient force associated with a pool of cold air over the ice shelf facilitates the formation of the katabatic jump.     

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M _w= 9.2) and 2005 March ( M _w= 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.     

Historic Landscape Reconstruction and Visualisation, West Oxfordshire, England

As part of the Wychwood Project in West Oxfordshire, England, aerial photography and historic maps are being used to reconstruct landscape change in the former medieval Royal Hunting Forest of Wychwood. An important objective of the project is to stimulate local interest in Wychwood and promote understanding of its evolution, thereby encouraging the support of stakeholders in the design of future landscapes. The parish of Leafield was selected to develop and test a method for reconstructing and visualising past and present landscapes to present to the public for evaluation and comparison. Leafield contains fine examples of both early and late enclosure landscapes and underwent large-scale disafforestation following the Enclosure Act of 1856. Ordnance Survey (O.S.) 1:10000 Landform Panorama contour data were used to create a digital elevation model (DEM) of the parish using ArcView 3D Analyst, enabling the user to create perspective views of the landscape from any height and viewing angle. Contemporary aerial photography was scanned and geo-rectified to create an overlay for 'draping' onto the DEM. Historical maps provided the reference for modifying the aerial photography using the PaintShop Pro graphics package. This paper presents and reviews the results and compares the method with other techniques, including VRML.     

Road Reduction Filtering for GPS-GIS Navigation

A novel method of map matching using the Global Positioning System (GPS) has been developed for civilian use, which uses digital mapping data to infer the systematic position errors of less than 100m which result largely from 'selective availability' (S/A) imposed by the US military. Selective availability was switched off on the 2nd of May 2000, and is to be replaced with 'regional denial capabilities in lieu of global degradation'. The system tracks a vehicle on all possible roads (road centre-lines) in a computed error region, then uses a method of rapidly detecting inappropriate road centre-lines from the set of all those possible. This is called the Road Reduction Filter (RRF) algorithm. Point positioning is computed using C/A code pseudorange measurements direct from a GPS receiver. The least squares estimation is performed in the software developed for the experiment described in this paper. Virtual differential GPS (VDGPS) corrections are computed and used from a vehicle's previous positions, thus providing an autonomous alternative to DGPS for in-car navigation and fleet management. Height aiding is used to augment the solution and reduce the number of satellites required for a position solution. Ordnance Survey (OS) digital map data was used for the experiment, i.e. OSCAR 1 m resolution road centre-line geometry and Land Form PANORAMA 1:50,000, 50 m-grid digital terrain model (DTM). Testing of the algorithm is reported and results are analysed. Vehicle positions provided by RRF are compared with the 'true' position determined using high precision (cm) GPS carrier phase techniques. It is shown that height aiding using a DTM and the RRF significantly improve the accuracy of position provided by inexpensive single frequency GPS receivers.     

Analysis And Simulation Of Surface-Layer Winds Using Multiplicative Cascade Models With Self-Similar Probability Densities

Statistical analysis techniques based on multiplicative cascades are investigated for use with surface-layer wind data sets collected in the atmospheric boundary layer over flat farm land. The data were found to exhibit multiscaling statistics, allowing the surface-layer winds to be simulated with the use of multiplicative random cascades. The study found evidence that, for the surface-layer at least, these cascade models (andhence the methods of multifractal analysis) should be applied in separate ways to the microscale inertial range, and the mesoscale. This is at odds with the view found in the existing literature, which proposes a `universal multifractal' model to replace the widely held view that there exists separate microscale, mesoscale and synoptic scales for which the processes governing each are different. At least two separate ranges of scaling are suggested for surface-layer wind data, corresponding to the microscale inertial range and the mesoscale. For the case of the mesoscale range, a self-similar distribution of weighting factors was found for the wind speed data themselves, rather than for an intermediate (dissipation) field, as is required for themicroscale data.     

Helium isotope geochemistry of mid-ocean ridge basalts from the South Atlantic

We report new helium isotope results for 49 basalt glass samples from the Mid-Atlantic Ridge between 1°N and 47°S.³He/⁴He in South Atlantic mid-ocean ridge basalts (MORB) varies between 6.5 and 9.0 R_A (R_A is the atmospheric ratio of1.39 × 10⁻⁶), encompassing the range of previously reported values for MORB erupted away from high³He/⁴He hotspots such as Iceland. He, Sr and Pb isotopes show systematic relationships along the ridge axis. The ridge axis is segmented with respect to geochemical variations, and local spike-like anomalies in³He/⁴He, Pb and Sr isotopes, and trace element ratios such as(La/Sm)_N are prevalent at the latitudes of the islands of St. Helena, Tristan da Cunha and Gough to the east of the ridge. The isotope systematics are consistent with injection beneath the ridge of mantle “blobs” enriched in radiogenic He, Pb and Sr, derived from off-axis hotspot sources. The variability in³He/⁴He along the ridge can be used to refine the hotspot source-migrating-ridge sink model.

MORB from the 2–7°S segment are systematically the least radiogenic samples found along the mid-ocean ridge system to date. Here the depleted mantle source is characterized by⁸⁷Sr/⁸⁶Sr of 0.7022, Pb isotopes close to the geochron and with²⁰⁶Pb/²⁰⁴Pb of 17.7, and³He/⁴He of 8.6–8.9 R_A. The “background contamination” of the subridge mantle, by radiogenic helium derived from off-ridge hotspots, displays a maximum between 20 and 24°S. The HePb and HeSr isotope relations along the ridge indicate that the³He/⁴He ratios are lower for the hotspot sources of St. Helena, Tristan da Cunha and Gough than for the MORB source, consistent with direct measurements of³He/⁴He ratios in the island lavas. Details of the HeSrPb isotope systematics between 12 and 22°S are consistent with early, widespread dispersion of the St. Helena plume into the asthenosphere, probably during flattening of the plume head beneath the thick lithosphere prior to continental breakup. The geographical variation in theHe/Pbratio deduced from the isotope systematics suggests only minor degassing of the plume during this stage. Subsequently, it appears that the plume component reaching the mid-Atlantic ridge was partially outgassed of He during off-ridge hotspot volcanism and related melting activity.

Overall, the similar behavior of He and Pb isotopes along the ridge indicates that the respective mantle sources have evolved under conditions which produced related He and Pb isotope variations.     

Experimental stream acidification — the influence of sediment and streambed vegetation

An acidification experiment was conducted on a small stream in the Loch Ard area of central Scotland. The stream was chosen because of its large, flow related, variation in pH (5.9-4.0). Two acid additions were made to approximately pH 3.5–3.7. The results indicated a strong correlation between labile aluminium and hydrogen, and a noticeable hysteresis in the response of calcium and hydrogen. It is hypothesised that divalent cation response is a result of ion-exchange mechanisms involving the streambed vegetation, with aluminium release resulting primarily from exchange reactions with streambed sediment stores. Data from a program of stream spot sampling have been analysed in an attempt to elucidate the contribution of different sources of aluminium under different flow conditions. Streambed sources of aluminium contribute significantly under low flow conditions; however, at high flow, additional sources of aluminium must contribute to match observed streamwater chemistry.