Estimation of the correlation structure of crustal velocity heterogeneity from seismic reflection data

Numerous sources of evidence point to the fact that heterogeneity within the Earth's deep crystalline crust is complex and hence may be best described through stochastic rather than deterministic approaches. As seismic reflection imaging arguably offers the best means of sampling deep crustal rocks in situ, much interest has been expressed in using such data to characterize the stochastic nature of crustal heterogeneity. Previous work on this problem has shown that the spatial statistics of seismic reflection data are indeed related to those of the underlying heterogeneous seismic velocity distribution. As of yet, however, the nature of this relationship has remained elusive due to the fact that most of the work was either strictly empirical or based on incorrect methodological approaches. Here, we introduce a conceptual model, based on the assumption of weak scattering, that allows us to quantitatively link the second-order statistics of a 2-D seismic velocity distribution with those of the corresponding processed and depth-migrated seismic reflection image. We then perform a sensitivity study in order to investigate what information regarding the stochastic model parameters describing crustal velocity heterogeneity might potentially be recovered from the statistics of a seismic reflection image using this model. Finally, we present a Monte Carlo inversion strategy to estimate these parameters and we show examples of its application at two different source frequencies and using two different sets of prior information. Our results indicate that the inverse problem is inherently non-unique and that many different combinations of the vertical and lateral correlation lengths describing the velocity heterogeneity can yield seismic images with the same 2-D autocorrelation structure. The ratio of all of these possible combinations of vertical and lateral correlation lengths, however, remains roughly constant which indicates that, without additional prior information, the aspect ratio is the only parameter describing the stochastic seismic velocity structure that can be reliably recovered.

A comparative long-memory Analysis between Spanish, Mexican and U.S. interest rates

Evidence exists that many natural facts are described better as a fractal. Although fractals are very useful for describing nature, it is also appropiate to review the concept of random fractal in finance. Due to the extraordinary importance of Brownian motion in physics, chemistry or biology, we will consider the generalization that supposes fractional Brownian motion introduced by Mandelbrot.The main goal of this work is to analyse the existence of long range dependence in instantaneous forward rates of different financial markets. Concretelly, we perform an empirical analysis on the Spanish, Mexican and U.S. interbanking interest rate. We work with three time series of daily data corresponding to 1 day operations from 28th March 1996 to 21st May 2002. From among all the existing tests on this matter we apply the methodology proposed in Taqqu, Teverovsky and Willinger (1995).

Reconstrucción grupal de cuentos: una instancia para la producción de narrativa

Este artículo presenta los resultados del estudio de las reconstrucciones grupales de los cuentos leídos en salas de jardín de infantes. Se analizó cualitativamente empleando el método comparativo constante, un corpus de 30 situaciones de lectura de cuentos, de las cuales se consideró el momento posterior a la lectura en el que se reconstruye el cuento grupalmente. Los resultados permitieron distinguir diferentes tipos de reconstrucciones de los cuentos a través de la interacción, de acuerdo a la cantidad y calidad de participación de la maestra y de los niños y de información del texto recuperada: reconstrucciones incompletas, reconstrucciones descriptivas, reconstrucciones-cuestionario, reconstrucciones-debate (focalizadas en la comprensión), reconstrucciones monopolizadas y reconstrucciones completas compartidas. Cada uno de estos tipos, que se describe en función de las operaciones cognitivas que demanda de sus participantes, da cuenta de modos diferentes de apoyar la producción de narrativa oral.

La simplificación impositiva en Cuba: Sistemas objetivos de determinación de bases y aglutinación de varios impuestos en uno sólo impuesto

El presente estudio analiza el Impuesto sobre ingresos brutos integrado en el sistema tributario cubano al que están sujetas las sociedades que tributan por el Impuesto sobre Utilidades. Consiste en un sistema que sigue un método de determinación de bases mediante la utilización de índices objetivos referidos a magnitudes físicas, coeficientes o rendimientos de operaciones reales al alcance del sujeto pasivo y de la Administración tributaria

Increased brain perfusion contrast with T2 -prepared intravoxel incoherent motion (T2prep IVIM) MRI.

