Seismic velocity and Poisson’s ratio tomography of the crust beneath East Anatolia

Eastern Anatolia is a region in the early stages of continent–continent collision and so provides a unique opportunity to study the early development of continental plateau. Located within the Alpine–Himalayan fold-thrust fault belt, the Anatolian plateau is geologically very complex, with over half of the surface area covered with late Cenozoic volcanics of diverse composition. The plateau is also seismically active and is dissected by numerous seismogenic faults predominantly of strike-slip motion. In this study, we determine 3-D tomographic images of the crust beneath eastern Anatolia by inverting a large number of arrival time data of P- and S-waves. From the obtained P- and S-wave velocity models, we estimated the Poisson’s ratio structures for a more reliable interpretation of the obtained velocity anomalies. Our tomographic results are generally consistent with the major tectonic features of the region. High P- and S-wave velocity anomalies are recognized near the surface, while at deeper crustal layers, low seismic wave velocities are widely distributed. Poisson’s ratio exhibits significant structural heterogeneities compared to the imaged velocity structure. The seismic activity is intense along highly heterogeneous zones and is closely associated with pre-existing faults in the central and western parts of the study area. Results of the checkerboard resolution test indicate that the imaged anomalies are reliable features down to a depth of about 40 km. The low-velocity/high Poisson’s ratio zones in the middle to lower crust are consistent with many geophysical observations such as strong Sn attenuation, low Pn and Sn velocity, and the absence of mantle lid, implying the presence of partial melt in the uppermost mantle.     

Crustal velocity and Vp/Vs structures beneath central Anatolia from local seismic tomography

We applied a seismic tomography technique to arrival time data generated by local crustal earthquakes in central Anatolia in order to study the three-dimensional velocity and Vp/Vs structures and their relation with the complex tectonic processes and seismic activity occurring in the study region. The relatively equal and large number of both P- and S-wave arrival times comprising a total of 51,650 arrivals and the relatively uniform distribution of the recording stations imply that the obtained velocity anomalies are reliable features down to a depth of 40 km. This is also evident from the results of the checkerboard resolution test, hit count, and the ray-path coverage. The inversion results indicate the existence of strong lateral heterogeneities in the crust and uppermost mantle beneath central Anatolia. Prominent low-velocity anomalies are clearly imaged at all layers especially beneath existing volcanoes and the active fault segments. Higher-than-average Vp/Vs ratios are widely distributed, indicating the possible existence of over-pressurized fluids that may be responsible for the triggering of the large crustal earthquakes along the north and east Anatolian fault zones. We noticed that the seismic activity occurs mainly at the low-velocity areas and to a lesser extent in some high-velocity zones, perhaps because of the complex tectonics and geological structures. These observations imply that all the zones with velocity anomalies—either low or high—are potential sites for strain energy accumulation and subsequent release. The obtained velocity and Vp/Vs models are consistent with previous geophysical measurements conducted beneath central Anatolia and give much deeper understanding of the current seismotectonic processes occurring in the region.     

First Ground-Based Photometry and Light Curve Analysis of the Recently Discovered Contact Binary Hx Uma

Photoelectric UBV light curves of the recently discovered eclipsing binary HX UMa were obtained and studied to determine the preliminary physical parameters of the system for the first time. The observations were taken at the TüBİTAK¹ – Turkish National Observatory (TUG) on three nights in April 2003. A simultaneous analysis of the light and radial velocity curves yields a typical A-type contact binary with a high degree of overcontact. The influence of the close visual companion to the total light of the system was taken into account during the analysis. ¹TüBİTAK: The Scientific and Technical Research Council of Turkey.     

The Period Variation of V839 Oph

In this study we present the results of the preliminary analysis of the period variation of V839 Oph based on the extensive series of minima times collected from the literature. The character of the (O–C) diagram can be approximated with a cyclic variation superimposed on a quadratic variation. The quadratic variation can be explained in terms of mass exchange/loss mechanism in the system, while the cyclic variation could be attributed to the light-time effect of a gravitationally bound third body to the system or the magnetic activity cycle of the primary component.     

