Comparison of Static,Kinematic and Dynamic Geodetic Deformation Models for Kutlugün Landslide in Northeastern Turkey

Being the cause of the loss of life and damaging property, landslide is an important natural hazard. Therefore, landslides have to be monitored and preventive measures taken accordingly. In Geodesy, landslides can be determined with static, kinematic and dynamic geodetic models. The aim of this study is to develop a dynamic geodetic model for landslides and compare it with static and kinematic geodetic models. A study area was selected in the Northeastern Black Sea Region of Turkey where landslides are the most effective natural hazard. Movements were determined with static, kinematic and dynamic geodetic models using geodetic, geologic and geophysical measurements made in the study area. Groundwater levels changes were regarded as causative forces in the formulation of the dynamic model. The dynamic model delivered more detailed information (direction, values, velocity, acceleration of movements) about landslide movements. It is possible to formulate more realistic strategies about prevention of landslides by using this information. As a result, it can be suggested that dynamic geodetic models are more useful in landslide studies.     

Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.     

Probabilistic seismic response and reliability assessment of isolated bridges

Bridge seismic isolation strategy is based on the reduction of shear forces transmitted from the superstructure to the piers by two means: shifting natural period and earthquake input energy reduction by dissipation concentrated in protection devices. In this paper, a stochastic analysis of a simple isolated bridge model for different bridge and device parameters is conducted to assess the efficiency of this seismic protection strategy. To achieve this aim, a simple nonlinear softening constitutive law is adopted to model a wide range of isolation devices, characterized by only three essential mechanical parameters. As a consequence of the random nature of seismic motion, a probabilistic analysis is carried out and the time modulated Kanai-Tajimi stochastic process is adopted to represent the seismic action. The response covariance in the state space is obtained by solving the Lyapunov equation for a stochastic linearized system. After a sensitivity analysis, the failure probability referred to extreme displacement and the mean value of dissipated energy are assessed by using the introduced stochastic indices of seismic bridge protection efficiency. A parametric analysis for protective devices with different mechanical parameters is developed for a proper selection of parameters of isolation devices under different situations.     

Centrifuge modeling of PGD response of buried pipe

A new centrifuge based method for determining the response of continuous buried pipe to PGD is presented. The physical characteristics of the RPI‘s 100 g-ton geotechnical centrifuge and the current lifeline experiment split-box are described: The split-box contains the model pipeline and surrounding soil and is manufactured such that half can be offset, in flight, simulating PGD. In addition, governing similitude relations which allow one to determine the physical characteristics, (diameter, wall thickness and material modulus of elasticity) of the model pipeline are presented. Finally, recorded strains induced in two buried pipes with prototype diameters of 0.63 m and 0.95 m (24 and 36 inch) subject to 0.6 and 2.0 meters (2 and 6 feet) of full scale fault offsets and presented and compared to corresponding FE results.     

Measurement of six degree-of-freedom ground motion by using eight accelerometers

A new integrated measuring system with eight force-balance accelerometers is proposed to obtain a direct measurement of six degree-of freedom （DOF） ground motions, including three rotational and three actual translational acceleration components without gyroscopes. In the proposed measuring system, the relationship between the output from eight force-balance accelerometer and the six DOF motion of the measuring system under an earthquake are described by differential equations. These equations are derived from the positions and directions of the eight force-balance accelerometers in the measuring system. The third-order Runge-Kutta algorithm is used to guarantee the accuracy of the numerical calculation. All the algorithms used to compute the six DOF components of the ground motion are implemented in a real-time in Digital Signal Processor （DSP）. The distortion of the measured results caused by position and direction errors of the accelerometers in the measuring system are reduced by multiplying a compensation coefficient C to the output and subtracting static zero drift from the measured results, respectively.     

Geochemical characteristics of crustal anatexis during the formation of migmatite at the Southern Sierra Nevada,California

