Changes in the photosynthetic coefficients α and Pm of Planktothrix rubescens grown on light-dark cycles

The relationship between photosynthesis and irradiance was investigated in a strain of the cyanobacterium Planktothrix rubescens isolated from Lake Zürich, using cultures grown on a 12:12-h light-dark cycle. From the photosynthesis-irradiance (P/I) curves, values of the light-affinity coefficient, f, and the maximum rate of photosynthesis, P_m, were determined: in different cultures f ranged from 0.02-0.08 7mol mg^-1 h^-1 (7mol m^-2 s^-1)^-1 and P_m from 0.6-2.6 7mol mg^-1 h^-1 (rates of O₂ production related to dry biomass). Comparisons made at different phases of the light-dark cycle in individual cultures showed that in the first hour of the light phase f rose by an average of 21 % and P_m by 7 %. In cultures that had been grown in many light-dark cycles, however, the rise occurred only if light was given at the beginning of the notional light phase; cultures kept in darkness showed no rise. It is concluded that the observed increases in f and P_m after dawn were in response to the exogenous irradiance rather than due to endogenous rhythms. These findings are discussed in relation to diel variations in the photosynthetic coefficients of natural populations of P. rubescens in Lake Zürich.     

Effects of stream acidification and habitat on fish populations of a North American river

Water quality, physical habitat, and fisheries at sixteen reaches in the Neversink River Basin were studied during 1991-95 to identify the effects of acidic precipitation on stream-water chemistry and on selected fish-species populations, and to test the hypothesis that the degree of stream acidification affected the spatial distribution of each fish-species population. Most sites on the East Branch Neversink were strongly to severely acidified, whereas most sites on the West Branch were minimally to moderately acidified. Mean density of fish populations ranged from 0 to 2.15 fish/m²; biomass ranged from 0 to 17.5 g/m². Where brook trout were present, their population density ranged from 0.04 to 1.09 fish/m², biomass ranged from 0.76 to 12.2 g/m², and condition (K) ranged from 0.94 to 1.07. Regression analyses revealed strong relations (r² - 0.41 to 0.99; p h 0.05) between characteristics of the two most common species (brook trout and slimy sculpin) populations and mean concentrations of inorganic monomeric aluminum (Al_im), pH, Si, K⁺, NO₃^-, NH₄⁺, DOC, Ca²⁺, and Na⁺; acid neutralizing capacity (ANC); and water temperature. Stream acidification may have adversely affected fish populations at most East Branch sites, but in other parts of the Neversink River Basin these effects were masked or mitigated by other physical habitat, geochemical, and biological factors.     

Thermal-mechanical Coupling in Shear Deformation of Viscoelastic Material as a Model of Frictional Constitutive Relations

-- We propose a thermal-mechanical model of shear deformation of a viscoelastic material to describe the temperature-dependence of friction law. We consider shear deformation of one-dimensional layer composed of a Maxwell linear viscoelastic material under a constant velocity V and temperature T_w at the boundary. The strain rate due to viscous deformation depends both on temperature and shear stress. The temperature inside the layer changes owing to frictional heating and conductive cooling. Steady-state calculations show that the sign of d†^ss/dV, where †^ss is steady-state stress, changes from positive to negative as V increases, and that the threshold velocity above which the sign of d†^ss/dV is negative increases with increasing T_w. These results are in accordance with the conjecture that the downdip limit of seismogenic zones is marked by the transition in the sign of d†^ss/dV due to temperature rise with depth. We also find that the response of steady state to a step change in V is quite similar to the response of frictional slip with constitutive laws which employ state variables. These findings suggest that by further improving the present model a model of constitutive relations along faults or plate boundaries can be developed which contains temperature-dependence in a physically-sound manner.     

Viscosity of hydrous Etna basalt: implications for Plinian-style basaltic eruptions

