Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.

Evidence of episodic brittle faulting in the cratonic part of the Peninsular India and its implications for seismic hazard in slow deforming regions

The region around Waclakkancheri, in the province of Kerala, India, which lies in the vicinity of Palghat-Cauvery ;hear zone (within the Precambrian crystalline terrain), has been a site of microseismic activity since 1989. Earlier studies had identified a prominent WNW-ESE structure overprinting on the E-W trending lineaments associated with Palghat-Cauvery shear zone. We have mapped this structure, located in a chamockite quarry near Desamangalam, Waclakkancheri, which we identify as a ca. 30 km-long south dipping reverse fault. This article presents the characteristics of this fault zone exposed on the exhumed crystalline basement and discusses its significance in understanding the earthquake potential of the region. This brittle deformation zone consists of fracture sets with small-scale displacement and slip planes with embedded fault gouges. The macroscopic as well as the microscopic studies of this fault zone indicate that it evolved through different episodes of faulting in the presence of fluids. The distinct zones within consolidated gouge and the cross cutting relationship of fractures indicate episodic fault activity. At least four faulting episodes can be recognized based on the sequential development of different structural elements in the fault rocks. The repeated ruptures are evident along this shear zone and the cyclic behavior of this fault consists of co-seismic ruptures alternating with inter-seismic periods, which is characterized by the sealed fractures and consolidated gouge. The fault zone shows a minimum accumulated dip/oblique slip of 2.1 m in the reverse direction with a possible characteristic slip of 52 cm (for each event). The ESR dating of fault gouge indicates that the deformation zone records a major event in the Middle Quaternary. The empirical relationships between fault length and slip show that this fault may generate events M >= 6. The above factors suggest that this fault may be characterized as potentially active. Our study offers some new pointers that can be used in other slow deforming cratonic hinterlands in exploring the discrete active faults.

Structural and clumped-isotope constraints on the mechanisms of displacement along low-angle detachments

Despite years of research on low-angle detachments, much about them remains enigmatic. This thesis addresses some of the uncertainty regarding two particular detachments, the Mormon Peak detachment in Nevada and the Heart Mountain detachment in Wyoming and Montana.

Constraints on the geometry and kinematics of emplacement of the Mormon Peak detachment are provided by detailed geologic mapping of the Meadow Valley Mountains, along with an analysis of structural data within the allochthon in the Mormon Mountains. Identifiable structures well suited to constrain the kinematics of the detachment include a newly mapped, Sevier-age monoclinal flexure in the hanging wall of the detachment. This flexure, including the syncline at its base and the anticline at its top, can be readily matched to the base and top of the frontal Sevier thrust ramp, which is exposed in the footwall of the detachment to the east in the Mormon Mountains and Tule Springs Hills. The ~12 km of offset of these structural markers precludes the radial sliding hypothesis for emplacement of the allochthon.

The role of fluids in the slip along faults is a widely investigated topic, but the use of carbonate clumped-isotope thermometry to investigate these fluids is new. Faults rocks from within ~1 m of the Mormon Peak detachment, including veins, breccias, gouges, and host rocks, were analyzed for carbon, oxygen, and clumped-isotope measurements. The data indicate that much of the carbonate breccia and gouge material along the detachment is comminuted host rock, as expected. Measurements in vein material indicate that the fluid system is dominated by meteoric water, whose temperature indicates circulation to substantial depths (c. 4 km) in the upper crust near the fault zone.

Slip along the subhorizontal Heart Mountain detachment is particularly enigmatic, and many different mechanisms for failure have been proposed, predominantly involving catastrophic failure. Textural evidence of multiple slip events is abundant, and include multiple brecciation events and cross-cutting clastic dikes. Footwall deformation is observed in numerous exposures of the detachment. Stylolitic surfaces and alteration textures within and around “banded grains” previously interpreted to be an indicator of high-temperature fluidization along the fault suggest their formation instead via low-temperature dissolution and alteration processes. There is abundant textural evidence of the significant role of fluids along the detachment via pressure solution. The process of pressure solution creep may be responsible for enabling multiple slip events on the low-angle detachment, via a local rotation of the stress field.

Clumped-isotope thermometry of fault rocks associated with the Heart Mountain detachment indicates that despite its location on the flanks of a volcano that was active during slip, the majority of carbonate along the Heart Mountain detachment does not record significant heating above ambient temperatures (c. 40-70°C). Instead, cold meteoric fluids infiltrated the detachment breccia, and carbonate precipitated under ambient temperatures controlled by structural depth. Locally, fault gouge does preserve hot temperatures (>200°C), as is observed in both the Mormon Peak detachment and Heart Mountain detachment areas. Samples with very hot temperatures attributable to frictional shear heating are present but rare. They appear to be best preserved in hanging wall structures related to the detachment, rather than along the main detachment.

