Inheritance of early Archaean Pb-isotope variability from long-lived Hadean protocrust

Comparison of initial Pb-isotope signatures of several early Archaean (3.65-3.82 Ga) lithologies (orthogneisses and metasediments) and minerals (feldspar and galena) documents the existence of substantial isotopic heterogeneity in the early Archaean, particularly in the Pb-207/Pb-204 ratio. The magnitude of isotopic variability at 3.82-3.65 Ga requires source separation between 4.3 and 4.1 Ga, depending on the extent of U/Pb fractionation possible in the early Earth. The isotopic heterogeneity could reflect the coexistence of enriched and depleted mantle domains or the separation of a terrestrial protocrust with a U-238/Pb-204 (mu) that was ca. 20-30% higher than coeval mantle. We prefer this latter explanation because the high-p signature is most evident in metasediments (that formed at the Earth's surface). This interpretation is strengthened by the fact that no straightforward mantle model can be constructed for these high-mu lithologies without violating bulk silicate Earth constraints. The Pb-isotope evidence for a long-lived protocrust complements similar Hf-isotope data from the Earth's oldest zircons, which also require an origin from an enriched (low Lu/Hf) environment. A model is developed in which greater than or equal to3.8-Ga tonalite and monzodiorite gneiss precursors (for one of which we provide zircon U-Pb data) are not mantle-derived but formed by remelting or differentiation of ancient (ca. 4.3 Ga) basaltic crust which had evolved with a higher U/Pb ratio than coeval mantle in the absence of the subduction process. With the initiation of terrestrial subduction at, we propose, ca. 3.75 Ga, most of the greater than or equal to3.8-Ga basaltic shell (and its differentiation products) was recycled into the mantle, because of the lack of a stabilising mantle lithosphere. We argue that the key event for preservation of all greater than or equal to3.8-Ga terrestrial crust was the intrusion of voluminous granitoids immediately after establishment of global subduction because of complementary creation of a lithospheric keel. Furthermore, we argue that preservation of !3.8-Ga material (in situ rocks and zircons) globally is restricted to cratons with a high U/Pb source character (North Atlantic, Slave, Zimbabwe, Yilgarn, and Wyoming), and that the Pb-isotope systematics of these provinces are ultimately explained by reworking of material that was derived from ca. 4.3 Ga (i.e. Hadean) basaltic crust.

Integrated Pb- and S-isotope investigation of sulphide minerals from the early Archaean of southwest Greenland

We present high-spatial resolution secondary ion mass spectrometry (SIMS) measurements of Pb and S isotopes in sulphides from early Archaean samples at two localities in southwest Greenland. Secondary pyrite from a 3.71 Ga sample of magnetite-quartz banded iron formation in the Isua Greenstone Belt, which has previously yielded unradiogenic Pb consistent with its ancient origin, contains sulphur with a mass independently fractionated (MIF) isotope signature (Delta(33)S =+3.3 parts per thousand). This reflects the secondary mineralization of remobilized sedimentary S carrying a component modified by photochemical reactions in the early Archaean atmosphere. It further represents one of the most extreme positive excursions so far known from the early Archaean rock record. Sulphides from a quartz-pyroxene rock and an ultramafic boudin from the island of Akilia, in the Godth (a) over circle bsfjord, have heterogeneous and generally radiogenic Pb isotopic compositions that we interpret to represent partial re-equilibration of Pb between the sulphides and whole rocks during tectonothermal events at 3.6, 2.7 and 1.6 Ga. Both these samples have Delta(33)S=0 (within analytical error) and therefore show no evidence for MIF sulphur. These data are consistent with previous interpretations that the rock cannot be proven to have a sedimentary origin. Our study illustrates that SIMS S-isotope measurements in ancient rocks can be used to elucidate early atmospheric parameters because of the ability to obtain combined S and Pb-isotope data, but caution must be applied when using such data to infer protolith. When information from geological context, petrography and chronology (i.e. by Pb isotopes) is combined and fully evaluated, Delta(33)S signatures from sulphides and their geological significance can be interpreted with a higher degree of confidence. (c) 2005 Elsevier B.V All rights reserved.

Archaean crustal evolution in eastern Finland : New geochronological and geochemical constraints from the Kuhmo terrain and the Nurmes belt

This thesis summarises the results of four original papers concerning U-Pb geochronology and geochemical evolution of Archaean rocks from the Kuhmo terrain and the Nurmes belt, eastern Finland. The study area belongs to a typical Archaean granite-greenstone terrain, composed of metavolcanic and metasedimentary rocks in generally N-S trending greenstone belts as well as a granitoid-gneiss complex with intervening gneissic and migmatised supracrustal and plutonic rocks. U-Pb data on migmatite mesosomes indicate that the crust surrounding the Tipasjärvi-Kuhmo-Suomussalmi greenstone belt is of varying age. The oldest protolith detected for a migmatite mesosome from the granitoid-gneiss complex is 2.94 Ga, whereas the other dated migmatites protoliths have ages of 2.84 2.79 Ga. The latter protoliths are syngenetic with the majority of volcanic rocks in the adjacent Tipasjärvi-Kuhmo-Suomussalmi greenstone belt. This suggests that the genesis of some of the volcanic rocks within the greenstone belt and surrounding migmatite protoliths could be linked. Metamorphic zircon overgrowths with ages of 2.84 2.81 Ga were also obtained. The non-migmatised plutonic rocks in the Kuhmo terrain and in the Nurmes belt record secular geochemical evolution, typical of Archaean cratons. The studied tonalitic rocks have ages of 2.83 2.75 Ga and they have geochemical characteristics similar to low-Al and high-Al TTD (tonalite-trondhjemite-dacite). The granodiorites, diorites, and gabbros with high Mg/Fe and LILE-enriched characteristics were mostly emplaced between 2.74 2.70 Ga and they exhibit geochemical characteristics typical of Archaean sanukitoid suites. The latest identified plutonic episode took place at 2.70 2.68 Ga, when compositionally heterogeneous leucocratic granitoid rocks, with a variable crustal component, were emplaced. U-Pb data on migmatite leucosomes suggest that leucosome generation may have been coeval with this latest plutonic event. On the basis of available U-Pb and Sm-Nd isotopic data it appears that the plutonic rocks of the Kuhmo terrain and the Nurmes belt do not contain any significant input from Palaeoarchaean sources. A characteristic feature of the Nurmes belt is the presence of migmatised paragneisses, locally preserving primary edimentary structures, with sporadic amphibolite intercalations. U-Pb studies on zircons indicate that the precursors of the Nurmes paragneisses were graywackes that were deposited between 2.71 Ga and 2.69 Ga and that they had a prominent 2.75 2.70 Ga source. Nd isotopic and whole-rock geochemical data for the intercalated amphibolites imply MORB sources. U-Pb data on zircons from the plutonic rocks and paragneisses reveal that metamorphic zircon growth took place at 2.72 2.63 Ga. This was the last tectonothermal event related to cratonisation of the Archaean crust of eastern Finland.

The Rondonian-San Ignacio Province in the SW Amazonian Craton: An overview

The Rondonian-San Ignacio Province (1.56-1.30 Ga) is a composite orogen created through successive accretion of arcs, ocean basin closure and final oblique microcontinent-continent collision. The effects of the collision are well preserved mostly in the Paragua Terrane (Bolivia and Mato Grosso regions) and in the Alto Guapore Belt and the Rio Negro-Juruena Province (Rondonia region), considering that the province was affected by later collision-related deformation and metamorphism during the Sunsas Orogeny (1.25-1.00 Ga). The Rondonian-San Ignacio Province comprises: (1) the Jauru Terrane (1.78-1.42 Ga) that hosts Paleoproterozoic basement (1.78-1.72 Ga), and the Cachoeirinha (1.56-1.52 Ga) and the Santa Helena (1.48-1.42 Ga) accretionary orogens, both developed in an Andean-type magmatic arc; (2) the Paragua Terrane (1.74-1.32 Ga) that hosts pre-San Ignacio units (>1640 Ma: Chiquitania Gneiss Complex, San Ignacio Schist Group and Lomas Manechis Granulitic Complex) and the Pensamiento Granitoid Complex (1.37-1.34 Ga) developed in an Andean-type magmatic arc; (3) the Rio Alegre Terrane (1.51-1.38 Ga) that includes units generated in a mid-ocean ridge and an intra-oceanic magmatic arc environments; and (4) the Alto Guapore Belt (<1.42-1.34 Ga) that hosts units developed in passive marginal basin and intra-oceanic arc settings. The collisional stage (1.34-1.32 Ga) is characterized by deformation, high-grade metamorphism, and partial melting during the metamorphic peak, which affected primarily the Chiquitania Gneiss Complex and Lomas Manechis Granulitic Complex in the Paragua Terrane, and the Colorado Complex and the Nova Mamore Metamorphic Suite in the Alto Guapore Belt. The Paragua Block is here considered as a crustal fragment probably displaced from its Rio Negro-Juruena crustal counterpart between 1.50 and 1.40 Ga. This period is characterized by extensive A-type and intra-plate granite magmatism represented by the Rio Crespo Intrusive Suite (ca. 1.50 Ga), Santo Antonio Intrusive Suite (1.40-1.36 Ga), and the Teotonio Intrusive Suite (1.38 Ga). Magmatism of these types also occur at the end of the Rondonian-San Ignacio Orogeny, and are represented by the Alto Candeias Intrusive Suite (1.34-1.36 Ga), and the Sao Lourenco-Caripunas Intrusive Suite (1.31-1.30 Ga). The cratonization of the province occurred between 1.30 and 1.25 Ga. (C) 2009 Elsevier Ltd. All rights reserved.

The Rondonian-San Ignacio Province in the SW Amazonian Craton: An overview

The Rondonian-San Ignacio Province (1.56-1.30 Ga) is a composite orogen created through successive accretion of arcs, ocean basin closure and final oblique microcontinent-continent collision. The effects of the collision are well preserved mostly in the Paragua Terrane (Bolivia and Mato Grosso regions) and in the Alto Guapore Belt and the Rio Negro-Juruena Province (Rondonia region), considering that the province was affected by later collision-related deformation and metamorphism during the Sunsas Orogeny (1.25-1.00 Ga). The Rondonian-San Ignacio Province comprises: (1) the Jauru Terrane (1.78-1.42 Ga) that hosts Paleoproterozoic basement (1.78-1.72 Ga), and the Cachoeirinha (1.56-1.52 Ga) and the Santa Helena (1.48-1.42 Ga) accretionary orogens, both developed in an Andean-type magmatic arc; (2) the Paragua Terrane (1.74-1.32 Ga) that hosts pre-San Ignacio units (>1640 Ma: Chiquitania Gneiss Complex, San Ignacio Schist Group and Lomas Manechis Granulitic Complex) and the Pensamiento Granitoid Complex (1.37-1.34 Ga) developed in an Andean-type magmatic arc; (3) the Rio Alegre Terrane (1.51-1.38 Ga) that includes units generated in a mid-ocean ridge and an intra-oceanic magmatic arc environments; and (4) the Alto Guapore Belt (<1.42-1.34 Ga) that hosts units developed in passive marginal basin and intra-oceanic arc settings. The collisional stage (1.34-1.32 Ga) is characterized by deformation, high-grade metamorphism, and partial melting during the metamorphic peak, which affected primarily the Chiquitania Gneiss Complex and Lomas Manechis Granulitic Complex in the Paragua Terrane, and the Colorado Complex and the Nova Mamore Metamorphic Suite in the Alto Guapore Belt. The Paragua Block is here considered as a crustal fragment probably displaced from its Rio Negro-Juruena crustal counterpart between 1.50 and 1.40 Ga. This period is characterized by extensive A-type and intra-plate granite magmatism represented by the Rio Crespo Intrusive Suite (ca. 1.50 Ga), Santo Antonio Intrusive Suite (1.40-1.36 Ga), and the Teotonio Intrusive Suite (1.38 Ga). Magmatism of these types also occur at the end of the Rondonian-San Ignacio Orogeny, and are represented by the Alto Candeias Intrusive Suite (1.34-1.36 Ga), and the Sao Lourenco-Caripunas Intrusive Suite (1.31-1.30 Ga). The cratonization of the province occurred between 1.30 and 1.25 Ga. (C) 2009 Elsevier Ltd. All rights reserved.

Caracterização estrutural do embasamento do Grupo Araí, na zona externa da faixa Brasília (GO)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Caracterização metamórfica do Greenstone Belt de Guarinos - GO

The study area comprises Guarinos Greenstone Belt which is located in the centralwest state of Goias and has an approximate size of 280 km2. The present work has as main objective to characterize the metamorphism of the greenstone belt and the relationship between the tectonic events and the metamorphism. For the execution of the work we used data collected in the course Geological Mapping of Crystalline Area of the geology course of UNESP, classes 2009 and 2010, the data collected by Alex Choupina Joaquim Andrade Silva in his master's thesis and the information collected in field in the study area. This compilation of data allows approximately 200 thin sections have been evaluated and allowed the construction of the metamorphic map in scale 1:125.000, also allowed the preparation of the geological map in scale 1:125.000 greestone belt, based on 1:25.000 scale maps performed by students of UNESP and the mapping conducted by Silva (2011). The study area was divided into two tectonic compartments depending on the Shear Zone Engenho Velho that longitudinally crosses study area, with direction approximately NW-SE. The west compartment is bordered by granite-gneiss complex and the east compartment, seems thicker to the south and narrows gradually towards the north and its structure is seen in general as a great monoclinal with NW-SE direction with trim to SW, this way your stack was considered as a package of rocks at the base and metaultramafic metamafic the top, with a second package in tectonic contact, which is composed of detrital and chemical metasedimentary rocks. The east compartment is represented by the intrusion of a large granitic body surrounded by amphibolites and metasediments, is bounded on the east by a granite-gneiss complexes through thrust fault, and is bounded on the west by the other compartment. The structure of the compartment is seen as a major axis oriented roughly... (Completo abstract click electronic access below)

Análise estrutural e microtectônica da inflexão do jordino e depósitos de ouro associados no greenstone belt de Pilar de Goiás

The study area is included in the geological context of the Goias Median Massif, a region where there are associations of Archean granite-gneiss complex (Block Moquém) and a Paleoproterozoic metavolcano-sedimentary sequence ( Pilar greenstone belt ). At the south of area, the greenstone sequence is partially overlain by Neopreoterozoic metasediments of the Araxá Group. The lithostratigraphic units of the Pilar greenstone belt define a shift from about N30W direction (north of the deposit) to N60 and 70W in the region south of the Jordino deposit, where are truncated by the Araxá Group rocks. Mineralogical associations described in this paper allow to indicate that the regional metamorphism that affected the rocks of the greenstone belt and Araxá Group, in the mapped area, reached the upper greenschist facies (garnet zone). Data obtained during mapping and by microtectonics analysis allow to indicate the existence of at least four deformational events that acted on the rocks of the Guarinos greenstone belt and Araxá Group, represented by the phases called Dn-1, Dn, Dn + 1 and Dn+2. It was observed that the pattern of sulphide porfitoblasts in mineralized levels is similar to garnet, biotite and muscovite porfiroblasts (tardi to post Dn) that marks the metamorphic peak of the area

A refined solution to Earth's hidden niobium: implications for evolution of continental crust and mode of core formation

New high-precision niobium (Nb) and tantalum (Ta) concentration data are presented for early Archaean metabasalts, metabasaltic komatiites and their erosion products (mafic metapelites) from SW Greenland and the Acasta gneiss complex, Canada. Individual datasets consistently show sub-chondritic Nb/Ta ratios averaging 15.1+/-11.6. This finding is discussed with regard to two competing models for the solution of the Nb-deficit that characterises the accessible Earth. Firstly, we test whether Nb could have sequestered into the core due to its slightly siderophile (or chalcophile) character under very reducing conditions, as recently proposed from experimental evidence. We demonstrate that troilite inclusions of the Canyon Diablo iron meteorite have Nb and V concentrations in excess of typical chondrites but that the metal phase of the Grant, Toluca and Canyon Diablo iron meteorites do not have significant concentrations of these lithophile elements. We find that if the entire accessible Earth Nb-deficit were explained by Nb in the core, only ca. 17% of the mantle could be depleted and that by 3.7 Ga, continental crust would have already achieved ca. 50% of its present mass. Nb/Ta systematics of late Archaean metabasalts compiled from the literature would further require that by 2.5 Ga, 90% of the present mass of continental crust was already in existence. As an alternative to this explanation, we propose that the average Nb/Ta ratio (15.1+/-11.6) of Earth's oldest mafic rocks is a valid approximation for bulk silicate Earth. This would require that ca. 13% of the terrestrial Nb resided in the Ta-free core. Since the partitioning of Nb between silicate and metal melts depends largely on oxygen fugacity and pressure, this finding could mean that metal/silicate segregation did not occur at the base of a deep magma ocean or that the early mantle was slightly less reducing than generally assumed. A bulk silicate Earth Nb/Ta ratio of 15.1 allows for depletion of up to 40% of the total mantle. This could indicate that in addition to the upper mantle, a portion of the lower mantle is depleted also, or if only the upper mantle were depleted, an additional hidden high Nb/Ta reservoir must exist. Comparison of Nb/Ta systematics between early and late Archaean metabasalts supports the latter idea and indicates deeply subducted high Nb/Ta eclogite slabs could reside in the mantle transition zone or the lower mantle. Accumulation of such slabs appears to have commenced between 2.5 and 2.0 Ga. Regardless of these complexities of terrestrial Nb/Ta systematics, it is shown that the depleted mantle Nb/Th ratio is a very robust proxy for the amount of extracted continental crust, because the temporal evolution of this ratio is dominated by Th-loss to the continents and not Nb-retention in the mantle. We present a new parameterisation of the continental crust volume versus age curve that specifically explores the possibility of lithophile element loss to the core and storage of eclogite slabs in the transition zone. (C) 2003 Elsevier Science B.V. All rights reserved.

Multiple garnet growth in garnet-kyanite-staurolite gneiss, Pangong metamorphic complex, Ladakh Himalaya: New constraints on tectonic setting

Garnet-kyanite-staurolite gneiss in the Pangong complex, Ladakh Himalaya, contains porphyroblastic euhedral garnets, blades of kyanite and resorbed staurolite surrounded by a fine-grained muscovite-biotite matrix associated with a leucogranite layer. Sillimanite is absent. The gneiss contains two generations of garnet in cores and rims that represent two stages of metamorphism. Garnet cores are extremely rich in Mn (X(Sps) = 0.35-038) and poor in Fe (X(Alm) = 0.40-0.45), whereas rims are relatively Mn-poor (X(Sps) =0.07-0.08), and rich in Fe (X(Alm), = 0.75-0.77). We suggest that garnet cores formed during prograde metamorphism in a subduction zone followed by abrupt exhumation, during early collision of the Ladakh arc and Karakoram block. The subsequent India-Asia continental collision subducted the metamorphic rocks to a mid-crustal level, where the garnet rims overgrew the Mn-rich cores at ca. 680 degrees C and ca. 8.5 kbar. PT calculations were estimated from phase diagrams calculated using a calculated bulk chemical composition in the Mn-NCKFMASHT system for the garnet-kyanite-staurolite-bearing assemblage. Muscovites from the metamorphic rocks and associated leucogranites have consistent K-Ar ages (ca. 10 Ma), closely related to activation of the Karakoram fault in the Pangong metamorphic complex. These ages indicate the contemporaneity of the exhumation of the metamorphic rocks and the cooling of the leucogranites. (C) 2011 Elsevier B.V. All rights reserved.

Micas in experimentally shocked gneiss

Biotites and muscovites from a gneiss have been experimentally shocked between 18 and 70 GPa using powder-propellant guns at NASA Johnson Space Center and at the California Institute of Technology. This study shows that shock in biotite and muscovite can produce homogeneous and devolatilized glasses within microseconds. Shock-deformed micas display fracturing, kinking, and complex extinction patterns over the entire pressure range investigated. However, these deformation features are not a sensitive pressure indicator. Localized melting of micas begins at 33 GPa and goes to completion at 70 GPa. Melted biotite and muscovite are optically opaque, but show extensive microvesiculation and flow when observed with the SEM. Electron diffraction confirms that biotite and muscovite have transformed to a glass. The distribution of vesicles in shock-vitrified mica shows escape of volatiles within the short duration of the shock experiment. Experimentally shocked biotite and muscovite undergo congruent melting. Compositions of the glasses are similar to the unshocked micas except for volatiles (H2O loss and K loss). These unusual glasses derived from mica may be quenched by rapid cooling conditions during the shock experiment. Based on these results, the extremely low H2O content of tektites may be reconciled with a terrestrial origin by impact. Release of volatiles in shock-melted micas affects the melting behavior of coexisting dry silicates during the short duration of the shock experiment. Transportation and escape of volatiles released from shock-melted micas may provide plausible mechanisms for the origin of protoatmospheres on terrestrial planets, hydrothermal activity on phyllosilicate-rich meteorite parent bodies, and fluid entrapment in meteorites.

K-Ar ages of biotite and muscovite from Pangong metamorphic complex, Shyok suture zone, India: Implications for the youngest post-collision metamorphic event in Ladakh Himalaya

The garnet-kyanite-staurolite and garnet-biotite-staurolite gneisses were collected from a locality within Lukung area that belongs to the Pangong metamorphic complex in Shyok valley, Ladakh Himalaya. The kyanite-free samples have garnet and staurolite in equilibrium, where garnets show euhedral texture and have flat compositional profile. On the other hand, the kyanite-bearing sample shows equilibrium assemblage of garnet-kyanite-staurolite along with muscovite and biotite. In this case, garnet has an inclusion rich core with a distinct grain boundary, which was later overgrown by inclusion free euhedral garnet. Garnet cores are rich in Mn and Ca, while the rims are poor in Mn and rich in Fe and Mg, suggesting two distinct generations of growth. However, the compositional profiles and textural signature of garnets suggests the same stage of P -T evolution for the formation of the inclusion free euhedral garnets in the kyanite-free gneisses and the inclusion free euhedral garnet rims in the kyanite-bearing gneiss. Muscovites from the four samples have consistent K-Ar ages, suggesting the cooling age (∼ 10 Ma) of the gneisses. These ages make a constraint on the timing of the youngest post-collision metamorphic event that may be closely related to an activation of the Karakoram fault in Pangong metamorphic complex.

The Palu Metamorphic Complex, NW Sulawesi, Indonesia: Origin and evolution of a young metamorphic terrane with links to Gondwana and Sundaland

The Palu Metamorphic Complex (PMC) is exposed in a late Cenozoic orogenic belt in NW Sulawesi, Indonesia. It is a composite terrane comprising a gneiss unit of Gondwana origin, a schist unit composed of meta-sediments deposited along the SE Sundaland margin in the Late Cretaceous and Early Tertiary, and one or more slivers of amphibolite with oceanic crust characteristics. The gneiss unit forms part of the West Sulawesi block underlying the northern and central sections of the Western Sulawesi Province. The presence of Late Triassic granitoids and recycled Proterozoic zircons in this unit combined with its isotopic signature suggests that the West Sulawesi block has its origin in the New Guinea margin from which it rifted in the late Mesozoic. It docked with Sundaland sometime during the Late Cretaceous. U–Th–Pb dating results for monazite suggest that another continental fragment may have collided with the Sundaland margin in the earliest Miocene. High-pressure (HP) and ultrahigh-pressure (UHP) rocks (granulite, peridotite, eclogite) are found as tectonic slices within the PMC, mostly along the Palu–Koro Fault Zone, a major strike-slip fault that cuts the complex. Mineralogical and textural features suggest that some of these rocks resided at depths of 60–120 km during a part of their histories. Thermochronological data (U–Th–Pb zircon and 40Ar/39Ar) from the metamorphic rocks indicate a latest Miocene to mid-Pliocene metamorphic event, which was accompanied by widespread granitoid magmatism and took place in an extensional tectonic setting. It caused recrystallization of, and new overgrowths on, pre-existing zircon crystals, and produced andalusite–cordierite–sillimanite–staurolite assemblages in pelitic protoliths, indicating HT–LP (Buchan-type) metamorphism. The PMC was exhumed as a core complex at moderate rates (c. 0.7–1.0 mm/yr) accompanied by rapid cooling in the Plio-Pleistocene. Some of the UHP rocks were transported to the surface at significantly higher rates (⩾16 mm/yr). The results of our study do not support recent plate tectonic reconstructions that propose a NW Australia margin origin for the West Sulawesi block (e.g. Hall et al., 2009).

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers. (C) 2013 Elsevier B.V. All rights reserved.