Coda <Emphasis Type="Italic">Q</Emphasis> in Qeshm Island,south of Iran,using aftershocks of the Qeshm earthquake of November 27, 2005

Coda wave attenuation is estimated for Qeshm Island which is located in the southeastern part of Zagros. For this purpose, the aftershocks of Qeshm earthquake in November 27, 2005, recorded within an epicentral distance less than 100 km, have been used. More than 829 earthquakes were recorded by a local temporary network consisting of 16 short period stations installed after a week after the main shock for ～10 weeks. The coda quality factor, Q _c, was estimated using the single-backscattering model in frequency bands of 0.5–24 Hz. In this research, lateral and vertical variations of coda Q in Qeshm Island are explored. In Qeshm Island, absence of significant lateral variation of coda Q is observed. To investigate the attenuation variation with depth, the coda Q value was calculated for coda time windows with different lengths (5, 10, 15, 20, 25, and 30 s). It is observed that coda Q increases with depth. However, in our study area, the rate of increase of coda Q with depth is not uniform. Beneath Qeshm Island, the rate of increase of coda Q is greater at depths less than ～40 km compared with those of larger depths. This is indicating the existence of a low attenuation anomalous structure under the ～40-km depth which may be correlated with the Moho depth in this region. The average frequency relation for this region is Q _c = 36 ± 1.2f ^{0.94 ± 0.039} at a 5 s-lapse time window length and Q _c = 110 ± 1.8f ^{0.88 ± 0.09} at a 30-s lapse time window length.     

Source parameters of some recent earthquakes in the Gulf of Aqaba, Egypt

Gulf of Aqaba is recognized as an active seismic zone where many destructive earthquakes have occurred. The estimation of source parameters and coda Q attenuation are the main target of this work. Fifty digital seismic events in eight short-period seismic stations with magnitude 2.5–5.2 are used. Most of these events occurred at hypocentral depths in the range of 7–20 km, indicating that the activity was restricted in the upper crust. Seismic moment, M _o, source radius, r, and stress drop, Δσ, are estimated from P- and S-wave spectra using the Brune’s seismic source model. The average seismic moment generated by the whole sequence of events was estimated to be 4.6E?+?22 dyne/cm. The earthquakes with higher stress drop occur at 10-km depth. The scaling relation between the seismic moment and the stress drop indicates a tendency of increasing seismic moment with stress drop. The seismic moment increases with increasing the source radius. Coda waves are sensitive to changes in the subsurface due to the wide scattering effects generating these waves. Single scattering model of local earthquakes is used to the coda Q calculation. The coda with lapse times 10, 20, and 30 s at six central frequencies 1.5, 3, 6, 12, 18, 24 Hz are calculated. The Q _c values are frequency dependent in the range 1–25 Hz, and are approximated by a least squares fit to the power law [ $ {Q_c}(f) = {Q_o}{(f/{f_o})^\eta } $ ]. The average of Q _c values increases from 53?±?10 at 1.5 Hz to 700?±?120 at 24 Hz. The average of Q _o values ranges from 13?±?1 at 1.5 Hz to 39?±?4 at 24 Hz. The frequency exponent parameter η ranges between 1.3?±?0.008 and 0.9?±?0.001.     

Attenuation studies based on local earthquake Coda waves in the southern Indian peninsular shield

The Dharwar craton, Cuddapah basin and the Godavari graben characterise three diverse geological and tectonic settings in the peninsular shield of India. Owing to their contiguous proximity, they offer a unique opportunity to document differences, if any, in their seismic wave attenuation characteristics that might have a bearing on the seismogenic nature of the crust in a craton, basin and a rift-like graben structure. An attempt is made here to bring out these differences using constraints from coda-Q. We considered local earthquakes with epicentral distances ranging from␣14 to 150 km recorded at the digital broadband stations at Dharwar (DHD), Cuddapah (CUD) and Kothagudem (KGD) regions to derive the frequency-dependent coda-Q relations. Using the single scattering method, we obtained the frequency-dependent Q _C relationship (Q _C = Q ₀ f ⁿ )for each of the three geological units separately: DHD: Q _C = (730.62 ± 0.09)f ^{(0.54 ± 0.01)}; CUD: Q _C = (535.06 ± 0.13)f ^{(0.59 ± 0.01)} and KGD: Q _C = (150.56 ± 0.08)f ^{(0.91 ± 0.01)}. The Q _C values obtained for all the three sub-regions show moderate to strong frequency dependence and essentially reflect the level of crustal heterogeneities to varying degrees.     

