Preliminary study on the occurrence of geothermal systems in the tectonic compressional regions: an example from the Derman geothermal field in the Biga Peninsula, Turkey

There are numerous hot springs with temperatures ranging from 30 to 100 °C in Biga peninsula and they occur throughout the peninsula. The result of this study shows that the region is under a tectonic compressional regime. The investigation of the faults and fractures in the region indicates that the region has been affected first by N–S and then E–W compression since the Middle Miocene. Opening fractures and antithetic and synthetic faults due to the compressional movements provide paths for the deep circulation of water. In addition, the tectonic movements, granitic intrusion and volcanic activity have also played important roles as heat sources for the geothermal systems.     

Charnockitic magmatism in southern India

Large charnockite massifs cover a substantial portion of the southern Indian granulite terrain. The older (late Archaean to early Proterozoic) charnockites occur in the northern part and the younger (late Proterozoic) charnockites occur in the southern part of this high-grade terrain. Among these, the older Biligirirangan hill, Shevroy hill and Nilgiri hill massifs are intermediate charnockites, with Pallavaram massif consisting dominantly of felsic charnockites. The charnockite massifs from northern Kerala and Cardamom hill show spatial association of intermediate and felsic charnockites, with the youngest Nagercoil massif consisting of felsic charnockites. Their igneous parentage is evident from a combination of features including field relations, mineralogy, petrography, thermobarometry, as well as distinct chemical features. The southern Indian charnockite massifs show similarity with high-Ba-Sr granitoids, with the tonalitic intermediate charnockites showing similarity with high-Ba-Sr granitoids with low K₂O/Na₂O ratios, and the felsic charnockites showing similarity with high-Ba-Sr granitoids with high K₂O/Na₂O ratios. A two-stage model is suggested for the formation of these charnockites. During the first stage there was a period of basalt underplating, with the ponding of alkaline mafic magmas. Partial melting of this mafic lower crust formed the charnockitic magmas. Here emplacement of basalt with low water content would lead to dehydration melting of the lower crust forming intermediate charnockites. Conversely, emplacement of hydrous basalt would result in melting at higher {ie565-01} favoring production of more siliceous felsic charnockites. This model is correlated with two crustal thickening phases in southern India, one related to the accretion of the older crustal blocks on to the Archaean craton to the north and the other probably related to the collision between crustal fragments of East and West Gondwana in a supercontinent framework.     

Crustal Growth by Magmatic Accretion Constrained by Metamorphic P-T Paths and Thermal Models of the Kohistan Arc, NW Himalayas

Magmatic accretion is potentially an important mechanism inthe growth of the continental crust and the formation of granulites.In this study, the thermal evolution of a magmatic arc in responseto magmatic accretion is modeled using numerical solutions ofthe one-dimensional heat conduction equation. The initial andboundary conditions used in the model are constrained by geologicalobservations made in the Kohistan area, NW Himalayas. Takingconsideration of the preferred intrusion locations for basalticmagmas, we consider two plausible modes of magmatic accretion:the first involves the repeated intrusion of basalt at mid-crustaldepths (‘intraplate model’), and the second evaluatesthe simultaneous intrusion of basalt and picrite at mid-crustaldepths and the base of the crust respectively (‘double-platemodel’). The results of the double-plate model accountfor both the inferred metamorphic P–T paths of the Kohistanmafic granulites and the continental geotherm determined frompeak P–T conditions observed for granulite terranes. Thedouble-plate model may be applicable as a key growth processfor the production of thick mafic lower crust in magmatic arcs. KEY WORDS: thermal model; magmatic underplating; P–T path; granulite; lower crust     

A quantitative, three-dimensional depositional model of gravelly braided rivers