Evidence is presented for the prevalence of relatively cold, meteoric fluids along both shallow crustal detachments studied, and for protracted histories of slip along both detachments. Frictional heating is evident from both areas, but is a minor component of the preserved fault rock record. Pressure solution is evident, and might play a role in initiating slip on the Heart Mountain fault, and possibly other low-angle detachments.

Une tradition technologique régionale de l’industrie de pierre polie dans la vallée de l’Outaouais au cours de l’Archaïque supérieur

La séquence culturelle régionale de l’Archaïque supérieur et terminal de la vallée de l’Outaouais n’est pas complétée, mais les variations stylistiques, technologiques et fonctionnelles des artéfacts sur le plan régional et local facilitent le découpage culturel. La reconstruction de la chaîne opératoire des haches, des herminettes et des gouges en pierre polie des sites Muldoon et Lamoureux, deux sites du sud-est ontarien, pourrait permettre de déceler une tradition technologique régionale à cheval entre l’Archaïque laurentien, l’Archaïque post-laurentien et d’autres courants des Grands Lacs. L’analyse des haches, des herminettes et des gouges des sites Muldoon et Lamoureux démontre l’utilisation massive de l’amphibolite. La moitié distale de ces outils est surtout polie et finement abrasée et la moitié proximale souvent éclatée, mais parfois abrasée. Ces éléments technologiques et l’industrie de pierre taillée démontrent des ressemblances et des différences avec les sites laurentiens de la région, et certains sites post-laurentiens du Québec et de l’Ontario méridional. Le matériel en pierre polie des sites Muldoon et Lamoureux démontre des liens technologiques vers l’Ouest tout en gardant un contact avec la sphère d’interaction postlaurentienne. La vallée de l’Outaouais semble alors prendre une place indépendante dans l’Archaïque supérieur, certaines continuités technologiques s’observent entre l’Archaïque laurentien et l’Archaïque post-laurentien.

Tectônica rúptil e sismicidade na área de inundação do açude do Castamhão (CE): implicações para o risco sísmico

This study describes brittle deformation and seismicity in the Castanhão Dam region, Ceará State, Brazil. This reservoir will include a hidroeletric power plant and will store about 6,7 billions m3 of water. Five main litostratigraphic unit were identified in the region: gneissic-migmatitic basement, metavolcanosedimentary sequence, granitoid plutons of Brasiliano age, Mesozoic basaltic dike swarm, and Cenozoic fluvial terraces of the Jaguaribe river. The region has experienced several faulting events that occurred at different crustal levels. Faults formed at depths less than about 12 km present left-lateral movement and are associated with epidote and quartz infillings. Faults formed at depths less than 7 km are mainly strike-slip present cataclastic rocks, fault breccia and gouge. Both fault groups form mainly NE-trendind lineaments and represent reactivation of ductile shear zones or new formed faults that cut across existing structures. Seismically-induced liquefaction fractures take place in Cenozoic terraces and indicate paleoearthquakes that may have reached at leat 6,8 MS. In short, this work indicate that the level of paleoseismicity is much greater than one observed in the instrumental record. Several faults are favourably oriented for reactivation and induced seismicity should be expected after the Castanhão Dam impoudment

Correlação entre contexto morfoestrutural e sismicidade nas regiões de João Câmara e São Rafael (RN)

This MSc thesis describes brittle deformation in two seismic zones located in north-eastern Brazil: João Câmara and São Rafael, Rio Grande do Norte State. Both areas show seismogenic faults, Samambaia and São Rafael, indicated by narrow zones of epicentres with a strike of 040o, a lenght of 30 km and 4 km, and a depth of 1-12 and 0,5-4 km, respectively. The first seismological and geological studies suggested blind faults or faults that were still in the beginning of the nucleation process. The region is under E-W-oriented compression and is underlain by Precambrian terrains, deformed by one or more orogenic cycles, which generated shear zones generally marked by strong pervasive foliation and sigmoidal shapes. The crystalline basement is capped by the Cretaceous Potiguar basin, which is also locally capped by Pliocene continental siliciclastic deposits (Barreiras Formation), and Quaternary alluvium. The main aim of this study was to map epicentral areas and find whether there are any surface geological or morphotectonic expression related to the seismogenic faults. A detailed geological map was carried out in both seismic areas in order to identify brittle structures and fault-related drainage/topographic features. Geological and morphotectonic evidence indicate that both seismogenic faults take place along dormant structures. They either cut Cenozoic rocks or show topographic expression, i.e., are related to topographic heights or depressions and straight river channels. Faults rocks in the Samambaia and São Rafael faults are cataclasite, fault breccia, fault gouge, pseudotachylyte, and quartz veins, which point to reactivation processes in different crustal levels. The age of the first Samambaia and the São Rafael faulting movement possibly ranges from late Precambrian to late Cretaceous. Both fault cut across Precambrian fabric. They also show evidence of brittle processes which took place between 4 and 12 km deep, which probably have not occurred in Cenozoic times. The findings are of great importance for regional seismic hazard. They indicate that fault zones are longer than previously suggested by seismogenic studies. According to the results, the methodology used during this thesis may also be useful in other neotectonic investigation in intraplate areas

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