Rethinking communication in risk interpretation and action

We have studied the attenuation characteristics of eastern Himalaya and southern Tibet by using local earthquake data set that consists of 123 well-located events, recorded by the Himalayan Nepal Tibet Seismic Experiment operated during 2001–2003. We have used single backscattering model to calculate frequency-dependent values of coda Q (\(Q_\mathrm{c}\)). The estimation of \(Q_\mathrm{c}\) is made at central frequencies 2, 4, 8 and 12 Hz through five lapse time windows from 10 to 50 s starting at double the travel time of the S-wave. The observed \(Q_\mathrm{c}\) is found to be strongly frequency-dependent and follows a similar trend as observed in other tectonically active parts of the Himalaya. The trend of variation of \(Q_\mathrm{c}\) with lapse time and the corresponding apparent depths is also studied. Increase in \(Q_\mathrm{c}\) values with the lapse time suggests that the deeper part of the study region is less heterogeneous than the shallower part. The observed values of \(Q_0\) (\(Q_\mathrm{c}\) at 1 Hz) and frequency parameter n indicate that the medium beneath the study area is highly heterogeneous and tectonically very active. A regionalization of the estimated \(Q_0\) is carried out, and a contour map is prepared for the whole region. Some segments of Lesser Himalaya and Sub-Himalaya exhibit very low \(Q_0\) , while the whole Tethyan Himalaya and some parts of Greater Himalaya are characterized by low \(Q_0\) values. Our results are comparable with those obtained from tectonically active regions in the world.     

Attenuation of coda waves in the Nubra-Siachen region,Himalaya, India

The attenuation properties have been estimated in the Nubra-Siachen region situated in the highly mountainous region of Himalayan belt. Coda wave quality factor (Q_c) has been determined for this virgin region by using the single backscattering method. A total of thirty earthquakes recorded in this region, which fall in the epicentral distance range of 3 to 115km have been used for the present work. A 30 sec window length of coda waves at different central frequencies 1.5, 3.0, 6.0, 9.0, 12.0, 18.0 and 24.0 Hz have been studied to determine Q_c at different recording stations. The frequency dependent coda wave quality factor relationships of the form Q_c(f) = Q_ofⁿ, have been computed at each recording stations separately: BASE: Q_c(f) = (137 ± 4.2) f^{(0.99 ± 0.12)}, CHALUNKA: Q_c(f) = (116±3.8)f^(1.0±0.05), PARTA: Q_c(f) = (122±3.0)f^(1.0±0.02), and SASOMA: Q_c(f) = (111±4.1)f^(1.0±0.03). A regional Q_c relation has been developed for the Nubra-Siachen region by using the average value of Q_c at different frequencies obtained at each recording station of the form Q_c(f) = (121±7.2)f^(1.0±0.04). The average Q_c values vary from 183 at 1.5 Hz to 3684 at 24 Hz central frequencies. The present regional relation developed for Nubra-Siachen region indicates heterogeneous and tectonically active region.     

Seismicity and source parameters of local earthquakes in Bilaspur region of Himachal Lesser Himalaya

