A selective phenomenology of the seismicity of Southern California.

Predictions of earthquakes that are based on observations of precursory seismicity cannot depend on the average properties of the seismicity, such as the Gutenberg-Richter (G-R) distribution. Instead it must depend on the fluctuations in seismicity. We summarize the observational data of the fluctuations of seismicity in space, in time, and in a coupled space-time regime over the past 60 yr in Southern California, to provide a basis for determining whether these fluctuations are correlated with the times and locations of future strong earthquakes in an appropriate time- and space-scale. The simple extrapolation of the G-R distribution must lead to an overestimate of the risk due to large earthquakes. There may be two classes of earthquakes: the small earthquakes that satisfy the G-R law and the larger and large ones. Most observations of fluctuations of seismicity are of the rate of occurrence of smaller earthquakes. Large earthquakes are observed to be preceded by significant quiescence on the faults on which they occur and by an intensification of activity at distance. It is likely that the fluctuations are due to the nature of fractures on individual faults of the network of faults. There are significant inhomogeneities on these faults, which we assume will have an important influence on the nature of self-organization of seismicity. The principal source of the inhomogeneity on the large scale is the influence of geometry--i.e., of the nonplanarity of faults and the system of faults.

Geochemistry of basalts from ODP Leg 49, North Atlantic Ocean

The geochemistry of basalts recovered from seven sites in the North Atlantic is described with particular reference to minor elements. Three sites (407, 408, and 409) along the same mantle flow line, transverse to the Reykjanes Ridge at about 63°N, provide information on the composition of basalts erupted over a 34-m.y. interval between 2.3 and 36 m.y. ago. At Site 410, at 45°N, penetration into 10 m.y.-old crust west of the ridge axis permits comparisons with young basalts dredged from the median valley at 45°N. Three sites in the FAMOUS area at about 36°N provided material from very young (1 m.y.) basaltic crust (Site 411), and material to test the geochemical coherence of basalts of different ages (1.5 and 3.5 m.y.) on either side of a fracture zone (Sites 412 and 413). These sites complement earlier data from dredged and drilled sites (Leg 37) in the FAMOUS area. At Site 407, four geochemically distinct basalt units occur, with different normative and rare-earth element (REE) characteristics, and there is a clear correlation with magnetic stratigraphy. Yet there is a remarkable consistency in incompatible element ratios between these units, indicating derivation from an essentially similar mantle source. The basalts from the younger sites, 408 and 409, show a similar range of normative and REE variation, but incompatible element ratios are identical to those at Site 407, indicating that basalts at all three sites were produced from a mantle source which was geochemically relatively uniform. Rare-earth differences between the basalts can be interpreted in terms of variations in the degree and depth of partial melting causing HREE (+Y) retention in the source, although there may be some inter-site differences with respect to REE. A similar picture is presented at 45°N. Apparently a range of tholeiitic, transitional, and alkalic basalts were being erupted 10 m.y. ago, which have almost identical geochemical characteristics to those recently erupted in the median valley at 45°N. Incompatible element ratios are markedly different from those recorded at the Reykjanes Ridge. Basalts recovered from the FAMOUS sites are geochemically similar to previous samples recovered from the FAMOUS area, and their incompatible element ratios are similar, but not identical, to those at 45°N. However, total trace element levels are consistently lower than in 45°N basalts, which might imply smaller degrees of partial melting and/or greater depths of magma generation at 45°N, or higher trace element levels in the mantle source at 45°N. Few of the basalts recovered on Leg 49 have the geochemical characteristics of typical "MORB" (e.g., Nazca Plate, Leg 34). The data strongly support models invoking geochemical inhomogeneity in the source regions of basalts produced at the Mid-Atlantic Ridge. However, the data also introduce an additional time factor into such models and demonstrate the uniformity of the mantle source at a particular ridge sector (over periods in excess of 30 m.y.), while emphasizing the marked differences along the ridge. Mixing models invoking "depleted" and "enriched" mantle sources would seem to be inadequate to account for the observed variations.

(Table 1) Lithium contents and isotopic composition in sediments from DSDP Holes 64-477 and 64-477A