Water dissolved in a silicate melt can strongly influence its physical properties and thus magma behavior during crystallization, degassing, foaming and fragmentation. Etna is a basaltic volcano whose activity is dominated by effusive eruptions which have long represented a threat to the densely populated, surrounding area. Recently, recognition of the products of a Plinian eruption (122 B.C.) has raised further issues for hazard assessment at Etna and other basaltic volcanoes. Constraining the behavior of Etna magma under conditions relevant to both effusive and explosive hazards requires viscosity data under conditions near the glass transition. Here we have investigated the viscosity of hydrous Etna lava in order to better understand eruptive processes which characterize this volcano. The experimental methods which have been used include piston cylinder synthesis of the hydrated melts, micropenetration viscometry for low-temperature viscosity measurements, and near-infrared spectroscopy for the evaluation of sample homogeneity and measurements of water content. Additionally, scanning calorimetric determinations were performed to check whether incipient crystallization had occurred. Sample compositions were determined using electron microprobe analysis and ⁵⁷Fe Mössbauer spectroscopy. Results from this study are compared with previous reports of trachytic, phonolitic and model calc-alkaline rhyolite (HPG8) compositions. The viscosity of the basaltic melt (dry and wet) has been parameterized in terms of temperature and water content via the non-Arrhenian equation: log₁₀-=-4.643+(5,812.44-427.042H₂O)/(T(K)-499.31+28.742ln(H₂O)) where - is the viscosity in Pa s, H₂O is the water content in wt%, and T is the temperature in Kelvin. We observe that the viscosity of alkali basalt (at more than 0.5 wt% H₂O) is similar to that of an alkaline trachyte (Agnano-Monte Spina eruption, Phlegrean Fields) and much higher than that of a peralkaline phonolite (Teide, Tenerife) at similar silica contents and NBO/T. For water contents above 1.5 wt%, the viscosity of the basalt is similar to that of rhyolitic melts with similar water contents. At temperatures ranging from 1,050 to 1,150 °C and with water contents between 0.5 and 2.3 wt% (eruptive conditions), the viscosities calculated using the equation defined in this study are (1) in reasonable agreement with those calculated using Shaw's model, and (2) much lower than those experimentally determined in a previous study. However, outside these temperature and water content ranges, the agreement with Shaw's model (1972) breaks down.     

Submarine silicic volcanism of the Healy caldera,southern Kermadec arc (SW Pacific): I – volcanology and eruption mechanisms

The submarine Healy volcano (southern Kermadec arc), with a 2-2.5 km wide caldera, is pervasively mantled with highly vesicular silicic pumice within a water depth of 1,150-1,800 m. Pumices comprise type 1 white-light grey pumice with ⢾ mm vesicles and weak-moderate foliation, type 2 grey pumice with millimetre-scale laminae, flow banded foliation, including stretched vesicles ⣗ mm in length, and a minor finely vesicular type 3 pumice. All types are sparsely porphyritic, with undevitrified glassy groundmass (68-70% SiO₂), which is microlite and lithic free. Coexisting pyroxenes yield magma temperatures of ~950 °C. Pumice density is А.5 g cm^-3 and vesicularity is 78-83%. Vesicle size distributions for types 1 and 2 pumice, range from ~20 µm to >20 mm, with a strong power-law relation (with d=-2.5ǂ.4) for vesicles <1-2 mm. Larger vesicles have variable size modes. The vesicle size distribution and packing indicates rapid magma decompression and ascent. Consideration of the pressure dependent, solubility of H₂O at a magma temperature of 𙧶 °C and water content of Ж wt%, with pumice petrography and vesicle granulometry, strongly suggests a pyroclastic eruption. Reconstructions of the submarine edifice between water depths of 1,000 and 550 m constrain the ambient hydrostatic pressure to ~6-9 MPa. Pressures >~9 MPa will limit vesicularity to less than the observed 78-83%, whereas pressure <~6 MPa require a more shallower reconstruction of the edifice and larger-volume syn-eruptive collapse. Uniformly high vesicularity is interpreted as evidence of insulation within an eruption column comprising steam and hot pyroclasts. Most pyroclasts cool, condensing and ingesting water into steam-inflated vesicles, and then sink. Progression into pyroclastic mode would expand the eruption column, displace ambient water, reduce the hydrostatic load, and further promote vesiculation and fragmentation. Pyroclasts within the column would quench at these reduced pressures. We argue that Healy eruptions deeper than ~1,000 m cannot be pyroclastic. Volumes for the lower and upper bounds of edifice size are 2.36 and 3.58 km³, respectively, but do not account for intra-caldera pumice fill. These volumes are considered to be predominantly primary eruption output, as shown by a dearth of accessory lithics in all pumice, yielding (at an average 81% vesicularity) eruptive pumice volumes of between 10 and 15 km³. Some pyroclasts may have risen to the sea surface and be a correlative of the sea-rafted Loisels pumice; the latter occurs in some New Zealand Holocene beach sequences and has a estimated age of 590ᇤ calendar years.     