The pattern of local seismicity (110 events) and the source parameters of 26 local events (1.0?≤?Mw?≤?2.5) that occurred during May 2008 to April 2009 in Bilaspur region of Himachal Lesser Himalaya were determined. The digital records available from one station have been used to compute the source parameters and f _max based on the Brune source model (1970) and a high-frequency diminution factor (Boore 1983) above f _max. The epicentral distribution of events within 30 km of local network is broadly divided into three clusters of seismic activity: (1) a cluster located to the south of the Jamthal (JAMT) station and falls to the north of the Main Boundary Thrust (MBT) which seems to reflect the contemporary local seismicity of the segment of the MBT, (2) an elongated zone of local seismicity NE–SW trending, delineated NE of JAMT station that falls in the Lesser Himalaya between the MBT and the Main Central Thrust, and (3) NE–SW trending zone of local seismic activity located at about 10 km east of NHRI station and about 15 km northeast of NERI station and extending over a distance of about 20 km. Majority of events occur at shallow depths up to 20 km, and the maximum number of events occurs in the focal depth range between 10 and 15 km. The entire seismic activity is confined to the crust between 5 and 45 km. The average values of these source parameters range from 3.29?×?10¹⁷ to 3.73?×?10¹⁹?dyne-cm for seismic moment, 0.1 to 9.7 bars for stress drops, and 111.78 to 558.92 m for source radii. The average value of f _max for these events varies from 7 to 18 Hz and seems to be source dependent.     

Seismicity and seismotectonics of the Gulf of Aqaba region

The historical seismicity of the last ten centuries and the instrumental data that occurred in the Gulf of Aqaba region during the period 1982–2008 are evaluated. It is found that 12 historical earthquakes have occurred with average recurrence periods 70–90 and 333–500 years for M?≥?6.0 and 7.0, respectively. Those with M?≤?6.5 appear to be incomplete and require further investigation. More than 98 % of the instrumental data has occurred in the form of swarms and sequences. The first have released about 32 % of the total energy and are most likely related to subsurface volcanic activities. Their epicentral distribution indicates that all regional faults of the gulf area are active in the present, but with clear concentration within the area bound by latitudes 28.2°–29.8° and longitudes 34.4°–35.2°. Regional strike-slip faults of the northern two basins appear to be as twice active as the normal, or more. An appreciable level of seismic hazard is envisaged as the “a” value is 6.0–6.2 while the “b” value shows a temporal variation, mostly in the range 0.8–1.05. More than 95 % of the seismic energy was released from earthquakes shallower than 22 km. This indicates a brittle upper crust and a ductile lower crust and upper mantle. Tectonic movements at depths?>?22 km appear to be aseismic. The epicentral distribution of the five swarms indicates that the lengths of the causative faults varied in the range 45–70 km. The maximum expected magnitude is Mw?=?6.8–7.2. This implies a seismic slip rate of about 0.54–0.8 Cm/year and some 20–30 % of aseismic tectonic movements. This and the sequence nature of the seismicity of this region result in a noticeable hazard reduction. Combining the seismicity data of the Gulf of Aqaba region with other geophysical, geological, tectonic, and environmental data, clearly indicate that the seismicity of this region is as old as the initiation of the gulf itself. No apparent southward or northward migration of activity is observed.     

Source component mixing controls the variability in Cu and Au endowment along the strike of the Eastern Andean Cordillera in Peru

Mississippian arc magmatic suites of the Au-rich Pataz and Cu-dominated Montañitas regions in Peru reveal distinct modes of magmatic-hydrothermal petro- and metallogenesis. The distinction is remarkable due to their broad contemporaneity (336–322 Ma), arc-parallel position, and close distance (<?50 km) to each other. In both arc regions, petrography, geochemistry, and the tectonic setting of magmatic suites suggest a rapid switch from syn-collisional/compressional to post-collisional/extensional (with ‘A₂-type’ signature) emplacement regime. Rocks of the Au-rich Pataz region originate from mixed sources with a contribution from the mantle (εHf?>?0 and δ¹⁸O of ~?5.3‰) and assimilated old crust (variously low εHf and δ¹⁸O >?5.3‰). The ultimate source of Au in the mineralised Pataz batholith was oxidised (fO₂ at FMQ buffer; based on zircon trace chemistry) and alkali-, LILE- and HFSE-enriched, most likely represented by the metasomatised mantle. The syn-extensional emplacement of the relatively reduced (ΔFMQ?<?0), but unmineralised, A₂-type suite involved assimilation of reduced crust. Associated, reduced, magmatic-hydrothermal fluids infiltrated the Au-bearing batholith suite and effectively mobilised and transported and concentrated Au. In the Montañitas region, rocks are oxidised (ΔFMQ?>?0) and dominantly mantle derived without significant incorporation of crustal material. Samples from the Cu-mineralised suites indicate the additional contribution of a δ¹⁸O <?5.3‰ source, potentially melted layer-2 gabbro. In addition, the elevated whole-rock La/Yb and Sr/Y ratios are compatible with minor addition of slab-derived material, which may have enhanced Cu endowment in this region. Late-magmatic, oxidised fluids derived from the younger A₂-type suite controlled Cu mobilisation and concentration, while Au behaved largely refractory. In general terms, it is postulated that source mixing in continental arcs is a first-order control of contrasting Cu and Au endowment and that sequential intrusion processes facilitate late-magmatic-hydrothermal mobilisation and concentration of specific metal assemblages.     