Lithium isotopic compositions of hydrothermally altered sediments of Deep Sea Drilling Project (DSDP) site 477/477A, as well as high temperature vent fluids of the Guaymas Basin, have been determined to gain an understanding of lithium exchange during fluid-sediment interaction at this sediment-covered spreading center. Unaltered turbidite of the basin has a d6Li value of -10%, 5-7% heavier than fresh oceanic basalts. Contact metamorphism induced by a shallow sill intrusion results in a decrease of the lithium content of the adjacent sediments and a lighter isotopic value (-8%). Below the sill, sediments altered by a deep-seated hydrothermal system show strong depletions in lithium, while lithium isotopic compositions vary greatly, ranging from -11 to +1%. The shift to lighter composition is the result of preferential retention of the lighter isotope in recrystallized phases after destruction of the primary minerals. The complexity of the isotope profile is attributed to inhomogeneity in mineral composition, the tortuous pathway of fluids and the temperature effect on isotopic fractionation. The range of lithium concentration and d6Li values for the vent fluids sampled in 1982 and 1985 overlaps with that of the sediment-free mid-ocean ridge systems. The lack of a distinct expression of sediment input is explained in terms of a flow-through system with continuous water recharge. The observations on the natural system agree well with the results of laboratory hydrothermal experiments. The experimental study demonstrates the importance of temperature, pressure, water/rock ratio, substrate composition and reaction time on the lithium isotopic composition of the reacted fluid. High temperature authigenic phases do not seem to constitute an important sink for lithium and sediments of a hydrothermal system such as Guaymas are a source of lithium to the ocean. The ready mobility of lithium in the sediment under elevated temperature and pressure conditions also has important implications for lithium cycling in subduction zones.

Phytoplankton and bacterial productivity in the Tugur Bay (Sea of Okhotsk) in July-August 1990

A study of quantitative characteristics: phytoplankton photosynthesis (Ph), bacterial assimilation of CO2 (BA), total abundance of bacteria (TAB) and organic matter destruction (D) was carried out in waters the Tugur Bay (tidal level fluctuations up to 7 m) in July-August 1990. Calculations were made of integral indices in some parts: Ph -10-630, BA - 8-29, D - 280-1015 and of total primary production (TPP) - 18-652 mg C/(m2 day). According to obtained data and TAB the ecosystem of the Tugur Bay can be regarded as oligotrophic-mesotrophic one. Dependence on spatial and temporal inhomogeneity of primary productional processes on tide-ebb cycles was found. The role of bacterial relations in the ecosystem of the bay was shown. Portion of ?A in TPP varied from 4 to 44% reaching the maximum in desalinated water during the tide-ebb.

A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging

A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings. especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging. (C) 2004 Elsevier Inc. All rights reserved.

A Lanczos-powered implementation of the Faber polynomial quantum time propagator for reaction probabilities

Complementing our recent work on subspace wavepacket propagation [Chem. Phys. Lett. 336 (2001) 149], we introduce a Lanczos-based implementation of the Faber polynomial quantum long-time propagator. The original version [J. Chem. Phys. 101 (1994) 10493] implicitly handles non-Hermitian Hamiltonians, that is, those perturbed by imaginary absorbing potentials to handle unwanted reflection effects. However, like many wavepacket propagation schemes, it encounters a bottleneck associated with dense matrix-vector multiplications. Our implementation seeks to reduce the quantity of such costly operations without sacrificing numerical accuracy. For some benchmark scattering problems, our approach compares favourably with the original. (C) 2004 Elsevier B.V. All rights reserved.

Iterative quantum computations of HO2 bound states and resonances for J=4 and 5

Bound and resonance states of HO2 have been calculated by both the complex Lanczos homogeneous filter diagonalisation (LHFD) method(1,2) and the real Chebyshev filter diagonalisation method(3,4) for non-zero total angular momentum J = 4 and 5. For bound states, the agreement between the two methods is quite satisfactory; for resonances while the energies are in good agreement, the widths are only in general agreement. The relative performances of the two iterative FD methods have also been discussed in terms of efficiency as well as convergence behaviour for such a computationally challenging problem. A helicity quantum number Ohm assignment (within the helicity conserving approximation) is performed and the results indicate that Coriolis coupling becomes more important as J increases and the helicity conserving approximation is not a good one for the HO2 resonance states.

Recent developments in the quantum dynamical characterization of unimolecular resonances

We give a selective review of quantum mechanical methods for calculating and characterizing resonances in small molecular systems, with an emphasis on recent progress in Chebyshev and Lanczos iterative methods. Two archetypal molecular systems are discussed: isolated resonances in HCO, which exhibit regular mode and state specificity, and overlapping resonances in strongly bound HO2, which exhibit irregular and chaotic behavior. Recent progresses for non-zero total angular momentum J calculations of resonances including parallel computing models are also included and future directions in this field are discussed.

Method for a detailed measurement of image intensity nonuniformity in magnetic resonance imaging