Geophysical Applications of Multidimensional Filtering with Wavelets

--We present imaging results in geophysics based on using multidimensional Gaussian wavelets as a filter in a 2-D Cartesian domain. Besides decomposing the field into various distinct lengthscales, we have also constructed the 2-D maps describing the spatial distributions of the maximum of the wavelet-transformed L₂-norm E_max (x,y) and its corresponding local wavenumber k_max (x,y), where x and y are the Cartesian coordinates. For geoid anomalies, using a wavelet filter extending to 90 degrees, we have discerned the distinct outlines of convergent and divergent tectonic zones and have conducted a quantitative comparison of the short-wavelength gravitational anomalies at those wavelengths between two different geographical locations. We have also compared the wavelet results with a nonlinear bandpass filter in the spectral domain where a Gaussian filter with the logarithm of the degree l acting as the argument has been employed. A wavelet solution, with a length-scale corresponding to 256 degrees, would need a filter with over 400 spherical harmonics centering around l=157 for an optimal spatial fit. The computational effort with the bandpass filter technique greatly exceeds those associated with wavelets. We have also shown the ability of the wavelets to analyze the vastly different scales present in high Rayleigh number convection and the mixing of passive heterogeneities driven by thermal convection. Wavelets will be a useful tool for rapid analyzing of the large multidimensional fields to be captured in many other geophysical endeavors, such as the upcoming gravity satellite missions and satellite radar interferometry images.     

Pyroclastic surges and flows from the 8–10 May 1997 explosive eruption of Bezymianny volcano, Kamchatka, Russia

The 8-10 May 1997 eruption of Bezymianny volcano began with extrusion of a crystallized plug from the vent in the upper part of the dome. Progressive gravitational collapses of the plug caused decompression of highly crystalline magma in the upper conduit, leading at 13:12 local time on 9 May to a powerful, vertical Vulcanian explosion. The dense pyroclastic mixture collapsed in boil-over style to generate a pyroclastic surge which was focused toward the southeast by the steep-walled, 1956 horseshoe-shaped crater. This surge, with a temperature <200 °C, covered an elliptical area >30 km² with deposits as much as 30 cm thick and extending 7 km from the vent. The surge deposits comprised massive to vaguely laminated, gravelly sand (Md -1.2 to 3.7J sorting 1.2 to 3J) of poorly vesiculated andesite (mean density 1.82 g cm^-3; vesicularity 30 vol%; SiO₂ content ~58.0 wt%). The deposits, with a volume of 5-15᎒⁶ m³, became finer grained and better sorted with distance; the maximal diameter of juvenile clasts decreased from 46 to 4 cm. The transport and deposition of the surge over a snowy landscape generated extensive lahars which traveled >30 km. Immediately following the surge, semi-vesiculated block-and-ash flows were emplaced as far as 4.7 km from the vent. Over time the juvenile lava in clasts of these flows became progressively less crystallized, apparently more silicic (59.0 to 59.9 wt% SiO₂) and more vesiculated (density 1.64 to 1.12 g cm^-3; vesicularity 37 to 57 vol%). At this stage the eruption showed transitional behavior, with mass divided between collapsing fountain and buoyant column. The youngest pumice-and-ash flows were accompanied by a sustained sub-Plinian eruption column ~14 km high, from which platy fallout clasts were deposited (~59.7% SiO₂; density 1.09 g cm^-3; vesicularity 58 vol%). The explosive activity lasted about 37 min and produced a total of ~0.026 km³ dense rock equivalent of magma, with an average discharge of ~1.2᎒⁴ m³ s^-1. A lava flow ~200 m long terminated the eruption. The evolutionary succession of different eruptive styles during the explosive eruption was caused by vertical gradients in crystallization and volatile content of the conduit magma, which produced significant changes in the properties of the erupting mixture.     

The solubility of iron sulphides in synthetic and natural waters at ambient temperature

