Acigöl rhyolite field,Central Anatolia (part 1): high-resolution dating of eruption episodes and zircon growth rates

Protracted pre-eruptive zircon residence is frequently detected in continental rhyolites and can conflict with thermal models, indicating briefer magma cooling durations if scaled to erupted volumes. Here, we present combined U-Th and (U-Th)/He zircon ages from the Acigöl rhyolite field (Central Anatolia, Turkey), which is part of a Quaternary bimodal volcanic complex. Unlike other geochronometers, this approach dates crystallization and eruption on the same crystals, allowing for internal consistency testing. Despite the overall longevity of Acigöl rhyolite volcanism and systematic trends of progressive depletion in compatible trace elements and decreasing zircon saturation temperatures, we find that zircon crystallized in two brief pulses corresponding to eruptions in the eastern and western part of the field during Middle and Late Pleistocene times, respectively. For Late Pleistocene zircon, resolvable differences exist between interior (average: 30.7 ± 0.9 ka; 1σ error) and rim (21.9 ± 1.3 ka) crystallization ages. These translate into radial crystal growth rates of ~10^?13 to 10^?14 cm/s, broadly consistent with those constrained by diffusion experiments. Rim crystallization and (U-Th)/He eruption ages (24.2 ± 0.4 ka) overlap within uncertainty. Evidence for brief zircon residence at Acigöl contrasts with many other rhyolite fields, suggesting that protracted zircon crystallization in, or recycling from, long-lived crystal mushes is not ubiquitous in continental silicic magma systems. Instead, the span of pre-eruptive zircon ages is consistent with autochthonous crystallization in individual small-volume magma batches that originated from basaltic precursors.     

The origin of oriented lakes in the Andean foreland, Parque Nacional Torres del Paine (Chilean Patagonia)

The Parque Nacional Torres Del Paine and surrounding area in the Magallanes foreland basin in Chilean Patagonia is the site for numerous lakes fed by glaciers and rivers in the Andean highlands to the west. The lakes are elongate and have conspicuously systematic orientations. We hypothesize that the origin of the oriented lakes lies in the fault system, composed of a right-lateral strike-slip fault set oriented 58° from north, a left-lateral strike-slip set oriented 87°, and a thrust fault set oriented 167°, that exists within the underlying rocks. To test this hypothesis quantitatively, we determined the shape and orientation of the lakes by fitting each lake with an ellipse of appropriate aspect ratio, and later with multiple ellipses consistent with the composite geometry of some lakes. We then examined the faults in the area in terms of their kinematics, orientation and distribution. The distribution of lake orientations showed three distinct groups which appear to correspond to the three main fault groups. For lakes fitted with multiple ellipses, the difference in means between the right-lateral, left-lateral, and thrust faults and their corresponding groups of lakes are 3.05°, 1.57°, and 5.17°. Using a Kolmogorov–Smirnov (K–S) statistical test to compare the orientations of faults with respect to the lakes suggests that there is not a strongly significant difference between the fault orientations and the corresponding lake groups. These results indicate that the faults have a profound control on the orientation, shape, and distribution of the lakes. We attribute this to faults and their damage zones being weaker and therefore prone to a faster rate of erosion, and to stress perturbations associated with discontinuous faults resulting in localized high density fracturing and surface subsidence. These results have implications for lake and drainage system morphologies in other foreland basins along the Andes and other similar settings.     

Geochemical and Sr–Nd–Pb isotopic compositions of the Eocene Dölek and Sariçiçek Plutons, Eastern Turkey: Implications for magma interaction in the genesis of high-K calc-alkaline granitoids in a post-collision extensional setting

