Accurate Location of Synthetic Acoustic Emissions and Location Sensitivity to Relocation Methods, Velocity Perturbations, and Seismic Anisotropy

Acoustic emission (AE) monitoring is a non-invasive method of monitoring fracturing both in situ, and in experimental rock deformation studies. Until recently, the major impediment for imaging brittle failure within a rock mass is the accuracy at which the hypocenters may be located. However, recent advances in the location of regional scale earthquakes have successfully reduced hypocentral uncertainties by an order of magnitude. The least-squares Geiger, master event relocation, and double difference methods have been considered in a series of synthetic experiments which investigate their ability to resolve AE hypocentral locations. The effect of AE hypocenter location accuracy due to seismic velocity perturbations, uncertainty in the first arrival pick, array geometry and the inversion of a seismically anisotropic structure with an isotropic velocity model were tested. Hypocenters determined using the Geiger procedure for a homogeneous, isotropic sample with a known velocity model gave a RMS error for the hypocenter locations of 2.6 mm; in contrast the double difference method is capable of reducing the location error of these hypocenters by an order of magnitude. We test uncertainties in velocity model of up to ±10% and show that the double difference method can attain the same RMS error as using the standard Geiger procedure with a known velocity model. The double difference method is also capable of precise locations even in a 40% anisotropic velocity structure using an isotropic model for location and attains a RMS mislocation error of 2.6 mm that is comparable to a RMS mislocation error produced with an isotropic known velocity model using the Geiger approach. We test the effect of sensor geometry on location accuracy and find that, even when sensors are missing, the double difference method is capable of a 1.43 mm total RMS mislocation compared to 4.58 mm for the Geiger method. The accuracy of automatic picking algorithms used for AE studies is ±0.5 μs (1 time sample when the sampling rate is 0.2 μs). We investigate how AE locations are effected by the accuracy of first arrival picking by randomly delaying the actual first arrival by up to 5 time samples. We find that even when noise levels are set to 5 time samples the double difference method successfully relocates the synthetic AE.     

Compaction of a Rock Fracture Moderated by Competing Roles of Stress Corrosion and Pressure Solution

Unusually rapid closure of stressed fractures, observed in the initial stages of loading and at low temperatures, is examined using models for subcritical crack growth and pressure solution. The model for stress corrosion examines tensile stress concentrations induced at the Hertzian contact of propping fracture asperities, and mediates fracture growth according to a kinetic rate law. Conversely, pressure solution is described by the rate-limiting process of dissolution, resulting from the elevated stresses realized at the propping asperity contact. Both models are capable of following the observed compaction of fractures in novaculite. However, closure rates predicted for stress corrosion cracking are orders of magnitudes faster than those predicted for pressure dissolution. For consistent kinetic parameters, predictions from stress corrosion better replicate experimental observations, especially in the short-term and at low temperature when mechanical effects are anticipated to dominate. Rates and magnitudes of both stress corrosion and pressure solution are dependent on stresses exerted over propping asperities. Rates of closure due to stress corrosion cracking are shown to be always higher than for pressure solution, except where stress corrosion ceases as contact areas grow, and local stresses drop below an activation threshold. A simple rate law is apparent for the progress of fracture closure, defined in terms of a constant and an exponent applied to the test duration. For current experimental observations, this rate law is shown to replicate early progress data, and shows promise to define the evolution of transport properties of fractures over extended durations.     

Patos Lagoon outflow within the Río de la Plata plume using an airborne salinity mapper: Observing an embedded plume

Major river systems discharging into continental shelf waters frequently form buoyant coastal currents that propagate along the continental shelf in the direction of coastal trapped wave propagation (with the coast on the right/left, in the northern/southern hemisphere). The combined flow of the Uruguay and Paraná Rivers, which discharges freshwater into the Río de la Plata estuary (Lat. ∼36°S), often gives rise to a buoyant coastal current (the ‘Plata plume’) that extends northward along the continental shelf off Uruguay and Southern Brazil. Depending upon the prevailing rainfall, wind and tidal conditions, the Patos/Mirim Lagoon complex (Lat. ∼32°S) may also produce a freshwater outflow plume that expands across the inner continental shelf. Under these circumstances the Patos outflow plume can be embedded in temperature, salinity and current fields that are strongly influenced by the larger Plata plume. The purpose of this paper is to present observations of such an embedded plume structure and to determine the dynamical characteristics of the ambient and embedded plumes.     

Cross-shelf transport into nearshore waters due to shoaling internal tides in San Pedro Bay, CA

