On the Generation of ULF Magnetic Variations by Conductivity Fluctuations in a Fault Zone

v--vPrior to the October 18, 1989 Loma Prieta M_s 7.1 earthquake, Fraser-Smith et al. (1990) recorded a 10-100 fold increase in ultra-low frequency (ULF) magnetic fields near the earthquake epicenter. Several mechanisms for generation of these ULF fields by fluid flow in the earth have been advanced, but all appear to require unrealistic fluid velocities or hydraulic permeabilities to match the observations. As an alternative explanation, Merzer and Klemperer (1998) proposed that the increase in ULF magnetic fields could result from induced electric currents flowing in a fault-zone made temporarily much more electrically conductive by stress-induced reorganization of pore geometry. Using a numerical model we show that while this mechanism could produce a significant increase in ULF variations, mutual induction between the fault zone and the surrounding crust would probably limit the amplitude increase to levels well below those observed at Loma Prieta. We consider a variant on this quasi-static conductive fault zone model in which low frequency telluric currents are modulated by small higher frequency variations of bulk fault zone conductivity. We show that because the spectrum of natural EM variations is red, substantially larger relative increases in ULF magnetic fields could be produced by this mechanism with even small conductivity fluctuations at these frequencies. These variations would be easy to detect with a well-designed experiment, if they occurred. In principle this mechanism could explain the Loma Prieta ULF observations, however the magnitude of conductivity fluctuations that would be required to match the very large reported amplification factors still appears to be too large to be physically plausible.     

直接消除电磁耦合的斩波去耦方法

应用激发极化法时经常遇到电磁耦合的干扰.基于耦合效应和激电效应在双频波测量波形上具有不同的表现形式,提出一种在野外采集激电数据时,直接剔除电磁耦合而又基本上不损失激电信息的“斩波”去耦方法.本文从“斩波去耦”的方法原理、斩波对正常双频波产生的畸变、以及存在EM和IP效应时的双频波斩波等方面进行了详细的推导和分析.数值计算结果表明,该方法去耦效果明显,可以在现场条件下直接去耦,值得在激发极化法中推广应用     

AN INTERPRETATION THEORY FOR INDUCED POLARIZATION VERTICAL SOUNDINGS (TIME-DOMAIN)*

In this paper it is shown how one may obtain a generalized Ohm's law which relates the induced polarization electric field to the steady-state current density through the introduction of a fictitious resistivity defined as the product of the chargeability and the resistivity of a given medium. The potential generated by the induced polarization is calculated at any point in a layered earth by the same procedure as used for calculating the potential due to a point source of direct current. On the basis of the definition of the apparent chargeability m_a, the expressions of m_a for different stratigraphie situations are obtained, provided the IP measurements are carried out on surface with an appropriate AMNB array. These expressions may be used to plot master curves for IP vertical soundings. Finally some field experiments over sedimentary formations and the quantitative interpretation procedure are reported.     

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.     

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.     

激发极化观测中电磁耦合的时间特性

为了克服激发极化观测中经常遇到的电磁耦合干扰,以往研究大多基于电磁耦合的频率特点.本文从电偶极子电磁耦合的基本公式出发,借助拉普拉斯变换和数值积分,着重研究了电磁耦合的时间历程及特点.研究表明,在足够长时间的矩形电流激励下,偶极装置电磁耦合的响应电位差表现为矩形激励电流开通时的上升斜坡和关断时的下降斜坡,并且使电磁信号到达的时间推迟;中间梯度装置的电磁耦合响应表现为矩形激励电流开通时的正尖脉冲和关断时的负尖脉冲,同时也使电磁信号到达的时间推迟.与激电法常用矩形电流的接通（关断）时间相比,不论采用何种观测装置,电磁感应耦合终归是一种“短命”的瞬变现象.不同装置的电磁耦合既有不同之处,又有相同之点.在固体矿产勘查常见的地电条件下,就激电法常用的中间梯度和偶极装置以及供电波形而言,电磁耦合的时间一般不超过10^-2 s,利用这些时间特点能够达到区分并克服电磁耦合的目的.     

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.     

电场差分法理论研究

对均匀大地和二层大地的频率域电场差分法理论曲线进行的计算和研究表明,电场差分法是一种以目标地质体和围岩在电阻率和极化率方面的差异为基础,根据测量到的能明显反映激电响应与电磁耦合响应的比值随频率变化的P参数曲线在性态上的差异来圈划油气藏的一种电法勘探方法. 因为它圈划油气藏的物理基础是油气藏上方与围岩上方测量到的电磁耦合响应和激电响应的差异,所以在圈划油气藏方面,该方法与复电阻率法类似,同样具有找油气的优良性能. 此外,该方法与复电阻率法相比,具有对目标地质体和围岩的电阻率和极化率差异反映灵敏、效率高、成本低等优点,所以它是一种值得推广的找油气的电法勘探方法.     