A quantitative, three‐dimensional depositional model of gravelly, braided rivers has been developed based largely on the deposits of the Sagavanirktok River in northern Alaska. These deposits were described using cores, wireline logs, trenches and ground‐penetrating radar profiles. The origin of the deposits was inferred from observations of: (1) channel and bar formation and migration and channel filling, interpreted from aerial photographs; (2) water flow during floods; and (3) the topography and texture of the river bed at low‐flow stage. This depositional model quantitatively represents the geometry of the different scales of strataset, the spatial relationships among them and their sediment texture distribution. Porosity and permeability in the model are related to sediment texture. The geometry of a particular type and scale of strataset is related to the geometry and migration of the bedform type (e.g. ripples, dunes, bedload sheets, bars) associated with deposition of the strataset. In particular, the length‐to‐thickness ratio of stratasets is similar to the wavelength‐to‐height ratio of associated bedforms. Furthermore, the wavelength and height of bedforms such as dunes and bars are related to channel depth and width. Therefore, the thickness of a particular scale of strataset (i.e. medium‐scale cross‐sets and large‐scale sets of inclined strata) will vary with river dimensions. These relationships between the dimensions of stratasets, bedforms and channels mean that this depositional model can be applied to other gravelly fluvial deposits. The depositional model can be used to interpret the origin of ancient gravelly fluvial deposits and to aid in the characterization of gravelly fluvial aquifers and hydrocarbon reservoirs.     

Gravity anomalies and associated tectonic features over the Indian Peninsular Shield and adjoining ocean basins

A combined gravity map over the Indian Peninsular Shield (IPS) and adjoining oceans brings out well the inter-relationships between the older tectonic features of the continent and the adjoining younger oceanic features. The NW–SE, NE–SW and N–S Precambrian trends of the IPS are reflected in the structural trends of the Arabian Sea and the Bay of Bengal suggesting their probable reactivation. The Simple Bouguer anomaly map shows consistent increase in gravity value from the continent to the deep ocean basins, which is attributed to isostatic compensation due to variations in the crustal thickness. A crustal density model computed along a profile across this region suggests a thick crust of 35–40 km under the continent, which reduces to 22/20–24 km under the Bay of Bengal with thick sediments of 8–10 km underlain by crustal layers of density 2720 and 2900/2840 kg/m³. Large crustal thickness and trends of the gravity anomalies may suggest a transitional crust in the Bay of Bengal up to 150–200 km from the east coast. The crustal thickness under the Laxmi ridge and east of it in the Arabian Sea is 20 and 14 km, respectively, with 5–6 km thick Tertiary and Mesozoic sediments separated by a thin layer of Deccan Trap. Crustal layers of densities 2750 and 2950 kg/m³ underlie sediments. The crustal density model in this part of the Arabian Sea (east of Laxmi ridge) and the structural trends similar to the Indian Peninsular Shield suggest a continent–ocean transitional crust (COTC). The COTC may represent down dropped and submerged parts of the Indian crust evolved at the time of break-up along the west coast of India and passage of Reunion hotspot over India during late Cretaceous. The crustal model under this part also shows an underplated lower crust and a low density upper mantle, extending over the continent across the west coast of India, which appears to be related to the Deccan volcanism. The crustal thickness under the western Arabian Sea (west of the Laxmi ridge) reduces to 8–9 km with crustal layers of densities 2650 and 2870 kg/m³ representing an oceanic crust.     

Regional moment tensor uncertainty due to mismodeling of the crust