In the summer of 2001, a coastal ocean measurement program in the southeastern portion of San Pedro Bay, CA, was designed and carried out. One aim of the program was to determine the strength and effectiveness of local cross-shelf transport processes. A particular objective was to assess the ability of semidiurnal internal tidal currents to move suspended material a net distance across the shelf. Hence, a dense array of moorings was deployed across the shelf to monitor the transport patterns associated with fluctuations in currents, temperature and salinity. An associated hydrographic program periodically monitored synoptic changes in the spatial patterns of temperature, salinity, nutrients and bacteria. This set of measurements show that a series of energetic internal tides can, but do not always, transport subthermocline water, dissolved and suspended material from the middle of the shelf into the surfzone. Effective cross-shelf transport occurs only when (1) internal tides at the shelf break are strong and (2) subtidal currents flow strongly downcoast. The subtidal downcoast flow causes isotherms to tilt upward toward the coast, which allows energetic, nonlinear internal tidal currents to carry subthermocline waters into the surfzone. During these events, which may last for several days, the transported water remains in the surfzone until the internal tidal current pulses and/or the downcoast subtidal currents disappear. This nonlinear internal tide cross-shelf transport process was capable of carrying water and the associated suspended or dissolved material from the mid-shelf into the surfzone, but there were no observation of transport from the shelf break into the surfzone. Dissolved nutrients and suspended particulates (such as phytoplankton) transported from the mid-shelf into the nearshore region by nonlinear internal tides may contribute to nearshore algal blooms, including harmful algal blooms that occur off local beaches.     

Storm surge computations in estuarine and near-coastal regions: the Mersey estuary and Irish Sea area

An unstructured grid storm surge model of the west coast of Britain, incorporating a high-resolution representation of the Mersey estuary is used to examine storm surge dynamics in the region. The focus of the study is the major surge that occurred during the period 11–14 November 1977, which has been investigated previously using uniform grid finite difference models and a finite element model of the west coast of Britain. However, none of these models included the Mersey estuary. Comparison of solutions in the eastern Irish Sea with those computed with these earlier models showed that, away from the Liverpool Bay region, the inclusion of the Mersey estuary had little effect. However, at the entrance to the Mersey, its inclusion did influence the solution. By including a detailed representation of the Mersey estuary within the model, it was possible to conduct a detailed study of storm surge propagation in the Mersey, which had never previously been performed. This detailed study showed for the first time that the surge’s temporal variability within the estuary is influenced by surge elevation at its entrance. This varies with time as a function of spatial and temporal variations of wind stress over the west coast of Britain. Within the Mersey, calculations show that the spatial variability is mainly determined by changes in bottom topography, which had not been included in earlier finite difference models. However, since water depth is influenced by variations in tidal elevation, this, together with tide surge interaction through bottom friction and momentum advection, influences the surge. The ability of the finite element model to vary the mesh in near-shore regions to such an extent that it can resolve the Mersey and hence the impact of the Mersey estuary upon the Liverpool Bay circulation shows that it has distinct advantages over earlier finite difference models. In the absence of detailed measurements within the Mersey at the time of the surge, it was not possible to validate predicted surge elevations within the Mersey. However, significant insight into physical processes influencing the surge propagation down the estuary, its reflection and spatial/temporal variability could be gained.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Diversity and spatial distribution of Opuntia spp. in the Mediterranean Basin

This work offers a methodology based on statistical tools and GIS applications in a study of the adaptability of Opuntia to the environment, applied on the Mediterranean island of Sardinia. Detailed surveys of habitats containing naturalized populations of Opuntia spp. were made in order to select study sites. The topographic, soil and climatic features of these are described. The local biotypes are ascribed mainly to the spineless Opuntia ficus–indica (L.) Mill. and to the spiny Opuntia amyclaea Ten. Although introduced at different times, for different uses, they have adapted to the Sardinian environment and established in the wild. The distributions of the two species differed. The FloraMap^® modelling approach was used to simulate Opuntia distribution as a function of the island's climate and produced individual distribution maps using the accession points for the two species. Our analysis showed that environmental conditions at elevations above 800 m are the main limiting factors for the presence and abundance of Opuntia in Sardinia. Spatial representation developed by GIS software indicated temperatures below 2 °C as the threshold impairing the species naturalized distribution. The statistic results also show that variation in rainfall does not affect Opuntia spatial distribution in Sardinia.     

Measuring proton energies and fluxes using EIT (SOHO) CCD areas outside the solar disk images

An indirect proton flux measuring tool based on discrimination of the energy deposited by protons in 128 × 128 pixel EIT CCD areas outside the solar disk images is presented. Single pixel intensity events are converted into proton incident energy flux using modeled energy deposition curves for angles of incidence ±60° in four EIT spatial areas with different proton stopping power. The extracted proton flux is corrected for both the loss of one‐pixel events in the range of angles of incidence as well as for the contribution to the single pixel events resulting from scattered middle‐energy protons (low‐energy or high‐energy particles are stopped by the EIT components or pass through them, accordingly). A simple geometrical approach was found and applied to correct for a non‐unique relation between the proton‐associated CCD output signal and the incident proton energy. With this geometrical approximation four unique proton incident energy ranges were determined as 45–49, 145–154, 297–335, and 390–440 MeV. The indirect proton flux measuring tool has been tested by comparing Solar Energetic Particles (SEP) flux temporal profiles extracted from the EIT CCD frames and downloaded from the GOES database for the Bastille Day (BD) of 2000 July 14 and the more recent 2005 January 20 events. The SEP flux temporal profiles and proton spectra extracted from the EIT in the relatively narrow energy ranges between 45 and 440 MeV reported here are consistent with the related GOES profiles. The four additional EIT extracted ranges provide higher energy resolution of the SEP data. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)