Shear wave attenuation characteristics over the Central India Tectonic Zone and its surroundings

The study area lies between latitude 18–26°N and longitude 73–83°E, and mainly covers the Central India Tectonic Zone (CITZ). The frequency-dependent shear wave quality factor (Q_s) has been estimated over the CITZ and its surroundings using Double Spectral Ratio (DSR) method. We have considered 25 local earthquakes with magnitude (M_L) varies from 3.0 to 4.7 recorded at 11 stations running under national seismic network. The Fast Fourier Transformed (FFT) spectra were computed from the recorded waveform having time-window from onset of S-phase to 1.0 s and for a frequency-band of 0.1–10 Hz. Three different shear wave velocities (i.e., 3.87, 3.39 and 3.96 km/s) were obtained over the study area based on a pair of earthquakes recorded at a pair of stations. The low Q_s values of 51–96 at 1 Hz (i.e., Q_s = 51f^0.49; Q_s = 90f^0.488 and Q_s = 96f^0.53) were found in the area covering the Son–Narmada–Tapti (SONATA) lineament, CITZ, eastern part of the Satpura fold belt, Vindhyan and Gondwana basins, Godavari and Mahanadi grabens, and southern part of Gangetic plain. Intermediate Q_s values of the order of 204–277 (i.e., Q_s = 204f^0.56 and Q_s = 277f^0.55) were noted in the cartonic areas, namely, Bundelkhand, Dharwar-Bhandara and Bastar. While the higher Q_s values of 391–628 at 1 Hz (i.e., Q_s = 391f^0.49, Q_s = 409f^0.48, Q_s = 417f^0.48, Q_s = 500f^0.66, Q_s = 585f^0.65 and Q_s = 628f^0.69) were found in the eastern part of the SONATA, CITZ, and the northeastern part of the Satpura fold belt. The low Q_s values might be attributing to the more heterogeneous SONATA rift system. Low Q_s values further may presumably be associated with lower-level of seismicity and apparently account for higher tectonic stress accumulation over long duration. The long-term accumulated stress is generally released through occasional triggering of moderate magnitude earthquakes in the SONATA zone. Surrounding the SONATA region, the higher Q_s values possibly accounts for a more homogeneous subsurface structure along the SONATA zone.     

Stress field in the central and northern parts of the Gulf of Suez area,Egypt from earthquake fault plane solutions

Earthquake focal mechanism solutions from 18 events in the central and northern parts of the Gulf of Suez with local magnitudes ranging from 2.8 to 5.2 and occurring between 1983 and 2004 are used to determine the type of motion and stress pattern of the region. Fault plane solutions show mostly normal component; pure normal faulting mechanisms and normal faulting with a strike-slip component. Only some mechanisms show pure strike-slip faulting. The fault planes strike in NW, WNW, NNE and ENE directions, in conformity with the geologically observed striking faults in the northern and central parts of the gulf. The principal stress orientation is also estimated by inverting the selected focal mechanism solutions. The results show that the northern part of the Gulf is subjected to NE–SW to NNE–SSW extension, with a horizontal σ₃ (plunge 3°) and subvertical σ₁ (plunge 80°). This means that the horizontal extensional stresses are still present in the central/northern Gulf of Suez.     

Triassic magmatism in the eastern part of the South China Block: Geochronological and petrogenetic constraints from Indosinian granites