The retrieval of earthquake moment tensor (MT) requires the response of the medium, in which seismic waves travel from the hypocenter to the stations, to be known. In inverting long-period (LP) seismic data (teleseismic and LP regional records), a gross earth model is sufficient; with decreasing periods, a more detailed model is needed. This is the case when waveforms of weak earthquakes at regional distances are to be inverted. Regional moment tensors (RMTs) of mostly Mediterranean earthquakes are determined on a routine basis by the Swiss Seismological Survey (SED) by using averaged models of the earth's crust. By inverting broad-band records of the M_w=4.8 earthquake near Udine, N Italy, on Feb. 14, 2002, we tested the sensitivity of the MT solution with respect to possible errors in the earth model used and in the location of the hypocenter depth. We perturbed the P and S velocities and the thickness in the 1-D earth model in the range from 3% to 30% of the parameter values and constructed estimates of confidence regions of the MT and error bars of the source time function (STF) and scalar moment in three frequency bands. Similarly, these error characteristics were determined assuming a mislocation in the hypocenter depth. We found that, in the band of periods from 25 to 50 s, the mechanism is resolved well (at the confidence level 95% at least) up to an earth model uncertainty of 30%, in the passband 10–25 s up to about 10%, but it is undetermined completely at periods of 5–10 s. An error in hypocenter depth of as much as double the value reported by the location procedure does not destroy the resolution of the mechanism at periods above 10 s. In the RMT catalog of the SED, earthquakes of M_w greater than about 3.5 are processed at periods above 30 s; thus, the solutions for these events are robust with respect to a possible uncertainty in the earth model used. Mechanisms of weaker earthquakes, retrieved from short periods, should be interpreted with caution.     

Analysis of geomagnetic field along seismic profile P4 of the International Project POLONAISE'97

This paper presents relative secular variations of the total intensity of the geomagnetic field against a background of results of magnetic anomaly interpretation along seismic profile P4. Profile P4 crosses a Variscan folding zone in the Paleozoic Platform (PLZ), the Trans-European Suture Zone (TESZ), and the Polish part of the East European Craton (EEC). Secular geomagnetic field variations measured in 1966–2000 along a line adjacent to seismic profile P4 were analysed. The study of secular variations, reduced to the base recordings at the Belsk Magnetic Observatory, showed that the growth of geomagnetic field at the East European Craton was slower than in the Trans-European Suture Zone and the Paleozoic Platform.A 2D crustal magnetic model was interpreted as a result of magnetic modelling, in which seismic, geological and geothermal data were also used. The modelling showed that there were significant differences in the magnetic model for geotectonic units, which had been earlier determined based on deep seismic survey data. It should be noted that a fundamental change of trend of the relative secular variations was observed at the slope of the Precambrian Platform. After analysing the geomagnetic field observed along profile P4, the hypothesis that the contact between Phanerozoic and Precambrian Europe lies in Poland's territory can be proven.     

Application of Remote Sensing in Flood Management with Special Reference to Monsoon Asia: A Review

The conventional means to record hydrological parameters of aflood often fail to record an extreme event. Remote sensingtechnology along with geographic information system (GIS)has become the key tool for flood monitoring in recent years.Development in this field has evolved from optical to radarremote sensing, which has provided all weather capabilitycompared to the optical sensors for the purpose of flood mapping.The central focus in this field revolves around delineation of floodzones and preparation of flood hazard maps for the vulnerable areas.In this exercise flood depth is considered crucial for flood hazardmapping and a digital elevation model (DEM) is considered to bethe most effective means to estimate flood depth from remotelysensed or hydrological data. In a flat terrain accuracy of floodestimation depends primarily on the resolution of the DEM. Riverflooding in the developing countries of monsoon Asia is very acutebecause of their heavy dependence on agriculture but any floodestimation or hazard mapping attempt in this region is handicappedby poor availability of high resolution DEMs. This paper presents areview of application of remote sensing and GIS in flood managementwith particular focus on the developing countries of Asia.     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.     

Simulation of Orissa Super Cyclone (1999) using PSU/NCAR Mesoscale Model

In this study a non-hydrostatic version of Penn State University (PSU) -- NationalCenter for Atmospheric Research (NCAR) mesoscale model is used to simulate thesuper cyclonic storm that crossed Orissa coast on 29 October 1999. The model isintegrated up to 123 h for producing 5-day forecast of the storm. Several importantfields including sea level pressure, horizontal wind and rainfall are compared with theverification analysis/observation to examine the performance of the model. The modelsimulated track of the cyclone is compared with the best-fit track obtained from IndiaMeteorological Department (IMD) and the track obtained from NCEP/NCAR reanalysis. The model is found to perform reasonably well in simulating the track and in particular, the intensity of the storm.