In magnetic resonance imaging (MRI), the MR signal intensity can vary spatially and this spatial variation is usually referred to as MR intensity nonuniformity. Although the main source of intensity nonuniformity arises from B, inhomogeneity of the coil acting as a receiver and/or transmitter, geometric distortion also alters the MR signal intensity. It is useful on some occasions to have these two different sources be separately measured and analyzed. In this paper, we present a practical method for a detailed measurement of the MR intensity nonuniformity. This method is based on the same three-dimensional geometric phantom that was recently developed for a complete measurement of the geometric distortion in MR systems. In this paper, the contribution to the intensity nonuniformity from the geometric distortion can be estimated and thus, it provides a mechanism for estimation of the intensity nonuniformity that reflects solely the spatial characteristics arising from B-1. Additionally, a comprehensive scheme for characterization of the intensity nonuniformity based on the new measurement method is proposed. To demonstrate the method, the intensity nonuniformity in a 1.5 T Sonata MR system was measured and is used to illustrate the main features of the method. (c) 2005 American Association of Physicists in Medicine.

Numerical modeling of 11.1 T MRI of a human head using MoM/FDTD methods

In a recent study, severe distortions in the proton images of an excised, fixed, human brain in an 11.1 Tesla/40 cm MR instrument have been observed, and the effect modeled on phantom images using a finite difference time domain (FDTD) model. in the present study, we extend these simulations to that of a complete human head, employing a hybrid FDTD and method of moments (MoM) approach, which provides a validated method for simulating biological samples in coil structures. The effect of fixative on the image distortions is explored. importantly, temperature distributions within the head are also simulated using a bioheat method based on parameters derived from the electromagnetic simulations. The MoM/FDTD simulations confirm that the transverse magnetic field (B,) from a ReCav resonator exhibits good homogeneity in air but strong inhomogeneity when loaded with the head with or without fixative. The fixative serves to increase the distortions, but they are still significant for the in vivo simulations. The simulated signal intensity (SI) distribution within the sample confirm the distortions in the experimental images are caused by the complex interactions of the incident electromagnetic fields with tissue, which is heterogeneous in terms of conductivity and permittivity. The temperature distribution is likewise heterogeneous, raising concerns regarding hot spot generation in the sample that may exceed acceptable levels in future in vivo studies. As human imaging at 11.1 T is some time away, simulations are important in terms of predicting potential safety issues as well as evaluating practical concerns about the quality of images. Simulation on a whole human head at 11.1 T implies the wave behavior presents significant engineering challenges for ultra-high-field (UHF) MRI. Novel strategies will have to be employed in imaging technique and resonator design for UHF MRI to achieve the theoretical signal-to-noise ratio (SNR) improvements it offers over lower field systems. (C) 2005 Wiley Periodicals, Inc.

Effect of smoke inhalation on viscoelastic properties and ventilation distribution in sheep

Smoke inhalation injuries are the leading cause of mortality from burn injury. Airway obstruction due to mucus plugging and bronchoconstriction can cause severe ventilation inhomogeneity and worsen hypoxia. Studies describing changes of viscoelastic characteristics of the lung after smoke inhalation are missing. We present results of a new smoke inhalation device in sheep and describe pathophysiological changes after smoke exposure. Fifteen female Merino ewes were anesthetized and intubated. Baseline data using electrical impedance tomography and multiple-breath inert-gas washout were obtained by measuring ventilation distribution, functional residual capacity, lung clearance index, dynamic compliance, and stress index. Ten sheep were exposed to standardized cotton smoke insufflations and five sheep to sham smoke insufflations. Measured carboxyhemoglobin before inhalation was 3.87 +/- 0.28% and 5 min after smoke was 61.5 +/- 2.1%, range 50-69.4% ( P < 0.001). Two hours after smoke functional residual capacity decreased from 1,773 +/- 226 to 1,006 +/- 129 ml and lung clearance index increased from 10.4 +/- 0.4 to 14.2 +/- 0.9. Dynamic compliance decreased from 56.6 +/- 5.5 to 32.8 +/- 3.2 ml/ cmH(2)O. Stress index increased from 0.994 +/- 0.009 to 1.081 +/- 0.011 ( P < 0.01) ( all means +/- SE, P < 0.05). Electrical impedance tomography showed a shift of ventilation from the dependent to the independent lung after smoke exposure. No significant change was seen in the sham group. Smoke inhalation caused immediate onset in pulmonary dysfunction and significant ventilation inhomogeneity. The smoke inhalation device as presented may be useful for interventional studies.

An inverse method for designing loaded RF coils in MRI