The granitic dykes in the Badu Group,Zhejiang Province,South China provide important insights on tectonic setting and crustal evolution of the South China Block(SCB) and the Indochina Block during Triassic.Here we report LA-ICP-MS U-Pb data of granitic rocks from the Hucun and Kengkou which show early Triassic ages of 242 ± 2 and 232 ± 3 Ma,respectively,representing their timing of emplacement.The dyke rocks are enriched in K,Al,LREE,Rb,Th.U,and Pb.and are depleted in Nb,Ta,Sr,and Ti.The rocks are characterized by highly fractionated REE patterns with(La/Yb)N ratios of 28.46-38.07 with strong negative Eu anomalies(Eu/Eu* = 0.65-0.73).In situ Hf isotopic analyses of zircons from the Hucun granite yielded ε_(Hf)(t) values of-13.9 to-6.4 and two-stage depleted mantle Hf model ages of 1.68-2.15 Ga,which indicate that the magma was formed by partial melting of the Paleoproterozoic metasedimentary protoliths in the Cathaysia Block.The zircons from the Kengkou granite have ε_(Hf)(t) values ranging from 40.7 to 31.5 and yield two-stage depleted mantle Hf model ages of 0.99-2.49 Ga,indicating magma origin from a mixed source.The Hucun and Kengkou dykes,together with the Triassic A-type granites in SE China were probably generated during magmatism associated with crust-mantle decoupling along the convergent plate boundary between SCB and the Indochina Block.     

A quantitative analysis of relative tectonic activity classes of Kali river basin, southwest coast of India

The Kali river basin of the southwest coast of India has been investigated to examine the influence of tectonic activity through an analysis of the geomorphic indices that were computed using geographical information system (GIS). Five geomorphic indices—stream length–gradient index, asymmetry factor, hypsometric integral, valley floor width-to-height ratio, and elongation ratio—were used for identification of evidences of tectonic activity. The results obtained from these indices were combined to develop an index of relative tectonic activity classes using GIS. The average of the five measured geomorphic indices was used to evaluate the distribution of relative tectonic activity in the study area. About 18 % of the study area (～726 km²) belongs to class?1; 37 % (～1506 km²) to class?2; and 44 % (～1,762 km²) to class 3. Sub-basins found at the Western Ghat region, relative tectonic activity classes show low and medium values which indicate the higher degree of tectonic activity compared to the upland plateau regions.     

Seismic attenuation of coda waves in the eastern region of Cuba

Cuba's seismic attenuation had never been studied in detail. In this paper we present the results of the research on the seismic attenuation of Cuba's eastern zone based upon the information collected by the seismological Cuban network from 1998 to 2003. 581 earthquakes were selected from the Cuban catalogue to make this study. All of them, recorded by at least three seismic stations, had their epicenters located in the eastern Cuban region (19.3–22 N, 79–73 W), epicentral distances between 15 km and 213 km, their coda duration magnitudes ranging from 2 and 4.1 and their focal depths reaching up to 30 km. The seismic wave attenuation was studied using coda waves. The single scattering method proposed by Sato in 1977 was applied, the attenuation and frequency dependency for different paths and the correlation of the results with the geotectonics of the region are presented in this paper.The mean Q_c value calculated was Q_c = (64 ± 2)f^{0.84 ± 0.01}. The relatively low Q₀ and the high frequency dependency agree with the values of a region characterized by a high tectonic activity. The Q_c values of seven subregions of eastern Cuba were calculated and correlated with the geology and tectonics of the area.     

Seismic source characteristics in Kachchh and Saurashtra regions of Western India: <Emphasis Type="Italic">b</Emphasis>-value and fractal dimension mapping of aftershock sequences

Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M _w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M _w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D _c) of the hypocenters, are estimated. The b-value and the D _c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D _c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D _c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D _c-values in both the tectonic regions.     

Mineralogy and microchemistry of the Egat gold deposit, South Eastern Desert of Egypt