We provide data on the geochemical and isotopic consequences of nonmodal partial melting of a thick Jurassic pelite unit at mid-crustal levels that produced a migmatite complex in conjunction with the intrusion of part of the southern Sierra Nevada batholith at ca. 100 Ma. Field relations suggest that this pelitic migmatite formed and then abruptly solidified prior to substantial mobilization and escape of its melt products. Hence, this area yields insights into potential mid-crustal level contributions of crustal components into Cordilleran-type batholiths. Major and trace-element analyses in addition to field and petrographic data demonstrate that leucosomes are products of partial melting of the pelitic protolith host. Compared with the metapelites, leucosomes have higher Sr and lower Sm concentrations and lower Rb/Sr ratios. The initial ⁸⁷Sr/⁸⁶Sr ratios of leucosomes range from 0.7124 to 0.7247, similar to those of the metapelite protoliths (0.7125–0.7221). However, the leucosomes have a much wider range of initial ε_Nd values, which range from −6.0 to −11.0, as compared to −8.7 to −11.3 for the metapelites. Sr and Nd isotopic compositions of the leucosomes, migmatites, and metapelites suggest disequilibrium partial melting of the metapelite protolith. Based on their Sr, Nd, and other trace-element characteristics, two groups of leucosomes have been identified. Group A leucosomes have relatively high Rb, Pb, Ba, and K₂O contents, Rb/Sr ratios (0.15<Rb/Sr<1.0), and initial ε_Nd values. Group B leucosomes have relatively low Rb, Pb, Ba, and K₂O contents, Rb/Sr ratios (<0.15), and initial ε_Nd values. The low Rb concentrations and Rb/Sr ratios of the group B leucosomes together suggest that partial melting was dominated by water-saturated or H₂O-fluxed melting of quartz + feldspar assemblage with minor involvement of muscovite. Breakdown of quartz and plagioclase with minor contributions from muscovite resulted in low Rb/Sr ratios characterizing both group A and group B leucosomes. In contrast, group A leucosomes have greater contributions from K-feldspar, which is suggested by: (1) their relatively high K concentrations, (2) positive or slightly negative Eu anomalies, and (3) correlation of their Pb and Ba concentrations with K₂O contents. It is also shown that accessory minerals have played a critical role in regulating the partitioning of key trace elements such as Sm, Nd, Nb, and V between melt products and residues during migmatization. The various degrees of parent/daughter fractionations in the Rb–Sr and Sm–Nd isotopic systems as a consequence of nonmodal crustal anatexis would render melt products with distinct isotopic signatures, which could profoundly influence the products of subsequent mixing events. This is not only important for geochemical patterns of intracrustal differentiation, but also a potentially important process in generating crustal-scale as well as individual pluton-scale isotopic heterogeneities.     

Chemical composition of rock-forming minerals in gold-related granitoid intrusions,southwestern New Brunswick,Canada: implications for crystallization conditions,volatile exsolution,and fluorine-chlorine activity

Chemical composition of rock-forming minerals in Appalachian Siluro-Devonian granitoid intrusions, southwestern New Brunswick, was systematically determined by electron microprobe. The mineral chemical data together with petrographic examination was used to test magmatic equilibration and to constrain crystallization conditions, volatile exsolution, and fluorine-chlorine activity of fluids associated with these intrusions. Mineralogical distinction between Late Silurian to Early Devonian granodioritic to monzogranitic series (GMS) and Late Devonian granitic series (GS) rocks is evident, although both are subsolvus I-type to evolved I-type granitoids. Oxidized to reduced GMS rocks consist of quartz, plagioclase (An>10), K-feldspar, biotite, apatite, titanite, zircon, monazite, ± hornblende, ± pyroxene, ± magnetite, ± ilmenite, and ± sulfide. GS rocks comprise quartz, K-feldspar, plagioclase (An<10), mica group minerals, zircon, monazite, apatite, sulfide, ± ilmente, ± magnetite, ± topaz, ± columbite, and ± xenotime. Inter-intrusion and intra-intrusion variations in mineral chemistry are interpreted to reflect petrogenetic processes (e.g., assimilation and fractional crystallization) during granitoid evolution. Although magmatic equilibration among rock-forming minerals are disturbed by subsolidus hydrothermal processes, GMS rocks appear to have higher magmatic temperatures, variable levels of emplacement, a range of (i.e., reduced intrusions 10^−16.7∼10^−13.4 and oxidized intrusions 10^−14.0∼10^−10.5 bars), and relatively low f _HF/f _HCl ratios (10^−3.0∼10^−1.0) in exsolved fluids, compared to GS rocks. Reduced GMS intrusions bear higher gold potential and thus may be prospective targets for intrusion-related gold systems. Electronic Supplementary Material Supplementary material is available for this article at     

Refining the timing of eclogite metamorphism: a geochemical,petrological, Sm-Nd and U-Pb case study from the Pohorje Mountains,Slovenia (Eastern Alps)