时间域航空电磁法激电效应对电磁扩散的影响

由于激发极化效应的影响,时间域航空电磁晚期道信号经常会出现变号现象.基于电阻率的传统反演方法无法对变号数据进行正确反演,因此通常在数据处理中予以剔除.为深入了解极化介质的电磁扩散特征,认识航空瞬变电磁负响应的产生机理,本文研究时间域航空电磁系统的电磁扩散特征.我们以均匀极化、非极化半空间及层状介质模型为例,通过直接积分的方法求解频率域电场响应,并由欧姆定律得到电流响应,再经过汉克尔变换得到时间域电流响应.通过研究电流随时间在地下极化介质中的传播特征研究电磁扩散过程;通过对比不同激电参数对电磁扩散的影响,研究极化介质中感应电流与极化电流的扩散规律,从而合理地解释极化介质中负响应的产生机理.基于本文研究和分析结果,可加深对时间域航空电磁法中激电效应的认识.     

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.     

Calibration of the Regional Crustal Waveguide and the Retrieval of Source Parameters Using Waveform Modeling

v--vRegional crustal waveguide calibration is essential to the retrieval of source parameters and the location of smaller (M < 4.8) seismic events. This path calibration of regional seismic phases is strongly dependent on the accuracy of hypocentral locations of calibration (or master) events. This information can be difficult to obtain, especially for smaller events. Generally, explosion or quarry blast generated travel-time data with known locations and origin times are useful for developing the path calibration parameters, but in many regions such data sets are scanty or do not exist. We present a method which is useful for regional path calibration independent of such data, i.e. with earthquakes, which is applicable for events down to M_w = 4 and which has successfully been applied in India, central Asia, western Mediterranean, North Africa, Tibet and the former Soviet Union. These studies suggest that reliably determining depth is essential to establishing accurate epicentral location and origin time for events. We find that the error in source depth does not necessarily trade-off only with the origin time for events with poor azimuthal coverage, but with the horizontal location as well, thus resulting in poor epicentral locations. For example, hypocenters for some events in central Asia were found to move from their fixed-depth locations by about 20 km. Such errors in location and depth will propagate into path calibration parameters, particularly with respect to travel times. The modeling of teleseismic depth phases (pP, sP) yields accurate depths for earthquakes down to magnitude M_w = 4.7. This M_w threshold can be lowered to four if regional seismograms are used in conjunction with a calibrated velocity structure model to determine depth, with the relative amplitude of the P_nl waves to the surface waves and the interaction of regional sPmP and pPmP phases being good indicators of event depths. We also found that for deep events a seismic phase which follows an S-wave path to the surface and becomes critical, developing a head wave by S to P conversion is also indicative of depth. The detailed characteristic of this phase is controlled by the crustal waveguide. The key to calibrating regionalized crustal velocity structure is to determine depths for a set of master events by applying the above methods and then by modeling characteristic features that are recorded on the regional waveforms. The regionalization scheme can also incorporate mixed-path crustal waveguide models for cases in which seismic waves traverse two or more distinctly different crustal structures. We also demonstrate that once depths are established, we need only two-stations travel-time data to obtain reliable epicentral locations using a new adaptive grid-search technique which yields locations similar to those determined using travel-time data from local seismic networks with better azimuthal coverage.     

两层大地中三维体的激发极化与电阻率响应的积分方程模拟

用散射、叠加方法推导出两层大地的并矢格林函数。使用这些函数,含三维异常体的二层大地的边值问题转变为积分方程。使用矩量法,可解此积分方程。 使用偶极-偶极装置计算了激发极化和电阻率响应的几个数值结果,并对这一算法作了几方面检验。     

MODEL STUDIES ON SOME ASPECTS OF RESISTIVITY AND MEMBRANE POLARIZATION BEHAVIOUR OVER A LAYERED EARTH*