The feasibility to measure brain perfusion using intravoxel incoherent motion (IVIM) MRI has been reported recently with currently clinically available technology. The method is intrinsically local and quantitative, but is contaminated by partial volume effects with cerebrospinal fluid (CSF). Signal from CSF can be suppressed by a 180° inversion recovery (180°-IR) magnetization preparation, but this also leads to strong suppression of blood and brain tissue signal. Here, we take advantage of the different T2 relaxations of blood and brain relative to CSF, and implement a T2 -prepared IVIM (T2prep IVIM) inversion recovery acquisition, which permits a recovery of between 43% and 57% of arterial and venous blood magnetization at excitation time compared with the theoretical recovery of between 27% and 30% with a standard 180°-IR. We acquired standard IVIM (IVIM), T2prep IVIM and dynamic susceptibility contrast (DSC) images at 3 T using a 32-multichannel receiver head coil in eight patients with known large high-grade brain tumors. We compared the contrast and contrast-to-noise ratio obtained in the corresponding cerebral blood volume images quantitatively, as well as subjectively by two neuroradiologists. Our findings suggest that quantitative cerebral blood volume contrast and contrast-to-noise ratio, as well as subjective lesion detection, contrast quality and diagnostic confidence, are increased with T2prep IVIM relative to IVIM and DSC.

MRI-based prostate brachytherapy seed localization.

A magnetic resonance imaging (MRI) pulse sequence and a corresponding image processing algorithm to localize prostate brachytherapy seeds during or after therapy are presented. Inversion-Recovery with ON-resonant water suppression (IRON) is an MRI methodology that generates positive contrast in regions of magnetic field susceptibility, as created by prostate brachytherapy seeds. Phantoms comprising of several materials found in brachytherapy seeds were created to assess the usability of the IRON pulse sequence for imaging seeds. Resulting images show that seed materials are clearly visible with high contrast using IRON, agreeing with theoretical predictions. A seed localization algorithm to process IRON images demonstrates the potential of this imaging technique for seed localization and dosimetry.

Juxtaposition of melt impregnation and high-temperature shear zones in the upper mantle; field and petrological constraints from the Lanzo peridotite (Northern Italy)

Results of a field and microstructural study between the northern and the central bodies of the Lanzo plagioclase peridotite massif (NW Italy) indicate that the spatial distribution of deformation is asymmetric across kilometre-scale mantle shear zones. The southwestern part of the shear zone (footwall) shows a gradually increasing degree of deformation from porphyroclastic peridotites to mylonite, whereas the northeastern part (hanging wall) quickly grades into weakly deformed peridotites. Discordant gabbroic and basaltic dykes are asymmetrically distributed and far more abundant in the footwall of the shear zone. The porphyroclastic peridotite displays porphyroclastic zones and domains of igneous crystallization whereas mylonites are characterized by elongated porphyroclasts, embedded between fine-grained, polycrystalline bands of olivine, plagioclase, clinopyroxene, orthopyroxene, spinel, rare titanian pargasite, and domains of recrystallized olivine. Two types of melt impregnation textures have been found: (1) clinopyroxene porphyroclasts incongruently reacted with migrating melt to form orthopyroxene plagioclase; (2) olivine porphyroclasts are partially replaced by interstitial orthopyroxene. The meltrock reaction textures tend to disappear in the mylonites, indicating that deformation in the mylonite continued under subsolidus conditions. The pyroxene chemistry is correlated with grain size. High-Al pyroxene cores indicate high temperatures (11001030C), whereas low-Al neoblasts display lower final equilibration temperatures (860C). The spinel Cr-number [molar Cr/(Cr Al)] and TiO2 concentrations show extreme variability covering almost the entire range known from abyssal peridotites. The spinel compositions of porphyroclastic peridotites from the central body are more variable than spinel from mylonite, mylonite with ultra-mylonite bands, and porphyroclastic rocks of the northern body. The spinel compositions probably indicate disequilibrium and would favour rapid cooling, and a faster exhumation of the central peridotite body, relative to the northern one. Our results indicate that melt migration and high-temperature deformation are juxtaposed both in time and space. Meltrock reaction may have caused grain-size reduction, which in turn led to localization of deformation. It is likely that melt-lubricated, actively deforming peridotites acted as melt focusing zones, with permeabilities higher than the surrounding, less deformed peridotites. Later, under subsolidus conditions, pinning in polycrystalline bands in the mylonites inhibited substantial grain growth and led to permanent weak zones in the upper mantle peridotite, with a permeability that is lower than in the weakly deformed peridotites. Such an inversion in permeability might explain why actively deforming, fine-grained peridotite mylonite acted as a permeability barrier and why ascending mafic melts might terminate and crystallize as gabbros along actively deforming shear zones. Melt-lubricated mantle shear zones provide a mechanism for explaining the discontinuous distribution of gabbros in oceancontinent transition zones, oceanic core complexes and ultraslow-spreading ridges.

Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petro-physical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine the parameters of a Gaussian mixture model that defines zonal geometries from joint-inversion tomograms. The resulting zones are used to estimate representative geophysical parameters of each zone, which are then used for field-scale petrophysical analysis. A synthetic study demonstrated how joint inversion of seismic and radar traveltimes and electrical resistance tomography (ERT) data greatly reduces misclassification of zones (down from 21.3% to 3.7%) and improves the accuracy of retrieved zonal parameters (from 1.8% to 0.3%) compared to individual inversions. We applied our scheme to a data set collected in northeastern Switzerland to delineate lithologic subunits within a gravel aquifer. The inversion models resolve three principal subhorizontal units along with some important 3D heterogeneity. Petro-physical analysis of the zonal parameters indicated approximately 30% variation in porosity within the gravel aquifer and an increasing fraction of finer sediments with depth.