Radio-frequency ( RF) coils are designed such that they induce homogeneous magnetic fields within some region of interest within a magnetic resonance imaging ( MRI) scanner. Loading the scanner with a patient disrupts the homogeneity of these fields and can lead to a considerable degradation of the quality of the acquired image. In this paper, an inverse method is presented for designing RF coils, in which the presence of a load ( patient) within the MRI scanner is accounted for in the model. To approximate the finite length of the coil, a Fourier series expansion is considered for the coil current density and for the induced fields. Regularization is used to solve this ill-conditioned inverse problem for the unknown Fourier coefficients. That is, the error between the induced and homogeneous target fields is minimized along with an additional constraint, chosen in this paper to represent the curvature of the coil windings. Smooth winding patterns are obtained for both unloaded and loaded coils. RF fields with a high level of homogeneity are obtained in the unloaded case and a limit to the level of homogeneity attainable is observed in the loaded case.

First-principles quantum simulations of dissociation of molecular condensates: Atom correlations in momentum space

We investigate the quantum many-body dynamics of dissociation of a Bose-Einstein condensate of molecular dimers into pairs of constituent bosonic atoms and analyze the resulting atom-atom correlations. The quantum fields of both the molecules and atoms are simulated from first principles in three dimensions using the positive-P representation method. This allows us to provide an exact treatment of the molecular field depletion and s-wave scattering interactions between the particles, as well as to extend the analysis to nonuniform systems. In the simplest uniform case, we find that the major source of atom-atom decorrelation is atom-atom recombination which produces molecules outside the initially occupied condensate mode. The unwanted molecules are formed from dissociated atom pairs with nonopposite momenta. The net effect of this process-which becomes increasingly significant for dissociation durations corresponding to more than about 40% conversion-is to reduce the atom-atom correlations. In addition, for nonuniform systems we find that mode mixing due to inhomogeneity can result in further degradation of the correlation signal. We characterize the correlation strength via the degree of squeezing of particle number-difference fluctuations in a certain momentum-space volume and show that the correlation strength can be increased if the signals are binned into larger counting volumes.

White and gray matter development in human fetal, newborn and pediatric brains

Brain anatomy is characterized by dramatic growth from the end of the second trimester through the neonatal stage. The characterization of normal axonal growth of the white matter tracts has not been well-documented to date and could provide important clues to understanding the extensive inhomogeneity of white matter injuries in cerebral palsy (CP) patients. However, anatomical studies of human brain development during this period are surprisingly scarce and histology-based atlases have become available only recently. Diffusion tensor magnetic resonance imaging (DTMRI) can reveal detailed anatomy of white matter. We acquired diffusion tensor images (DTI) of postmortem fetal brain samples and in vivo neonates and children. Neural structures were annotated in two-dimensional (2D) slices, segmented, measured, and reconstructed three-dimensionally (3D). The growth status of various white matter tracts was evaluated on cross-sections at 19-20 gestational weeks, and compared with 0-month-old neonates and 5- to 6-year-old children. Limbic, commissural, association, and projection white matter tracts and gray matter structures were illustrated in 3D and quantitatively characterized to assess their dynamic changes. The overall pattern of the time courses for the development of different white matter is that limbic fibers develop first and association fibers last and commissural and projection fibers are forming from anterior to posterior part of the brain. The resultant DTNIRI-based 3D human brain data will be a valuable resource for human brain developmental study and will provide reference standards for diagnostic radiology of premature newborns. (c) 2006 Elsevier Inc. All rights reserved.

High-field magnetic resonance imaging with reduced field/tissue RF artefacts - A modeling study using hybrid MoM/FEM and FDTD technique

Due to complex field/tissue interactions, high-field magnetic resonance (MR) images suffer significant image distortions that result in compromised diagnostic quality. A new method that attempts to remove these distortions is proposed in this paper and is based on the use of transceiver-phased arrays. The proposed system uses, in the examples presented herein, a shielded four-element transceive-phased array head coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both the images together, the image distortion can be reduced several fold. A combined hybrid method of moments (MoM)/finite element method (FEM) and finite-difference time-domain (FDTD) technique is proposed and used to elucidate the concept of the new method and to accurately evaluate the electromagnetic field (EMF) in a human head model. In addition, the proposed method is used in conjunction with the generalized auto-calibrating partially parallel acquisitions (GRAPPA) reconstruction technique to enable rapid imaging of the two scans. Simulation results reported herein for 11-T (470-MHz) brain imaging applications show that the new method with GRAPPA reconstruction theoretically results in improved image quality and that the proposed combined hybrid MoM/FEM and FDTD technique is. suitable for high-field magnetic resonance imaging (MRI) numerical analysis.