Gold-bearing quartz lodes from the Egat gold mine, South Eastern Desert of Egypt, are associated with pervasively silicified, highly sheared ophiolitic metagabbro and island-arc metavolcanic rocks. The mineralized quartz veins and related alteration haloes are controlled by NNW-trending shear/fault zones. Microscopic and electron probe microanalyses (EPMA) data of the ore and gangue minerals reveal that fine-grained auriferous sulfarsenides represent early high-temperature (355–382 °C) phases, with formation conditions as (fS₂?=??10, and fO₂ around ?31). A late, low-temperature (302–333 °C) assemblage includes coarse pyrite, arsenopyrite, and free-milling gold grains (88–91 wt.% Au), with formation conditions as (fS₂?=??8 and fO₂ around ?30). Gold was impounded within early sulfarsenides, while free-milling gold blebs occur along microfractures in quartz veins and as inclusions in late sulfides. Infiltration of hydrothermal fluids under brittle–ductile shear conditions led to mobilization of refractory Au from early sulfarsenide phases and reprecipitated free gold, simultaneous with silicification of the host rocks. The positive correlation between Au and As favors and verifies the use of As as the best pathfinder for gold targets, along the NNW-trending shear zones.     

CODA Q estimates for Kumaun Himalaya

CodaQ (Q _c) estimates for the Kumaun Himalaya region have been obtained in high frequency range. Local earthquakes, recorded by a digital seismic network in the region, which fall in the epicentral distances range of 10 to 80 km and with a local magnitude range of 1.4 to 2.8, have been used. The coda waves of 30 sec window length, filtered at seven frequency bands centered at 1.5, 3, 6, 9, 12, 18 and 24Hz, have been analysed using the single backscattering model. The values ofQ _c estimates vary from 65 to 283 at 1.5 Hz to 2119 to 3279 at 24.0 Hz which showed thatQ _c is frequency dependent and its value increases as frequency increases. A frequency-dependentQ _c relationship,Q _c = (92 ± 4.73)f^(1.07±.023), is obtained for the region representing the average attenuation characteristics of seismic waves for Kumaun Himalaya region.     

Application of time- and magnitude-predictable model for long-term earthquake prediction in Iran

The regional time- and magnitude-predictable model has been applied successfully in diverse regions of the world to describe the occurrence of main shocks. In the current study, the model has been calibrated against the historical and instrumental catalog of Iranian earthquakes. The Iranian plateau is divided into 15 seismogenic provinces; then, the interevent times for strong main shocks have been determined for each one. The empirical relations reported by Papazachos et al. (Tectonophysics 271:295–323, 1997a) for the Alpine–Himalayan belt (including Iran) were adopted except for the constant terms that were calculated separately for every seismotectonic area. By using the calibrated equations developed for the study area and taking into account the occurrence time and magnitude of the last main shocks in each seismogenic source, the time-dependent conditional probabilities of occurrence P(?t) of the next main shocks during next 10, 20, 30, 40 and 50 years as well as the magnitude of the expected main shocks (M _f) have been estimated. The immediate probability (within next 10 years) of a large main shock is estimated to be high and moderate (>35 %) in all regions except zones 9 (M _f = 5.8) and 15 (M _f = 6.1). However, it should be noted that the probabilities have been estimated for different M _f values in 15 regions. Comparing the model predictions with the actual earthquake occurrence rates shows the good performance of the model for Iranian plateau.     

Temporal and genetic relationships between the Kidston gold-bearing Breccia Pipe and the Lochaber Ring Dyke Complex, North Queensland, Australia: insights from in situ U–Pb and Hf-isotope analysis of zircon