Modelled effect of changes in the CO2 concentration on the middle and upper atmosphere: Sensitivity to gravity wave parameterization

We have studied the temperature response to changes in the CO₂ concentration in the middle and upper atmosphere using the Coupled Middle Atmosphere–Thermosphere Model 2 (CMAT2). We have performed simulations with a range of CO₂ concentrations and three different ways of accounting for the effects of gravity waves, to allow for comparison with previous studies and sensitivity analyses. We initially find that the response of the model to the changes in CO₂ concentration which took place between 1965 and 1995 (320–360 ppm) is strongly dependent on the gravity wave parameterization that is used, but this is to a large degree due to steps or kinks in an otherwise nearly linear curve describing the temperature as a function of CO₂ concentration. We have not been able to identify the cause of these steps as part of the present study, which is a limitation and must be studied in future work. Here we treated the steps as model noise and rather focused on correcting for their effects by fitting straight lines to the temperature–CO₂ curves to estimate the overall slope of the curves. From these slopes we were able to obtain more robust trend estimates than can be obtained by comparing only two model simulations, as is normally done in other, similar studies. The corrected temperature responses to a 40 ppm change in CO₂ concentration still show up to 15–17% sensitivity to the gravity wave parameterization in the mesosphere and thermosphere. This remaining sensitivity is likely to be related to the fundamental differences in the way a change in temperature modifies the propagation and dissipation characteristics of gravity waves in each parameterization, which is particularly different for linear and non-linear schemes. The corrected trends we find are largely in agreement with other modelling studies, and therefore do not fully explain observed trends, which are typically larger than those predicted by modelling studies. However, modelling results could be similarly sensitive to other model parameters and settings, for example to gravity wave characteristics or solar activity level, and this should be further investigated as well.     

Simple risk calculations in dependent hydrological series

Abstract

Various theoretical expressions are presented for determining risk in water resources systems design based on floods modelled by dependent processes. The effect of dependence on the simple risks involved in any engineering design is investigated on the basis of the lag-one Markov process. Asymptotically dependent process formulations reduce to independent case solutions that are already available in the literature. It is shown that a design risk value can be determined if the expected project life, serial correlation coefficient and simple risk are known. Necessary procedures and tables are presented for risk calculations when the variables concerned are serially correlated.     

Spectroscopic and photometric study of the contact binary BO CVn

We present the results of the study of the contact binary system BO CVn. We have obtained physical parameters of the components based on combined analysis of new, multi-color light curves and spectroscopic mass ratio. This is the first time the latter has been determined for this object. We derived the contact configuration for the system with a very high filling factor of about 88%. We were able to reproduce the observed light curve, namely the flat bottom of the secondary minimum, only if a third light has been added into the list of free parameters. The resulting third light contribution is significant, about 20–24%, while the absolute parameters of components are: M₁ = 1.16, M₂ = 0.39, R₁ = 1.62 and R₂ = 1.00 (in solar units).The O-C diagram shows an upward parabola which, under the conservative mass transfer assumption, would correspond to a mass transfer rate of dM/dt = 6.3 × 10^?8M_⊙/yr, matter being transferred from the less massive component to the more massive one. No cyclic, short-period variations have been found in the O-C diagram (but longer-term variations remain a possibility).     

An empirical method for design of grouted bolts in rock tunnels based on the Geological Strength Index (GSI)