Direct in vitro comparison of six three-dimensional positive contrast methods for susceptibility marker imaging.

PURPOSE: To compare different techniques for positive contrast imaging of susceptibility markers with MRI for three-dimensional visualization. As several different techniques have been reported, the choice of the suitable method depends on its properties with regard to the amount of positive contrast and the desired background suppression, as well as other imaging constraints needed for a specific application. MATERIALS AND METHODS: Six different positive contrast techniques are investigated for their ability to image at 3 Tesla a single susceptibility marker in vitro. The white marker method (WM), susceptibility gradient mapping (SGM), inversion recovery with on-resonant water suppression (IRON), frequency selective excitation (FSX), fast low flip-angle positive contrast SSFP (FLAPS), and iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL) were implemented and investigated. RESULTS: The different methods were compared with respect to the volume of positive contrast, the product of volume and signal intensity, imaging time, and the level of background suppression. Quantitative results are provided, and strengths and weaknesses of the different approaches are discussed. CONCLUSION: The appropriate choice of positive contrast imaging technique depends on the desired level of background suppression, acquisition speed, and robustness against artifacts, for which in vitro comparative data are now available.

Flow targeted 3D steady-state free-precession coronary MR angiography: comparison of three different imaging approaches.

PURPOSE: To compare 3 different flow targeted magnetization preparation strategies for coronary MR angiography (cMRA), which allow selective visualization of the vessel lumen. MATERIAL AND METHODS: The right coronary artery of 10 healthy subjects was investigated on a 1.5 Tesla MR system (Gyroscan ACS-NT, Philips Healthcare, Best, NL). A navigator-gated and ECG-triggered 3D radial steady-state free-precession (SSFP) cMRA sequence with 3 different magnetization preparation schemes was performed referred to as projection SSFP (selective labeling of the aorta, subtraction of 2 data sets), LoReIn SSFP (double-inversion preparation, selective labeling of the aorta, 1 data set), and inflow SSFP (inversion preparation, selective labeling of the coronary artery, 1 data set). Signal-to-noise ratio (SNR) of the coronary artery and aorta, contrast-to-noise ratio (CNR) between the coronary artery and epicardial fat, vessel length and vessel sharpness were analyzed. RESULTS: All cMRA sequences were successfully obtained in all subjects. Both projection SSFP and LoReIn SSFP allowed for selective visualization of the coronary arteries with excellent background suppression. Scan time was doubled in projection SSFP because of the need for subtraction of 2 data sets. In inflow SSFP, background suppression was limited to the tissue included in the inversion volume. Projection SSFP (SNR(coro): 25.6 +/- 12.1; SNR(ao): 26.1 +/- 16.8; CNR(coro-fat): 22.0 +/- 11.7) and inflow SSFP (SNR(coro): 27.9 +/- 5.4; SNR(ao): 37.4 +/- 9.2; CNR(coro-fat): 24.9 +/- 4.8) yielded significantly increased SNR and CNR compared with LoReIn SSFP (SNR(coro): 12.3 +/- 5.4; SNR(ao): 11.8 +/- 5.8; CNR(coro-fat): 9.8 +/- 5.5; P < 0.05 for both). Longest visible vessel length was found with projection SSFP (79.5 mm +/- 18.9; P < 0.05 vs. LoReIn) whereas vessel sharpness was best in inflow SSFP (68.2% +/- 4.5%; P < 0.05 vs. LoReIn). Consistently good image quality was achieved using inflow SSFP likely because of the simple planning procedure and short scanning time. CONCLUSION: Three flow targeted cMRA approaches are presented, which provide selective visualization of the coronary vessel lumen and in addition blood flow information without the need of contrast agent administration. Inflow SSFP yielded highest SNR, CNR and vessel sharpness and may prove useful as a fast and efficient approach for assessing proximal and mid vessel coronary blood flow, whereas requiring less planning skills than projection SSFP or LoReIn SSFP.

