Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

频率域电磁法中视电阻率全区定义

视电阻率定义在电法勘探中得到广泛的应用。目前在频率域电磁法中常常采用远区视电阻率定义。由于在频率域中收发距并非足够大,以致在低频段远区定义视电阻率往往不能反映地下电阻率值。我们定义全区电阻率及引入校正系数K,经二层三层模型及野外实测资料的计算,证明了方法的有效性。     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Crustal structure of the Eastern Alps along the TRANSALP profile from wide-angle seismic tomography

The objective of the TRANSALP project is an investigation of the Eastern Alps with regard to their deep structure and dynamic evolution. The core of the project is a 340-km-long seismic profile at 12°E between Munich and Venice. This paper deals with the P-wave velocity distribution as derived from active source travel time tomography. Our database consists of Vibroseis and explosion seismic travel times recorded at up to 100 seismological stations distributed in a 30-km-wide corridor along the profile. In order to derive a velocity and reflector model, we simultaneously inverted refractions and reflections using a derivative of a damped least squares approach for local earthquake tomography. 8000 travel time picks from dense Vibroseis recordings provide the basis for high resolution in the upper crust. Explosion seismic wide-angle reflection travel times constrain both deeper crustal velocities and structure of the crust–mantle boundary with low resolution. In the resulting model, the Adriatic crust shows significantly higher P-wave velocities than the European crust. The European Moho is dipping south at an angle of 7°. The Adriatic Moho dips north with a gentle inclination at shallower depths. This geometry suggests S-directed subduction. Azimuthal variations of the first-break velocities as well as observations of shear wave splitting reveal strong anisotropy in the Tauern Window. We explain this finding by foliations and laminations generated by lateral extrusion. Based on the P-wave model we also localized almost 100 local earthquakes recorded during the 2-month acquisition campaign in 1999. Seismicity patterns in the North seem related to the Inn valley shear zone, and to thrusting of Austroalpine units over European basement. The alignment of deep seismicity in the Trento-Vicenza region with the top of the Adriatic lower crust corroborates the suggestion of a deep thrust fault in the Southern Alps.