The major and trace elements and Sr–Nd–Pb isotopes of the host rocks and the mafic microgranular enclaves (MME) gathered from the Dölek and Sariçiçek plutons, Eastern Turkey, were studied to understand the underlying petrogenesis and geodynamic setting. The plutons were emplaced at  43 Ma at shallow depths ( 5 to 9 km) as estimated from Al-in hornblende geobarometry. The host rocks consist of a variety of rock types ranging from diorite to granite (SiO₂ = 56.98–72.67 wt.%; Mg# = 36.8–50.0) populated by MMEs of gabbroic diorite to monzodiorite in composition (SiO₂ = 53.21–60.94 wt.%; Mg# = 44.4–53.5). All the rocks show a high-K calc-alkaline differentiation trend. Chondrite-normalized REE patterns are moderately fractionated and relatively flat [(La/Yb)_N = 5.11 to 8.51]. They display small negative Eu anomalies (Eu/Eu = 0.62 to 0.88), with enrichment of LILE and depletion of HFSE. Initial Nd–Sr isotopic compositions for the host rocks are ε_Nd(43 Ma) = − 0.6 to 0.8 and mostly I_Sr = 0.70482–0.70548. The Nd model ages (T_DM) vary from 0.84 to 0.99 Ga. The Pb isotopic ratios are (²⁰⁶Pb/²⁰⁴Pb) = 18.60–18.65, (²⁰⁷Pb/²⁰⁴Pb) = 15.61–15.66 and (²⁰⁸Pb/²⁰⁴Pb) = 38.69–38.85. Compared with the host rocks, the MMEs are relatively homogeneous in isotopic composition, with I_Sr ranging from 0.70485 to 0.70517, ε_Nd(43 Ma) − 0.1 to 0.8 and with Pb isotopic ratios of (²⁰⁶Pb/²⁰⁴Pb) = 18.58–18.64, (²⁰⁷Pb/²⁰⁴Pb) = 15.60–15.66 and (²⁰⁸Pb/²⁰⁴Pb) = 38.64–38.77. The MMEs have T_DM ranging from 0.86 to 1.36 Ga. The geochemical and isotopic similarities between the MMEs and their host rocks indicate that the enclaves are of mixed origin and are most probably formed by the interaction between the lower crust- and mantle-derived magmas. All the geochemical data, in conjunction with the geodynamic evidence, suggest that a basic magma derived from an enriched subcontinental lithospheric mantle, probably triggered by the upwelling of the asthenophere, and interacted with a crustal melt that originated from the dehydration melting of the mafic lower crust at deep crustal levels. Modeling based on the Sr–Nd isotope data indicates that  77–83% of the subcontinental lithospheric mantle involved in the genesis. Consequently, the interaction process played an important role in the genesis of the hybrid granitoid bodies, which subsequently underwent a fractional crystallization process along with minor amounts of crustal assimilation, en route to the upper crustal levels generating a wide variety of rock types ranging from diorite to granite in an extensional regime.     

The Effect of Regional Borders when Using the Gutenberg-Richter Model, Case Study: Western Anatolia

A probabilistic seismic hazard analysis (PSHA) for Western Anatolia is presented using the Gutenberg-Richter (G-R) frequency-magnitude relation. Since the modeling is sensitive to the location of seismotectonics boundaries, to use the information content of the observed earthquake data, as a general rule the borders of the affected area are extended. In this study, the effect of the region’s definition on the G-R model is debated on the Western Anatolian region, which is one of the most seismically active and rapidly deforming regions of the world throughout the ages. Calculations are carried out for two subregions and one combined region as a whole using the seismic catalog from 1900 to 2005. The data sets are determined by the region’s borders with the parameters computed according to these data sets by the least-squares and maximum likelihood methods, and then future predictions are estimated via these parameters. Comparing the results with historical earthquake records, most appropriate regional borders for Western Anatolia are defined, and for this new region G-R model parameters are obtained.     

Water Quality and Phosphorus Budget of Mogan Lake,Turkey

Mogan Lake is the largest recreational area near Ankara, which is the capital city of Turkey. Increased macrophyte growth in the water body and the present levels of urban development within the catchment are reflected in declining water quality and aesthetic deterioation. A study of the water quality of the lake was undertaken to quantify the variation of phosphorus, nitrogen, and chlorophyll-a concentrations during ice-free seasons from 1992 to 1994. Its total phosphorus and chlorophyll-a concentrations indicate a meso-eutrophic status. The total phosphorus budget of Mogan Lake was measured for a period of 22 months. The estimation of nutrient loading using Dillon-Rigler nutrient budget shows that an artificial load caused too high phosphorus inputs. The management implications of phosphorus loading and budget are discussed.     

地球介质的非定常参数黏弹模型

用开尔文黏弹模型对实验蠕变曲线和冰后期地面抬升观测曲线进行反演,得到相关的黏弹参数,发现所有的黏弹参数都随实验（或观测）时间长度发生有规律的变化, 表明黏弹参数是实验(或观测)时间的函数.在此基础上提出非定常参数的黏弹模型,将通常采用的剪切模量μ2、弛豫时间τ 和黏度η 相应地从常数修正为时间的函数μ2(t)、τ(t) 和η（t）,检验表明修正后的开尔文模型可以更好地符合实测蠕变曲线.用修正后的μ2(t)对软流层的黏度进行估算,其值约为4.6×1020 Pa·s.     

Reservoir operation using system dynamics under climate change impacts: a case study of Yamchi reservoir,Iran

This paper proposes a decision support system for Yamchi reservoir operation in semi-arid region of Iran. The paper consists of the following steps: Firstly, the potential impacts of climate change on the streamflow are predicted. The study then presents the projections of future changes in temperature and precipitation under A2 scenario using the LARS-WG downscaling model and under RCP2.6, RCP4.5, and RCP8.5 using the statistical downscaling model (SDSM) in the northwestern of Iran. To do so, a general circulation model of HadCM3 is downscaled by using the LARS-WG model. As a result, the average temperature, for the horizon 2030 (2011–2030), will increase by 0.77 °C and precipitation will decrease by 11 mm. Secondly, the downscaled variables are used as input to the artificial neural network to investigate the possible impact of climate change on the runoffs. Thirdly, the system dynamics model is employed to model different scenarios for reservoir operation using the Vensim software. System dynamics is an effective approach for understanding the behavior of complex systems. Simulation results demonstrate that the water shortage in different sectors (including agriculture, domestic, industry, and environmental users) will be enormously increased in the case of business-as-usual strategy. In this research, by providing innovative management strategies, including deficit irrigation, the vulnerability of reservoir operation is reduced. The methodology is evaluated by using different modeling tests which then motivates using the methodology for other arid/semi-arid regions.     

Spatial Statistical Properties of Pressure Solution Seams in Clastic Rocks in Southwest Ireland

Length, height, thickness and spacing measurements of pressure solution seams at outcrop, hand sample and thin section scale were taken from clastic rocks located in the southwest of Ireland. The lengths and spacings of pressure solution seams have similarly shaped (approximately log-hyperbolic) distributions at the observed scales suggesting that length and spacing distributions are scale-independent over the scales studied with a fractal dimension in the range of 1.4 to 1.6. Pressure solution seam lengths and thicknesses are related by a power-law and their spacings have a linear relationship to bed thickness. Although pressure solution seams are often considered as anticracks (forming under the same remote stresses as joints, but with opposite sign) we describe how the mechanism of pressure solution differs substantially from that of jointing. We use an existing mechanical model to show that stresses around pressure solution seam tips are much lower than those for joints under equal but opposite loading conditions. Pressure solution seams also have a decreasing tendency to lengthen as they grow, which is reflected in their length distributions. We propose that pressure solution seams, unlike joints, do not reach fracture saturation spacing because of transverse coalescence.     

Analysis of faulting in porous sandstones

Faults in porous sandstones occur in three forms: deformation bands about 1-mm thick and tens of m long and across which offsets are a few mm; zones of deformation bands constituted of many closely spaced deformation bands across which offsets are a few cm or dm; and slip surfaces, that is, distinct surfaces within zones of deformation bands across which offsets are a few m to a few tens of m. Deformation bands represent highly localized deformation; analogous localization within a field of homogeneous deformation is theoretically possible in inelastic materials with certain ranges of constitutive parameters. Crushing and consolidation of sandstone within a band cause the material there to become stiffer than the surrounding porous sandstone. A zone of deformation bands behaves mesoscopically much as a stiff inclusion in a soft matrix. According to the constitutive model assumed to investigate the formation of deformation bands, an instability can develop, and strain increments within the zone of deformation bands can become boundlessly large when the far-field stresses reach critical values. This instability is here associated with the formation of slip surfaces.