Using spatially distributed parameters and multi-response objective functions to solve parameterization of complex applications of semi-distributed hydrological models

Application of semi-distributed hydrological models to large, heterogeneous watersheds deals with several problems. On one hand, the spatial and temporal variability in catchment features should be adequately represented in the model parameterization, while maintaining the model complexity in an acceptable level to take advantage of state-of-the-art calibration techniques. On the other hand, model complexity enhances uncertainty in adjusted model parameter values, therefore increasing uncertainty in the water routing across the watershed. This is critical for water quality applications, where not only streamflow, but also a reliable estimation of the surface versus subsurface contributions to the runoff is needed. In this study, we show how a regularized inversion procedure combined with a multiobjective function calibration strategy successfully solves the parameterization of a complex application of a water quality-oriented hydrological model. The final value of several optimized parameters showed significant and consistentdifferences across geological and landscape features. Although the number of optimized parameters was significantly increased by the spatial and temporal discretization of adjustable parameters, the uncertainty in water routing results remained at reasonable values. In addition, a stepwise numerical analysis showed that the effects on calibration performance due to inclusion of different data types in the objective function could be inextricably linked. Thus caution should be taken when adding or removing data from an aggregated objective function.

New magnetotelluric data trough the boundary between the Ossa Morena and Centroiberian Zones

The south-western part of the Iberian Peninsula, including the southern branch of the Iberian Massif, has recently been the subject of several magnetotelluric (MT) studies. This area is made up of three different tectonic terranes: the South Portuguese Zone (SPZ), the Ossa Morena Zone (OMZ) and the Central Iberian Zone (CIZ). The boundaries between these zones are considered to be sutures, which appear as high electrical conductivity anomalies in the MT surveys. The OMZ is characterised by a conductive layer at middle-lower crustal levels. To investigate the continuity of this conductive layer into the CIZ, a new MT profile was carried out. This 75-km long ENE profile goes through the boundary between the OMZ and the CIZ. The results of a two-dimensional magnetotelluric inversion revealed a high-conductivity anomaly in the transition OMZ/CIZ (the so-called Central Unit), which is interpreted as due to interconnected graphite along shear planes. High-conductivity anomalies appeared in the middle crust of the CIZ, whose geometry and location are consistent with the conductive layer previously found in the OMZ, thus confirming the prolongation of the conductive layer into the CIZ. The top of this layer correlated spatially with a broad reflector detected by a seismic profile previously acquired in the same area. This, together with other geological and petrological evidence, points to a common origin for both features.

Deep crustal electromagnetic structure of Bhuj earthquake region (India) and its implications

The existence of fluids and partial melt in the lower crust of the seismically active Kutch rift basin (on the western continental margin of India) owing to underplating has been proposed in previous geological and geophysical studies. This hypothesis is examined using magnetotelluric (MT) data acquired at 23 stations along two profiles across Kutch Mainland Uplift and Wagad Uplift. A detailed upper crustal structure is also presented using twodimensional inversion of MT data in the Bhuj earthquake (2001) area. The prominent boundaries of reflection in the upper crust at 5, 10 and 20 km obtained in previous seismic reflection profiles correlate with conductive structures in our models. The MT study reveals 1-2 km thick Mesozoic sediments under the Deccan trap cover. The Deccan trap thickness in this region varies from a few meters to 1.5 km. The basement is shallow on the northern side compared to the south and is in good agreement with geological models as well as drilling information. The models for these profiles indicate that the thickness of sediments would further increase southwards into the Gulf of Kutch. Significant findings of the present study indicate 1) the hypocentre region of the earthquake is devoid of fluids, 2) absence of melt (that is emplaced during rifting as suggested from the passive seismological studies) in the lower crust and 3) a low resistive zone in the depth range of 5-20 km. The present MT study rules out fluidsand melt (magma) as the causative factors that triggered the Bhuj earthquake. The estimated porosity value of 0.02% will explain 100-500 ohm·m resistivity values observed in the lower crust. Based on the seismic velocities and geochemical studies, presence of garnet is inferred. The lower crust consists of basalts - probably generated by partial melting of metasomatised garnet peridotite at deeper depths in the lithosphere - and their composition might be modified by reaction with the spinel peridotites.