(Table 2) Trace element composition of seleceted clast and matrix materials from ODP Site 125-786

Shipboard examination of volcanic and sedimentary strata at Site 786 suggested that at least four types of breccias are present: flow-top breccias, associated with cooling and breakup on the upper surface of lava flows; autobreccias, formed by in-situ alteration at the base of flows; fault-gouge breccias; and true sedimentary breccias derived from weathering and erosion of underlying flows. It is virtually impossible to assess the origin of breccia matrix by textural and mineralogical analyses alone. However, it is fundamental for our understanding of breccia provenance to determine the source component of the matrix material. Whether the matrix is uniquely clastderived can be determined by geochemical fingerprinting. Trace elements that are immobile during weathering and alteration do not change their relative abundances. A contribution to the matrix from any source with an immobile trace element signature different from that of the clasts would appear as a perturbation of the trace element signature of the matrix. Trace element analysis of bulk samples from clasts and matrix material in individual breccia units was undertaken in a fashion similar to that used by Brimhall and Dietrich (1987, doi:10.1016/0016-7037(87)90070-6) in analyzing soil provenance: (1) to help distinguish between sedimentary and volcanic breccias, (2) to determine the degree of mixing and depth of erosion in sedimentary breccias, and (3) to analyze the local provenance of the individual breccia components (matrix and clasts). The following elements were analyzed by X-ray fluorescence (XRF): Rb, Sr, Ba, U, Zr, Cu, Zn, Ti, Cr, and V. Of these elements, Zr and Ti probably exhibit truly immobile behavior (Humphris and Thompson, 1978, doi:10.1016/0016-7037(78)90222-3 ). The remaining elements are useful as a reference for the extent of compositional change during the formation of matrix material (Brimhall and Dietrich, 1987, doi:10.1016/0016-7037(87)90070-6).

Documentation, profile descriptions and analytical results of loes profiles from Bavaria, southern Germany

This study investigates the landscape evolution and soil development in the loess area near Regensburg between approximately 6000-2000 yr BP (radiocarbon years), Eastern Bavaria. The focus is on the question how man and climate influenced landscape evolution and what their relative significance was. The theoretical background concerning the factors that controlled prehistoric soil erosion in Middle Europe is summarized with respect to rainfall intensity and distribution, pedogenesis, Pleistocene relief, and prehistoric farming. Colluvial deposits , flood loams, and soils were studied at ten different and representative sites that served as archives of their respective palaeoenvironments. Geomorphological, sedimentological, and pedological methods were applied. According to the findings presented here, there was a high asynchronity of landscape evolution in the investigation area, which was due to prehistoric land-use patterns. Prehistoric land use and settlement caused highly difIerenciated phases of morphodynamic activity and stability in time and space. These are documented at the single catenas ofeach site. In general, Pleistocene relief was substantially lowered. At the same time smaller landforms such as dells and minor asymmetric valleys filled up and strongly transformed. However, there were short phases at many sites, forming short lived linear erosion features ('Runsen'), resulting from exceptional rainfalls. These forms are results of single events without showing regional trends. Generally, the onset of the sedimentation of colluvial deposits took place much earlier (usually 3500 yr BP (radiocarbon) and younger) than the formation of flood loams. Thus, the deposition of flood loams in the Kleine Laaber river valley started mainly as a consequence of iron age farming only at around 2500 yr BP (radiocarbon). A cascade system explains the different ages of colluvial deposits and flood loams: as a result of prehistoric land use, dells and other minor Pleistocene landforms were filled with colluvial sediments. After the filling of these primary sediment traps , eroded material was transported into flood plains, thus forming flood loams. But at the moment we cannot quantify the extent ofprehistoric soil erosion in the investigation area. The three factors that controlled the prehistoric Iandscapc evolution in the Ioess area near Regensburg are as follows: 1. The transformation from a natural to a prehistoric cultural landscape was the most important factor: A landscape with stable relief was changed into a highly morphodynamic one with soil erosion as the dominant process of this change. 2. The sediment traps of the pre-anthropogenic relief determined where the material originated from soil erosion was deposited: either sedimentation took place on the slopes or the filled sediment traps of the slopes rendered flood loam formation possible. Climatic influence of any importance can only be documented as the result of land use in connection with singular and/or statistic events of heavy rainfalls. Without human impact, no significant change in the Holocene landscape would have been possible.