The procedure presented in this paper has been developed for the design of grouted rock bolts in rock tunnels during preliminary design stage. The proposed approach provides a step-by-step procedure to set up a series of practical guidelines for optimum pattern of rock bolting in a variety of rock mass qualities. For this purpose, a new formula for the estimation of the rock load (support pressure) is recommended. Due to its wide-spread acceptance in the field of rock engineering, the Geological Strength Index (GSI) is adopted in support pressure equation. For poor and very poor rock mass where the GSI < 27, the use of Modified-GSI is, instead, recommended. The supporting action is assumed to be provided by rock bolts carrying a total load defined by the rock load height. The mechanism of bolting is assumed to rely on roof arch forming and suspension principle. Integrated with support pressure function, the bolt density parameter is modified in order to provide an optimized bolt pattern for any shape of tunnel. The modified bolt density can also be used in analysis of a reinforced tunnel in terms of Ground Reaction Curve (GRC) in such a way as to evaluate the reinforced rock mass and the tunnel convergence. By doing so, the effectiveness of the bolting pattern is well evaluated. The proposed approach based on GSI is believed to overcome constrains and limitations of existing empirical bolt design methods based on RMR or Q-system, which are doubtful in poor rock mass usage. The applicability of the proposed method is illustrated by the stability analysis and bolt design of a rail-road tunnel in Turkey.     

Geochemical characteristics of ophiolitic rocks from the southern margin of the Sivas basin and their implications for the Inner Tauride Ocean,Central-Eastern Turkey

Number of dismembered ophiolite bodies crop out between Sivas and Malatya on the top of the Eastern Tauride platform in the central-eastern Turkey. One of which at the southern margin of the Sivas basin in the Tecer Mountain area comprises melange and the lower part of an oceanic lithospheric section on top of the Tauride platform. The mantle tectonites are characterized by variably serpentinized harzburgites and dunites, and are intruded by numerous isolated dykes. The gabbroic cumulates consist of olivine gabbro, gabbro and gabbronorite. The major and trace element geochemistry of the mafic cumulate rocks suggests that the primary magma was compositionally similar to those observed in modern island-arc tholeiitic sequences. The isolated dykes are exclusively basaltic in composition and display geochemically two distinct subgroups: Group I is represented by high TiO₂ (.87–1.47 wt.%) and other incompatible elements, whereas Group II is characterized by low TiO₂ (.36–.66 wt.%) and other incompatible elements. The Group I isolated diabase dykes have flat to slightly LREE-depleted profiles (La/Yb_N = .32–.79), whereas the Group II isolated diabase dykes are more depleted in general and have a LREE-depleted character (La/Yb_N = .19–.49). This suggests that the isolated dykes were derived from an island arc tholeiitic magma (Nb/Y = .02–.05) with different degrees of partial melting (Group II > Group I) and relatively high oxygen fugacity in intra-oceanic subduction zone. The ophiolitic rocks in the study area may well be compared with the Divri?i ophiolite to the southeast. All the evidence suggests that the isolated dykes in the Tecer Mountain area differ from the alkaline isolated dykes cutting the Divri?i ophiolite. Since the late stage dykes (~76 Ma) in the Divri?i area are alkaline, the tholeiitic isolated dykes in the present study should have been emplaced prior to the alkaline dykes during Late Cretaceous SSZ-spreading (~90 Ma) within the Inner Tauride Ocean.     

Effect of organic solvent type on the removal of nitrate ion using liquid membrane technique

In this study, the removal of nitrate ions from aqueous solutions with liquid membrane technique has been investigated for different organic solvent types in which solubilized tetradecyl trimethyl ammonium bromide (TDTMABr) as carrier. n-butyl alcohol, chloroform, and mixture of chloroform + n-hexane (n-hexane 85% + chloroform 15%) were used as organic solvent. Kinetic parameters (k _1d, k _2m, k _2a, t _max, R _m^max, J _m^max, J _a^max) were calculated from obtained data. time R _a values of mixture, butyl alcohol, and chloroform are 0.81, 0.78, and 0.55, respectively. Similarly R _d, R _m, and t _max values of the mixture equal to 0.14, 0.04, and 87.92 min, respectively. This behavior of the system shows the organic solvent type is an effective parameter on separation yield. It can be concluded that the mixture is the most effective organic solvent type among the investigated ones, because liquid membrane systems should be operated within the range of having the R _m, R _d, and t _max values are minimum while R _a values are maximum.