Exhumation Paths of High-Pressure-Low-Temperature Metamorphic Rocks from the Lycian Nappes and the Menderes Massif (SW Turkey): a Multi-Equilibrium Approach

The Menderes Massif and the overlying Lycian Nappes occupy anextensive area of SW Turkey where high-pressure–low-temperaturemetamorphic rocks occur. Precise retrograde P–T pathsreflecting the tectonic mechanisms responsible for the exhumationof these high-pressure–low-temperature rocks can be constrainedwith multi-equilibrium P–T estimates relying on localequilibria. Whereas a simple isothermal decompression is documentedfor the exhumation of high-pressure parageneses from the southernMenderes Massif, various P–T paths are observed in theoverlying Karaova Formation of the Lycian Nappes. In the uppermostlevels of this unit, far from the contact with the MenderesMassif, all P–T estimates depict cooling decompressionpaths. These high-pressure cooling paths are associated withtop-to-the-NNE movements related to the Akçakaya shearzone, located at the top of the Karaova Formation. This zoneof strain localization is a local intra-nappe contact that wasactive in the early stages of exhumation of the high-pressurerocks. In contrast, at the base of the Karaova Formation, alongthe contact with the Menderes Massif, P–T calculationsshow decompressional heating exhumation paths. These paths areassociated with severe deformation characterized by top-to-the-eastshearing related to a major shear zone (the Gerit shear zone)that reflects late exhumation of high-pressure parageneses underwarmer conditions. KEY WORDS: exhumation; high-pressure–low-temperature metamorphism; multi-equilibrium P–T estimates; Lycian Nappes; Menderes Massif     

Phylogeography of the pelagic fish Seriola lalandi at different scales: confirmation of inter-ocean population structure and evaluation of southern African genetic diversity

The study investigated the global and regional phylogeography of the yellowtail kingfish Seriola lalandi by examining genetic diversity and population genetic structure of this species at inter-and intra-ocean level and on a regional scale. DNA fragments of two mitochondrial genes, cytochrome b (Cytb) and cytochrome c oxidase subunit I (COI) and one nuclear gene, recombination activating gene 1 (RAG1), were sequenced to investigate the global-scale phylogeography of this species. The population genetic structure within the South Pacific, as well as along the South African coastline, was examined further using six microsatellite markers. Three distinct clades were identified for S. lalandi, which correspond with previously described subspecies of the North-East Pacific, North-West Pacific and the Southern Hemisphere. Within the latter, additional divergence was observed between the South Pacific and the South-East Atlantic regions. Divergence estimates were indicative of a Pacific origin for S. lalandi populations, because of Pleistocene vicariant events. Microsatellite analyses revealed overall significant genetic differentiation between South African and South Pacific samples. This corroborates recent findings on the global phylogeography of the species. No population differentiation was observed within South Africa, indicating high levels of gene flow.     

Brittleness index predictions from Lower Barnett Shale well-log data applying an optimized data matching algorithm at various sampling densities

The capability of accurately predicting mineralogical brittleness index(BI)from basic suites of well logs is desir-able as it provides a useful indicator of the fracability of tight formations.Measuring mineralogical components in rocks is expensive and time consuming.However,the basic well log curves are not well correlated with BI so correlation-based,machine-learning methods are not able to derive highly accurate BI predictions using such data.A correlation-free,optimized data-matching algorithm is configured to predict BI on a supervised basis from well log and core data available from two published wells in the Lower Barnett Shale Formation(Texas).This transparent open box(TOB)algorithm matches data records by calculating the sum of squared errors be-tween their variables and selecting the best matches as those with the minimum squared errors.It then applies optimizers to adjust weights applied to individual variable errors to minimize the root mean square error(RMSE)between calculated and predicted(BI).The prediction accuracy achieved by TOB using just five well logs(Gr,pb,Ns,Rs,Dt)to predict BI is dependent on the density of data records sampled.At a sampling density of about one sample per 0.5 ft BI is predicted with RMSE～0.056 and R2～0.790.At a sampling density of about one sample per 0.1 ft BI is predicted with RMSE～0.008 and R2～0.995.Adding a stratigraphic height index as an additional(sixth)input variable method improves BI prediction accuracy to RMSE～0.003 and R2～0.999 for the two wells with only 1 record in 10,000 yielding a BI prediction error of＞±0.1.The model has the potential to be applied in an unsupervised basis to predict BI from basic well log data in surrounding wells lacking mineralogical measure-ments but with similar lithofacies and burial histories.The method could also be extended to predict elastic rock properties in and seismic attributes from wells and seismic data to improve the precision of brittleness index and fracability mapping spatially.     

海河流域近60 a降水极值的频率分析及时空分布特征