Combined use of the GGSFT data base and on Board Marine Collected Data to Model the Moho Beneath the Powell Basin, Antarctica

The Powell Basin is a small oceanic basin located at the NE end of the Antarctic Peninsula developed during the Early Miocene and mostly surrounded by the continental crusts of the South Orkney Microcontinent, South Scotia Ridge and Antarctic Peninsula margins. Gravity data from the SCAN 97 cruise obtained with the R/V Hespérides and data from the Global Gravity Grid and Sea Floor Topography (GGSFT) database (Sandwell and Smith, 1997) are used to determine the 3D geometry of the crustal-mantle interface (CMI) by numerical inversion methods. Water layer contribution and sedimentary effects were eliminated from the Free Air anomaly to obtain the total anomaly. Sedimentary effects were obtained from the analysis of existing and new SCAN 97 multichannel seismic profiles (MCS). The regional anomaly was obtained after spectral and filtering processes. The smooth 3D geometry of the crustal mantle interface obtained after inversion of the regional anomaly shows an increase in the thickness of the crust towards the continental margins and a NW-SE oriented axis of symmetry coinciding with the position of an older oceanic spreading axis. This interface shows a moderate uplift towards the western part and depicts two main uplifts to the northern and eastern sectors.

Site-specific coupling between vascular wall thickness and function: an observational MRI study of vessel wall thickening and stiffening in hypertension.

OBJECTIVES: The objective of this study was to evaluate associations between aortic pulse wave velocity (PWV) and aortic and carotid vessel wall thickness (VWT) using cardiovascular magnetic resonance imaging (MRI) in patients with hypertension as compared with healthy adult volunteers. MATERIALS AND METHODS: Local medical ethics approval was obtained and the participants gave informed consent. Fifteen patients with hypertension (5 men and 10 women; mean [SD] age, 49 [14] years) and 15 age- and sex-matched healthy volunteers were prospectively included and compared. All participants underwent MRI examination for measuring aortic and carotid VWT and aortic PWV with well-validated MRI techniques at 1.5- and 3-T MRI systems: PWV was assessed from velocity-encoded MRI and VWT was assessed by using dual-inversion black-blood gradient-echo imaging techniques. Paired t tests were used for testing differences between the volunteers and the patients and Pearson correlation (r) and univariable and multivariable stepwise linear regression analyses were used to test associations between aortic and carotid arterial wall thickness and stiffness. RESULTS: Mean values for aortic PWV and aortic and carotid VWT (indexed for body surface area [BSA]) were all significantly higher in patients with hypertension as compared with the healthy volunteers (ie, aortic PWV, 7.0 ± 1.4 m/s vs 5.7 ± 1.3 m/s; aortic VWT/BSA, 0.12 ± 0.03 mL/m vs 0.10 ± 0.03 mL/m; carotid VWT/BSA, 0.04 ± 0.01 mL/m vs 0.03 ± 0.01 mL/m; all P < 0.01). Aortic PWV was highly correlated with aortic VWT/BSA (r = 0.76 and P = 0.002 in the patients vs r = 0.63 and P = 0.02 in the volunteers), and in the patients, aortic PWV was moderately correlated with carotid VWT/BSA (r = 0.50; P = 0.04). In the volunteers, correlation between aortic PWV and carotid VWT/BSA was not significant (r = 0.40; P = 0.13). In addition, aortic VWT/BSA was significantly correlated with carotid VWT/BSA, in both the patients (r = 0.60; P = 0.005) and volunteers (r = 0.57; P = 0.007). CONCLUSIONS: In the patients with hypertension and the healthy volunteers, the aortic PWV is associated more strongly with aortic wall thickness than with carotid wall thickness, reflecting site-specific coupling between vascular wall thickness and function.