High-pressure metamorphism in the Pohorje Mountains of Slovenia (Austroalpine unit, Eastern Alps) affected N-MORB type metabasic and metapelitic lithologies. Thermodynamic calculations and equilibrium phase diagrams of kyanite–phengite-bearing eclogites reveal PT conditions of >2.1 GPa at T<750°C, but within the stability field of quartz. Metapelitic eclogite country rocks contain the assemblage garnet + phengite + kyanite + quartz, for which calculated peak pressure conditions are in good agreement with results obtained from eclogite samples. The eclogites contain a single population of spherical zircon with a low Th/U ratio. Combined constraints on the age of metamorphism come from U/Pb zircon as well as garnet–whole rock and mineral–mineral Sm-Nd analyses from eclogites. A coherent cluster of single zircon analyses yields a ²⁰⁶Pb/²³⁸U age of 90.7±1.0 Ma that is in good agreement with results from Sm-Nd garnet–whole rock regression of 90.7±3.9 and 90.1±2.0 Ma (εNd: +8) for two eclogite samples. The agreement between U-Pb and Sm-Nd age data strongly suggests an age of approximately 90 Ma for the pressure peak of the eclogites in the Pohorje Mountains. The presence of garnet, omphacite and quartz inclusions in unfractured zircon indicates high-pressure rather than ultrahigh pressure conditions. The analysed metapelite sample yields a Sm-Nd garnet–whole rock scatterchron age of 97±15 Ma. These data probably support a single P-T loop for mafic and pelitic lithologies of the Pohorje area and a late Cretaceous high-pressure event that affected the entire easternmost Austroalpine basement including the Koralpe and Saualpe eclogite type locality in the course of the complex collision of the Apulian microplate and Europe.     

A Mössbauer study of pseudotachylytes: redox conditions during seismogenic faulting

The redox conditions during frictional melting provide information on the physical and chemical conditions during seismic slip in the crust. Here we examine frictional melts from five localities by analyzing host rocks and corresponding pseudotachylytes using Mössbauer spectroscopy. The faults examined are located at South Mountain, Arizona; Fort Foster, Maine (two localities); Long Ridge fault, North Carolina; and the Homestake shear zone, Colorado. The main iron-bearing phases in the pseudotachylytes are phyllosilicates (biotite, muscovite and clays) and iron oxides (magnetite and hematite) and minor pyrite. The ferrous/ferric ratios of the phyllosilicates in the host rocks are the same as those in the pseudotachylytes, except for the hematite-bearing pseudotachylyte from the Long Ridge fault, which is more oxidized. The magnetites in the host rocks and the corresponding pseudotachylytes have different ferric and ferrous iron distributions, which is attributed to different cation chemistry, rather than redox conditions. With the exception of the South Mountain locality, the ferric/ferrous ratios of the micas are interpreted to record the primary redox state of the pseudotachylyte melt as the calculated oxygen fugacities are consistent with magnetite and hematite equilibria. Pyrite-bearing pseudotachylytes plot ~0–1 log₁₀ units above the fayalite-magnetite-quartz (FMQ) buffer. Magnetite-bearing pseudotachylytes plot ~2–4 log₁₀ units above the FMQ buffer, and hematite-bearing pseudotachylytes plot 3.5 log₁₀ units above the hematite-magnetite (HM) buffer. Hematite-bearing pseudotachylytes, together with previous oxygen isotope data, are inferred to represent melting in the presence of externally derived pressurized water. Other localities are inferred to represent melting under rock-buffered, closed system, conditions. If the localities studied are representative of seismogenic faulting, the calculated oxygen fugacities indicate that, in the system C–O–H–S, H₂O and CO₂ should be the dominant fluid species. This is the first detailed study of the redox state of pseudotachylytes.     

Contrasting interactions of sodium and potassium with H<Subscript>2</Subscript>O in haplogranitic liquids and glasses at 200 MPa from hydration–diffusion experiments

This study examines hydration–diffusion in the metaluminous haplogranite system at 200 MPa H₂O and 800–300°C. At 800°C hydration is accompanied by melting and uphill diffusion of sodium from anhydrous glass toward the region of hydration and melting, whereas potassium diffuses away from the hydration front and into anhydrous glass. Silicon and aluminum are simply diluted upon hydration. There is no change in molecular Al/(Na + K) throughout the entire hydration-diffusion aureole and, therefore, (1) there is no loss of alkalis to the vapor, and (2) K migrates to replace Na in order to maintain local charge balance required by ^IVAl. Alkali diffusion occurs over a viscosity contrast from 10^4.1 Pa s in hydrous liquid to 10^11.8–10^13.5 Pa s in anhydrous glass. From these results, we interpret that: (1) Na is structurally or energetically favored over K as a charge-balancing cation for ^IVAl in hydrous granitic liquids, whereas the opposite behavior has been observed for anhydrous melts, and (2) the diffusion of alkalis through silicate melts is largely independent of viscosity. Results from 600°C are similar to those at 800°C, but hydration at 300°C involves a loss of Na and concomitant increase in molar Al/(Na + K) in the hydration zone due to hydrogen-alkali exchange between fluid and glass. Hydration behavior at 400°C is transitional between those at 300°C and 600°C, suggesting that the change in hydration mechanism occurs near the glass transition.