以海河流域为研究区域,利用MK和F检验对站点的年最大日降雨序列进行趋势、突变和跳跃分析。基于水文气象分区线性矩法进行一致区的划分、最优分布的选择和降水极值的频率估计值计算,并分析其空间分布特征。结果表明:MK趋势和突变检验显示只有8个站点(22.9%)呈现显著下降趋势,通过了信度为0.1的显著性水平检验,特别是京津唐区域具有显著的下降趋势,突变时段主要发生在1980-1990年。均值和方差的跳跃性显示大部分站点都呈现出显著的向下跳跃,主要分布在流域的滦河子流域和北三河水系;趋势和跳跃的综合分析能够对降水极值有可能引起的旱涝灾害进行更全面合理的判断。不同重现期下的频率估计值的空间分布总体趋势一致,从东南到西北、从沿海到内陆逐渐减少,并与地形表现出很好的一致性;降水极值空间分布的中心主要集中在滦河子流域的遵化和青龙附近;大部分站点50 a一遇的估计值和AMP序列的最大观测值能够保持较好的一致性,间接反映了估计值的合理准确性。     

Applications of satellite data in the study of monsoon variability

This paper discusses the use of satellite data for studying climate change, with particular emphasis on the inter-annual variability of the Indian southwest monsoon. Precipitation estimates made from INSAT-1B radiance data are shown to bring out the variations that occurred in the monsoon rainfall of 1987 and 1988. Outgoing Longwave Radiation derived from INSAT-1B shows good correspondence with precipitation patterns.     

Thermal evolution of the Tertiary Shimanto Belt, Muroto Peninsula, Shikoku, Japan

Abstract The Shimanto accretionary complex on the Muroto Peninsula of Shikoku comprises two major units of Tertiary strata: the Murotohanto Sub-belt (Eocene-Oligocene) and the Nabae Sub-belt (Oligocene-Miocene). Both sub-belts have been affected by thermal overprints following the peak of accretion-related deformation. Palaeotemperatures for the entire Tertiary section range from ～ 140 to 315°C, based upon mean vitrinite reflectance values of 0.9–5.0%R_m. Values of illite crystallinity index are consistent with conditions of advanced diagenesis and anchimetamorphism. Illite/mica b₀ lattice dimensions indicate that burial pressures were probably no greater than 2.5kbar. In general, levels of thermal maturity are higher for the Murotohanto Sub-belt than for the Nabae Sub-belt. The Eocene-Oligocene strata also display a spatial decrease in thermal maturity from south to north and this pattern probably was caused by regional-scale differential uplift following peak heating. Conversely, the palaeothermal structure within the Nabae Sub-belt is fairly uniform, except for the local effects of mafic intrusions at the tip of Cape Muroto. There is a paleotemperature difference of ～ 90°C across the boundary between the Murotohanto and Nabae Sub-belts (Shiina-Narashi fault), and this contrast is consistent with approximately 1200 m of post-metamorphic vertical offset. Subduction prior to Middle Miocene probably involved the Kula or fused Kula-Pacific plate and the background geothermal gradient during the Eocene-Oligocene phase of accretion was ～ 30–35°C/km. Rapid heating of the Shimanto Belt evidently occurred immediately after a Middle Miocene reorganization of the subduction boundary. Hot oceanic lithosphere from the Shikoku Basin first entered the subduction zone at ～ 15 Ma; this event also coincided with the opening of the Sea of Japan and the rapid clockwise rotation of southwest Japan. The background geothermal gradient at that time was ～ 70°C/km. Whether or not all portions of the inherited (Eocene-Oligocene) palaeothermal structure were overprinted during the Middle Miocene remains controversial.     

A computational method for approximating a Darcy–Stokes system governing a vuggy porous medium

We develop and analyze a mixed finite element method for the solution of an elliptic system modeling a porous medium with large cavities, called vugs. It consists of a second-order elliptic (i.e., Darcy) equation on part of the domain coupled to a Stokes equation on the rest of the domain, and a slip boundary condition (due to Beavers–Joseph–Saffman) on the interface between them. The tangential velocity is not continuous on the interface. We consider a 2-D vuggy porous medium with many small cavities throughout its extent, so the interface is not isolated. We use a certain conforming Stokes element on rectangles, slightly modified near the interface to account for the tangential discontinuity. This gives a mixed finite element method for the entire Darcy–Stokes system with a regular sparsity pattern that is easy to implement, independent of the vug geometry, as long as it aligns with the grid. We prove optimal global first-order L ² convergence of the velocity and pressure, as well as the velocity gradient in the Stokes domain. Numerical results verify these rates of convergence and even suggest somewhat better convergence in certain situations. Finally, we present a lower dimensional space that uses Raviart–Thomas elements in the Darcy domain and uses our new modified elements near the interface in transition to the Stokes elements.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: the discrete-in-time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.