Evaluation of the temporal variations of groundwater storage and its interactions with climatic variables using GRACE data and hydrological models: A study from Turkey

Monitoring of the fluctuations of groundwater storage is particularly important in arid and semi-arid regions where water scarcity brings about various challenges. Remote sensing data and techniques play a preponderant role in developing solutions to environmental problems. The launch of Gravity Recovery and Climate Experiment (GRACE) satellites has eased the remote monitoring and evaluation of groundwater resources with an unprecedented precision over large scales. Within the scope of the current study, the latest release (RL06) of GRACE mass concentrations (Mascons) from Jet Propulsion Laboratory (JPL) dataset as well as Global Land Data Assimilation System (GLDAS) models of Noah and Catchment Land Surface Model (CLSM) were used to provide Groundwater Storage Anomalies (GWSA) over Turkey. The temporal interactions of the estimated GWSA with the climatic variables of precipitation and temperature (derived from the reanalysis datasets of CHELSA [Climatologies at High resolution for the Earth's Land Surface Areas] and FLDAS [the Famine Early Warning Systems Network Land Data Assimilation System], respectively) were investigated statistically. The results suggest that there is a descending trend (from 2003 to 2016) for Terrestrial Water Storage Anomalies (TWSA) and GWSA over Turkey with a total loss of 11 and 6 cm of water, respectively. The statistical analysis results also indicate that the monthly variations of GWSA over Turkey are highly correlated with precipitation and temperature at 2-month lag. The analysis of the climatology (long-term) values of monthly GWSA, precipitation and temperature also revealed high agreement between the variables.     

Optimal Compensation and Implementation for Adaptive Optics Systems

This paper develops a compensation algorithm based on Linear–Quadratic–Gaussian (LQG) control system design whose parameters are determined (in part) by a model of the atmosphere. The model for the atmosphere is based on the open-loop statistics of the atmosphere as observed by the wavefront sensor, and is identified from these using an auto-regressive, moving average (ARMA) model. The (LQG) control design is compared with an existing compensation algorithm for a simulation developed at ESO that represents the operation of MACAO adaptive optics system on the 8.2 m telescopes at Paranal, Chile. This revised version was published online in July 2006 with corrections to the Cover Date.     

Rate of strike-slip motion on the Amanos Fault (Karasu Valley, southern Turkey) constrained by K–Ar dating and geochemical analysis of Quaternary basalts

The left-lateral Amanos Fault follows a 200-km-long and up to 2-km-high escarpment that bounds the eastern margin of the Amanos mountain range and the western margin of the Karasu Valley in southern Turkey, just east of the northeastern corner of the Mediterranean Sea. Regional kinematic models have reached diverse conclusions as to the role of this fault in accommodating relative motion between either the African and Arabian, Turkish and African, or Turkish and Arabian plates. Local studies have tried to estimate its slip rate by K–Ar dating Quaternary basalts that erupted within the Amanos Mountains, flowed across it into the Karasu Valley, and have since become offset. However, these studies have yielded a wide range of results, ranging from 0.3 to 15 mm a⁻¹, which do not allow the overall role and significance of this fault in accommodating crustal deformation to be determined. We have used the Cassignol K–Ar method to date nine Quaternary basalt samples from the vicinity of the southern part of the Amanos Fault. These basalts exhibit a diverse chemistry, which we interpret as a consequence varying degrees of partial melting of their source combined with variable crustal contamination. This dating allows us to constrain the Quaternary slip rate on the Amanos fault to 1.0 to 1.6 mm a⁻¹. The dramatic discrepancies between past estimates of this slip rate are partly due to technical difficulties in K–Ar dating of young basalts by isotope dilution. In addition, previous studies at the key locality of Hacılar have unwittingly dated different, chemically distinct, flow units of different ages that are juxtaposed. This low slip rate indicates that, at present, the Amanos Fault takes up a small proportion of the relative motion between the African and Arabian plates, which is transferred southward to the Dead Sea Fault Zone. It also provides strong evidence against the long-standing view that its slip continues offshore to the southwest along a hypothetical left-lateral fault zone located south of Cyprus.     

Estimation of fundamental periods of shear‐wall dominant building structures

Shear‐wall dominant multistorey reinforced concrete structures, constructed by using a special tunnel form technique are commonly built in countries facing a substantial seismic risk, such as Chile, Japan, Italy and Turkey. In spite of their high resistance to earthquake excitations, current seismic code provisions including the Uniform Building Code (International Conference of Building Officials, Whittier, CA, 1997) and the Turkish Seismic Code (Specification for Structures to be Built in Disaster Areas, Ankara, Turkey, 1998) present limited information for their design criteria. In this study, consistency of equations in those seismic codes related to their dynamic properties are investigated and it is observed that the given empirical equations for prediction of fundamental periods of this specific type of structures yield inaccurate results. For that reason, a total of 80 different building configurations were analysed by using three‐dimensional finite‐element modelling and a set of new empirical equations was proposed. The results of the analyses demonstrate that given formulas including new parameters provide accurate predictions for the broad range of different architectural configurations, roof heights and shear‐wall distributions, and may be used as an efficient tool for the implicit design of these structures. Copyright © 2003 John Wiley & Sons, Ltd.     

A finite element model for the diffusion—convection equation with applications to air pollution problems

A method of solution for the diffusion-convection equation in one spatial dimension has been developed previously. The generalization of this method into two-spatial dimensions will be presented. The method employs space-time volume elements with edges joining the nodes at subsequent time levels oriented along the characteristics of the associated pure convection problem. The accuracy and utility of the method are demonstrated by solving several examples and results are compared with the exact solutions.     

Impact of different sampling rates on precise point positioning performance using online processing service

In this study, the effect of different sampling rates (i.e. observation recording interval) on the Precise Point Positioning (PPP) solutions in terms of accuracy was investigated. For this purpose, a field test was carried out in ?orum province, Turkey, on 11 September 2019. Within this context, a Geodetic Point (GP) was established and precisely coordinated. A static GNSS measurement was occupied on the GP for about 4-hour time at 0.10 second (s)/10 Hz measurement intervals with the Trimble R10 geodetic grade GNSS receiver. The original observation file was converted to RINEX format and then decimated into the different data sampling rates as 0.2 s, 0.5 s, 1 s, 5 s, 10 s, 30 s, 60 s, and 120 s. All these RINEX observation files were submitted to the Canadian Spatial Reference System-Precise Point Positioning (CSRS-PPP) online processing service the day after the data collection date by choosing both static and kinematic processing options. In this way, PPP-derived static coordinates, and the kinematic coordinates of each measurement epoch were calculated. The PPP-derived coordinates obtained from each decimated sampling intervals were compared to known coordinates of the GP for northing, easting, 2D position, and height components. According to the static and kinematic processing results, high data sampling rates did not change the PPP solutions in terms of accuracy when compared to the results obtained using lower sampling rates. The results of this study imply that it was not necessary to collect GNSS data with high-rate intervals for many surveying projects requiring cm-level accuracy.     

Relative contributions of crust and mantle to generation of Campanian high-K calc-alkaline I-type granitoids in a subduction setting, with special reference to the Harşit Pluton, Eastern Turkey

We present elemental and Sr–Nd–Pb isotopic data for the magmatic suite (~79 Ma) of the Harşit pluton, from the Eastern Pontides (NE Turkey), with the aim of determining its magma source and geodynamic evolution. The pluton comprises granite, granodiorite, tonalite and minor diorite (SiO₂ = 59.43–76.95 wt%), with only minor gabbroic diorite mafic microgranular enclaves in composition (SiO₂ = 54.95–56.32 wt%), and exhibits low Mg# (<46). All samples show a high-K calc-alkaline differentiation trend and I-type features. The chondrite-normalized REE patterns are fractionated [(La/Yb) _n  = 2.40–12.44] and display weak Eu anomalies (Eu/Eu* = 0.30–0.76). The rocks are characterized by enrichment of LILE and depletion of HFSE. The Harşit host rocks have weak concave-upward REE patterns, suggesting that amphibole and garnet played a significant role in their generation during magma segregation. The host rocks and their enclaves are isotopically indistinguishable. Sr–Nd isotopic data for all of the samples display I _Sr = 0.70676–0.70708, ε _Nd(79 Ma) = −4.4 to −3.3, with T _DM = 1.09–1.36 Ga. The lead isotopic ratios are (²⁰⁶Pb/²⁰⁴Pb) = 18.79–18.87, (²⁰⁷Pb/²⁰⁴Pb) = 15.59–15.61 and (²⁰⁸Pb/²⁰⁴Pb) = 38.71–38.83. These geochemical data rule out pure crustal-derived magma genesis in a post-collision extensional stage and suggest mixed-origin magma generation in a subduction setting. The melting that generated these high-K granitoidic rocks may have resulted from the upper Cretaceous subduction of the Izmir–Ankara–Erzincan oceanic slab beneath the Eurasian block in the region. The back-arc extensional events would have caused melting of the enriched subcontinental lithospheric mantle and formed mafic magma. The underplating of the lower crust by mafic magmas would have played a significant role in the generation of high-K magma. Thus, a thermal anomaly induced by underplated basic magma into a hot crust would have caused partial melting in the lower part of the crust. In this scenario, the lithospheric mantle-derived basaltic melt first mixed with granitic magma of crustal origin at depth. Then, the melts, which subsequently underwent a fractional crystallization and crustal assimilation processes, could ascend to shallower crustal levels to generate a variety of rock types ranging from diorite to granite. Sr–Nd isotope modeling shows that the generation of these magmas involved ~65–75% of the lower crustal-derived melt and ~25–35% of subcontinental lithospheric mantle. Further, geochemical data and the Ar–Ar plateau age on hornblende, combined with regional studies, imply that the Harşit pluton formed in a subduction setting and that the back-arc extensional period started by least ~79 Ma in the Eastern Pontides.     

Sources and petrogenesis of Jurassic granitoids in the Yusufeli area,Northeastern Turkey: Implications for pre- and post-collisional lithospheric thinning of the eastern Pontides

Small granitoids emplaced into the early Jurassic volcani-clastic succession in the Yusufeli area, northeastern Turkey, can be temporally and geochemically classified into two groups: early Jurassic low-K and late Jurassic high-K. ⁴⁰Ar–³⁹Ar hornblende analyses yielded 188.0 ± 4.3 Ma for the Dutlup?nar intrusion, dating the subduction related rifting in the region. It comprises metaluminous to weakly peraluminous (ASI = 0.94–1.11) granodiorite and, to a lesser extent, tonalite whose K₂O-poor (< 2.04 wt.%) nature and weak negative Eu anomalies (Eu/Eu? = 0.9–0.7) preclude derivation by fractional crystallization from a K-rich melt. Sr, Nd and Pb isotopic data reveal derivation by partial melting from an already cooled tholeiitic basic rocks which had mantle-like isotope signature. The Sumbated intrusion formed in the late Jurassic (153.0 ± 3.4 Ma) and consists chiefly of metaluminous (ASI = 0.84–0.99) quartz monzodiorite. Medium to high-K₂O, relatively high MgO and Sr contents, flat HREE patterns without prominent Eu anomalies, slightly positive ε_Nd(t) values (+ 1.5 to + 2.5) and low I_Sr ratios (0.7046–0.7056) are consistent with an origin by dehydration melting of a juvenile source, above the garnet stability field, dominated by likely K-amphibole bearing calc-alkaline mafic rocks. Geochemical data show that fractional crystallization from a Sumbated-like quartz monzodioritic magma is the fundamental process responsible for the evolved compositional range of the Keçikaya intrusion. The geochemical and geochronological data presented here indicate that the late Jurassic magmatism occurred in a post-collisional setting. Slab-breakoff, which was followed by shortly after collision, seems to be the most plausible mechanism for the generation of medium to high-K calc-alkaline rocks of the Sumbated and the Keçikaya intrusions, indicating a switch in the geodynamic setting, e.g., from pre-collision to post-collision in the middle Jurassic in the eastern Pontides.