A critical evaluation of literature values for the solubility products, K _sp ^NBS = [Fe²⁺][HS^–] Fe²⁺ HS^– (H _NBS ⁺ )^–1, of various iron sulphide phases results in consensus values for the pKs of 2.95 ± 0.1 for amorphous ferrous sulphide, 3.6 ± 0.2 for mackinawite, 4.4 ± 0.1 for greigite, 5.1 ± 0.1 for pyrrhotite, 5.25 ± 0.2 for troilite and 16.4 ± 1.2 for pyrite.Where the analogous ion activity products have been measured in anoxic freshwaters in which there is evidence for the presence of solid phase FeS, the values lie within the range of 2.6–3.22, indicating that amorphous iron sulphide is the controlling phase. The single value for a groundwater of 2.65 (2.98 considering carbonate complexation) agrees. In seawater four values range between 3.85 to 4.2, indicating that mackinawite or greigite may be the controlling phase. The single low value of 2.94 is in a situation where particularly high fluxes of Fe (II) and S (–II) may result in the preferential precipitation of amorphous iron sulphide. Formation of framboidal pyrite in these sulphidic environments may occur in micro-niches and does not appear to influence bulk concentrations. Calculations show that the formation of Fe₂S₂ species probably accounts for very little of the iron or sulphide in most natural waters. Previously reported stability constants for the formation of Fe (HS)₂ and (Fe (HS)₃)^– are shown to be suspect, and these species are also thought to be negligible in natural waters. In completely anoxic pore waters polysulphides also have a negligible effect on speciation, but in tidal sediments they may reach appreciable concentrations and lead to the direct formation of pyrite. Concentrations of iron and sulphide in pore waters can be controlled by the more soluble iron sulphide phase. The change in the IAP with depth within the sediment may reflect ageing of the solid phase or a greater flux of Fe (II) and S (–II) nearer the sediment surface. This possible kinetic influence on the value of IAPs has implications for their use in geochemical studies involving phase formation.     

The Problem of Regionalization in the Study of Magnitude-intensity Correspondences in the Mediterranean Area: A Distribution-free Approach

--A distribution-free statistical approach is proposed for tackling the problem of regionalization devoted to the study of magnitude-intensity relationships in the Mediterranean area. The training-set constitutes a compilation of more than 2000 earthquakes which occurred in the Mediterranean region since the end of the nineteenth century, when instrumental data became available, with epicentral or maximum intensity S VI and M_S magnitude values. Firstly an empirical magnitude-intensity correspondence has been computed for each intensity class by using the entire data set in the assumption of homogeneity at the regional scale. Residuals of such relation have been analyzed by a distribution-free statistical approach in order to evaluate the opportunity of a regionalization able to locally improve the performances of magnitude-intensity relationships. The analysis indicates that data concerning larger earthquakes (intensity S VII) do not suggest the opportunity of zonation, and that unbiased estimates of macroseismic magnitude can also be obtained in the assumption that magnitude-intensity correspondences are uniform over the entire Mediterranean area. Therefore, better constrained relations determined for the entire Mediterranean region should be preferred to ill-defined local ones. As concerns smaller events (intensity VI), the procedure suggests that medium/small-scale lateral variations (on a wavelength lower than 10² km) should be taken into account if an efficient estimate of magnitudes from maximum observed intensity is pursued, but that data presently available are not sufficient to suggest any reliable zonation of the area under study.     

Characteristics of Cloud Drop Spectra at Different Levels and with Respect to Cloud Thickness

--Measurements on drop size spectra were made in cumulus clouds over Pune (inland) region on many days during the summer monsoon seasons. In this paper, the measurements in non-raining cumulus clouds made in the years 1984, 1985 and 1986 at different levels and for different cloud thickness have been studied. In general, the drop size spectra broadened with height and the concentration of drops with diameter > 50 wm (N_L), mean volume diameter (MVD), liquid water content (LWC) and dispersion increased with height while the concentration of drops with diameter < 20 wm (N_S) and the total concentration of drops (N_T) decreased with height. The average drop size distributions were unimodal at the lower levels while they were bimodal at the higher levels. High water contents were confined to drops in the size range 5-25 wm at both higher and lower levels. The average drop size spectra were broader and N_L, LWC, MVD and dispersion greater while N_T and N_S smaller for thicker clouds (range of vertical extent 1.1-2.1 km) as compared to those for thinner clouds (range of vertical extent 0.3-1.1 km). Water contents for the drops > 28 wm were higher while those for the drops > 28 wm lower in thicker clouds than in thinner clouds. The average drop size distributions were bimodal in the former case, while they were unimodal in the other case.     