Electrolytic model tank experiments to study resistivity and time domain induced polarization (IP) response over layered earth models were initiated primarily to facilitate the understanding of field results. Alternate layers of clay and sand (or clay-coated sand) with, in some cases, a surficial layer of water were assembled in the tank and resistivity and IP measurements made for a range of electrode spacings using the Wenner configuration. Graphite and silver-silver chloride electrodes were used as current and potential electrodes respectively. Clay-coated (3% by weight) sand was found to generate stronger polarization than either clay or sand alone. Apparent chargeability m_a was observed to be positive for a nonpolarizable surface layer. For a polarizable surface layer, the sign of IP was controlled by the polarizability, the thickness of the second layer, and the spacing of the electrode spreads. The apparent chargeability m_a can theoretically change sign from positive to negative and vice versa with a gradual increase in electrode spacing, and such negative IP effects were obtained in a few observations. A simultaneous decrease in IP and an increase in resistivity, which is a qualitative diagnostic feature for the occurrence of clean freshwater sand aquifers, could also be generated in the model tank experiment. Combined resistivity and IP soundings were carried out near Fredericton Junction and Tracy, New Brunswick, Canada. Field curves are presented along with the model curves for qualitative comparison and understanding of IP behaviour over a layered earth. Twenty-five out of twenty-seven soundings show only positive apparent chargeabilities, whereas two show chargeability sign changes (positive/negative/positive). The model study gives reason to believe that surface soils and Quaternary gravel boulder deposits near Fredericton Junction are relatively non-polarizable. As an auxiliary experiment, sand and clay were taken in different proportions by weight and mixed thoroughly with water in a cement mixer. The mixtures were then compressed with a suitable die and plunger under 3.6 Pa pressure to prepare cylindrical samples of height 18 cm and diameter 15.5 cm. IP measurements were done on the flat faces using the Wenner configuration with a= 2 cm. Chargeability was found to be negative for 100 and 90% clay mixtures. It reached a positive maximum for an 80% clay-20% sand mixture and then decreased gradually with increasing sand and decreasing clay content.     

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.     

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.     

1D and 2D Cole-Cole-inversion of time-domain induced-polarization data

A new method for the 2D inversion of induced polarization (IP) data in the time domain has been developed. The entire IP transients were observed and inverted into 2D Cole-Cole earth models, including resistivity, chargeability, relaxation time and the frequency constant. Firstly, a modified 1D time-domain electromagnetic algorithm was used to calculate the response of a layered polarizable ground. The transient signals were then inverted using the Marquardt method to derive the Cole-Cole parameters of each layer. However, model calculations showed that the EM effects could be neglected for the time range (>1 ms) and for the transmitter–receiver distances (<50 m) used in this study. Therefore, the induction effects were not considered for the solution of the 2D inverse problem and a DC solution was applied. An approximative forward algorithm was introduced in order to calculate the IP transients directly in the time domain and in order to speed up the inverse procedure. The approximation is highly accurate, and this is demonstrated by comparing the approximations with their exact solutions up to 3D. The inverse algorithm presented consists of two steps. The transient voltages of an array data set were inverted separately into a two-dimensional resistivity model for each time channel. The time-dependent resistivity of each cell was then interpreted as the response of a homogeneous half-space. In the 2D inversion algorithm, a 3D DC algorithm was used as a forward operator. The method only requires a standard 2D DC inversion and a homogenous half-space Cole-Cole inversion. The developed algorithm has been successfully applied to synthetic data sets and to a field data set obtained from a waste site situated close to Düren in Germany.     

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.     

INDUCED POLARIZATION RESPONSE OF A HORIZONTALLY MULTILAYERED EARTH WITH NO RESISTIVITY CONTRAST*

The induced polarization response of a horizontally multilayered earth with no resistivity contrast can rapidly be calculated on a desk calculator or minicomputer for any electrode array. The formulation is a simple series summation of the products of weighting coefficients and the true induced polarization responses for each of the layers. The coefficients are directly derivable from the corresponding resistivity model. This series approach to IP formulation was originally described by Seigel but has not been treated extensively in the present-day geophysical literature. This method can be applied to either time or frequency domain induced polarization measurements. Once the coefficients are known, apparent induced polarization response can readily be obtained by judicious substitution of known, suspected, or assumed values of the true induced polarization of each layer. Basic formulation is presented for the IP potential coefficients (pole-pole or two array) with no resistivity contrast between the layers. From these coefficients, response of any number of layers for any electrode array can be obtained by suitable differentiation. Some examples of Wenner array for a three-layered earth and dipole-dipole array for a four-layered earth are used to illustrate the application. The results of this technique are valid for many natural situations of modest resistivity contrast. However, they definitely cannot be used if there are highly contrasting resistivity layers present. Such an approach is conceptually simple and is useful for survey planning, checking or setting the “depth-of-penetration”of a given array. For field induced polarization data that fits reasonably well to the no-resistivity-contrast model, this simple approach facilitates quantitative interpretation.     

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.     

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.