An existing model for the temporal and genetic relationships between the Kidston gold-bearing Breccia Pipe and the nearby Lochaber Ring Dyke Complex has been evaluated using in situ U–Pb and Hf-isotope analyses of zircon grains. The Oak River Granodiorite, the host rock to the Kidston Breccia Pipe, has 1,551?±?6 Ma old zircon cores overgrown by 417.7?±?2.2 Ma rims. The Black Cap Diorite and Lochaber Granite within the Lochaber Ring Dyke Complex have crystallisation ages of 350.7?±?1.3 and 337.9?±?2.6 Ma respectively. The gold-rich Median Dyke within the Kidston Breccia Pipe has a crystallisation age of 335.7?±?4.2 Ma, and thus is temporally related to the Lochaber Granite. However, zircon grains from the Median Dyke have less radiogenic Hf-isotope compositions (? _Hf from ?7.8 to ?15.8) than those from the Black Cap Diorite ?_Hf?=?0.4 to ?7.2) and the Lochaber Granite (? _Hf?=??1.0 to ?7.5), but within the range defined by zircons from the Oak River Granodiorite ? _Hf?=??8.0 to ?29.2). The Hf-isotope data thus rule out the proposed fractional crystallisation relationship between the Kidston gold-bearing rocks and the Lochaber Ring Dyke Complex. The Kidston Median Dyke may have been produced by mixing between Lochaber Granite magmas and magmas derived by remelting of the Oak River Granodiorite, which was itself derived from Proterozoic crust. There is no evidence for a juvenile component in the Lochaber Ring Dyke Complex or the Median Dyke. The gold enrichment in the Kidston rocks thus may reflect the multi-stage reworking of the Proterozoic crust, which ultimately produced the Carboniferous felsic magmas.     

Petrogenesis and tectonic significance of an early Palaeozoic mafic-intermediate suite of rocks from the Central Tianshan,northwest China

We report petrologic, geochemical, and zircon U?Pb and Hf isotope data from an early Palaeozoic mafic-intermediate suite of rocks in Baluntai, heartland of the Central Tianshan. Analysed major elements of the intermediate rocks show a close affinity to metaluminous I-type granitoids, resembling arc plutons. The mafic rocks display depletions of Nb, Ta, and Ti typical of volcanic arc basalts (VAB). All of the samples fall within the subduction-related field in tectonic discrimination diagrams. Zircons from a granodiorite and a hornblende diorite yield indistinguishable crystallization ages of 475.1 ± 1.7 million years and 473.7 ± 1.6 Ma, indicating an early Palaeozoic magmatic event. Zircons in the mafic rocks yield younger crystallization ages of 427 ± 1 Ma and 426.5 ± 1.4 Ma. Combined with previous published results, our data suggest that the southward subduction of the Palaeo-Tianshan oceanic crust beneath the northern margin of the Tarim block probably began no later than Early Ordovician time and did not end prior to the middle Silurian. Zircons from the granodiorite and hornblende diorite show ?_Hf(t) values of??11.19 to??5.98 and??12.85 to??6.01, respectively, suggesting the reworking of ancient crust. Correspondingly, these zircons yield two T _DM2 age ranges: 2140–1812 Ma and 2241–1812 Ma, probably representing a significant episode of juvenile addition during the assembly of Columbia. In contrast, zircons from the mafic rocks yield ?_Hf(t) values of 3.12–8.91 and 3.19–8.76, corresponding to the T _DM1 ages ranging from 911 to 685 Ma and from 905 to 688 Ma, respectively, suggesting crustal growth in the Central Tianshan microcontinent associated with the initial break-up of Rodinia.     

Relationship between hydrogeological parameters for data-scarce regions: the case of the Araripe sedimentary basin, Brazil

This paper applies and validates a method for generating spatially distributed hydraulic conductivity (k) based on the specific capacity (Q _s) for data-scarce regions. This method has been applied to the Araripe sedimentary basin, Brazil, and consists of four steps: (1) selection of (32) wells for which both k and Q _s data are available; (2) estimation of k as a function of Q _s for the (128) wells for which only specific capacity data are available; (3) spatial distribution of k using the kriging geostatistical tool; (4) validation of the method, using (17) representative wells with k measured data. The equation relating k and Q _s showed a statistically significant linear relationship (R = 0.93), from which a database has been generated using kriging with the spherical model. The results showed a calibration coefficient of Nash and Sutcliffe (NS) of 0.54 and moderate spatial dependence ratio of 69 %. The validation process provided only a moderate efficiency (NS = 0.22), possibly due to the geological complexity of the focus system. Despite its limitations, the method indicates the possibility of application of ordinary kriging to generate reliable data from auxiliary variables, especially for the water management of data-scarce areas.