Geomechanical Evaluation of Derinkuyu Antique Underground City and its Implications in Geoengineering

Derinkuyu Underground City, located in the Cappadocia Region of Turkey, is an important structure not only for its antique and archaeological characteristics, but also as a structure in terms of the long-term stability of underground rock structures excavated by mankind. The authors carried out some observational, experimental and theoretical rock mechanics studies in the region from 1996 in the context of a research project supported by the Ministry of Education, Culture, Sports, Science and Technology of Japan for the assessment of the long-term behaviour of Derinkuyu Underground City, and these studies are still continuing. In addition to the monitoring of the environmental conditions such as temperature, moisture and air pressure, they also installed acoustic emission (AE) and electrical potential (EP) measurement systems to monitor the behaviour and response of the surrounding rock at the fifth and seventh floors of the underground city. In this article, the geology, seismicity and state of stress of the Cappadocia Region, climatic conditions in the underground city and its vicinity, short- and long-term behaviours of the surrounding rock, its index and mechanical properties, and effects of water content and freezing–thawing processes were investigated. The stability of Derinkuyu Underground City was also evaluated using theoretical and numerical methods, and the results were presented. Furthermore, its implications in modern geoengineering are also discussed.     

A new moment method for continuum radiative transfer in cosmological re-ionization

We introduce a new code for computing time-dependent continuum radiative transfer and non-equilibrium ionization states in static density fields with periodic boundaries. Our code solves the moments of the radiative transfer equation, closed by an Eddington tensor computed using a long characteristics (LC) method. We show that traditional short characteristics and the optically thin approximation are inappropriate for computing Eddington factors for the problem of cosmological re-ionization. We evolve the non-equilibrium ionization field via an efficient and accurate (errors <1 per cent) technique that switches between fully implicit or explicit finite differencing depending on whether the local time-scales are long or short compared to the time-step. We tailor our code for the problem of cosmological re-ionization. In tests, the code conserves photons, accurately treats cosmological effects and reproduces analytic Strömgren sphere solutions. Its chief weakness is that the computation time for the LC calculation scales relatively poorly compared to other techniques  ( t _LC∝ N ^∼1.5_cells)  ; however, we mitigate this by only recomputing the Eddington tensor when the radiation field changes substantially. Our technique makes almost no physical approximations, so it provides a way to benchmark faster but more approximate techniques. It can readily be extended to evolve multiple frequencies, though we do not do so here. Finally, we note that our method is generally applicable to any problem involving the transfer of continuum radiation through a periodic volume.     

The 150d Modulation of the X-ray Emission of CYG X-1

The X-ray observations of Cyg X-1, covering a period of 3.8 years, by the All Sky Monitor (ASM) on board the Rossi X-ray Timing Explorer (RXTE) were analysed to search for long periodicities. The periodicities in the original data were seen to be distorted by the high/soft state and the failed state counts of the system, and thus the power spectra showed two unrealistic peaks at 1608and 406 days. The ∼ 300 day modulation reported previously(e.g. Priedhorsky et al. 1983) is not confirmed by the present ASM data. Truncation of the high counts above 30 cnts/s partly eliminates the distortion and reveals the true period of 150^d in the power spectra. It is interesting that the high/soft state (HS) and the failed state transitions (FST) always occur on the maximum phase of the150^d cycle. This revised version was published online in July 2006 with corrections to the Cover Date.     

78-Day periodicity in the solar 1–8 Å X-ray index during the ascending branch of cycle 21

In the present paper, the time series of 1–8 Å X-ray index for the time interval between January 1977 to September 1981 has been analyzed by carrying out power spectrum using Discrete Fourier Transform. The 78-day periodicity which is one-half of the period found from various solar parameters of cycles 20 and 21 by various authors, was found. This periodicity is thought to be related to the 152-day periodicity.     

Back-analysis of a seismically induced highway embankment failure during the 1999 Düzce earthquake

An earthquake of magnitude 7.2 on the Richter scale occurred on 12 November 1999 in the Düzce-Bolu region of Turkey. The region was also hit approximately 3 months before during the devastating Kocaeli earthquake of 17 August 1999. Besides high casualties and damage to various engineering structures and buildings, slope and embankment failures on the highway and the country roads occurred. In this study, the authors are concerned with the back-analysis of an embankment failure that occurred on the four lanes E-5 highway connecting Ankara to Gstanbul at Bakacak in the Bolu Province during the 1999 Düzce earthquake. Both pseudo-dynamic and dynamic limiting equilibrium methods are used to back analyse the conditions for the initiation of failure and also displacement responses of the embankment during the earthquake shaking. After having given a brief summary of the investigations on the post-failure state, the geology and geotechnical characteristics of the site, and the dynamic limiting equilibrium method developed and used for analyses are described. The results of the back-analysis based on a pseudo-dynamic approach, revealed that the maximum ground acceleration to initiate the failure of the embankment was probably about 0.125 g. On the other hand, the application of a method for computing the displacement and velocity responses of the failed embankment showed that the failure was possible when the embankment was subjected to ground accelerations similar to that recorded at the Düzce station. In addition, the computations also revealed that the time history of accelerations could be very important for the failure of slopes and embankments.     

Stratigraphy and palaeoenvironment of the Lower-“middle” Oligocene units in the northern part of the Western Taurides (İncesu area,Isparta, Turkey)