Relocation of Swarm Events in the June 1987 Earthquake Sequence from the Lake Aswan Area in Egypt

--We have examined the digital waveform data and relocated a number of events within the June 1987 earthquake swarm, which occurred beneath the northern part of Lake Aswan, 70 km southwest of the Aswan High Dam in Egypt. This swarm occurred between June 17th and 19th with a maximum magnitude event of "M_L"=3.5.¶Cross correlation between a chosen master and the analyzed events has been carried out on seismograms from stations of the Aswan network. The cross correlation demonstrates the presence of a difference in both the P wave ((t_p) and the S wave ((t_s) arrival times at each station in the network relative to the arrival times of the master event at the same stations. (t_p ranges between т.15 and 0.11 second, while (t_s ranges between т.17 and 0.11 second.¶The primary interpretation is that the se time differences represent an error in the manually picking arrival times. Then, (t_p and (t_s values for each event result from a change in the hypocentral parameters from those of the master event, assuming the P- and S-wave velocity distribution remains constant during the swarm activity. This interpretation enables us to determine the relative distribution of hypocenters with respect to the hypocentral location of the master event. We present the results from a swarm of 9 events demonstrating they originate from a nearly unique location, rather than the zone identified from the preliminary locations which used manually picked onset times.     

Biochemical characteristics of alga-bacterial mats and invertebrates from shallow-water hydrothermal fields of the West Pacific Ocean

In this study we assessed the biochemical adaptations of biotic specimens from the hydrothermal vent fields from different geographical zones. Algabacterial and bacterial mats and specimens of bivalves (9 species) and sea urchins (2 species) from different coastal shallow-water fields of hydrothermal activity along the West Rift zones of the Pacific Ocean were sampled to compare biochemical activities to volcanogenic chemical characteristics. The algabacterial and bacterial mat cells and subcellular fractions of bivalve gills and sea urchin guts were measured for 5-aminolevulinic acid (ALA), cytochromes b, (c + c₁), (a + a₃) and P450 (CYP), and activities of ALA synthetase and CYP-dependent aminopyrine N-demethylase (APND). The algabacterial mats from the Kuril Islands and New Zealand displayed comparable ALA and cytochrome levels and enzyme activities. Similar ALA synthesis, APND activity, and levels of cytochromes b and (c + c₁) were observed in bacterial mats from hydrothermal fields from Kuriles, New Britain Island and New Zealand. More active ALA synthesis, APND activity and enhanced content of cytochromes b and (c + c₁) and lower levels of CYP were found in bacterial mats than in algabacterial mats. Further, all bivalves influenced by hydrothermal discharges had elevated CYP levels and APND activity as well as a significant increase in ALA content and ALA synthesis (P < 0.05). Statistically sufficient (P < 0.05) alterations in ALA synthesis, APND activity and total CYP level were found in sea urchins under hydrothermal influence. Our findings contribute to the global assessment of hydrothermal effluents on biota and indicate that living conditions near hydrothermal vents accelerate metabolism of bivalves and sea urchins in examined hydrothermal zones in the northern and southern West Pacific Ocean.     

Observation of Double-mixing Line Structure in the Convective Boundary Layer During the Summer Monsoon Season

v--vThe special aerological observations carried out as a part of Land Surface Processes Experiment (LASPEX) were used to investigate the thermodynamic structure of the convective boundary layer during the summer monsoon of 1997. The analysis suggested that the convective boundary layer top was found at 700 hPa which was associated with Š_e minimum and Š_es maximum values. Double-mixing line structure was noticed in the conserved variable diagrams which was possibly attributed to the radiative warming/evaporation of falling precipitation.     

Factors controlling permeability–porosity relationships in magma

A model has been developed to investigate the sensitivity of magma permeability, k, to various parameters. Power-law relationships between k and porosity J are revealed, in agreement with previous experimental and theoretical studies. These relationships take the form % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavTnhis1MBaeXatLxBI9gBae % rbd9wDYLwzYbWexLMBbXgBcf2CPn2qVrwzqf2zLnharyavP1wzZbIt % LDhis9wBH5garqqtubsr4rNCHbGeaGqiVCI8FfYJH8sipiYdHaVhbb % f9v8qqaqFr0xc9pk0xbba9q8WqFfeaY-biLkVcLq-JHqpepeea0-as % 0Fb9pgeaYRXxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaabau % aaaOqaaiqbdUgaRzaajaGaeyypa0Jaem4AaSMaei4la8IaemOCai3a % aWbaaSqabeaacqaIYaGmaaGccqGH9aqpcqWGHbqycqGGOaakcqaHgp % GzcqGHsislcqaHgpGzdaWgaaWcbaGaem4yamMaemOCaihabeaakiab % cMcaPmaaCaaaleqabaGaemOyaigaaaaa!4CE4! [^(k)] = k/r² = a(f- f_cr )^b\hat k = k/r^2 = a(\phi - \phi _{cr} )^b where r is the mean bubble radius, J_cr is the percolation threshold below which permeability is zero, and a and b are constants. It is discovered that % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavTnhis1MBaeXatLxBI9gBae % rbd9wDYLwzYbWexLMBbXgBcf2CPn2qVrwzqf2zLnharyavP1wzZbIt % LDhis9wBH5garqqtubsr4rNCHbGeaGqiVCI8FfYJH8sipiYdHaVhbb % f9v8qqaqFr0xc9pk0xbba9q8WqFfeaY-biLkVcLq-JHqpepeea0-as % 0Fb9pgeaYRXxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaabau % aaaOqaaiqbdUgaRzaajaGaeyOeI0IaeqOXdygaaa!3CDB! [^(k)] - f\hat k - \phi relationships are independent of bubble size. The percolation threshold was found to lie at around 30% porosity. Polydisperse bubble-size distributions (BSDs) give permeabilities around an order of magnitude greater than monodisperse distributions at the same porosity. If bubbles are elongated in a preferred direction then permeability in this direction is increased, but, perpendicular to this direction, permeability is unaffected. In crystal-free melts the greatest control on permeability is the ease of bubble coalescence. In viscous magmas, or when the cooling time-scale is short, bubble coalescence is impeded and permeability is much reduced. This last effect can cause variations in permeability of several orders of magnitude.     

The Skin Effect of the Anelastic Earth

--Assuming that the rheology of the earth's crust is defined by stress-strain relations containing a memory mechanism, it is seen that the stress field and the maximum shear stress (mss), caused by elongated mountain ranges, are limited mostly to the crust. In a geological, relatively short time, where relaxation is reached, the mss assumes the depth distribution (2~x)exp(-2~x) where x is the depth h measured in units of wavelength w of the mountain range; the mss reaches its maximum value, of about 1/3 the surface load, at the depth x = 1/2~ (x = h/w). Contrary to the case of the perfectly elastic earth model, where the mss at the surface is about 0.25 the maximum surface load, in the case of the anelastic earth model the mss is asymptotically nil at the surface and for depths larger than w, the mss is less than 1.2% the maximum surface load. It is seen that the relaxation of the medium causes a rotation of the planes of the mss in a layer near the surface and that, in the case when the surface load is applied periodically, there is a phase shift between the surface load and the stress at depth which is less than four days when the period is one year; it is smaller for shorter periods. It is also seen that the creation of reservoirs lasting a finite time causes a small variation of the mss in the crust which, when the reservoir is removed, leaves a residual mss with maximum at the surface and decreasing exponentially with depth; for a load lasting one year it takes 250 days to reduce to 40% the stress at depth. As a consequence of the residual mss, caused by the load of the topography, after earthquakes the local mss will not drop to zero; this decreases the time required by tectonic forces to release the next earthquake and increases the rate of seismicity.     

Earthquake Hazard Parameters in Crete Island and its Surrounding Area Inferred from Bayes Statistics: An Integration of Morphology of the Seismically Active Structures and Seismological Data

-- The study addresses the evaluation of earthquake hazard parameters such as maximum regional magnitude (M_max) and the slope of Gutenberg-Richter law # (where b=# log e) for the Hellenic Wadati-Benioff zone and the overriding lithospheric plate in the area of Crete and its surroundings. The seismicity of the area is divided in a cellular (1.0° 2 1.0°) manner allowing analysis of the localized earthquake hazard parameters and graphical representation of their spatial variation. Our approach incorporates the recently updated earthquake catalogue for Greece and the adjacent areas, the consideration of the morphology of the deep seismically active structures in the studied area and use of a probabilistic procedure for estimating the earthquake hazard parameters.¶One of the main inconsistencies in the earthquake hazard assessment is the estimation of the maximum magnitude and the related uncertaint y. The Bayesian approach, applied in the present, is a straightforward technique for evaluating the earthquake hazard parameters and is based on the following assumptions: Poissonian character of seismic events flow, a frequency-magnitude law of Gutenberg-Richter's type with cutoff maximal value for estimated parameter and a seismic catalogue, having a rather sizeable number of events (i.e., 50 events at least per cell). For five cells in which the number of events is less than 50, an effort is made to produce synthetic data. The re-assessed parameters obtained from the synthetic data show no significant difference and the real data (of the five cells) are finally taken into account although the estimated uncertainty is high.¶For four random cells we constructed hazard curves showing the probabilities that a certain magnitude M will be exceeded in one year and the return periods (in years) that are expected for a given magnitude. These are particularly useful for the mapping of earthquake hazard in regions of either low or high seismic activity, as is Crete and the adjacent area.¶The obtained results show that the W and E parts of both subducting and overriding plates differ in the spatial distribution of all the estimated earthquake hazard parameters. The M_max distribution indicates strong coupling between the western portions of the interacting plates (M_max > 6.3) to the south of 36°N. The smaller values of M_max (M_max < 6.3) estimated in the SE part of the studied area indicate weak coupling between the eastern portions of the subducting and overriding plates.¶Values of b > 1.0 are found to the south and east of Crete for the Wadati-Benioff zone, and over the central part of the island and the area to the northeast of it (cell 11) for the continental wedge, which suggests nonuniform stress field and/or heterogeneous material.     

Optimized Seismic Threshold Monitoring – Part 2: Teleseismic Processing

v--vThis second paper (Part 2) pertaining to optimized site-specific threshold monitoring addresses the application of the method to regions covered by a teleseismic or a combined regional-teleseismic network. In the first paper (Part 1) we developed the method for the general case, and demonstrated its application to an area well-covered by a regional network (the Novaya Zemlya nuclear test site). In the present paper, we apply the method to the Indian and Pakistani nuclear test sites, and show results during the periods of nuclear testing by these two countries in May 1998. Since the coverage by regional stations in these areas is poor, an optimized approach requires the use of selected, high-quality stations at teleseismic distances.¶To optimize the threshold monitoring of these test sites, we use as calibration events either one of the nuclear explosions or a nearby earthquake. From analysis of the calibration events we derive values for array beamforming steering delays, filter bands, short-term averages (STA) lengths, phase travel times (P waves), and amplitude-magnitude relationships for each station. By applying these parameters, we obtain a monitoring capability of both test sites ranging from m_b 2.8-3.0 using teleseismic stations only. When including the nearby Nilore station to monitor the Indian tests, we show that the threshold can be reduced by about 0.4 magnitude units. In particular, we demonstrate that the Indian tests on 13 May, 1998, which were not detected by any known seismic station, must have corresponded to a magnitude (m_b) of less than 2.4.¶We also discuss the effect of a nearby aftershock sequence on the monitoring capability for the Pakistani test sites. Such an aftershock sequence occurred in fact on the day of the last Pakistani test (30 May, 1998), following a large (m_b 5.5) earthquake in Afghanistan located about 1100 km from the test site. We show that the threshold monitoring technique has sufficient resolution to suppress the signals from these interfering aftershocks without significantly affecting the true peak of the nuclear explosion on the threshold trace.     

Optimized Seismic Threshold Monitoring – Part 1: Regional Processing

v--v Continuous seismic threshold monitoring is a technique that has been developed over the past several years to assess the upper magnitude limit of possible seismic events that might have occurred in a geographical target area. The method provides continuous time monitoring at a given confidence level, and can be applied in a site-specific, regional or global context.¶In this paper (Part 1) and a companion paper (Part 2) we address the problem of optimizing the site-specific approach in order to achieve the highest possible automatic monitoring capability of particularly interesting areas. The present paper addresses the application of the method to cases where a regional monitoring network is available. We have in particular analyzed events from the region around the Novaya Zemlya nuclear test site to develop a set of optimized processing parameters for the arrays SPITS, ARCES, FINES, and NORES. From analysis of the calibration events we have derived values for beam-forming steering delays, filter bands, short-term average (STA) lengths, phase travel times (P and S waves), and amplitude-magnitude relationships for each array. By using these parameters for threshold monitoring of the Novaya Zemlya testing area, we obtain a monitoring capability varying between m_b 2.0 and 2.5 during normal noise conditions.¶The advantage of using a network, rather than a single station or array, for monitoring purposes becomes particularly evident during intervals with high global seismic activity (aftershock sequences), high seismic noise levels (wind, water waves, ice cracks) or station outages. For the time period November-December 1997, all time intervals with network magnitude thresholds exceeding m_b 2.5 were visually analyzed, and we found that all of these threshold peaks could be explained by teleseismic, regional, or local signals from events outside the Novaya Zemlya testing area. We could therefore conclude within the confidence level provided by the method, that no seismic event of magnitude exceeding 2.5 occurred at the Novaya Zemlya test site during this two-month time interval.¶As an example of particular interest in a monitoring context, we apply optimized threshold processing of the SPITS array for a time interval around 16 August 1997 m_b 3.5 event in the Kara Sea. We show that this processing enables us to detect a second, smaller event from the same site (m_b 2.6), occurring about 4 hours later. This second event was not defined automatically by standard processing.     

Application of Regional Phase Amplitude Tomography to Seismic Verification

v--vWe have applied tomographic techniques to amplitude data to quantify regional phase path effects for use in source discrimination studies. Tomography complements interpolation (kriging) methods by extending our ability to apply path corrections into regions devoid of calibration events, as well as raising levels of confidence in the corrections because of their more physical basis. Our tomography technique solves for resolvable combinations of attenuation, source-generation, site and spreading terms. First difference regularization is used to remove singularities and reduce noise effects.¶In initial tests the technique was applied to a data set of 1488, 1.0 Hz, P_g/L_g amplitude ratios from 13 stations for paths inside a 30° by 40° box covering western China and surrounding regions. Tomography reduced variance 60%, relative to the power-law distance correction traditionally applied to amplitude ratios. Relative P_g/L_g attenuation varied with geologic region, with low values in Tibet, intermediate values in basins and high values for platforms and older crust. Spatial patterns were consistent with previous path effect studies in Asia, especially local earthquake coda-Q. Relative spreading was consistent with expected values for P_g and L_g. Relative site terms were similar to one another, yet some tradeoff with attenuation was evident.¶Tomography residuals followed systematic trends with distance, which may result from the evolution from direct to coda phases, focusing, model tradeoff or data windowing effects. Examination of the residuals using a kriging interpolator showed coherent geographical variations, indicating unmodeled path effects. The residual patterns often follow geological boundaries, which could result from attenuating zones or minor blockages that are too thin to be resolved, or that have anisotropic effect on regional phases. These results will guide efforts to reparameterize tomography models to more effectively represent regional wave attenuation and blockage. The interpolated residuals also can be combined with predictions of the tomographic model to account for path effects in discrimination studies on a station by station basis.     

Critical loads of acidity for Irish lakes

Critical loads of acidity have been estimated for Irish lakes (277) using survey data from the Republic of Ireland and Northern Ireland, and compatible inputs and methods across both regions. This is the first trans-region application of the Steady-State Water Chemistry model for the island of Ireland. Results indicate that a small percentage of lakes (10.8%, 30 lakes) have low critical loads of acidity (<50.0 mmol_c m^-2 yr^-1). Based on bulk sulphate deposition, exceedance of critical load was observed in 6.9% of the sampled lakes (19 lakes). Inferences relating to the whole lake population cannot be made, as statistical lake surveys were not carried out. However, as the surveys were weighted towards acid-sensitive lakes, the current study suggests that Irish lakes are not under a significant threat from anthropogenic acidification. The acidifying impact of nitrogen was incorporated into the exceedance calcula tion by superimposing the present leaching of nitrogen onto sulphate deposition. This increased exceedance by ഌ% (i.e., from 11 to 19 lakes). Despite the potential uncertainties associated with using default values and relationships in the Steady-State Water Chemistry model, uncertainty analysis indicated that the critical loads were relatively stable. The critical acid neutralising capacity limit (specified water quality criterion) had the greatest effect on the critical load. In the current study, a catchment dependent limit was employed (ranged: 0-50 7mol_c L^-1; mean: 35 7mol_c L^-1) rather than a single fixed value.