This study describes the stratigraphic and palaeoenvironmental significance of the Lower-“middle” Oligocene sediments based on the fauna from the Delikarkas? Formation and the microflora from the ?ncesu Formation of the ?ncesu area (northern part of the western Taurides, Isparta province, Turkey). In the area, the Oligocene sediments show a regressive succession, which begins with the limestones of the Delikarkas? Formation indicating marine conditions followed by conglomerates, sandstones and coaly mudstones of alluvial and fluvial (shallow marine-continental) origin. A well preserved foraminiferal assemblage including Nummulites intermedius, Nummulites vascus and Halkyardia maxima proves an Early Oligocene age for the Delikarkas? Formation. Due to palynological markers such as Boehlensipollis hohli, Slowakipollis hippophaëoides, Dicolpopollis kockelii, Magnolipollis neogenicus ssp. minor, Plicapollis pseudoexcelsus, Caryapollenites simplex and Intratriporopollenites instructus the ?ncesu Formation, which concordantly rests on the Delikarkas? Formation, may be assigned to the Early-“middle” Oligocene. From the palynomorph assemblage, three zones have been recognised according to abundance of species. Zone 1 is characterized by predominance of C. simplex and Momipites punctatus and rarely presence of tricolpate and tricolporate pollen. Zone 2 consists mainly of Inaperturopollenites dubius, Leiotriletes maxoides ssp. maximus, Verrucatosporites favus, Verrucatosporites alienus and infrequently marine dinoflagellate cysts. Zone 3 is characterized by a high percentage of ferns such as Echinatisporis?chattensis and Polypodiaceoisporites saxonicus. The presence of marine dinoflagellate cysts like Apectodinium sp. and Cleistosphaeridium sp., back-mangrove elements such as Acrostichum aureum and lepidocaryoid palms (e.g. Longapertites discordis, Longapertites punctatus and Longapertites psilatus) in the sediments of the ?ncesu Formation imply coastal or near-coastal conditions. Terrestrial palynomorphs in more inland settings were transported by running water towards the sea. Conifers are represented by poorly preserved and rare pollen grains of Pityosporites, Cathayapollis and Piceapollis which may have been transported by wind. In this study, the terrestrial climate of the ?ncesu Formation is also discussed on the basis of the Coexistence Approach method. The climate was warm at the coast (over 20 °C), as evidenced by A. aureum and lepidocaryoid palms, whereas there was a mean annual temperature of 17.2–17.4 °C must be assumed for the upland environment(s).     

Comparison of 2D versus 3D modeling approaches for the analysis of the concrete faced rock‐fill Cokal Dam

This paper's primary purpose is to compare the 2D and 3D analysis methodologies in investigating the performance of a concrete faced rock‐fill dams under dynamic loading conditions. The state of stress on the face plate was obtained in both cases using a total strain based crack model to predict the spreading of cracks on the plate and the corresponding crack widths. Results of the 2D and 3D analyses agree well. Although significantly more demanding, 3D analyses have the advantage of predicting the following: (i) the opening of the vertical construction joints; (ii) the cracking at the valley sides; and (iii) the crushing of the plate during the seismic action. During the earthquake loading, the cracking predicted at the base of the face plate after the impounding spread significantly towards the crest of the dam; however, the crack widths are obtained relatively small. Crushing of the face plates at the construction joints is the primary cause of the large scale cracking predicted on the face plate for high intensity earthquakes. Earthquake induced permanent deformation of the fill increases the compressive stresses on the face plate, thereby reducing the minute cracking on the plate. However, this effect also leads to significant increases in the residual openings at the construction joints and at the plinth level. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of Aftershocks on Earthquake Hazard Estimation: An Example from the North Anatolian Fault Zone

In order to investigate the effect of aftershocks on earthquake hazard estimation, earthquake hazard parameters (_m, b and M_max) have been estimated by the maximum likelihood method from the main shocks catalogue and the raw earthquakes catalogue for the North Anatolian Fault Zone (NAFZ). The main shocks catalogue has been compiled from the raw earthquake catalogue by eliminating the aftershocks using the window method. The raw earthquake catalogue consisted of instrumentally detected earthquakes between 1900 and 1992, and historical earthquakes that occurred between 1000–1900. For the events of the mainshock catalogue the Poisson process is valid and for the raw earthquake catalogue it does not fit. The paper demonstrates differences in the hazard outputs if on one hand the main catalogues and on the other hand the raw catalogue is used. The maximum likelihood method which allows the use of the mixed earthquake catalogue containing incomplete (historical) and complete (instrumental) earthquake data is used to determine the earthquake hazard parameters. The maximum regional magnitude (M_max, the seismic activity rate (_m), the mean return period (R) and the b value of the magnitude-frequency relation have been estimated for the 24°–31° E, 31°–41° E, 41°–45° E sections of the North Anatolian Fault Zone from the raw earthquake catalogue and the main shocks catalogue. Our results indicate that inclusion of aftershocks changes the b value and the seismic activity rate _m depending on the proportion of aftershocks in a region while it does not significantly effect the value of the maximum regional magnitude since it is related to the maximum observed magnitude. These changes in the earthquake hazard parameters caused the return periods to be over- and underestimated for smaller and larger events, respectively.