This research project on volcanic debris avalanches aims to provide a better understanding of the influence of the volume, fall height and topography on the deposit location and morphology. This will enable improvements in delineation of the areas at risk from volcanic debris avalanches, and improvements in management of a disaster should it occur. Undertaken to fulfil the requirements for a double degree (Geological Engineering and MSc in Hazard and Disaster Management) this work is the result of a collaboration between Polytechnic Institute LaSalle-Beauvais in France and the University of Canterbury in New Zealand. Following a brief introduction to the topic, statistical analyses of volcanic debris avalanche deposits are undertaken. Multiple variables analyses (Principal Components Analyses and Regressions) were carried out using a database of 298 volcanic debris avalanches derived from modification of Dufresne’s recent database. It was found that the volume has the main influence on the deposits rather than the fall height; the latter seems to have greater effect on avalanches of small volume. The topography into which the deposit is emplaced mainly determines its geometrical characteristics. These statistical results were compared with the results of laboratory-scale analogue modelling. A model similar to that used by Shea in 2005 provided data indicating similar trends of the influence of volume, fall height and topography on mass movement deposits at all scales. The final aspect of this project was a numerical simulation of a large debris avalanche from the north flank of the Taranaki volcano in the direction of the city of New Plymouth. The numerical code VolcFlow developed by Kelfoun in 2005 was used, after being tested against the laboratory experiments to verity its accuracy. The simulations showed that the Pouaki range protects the city of New Plymouth form major impacts from Taranaki collapses, but also indicated some potential problems with the hazard zoning and evacuation zones presently in place.

Sasol Limited experiences extremely high particulate fouling rates inside shell-and-tube heat exchangers that utilize process cooling water. The water and foulants are obtained from various natural and process sources and have irregular fluid properties. The fouling eventually obstructs flow on the shell side of the heat exchanger to such an extent that the tube bundles have to be replaced every nine months. Sasol requested that certain aspects of this issue be addressed. To better understand the problem, the effects of various tube and baffle configurations on the sedimentation rate in a shell-and-tube heat exchanger were numerically investigated. Single-segmental, double-segmental and disc-and-doughnut baffle configurations, in combination with square and rotated triangular tube configurations, were simulated by using the CFD software package, STAR-CCM+. In total, six configurations were investigated. The solution methodology was divided into two parts. Firstly, steady-state solutions of the six configurations were used to identify the best performing model in terms of large areas with high velocity flow. The results identified both single-segmental baffle configurations to have the best performance. Secondly, transient multiphase simulations were conducted to investigate the sedimentation characteristics of the two single-segmental baffle configurations. It was established that the current state of available technology cannot adequately solve the detailed simulations in a reasonable amount of time and results could only be obtained for a time period of a few seconds. By simulating the flow fields for various geometries in steady-state conditions, many of the observations and findings of literature were verified. The single-segmental baffle configurations have higher pressure drops than double-segmental and disc-and-doughnut configurations. In similar fashion, the rotated triangular tube configuration has a higher pressure drop than the square arrangement. The single-segmental configurations have on average higher flow velocities and reduced cross-flow mass flow fractions. It was concluded from this study that the single-segmental baffle with rotated triangular tube configuration had the best steady-state performance. Some results were extracted from the transient multiphase simulations. The transient multiphase flow simulation of the single-segmental baffle configurations showed larger concentrations of stagnant sediment for the rotated triangular tube configuration versus larger concentrations of suspended/flowing sediment in the square tube configuration. This result was offset by the observation that the downstream movement of sediment was quicker for the rotated triangular tube configuration. No definitive results could be obtained, but from the available results, it can be concluded that the configuration currently implemented at Sasol is best suited to handle sedimentation. This needs to be verified in future studies by using advanced computational resources and experimental results. / Thesis (MIng (Mechanical Engineering))--North-West University, Potchefstroom Campus, 2013

Sasol Limited experiences extremely high particulate fouling rates inside shell-and-tube heat exchangers that utilize process cooling water. The water and foulants are obtained from various natural and process sources and have irregular fluid properties. The fouling eventually obstructs flow on the shell side of the heat exchanger to such an extent that the tube bundles have to be replaced every nine months. Sasol requested that certain aspects of this issue be addressed. To better understand the problem, the effects of various tube and baffle configurations on the sedimentation rate in a shell-and-tube heat exchanger were numerically investigated. Single-segmental, double-segmental and disc-and-doughnut baffle configurations, in combination with square and rotated triangular tube configurations, were simulated by using the CFD software package, STAR-CCM+. In total, six configurations were investigated. The solution methodology was divided into two parts. Firstly, steady-state solutions of the six configurations were used to identify the best performing model in terms of large areas with high velocity flow. The results identified both single-segmental baffle configurations to have the best performance. Secondly, transient multiphase simulations were conducted to investigate the sedimentation characteristics of the two single-segmental baffle configurations. It was established that the current state of available technology cannot adequately solve the detailed simulations in a reasonable amount of time and results could only be obtained for a time period of a few seconds. By simulating the flow fields for various geometries in steady-state conditions, many of the observations and findings of literature were verified. The single-segmental baffle configurations have higher pressure drops than double-segmental and disc-and-doughnut configurations. In similar fashion, the rotated triangular tube configuration has a higher pressure drop than the square arrangement. The single-segmental configurations have on average higher flow velocities and reduced cross-flow mass flow fractions. It was concluded from this study that the single-segmental baffle with rotated triangular tube configuration had the best steady-state performance. Some results were extracted from the transient multiphase simulations. The transient multiphase flow simulation of the single-segmental baffle configurations showed larger concentrations of stagnant sediment for the rotated triangular tube configuration versus larger concentrations of suspended/flowing sediment in the square tube configuration. This result was offset by the observation that the downstream movement of sediment was quicker for the rotated triangular tube configuration. No definitive results could be obtained, but from the available results, it can be concluded that the configuration currently implemented at Sasol is best suited to handle sedimentation. This needs to be verified in future studies by using advanced computational resources and experimental results. / Thesis (MIng (Mechanical Engineering))--North-West University, Potchefstroom Campus, 2013

The objective of this thesis is to further understand the austenite Abnormal Grain Growth (AGG) phenomenon in relation with precipitation state in a low alloy steel. The abnormal grain growth is addressed from both experimental and numerical modeling point of view. Prior austenite grain size distribution, precipitation volume fraction and size distribution evolution of the different heat treated states are experimentally determined for two different industrial alloys (steel-A and steel-B) in different heat treated states and experimental results are compared with model predictions. A two-step modeling technique is adopted in this study: precipitation modeling and abnormal/normal grain growth modeling. The abnormal/normal grain growth modeling is done using a simplified analytical model where the grain growth is assumed to be driven by the decrease in interfacial energy. Both the conventional Zener pinning and corner pinning by precipitate is considered as boundary movement retarding forces. The precipitation model is based on the Classical Nucleation and Growth Theories. The assumption of homogeneous precipitate nucleation and growth gave a good prediction of volume fraction, mean radius and size distribution in comparison with the experimental results. Two coupled modeling approaches of abnormal grain growth and precipitation model: ①Soft coupling and ②Dynamic coupling; shed light on the different physical parameters controlling the grain growth condition in a particular material's state. A reasonable prediction of AGG and NGG is obtained from both approaches. The dynamic coupled modeling enabled us to paint a comprehensive time-temperature mechanism map of grain growth conditions. It is found that AGG in the austenitic state depends strongly on the initial grain size distribution and precipitation state. The modeling and the experimental results showed that the precipitation state evolution (increasing or decreasing volume fraction) also impact normal/abnormal grain growth. Plausible explanations in relation with the mean austenite grain size and the precipitation state are derived for the AGG phenomenon from the present work.

A quasi-two dimensional model for the southern Vancouver Island shelf was developed with the Regional Ocean Modelling System (ROMS) to study coupling of the carbon, oxygen, and nitrogen cycles in a summer wind-driven upwelling region. The physical model is coupled to an ecosystem module that includes a simple representation of a sediment layer and considers non-fixed C:N ratios for detritus and dissolved organic matter (i.e., explicitly modelled pools of carbon and nitrogen for those variables). The model accounts for denitrification within the sediments as well as within the water column when oxygen concentrations are low (below 5 mmol-O2 m-3). The objective is to identify the dominant processes controlling the cycles, their coupling, and their sensitivity to changes in environmental forcing. Results demonstrate how low oxygen and low pH events are tightly coupled in the coastal study region, especially through local ecosystem processes. In particular, exchange with the sediments plays a dominant role in consuming oxygen from and releasing inorganic carbon to the bottom waters on the shelf. Two key features distinguish the southern Vancouver Island shelf from other coastal regions in the California Current System and protect inner shelf waters from severe hypoxia and corrosive (i.e., undersaturated in aragonite) conditions. First, the greater width of the shelf reduces the penetration of subsurface offshore high-carbon and low-oxygen waters into shallower waters; and second, the relatively fresh Vancouver Island Coastal Current (VICC) brings oxygen-rich and carbon-poor waters to the bottom layer over the inner shelf. Sensitivity experiments show that carbon and oxygen cycles on the southern Vancouver Island shelf may be significantly affected by an altered upwelling season, a shallower offshore Oxygen Minimum Zone, a warmer ocean, and a carbon-enriched environment. Combinations of these scenarios suggest a potential increasing risk for the development of coastal hypoxia and corrosive conditions in the future. Further sensitivity simulations indicate that sedimentary denitrification provides an additional coupling between the carbon, oxygen, and nitrogen cycles. Total alkalinity generated by sediment denitrification has the potential to buffer anthropogenic ocean acidification. However, this alkalinity effect over the Vancouver Island shelf in late spring and summer simulations is small compared with studies for other locations at annual scales. Longer time scales need to be examined in this region to confirm whether the role of alkalinity generation in the sediments is significant. In conclusion, this dissertation not only demonstrates the coupled nature of biogeochemical cycles in the coastal ocean, but also the importance of this coupling as we try to estimate how coastal ecosystems will respond to human modifications of shelf waters and the climate.

Disertacijos tyrimo objektu buvo Pietryčių Lietuvos kvartero komplekso daugiasluoksnės uolienų storymės požeminio vandens ištekliai (balansas, pasiskirstymas ir formavimas). Disertacinio darbo metu pagrįsti ir suformuluoti kvartero daugiasluoksnės uolienų storymės schematizavimo principai – patikslinta hidrogeodinaminės sistemos sandara remiantis litologiniais duomenimis ir vandeningų sluoksnių tarpusavio ryšio sąlygomis, naudojant hidrogeocheminius, izotopinius ir hidraulinius tyrimus. Išskirti vandeningieji ir pusiau laidūs sluoksniai bei sudarytas kvartero vandeningojo komplekso matematinis modelis. Atlikta ekspertinė modeliavimui naudotų hidrogeologinių parametrų analizė – jų patikimumas tikslintas kalibruojant modelį hidrogeologiniais, infiltracinės mitybos, požeminio nuotėkio balansiniais paskaičiavimais bei požeminio vandens cheminės sudėties ir izotopiniais tyrimais. Parinktame eksperimentiniame hidrogeologiniame pjūvyje – vandens filtracijos kelyje įvertinta požeminio vandens cheminės-izotopinės sudėties transformacija. Remiantis tričio tyrimo ir matematinio modelio duomenimis regioniniame pjūvyje atlikta filtracijos schemos bei filtracijos greičių korekcija. Pirmą kartą, naudojant matematinio modeliavimo metodus, įvertintos Pietryčių Lietuvos kvartero vandeningojo komplekso gamtinių ir eksploatacinių prognozinių gėlo požeminio vandens išteklių formavimosi sąlygos, nustatytas galimas požeminio vandens gavybos poveikis aplinkai. / The aim of the research was accomplish an assessment of the groundwater resources (budget, distribution and formation) and reveal their formation peculiarities in the Quaternary aquifer system of the southeaster part of Lithuania. Substantiation and formulation of schematization principles of the Quaternary multi-aquifer system for conceptualization of groundwater flow model based on the complex of lithologic, hydrogeochemistry, isotopic and hydraulic data have been developed by this dissertation. Performing data analysis of hydraulically tested aquifers, the hydrogeological parameters which have been used in simulation and their reliability have been adjusted during model calibration based on hydrogeological–hydrological water balance methods with coupled groundwater isotopic-chemistry investigations. Transformation of groundwater isotope-hydrogeochemical composition in the direction from the recharge area toward the discharge places have been estimated on a selected experimental hydrogeological cross–section. According to tritium and simulation modelling data, groundwater flow scheme and average velocity comparison have been made from the highland to lowland. The pre–development groundwater flow and prognostic assessment of groundwater exploitation resources of Quaternary aquifer system with extraction impact on the geoenvironment have been evaluated.

The aim of the research was accomplish an assessment of the groundwater resources (budget, distribution and formation) and reveal their formation peculiarities in the Quaternary aquifer system of the southeaster part of Lithuania. Substantiation and formulation of schematization principles of the Quaternary multi-aquifer system for conceptualization of groundwater flow model based on the complex of lithologic, hydrogeochemistry, isotopic and hydraulic data have been developed by this dissertation. Performing data analysis of hydraulically tested aquifers, the hydrogeological parameters which have been used in simulation and their reliability have been adjusted during model calibration based on hydrogeological–hydrological water balance methods with coupled groundwater isotopic-chemistry investigations. Transformation of groundwater isotope-hydrogeochemical composition in the direction from the recharge area toward the discharge places have been estimated on a selected experimental hydrogeological cross–section. According to tritium and simulation modelling data, groundwater flow scheme and average velocity comparison have been made from the highland to lowland. The pre–development groundwater flow and prognostic assessment of groundwater exploitation resources of Quaternary aquifer system with extraction impact on the geoenvironment have been evaluated. / Disertacijos tyrimo objektu buvo Pietryčių Lietuvos kvartero komplekso daugiasluoksnės uolienų storymės požeminio vandens ištekliai (balansas, pasiskirstymas ir formavimas). Disertacinio darbo metu pagrįsti ir suformuluoti kvartero daugiasluoksnės uolienų storymės schematizavimo principai – patikslinta hidrogeodinaminės sistemos sandara remiantis litologiniais duomenimis ir vandeningų sluoksnių tarpusavio ryšio sąlygomis, naudojant hidrogeocheminius, izotopinius ir hidraulinius tyrimus. Išskirti vandeningieji ir pusiau laidūs sluoksniai bei sudarytas kvartero vandeningojo komplekso matematinis modelis. Atlikta ekspertinė modeliavimui naudotų hidrogeologinių parametrų analizė – jų patikimumas tikslintas kalibruojant modelį hidrogeologiniais, infiltracinės mitybos, požeminio nuotėkio balansiniais paskaičiavimais bei požeminio vandens cheminės sudėties ir izotopiniais tyrimais. Parinktame eksperimentiniame hidrogeologiniame pjūvyje – vandens filtracijos kelyje įvertinta požeminio vandens cheminės-izotopinės sudėties transformacija. Remiantis tričio tyrimo ir matematinio modelio duomenimis regioniniame pjūvyje atlikta filtracijos schemos bei filtracijos greičių korekcija. Pirmą kartą, naudojant matematinio modeliavimo metodus, įvertintos Pietryčių Lietuvos kvartero vandeningojo komplekso gamtinių ir eksploatacinių prognozinių gėlo požeminio vandens išteklių formavimosi sąlygos, nustatytas galimas požeminio vandens gavybos poveikis aplinkai.

This Thesis reports on the integration of different modelling techniques to construct a unified and better constrained conceptual model of structural evolution of the Pico del Águila anticline (External Sierras, Southern Pyrenees, Spain). The structure is a well-known example of detachment fold, which exhibits a N-S structural trend, parallel to the direction of tectonic transport in the Southern Pyrenees. Based on field observations of an unevenly distributed Triassic décollement, analogue modelling show how to generate orogen-perpendicular structures which may result in transverse anticlines. Numerical models investigate the effect of a complex mechanical stratigraphy, characterized by an interlayering of competent and incompetent layers, plus syn-kinematic sedimentation in the fold growth. Based on field data and seismic interpretations, a 3D reconstruction and sequential geomechanical restoration of the Pico del Águila anticline reports the coexistence of multiple folding mechanisms occurring simultaneously in different units and structural domains of the fold, leading to a complex strain pattern that can not be assessed by simplistic kinematic 2D approaches. By integrating the presented models with the previous data in the region, we discuss about the benefits and drawbacks of each modelling technique and present an integrated model of structural evolution for the Pico del Águila anticline. This brings a better comprehension of the structure as well as the processes that drove the evolution of the N-S detachment anticlines in the External Sierras of the Southern Pyrenees. / L’anticlinal del Pico del Águila és un conegut exemple de plec de desenganxament amb sedimentació marina a fluvio-deltaica associada, amb una tendència estructural N-S, paral•lela a la direcció de transport tectònic dels Pirineus Meridionals. Basat en observacions de camp que indiquen una distribució heterogènia del nivell de desenganxament Triàsic, els models analògics mostren el procés de generació d’estructures perpendiculars al sistema orogènic que poden donar lloc als anticlinals N-S descrits a les Sierras Exteriores Aragonesas. Els models numèrics investiguen l’efecte d’una estratigrafia mecànica complexa, caracteritzada per una intercalació d’unitats competents i incomptetents (amb marcades diferències en el grau de competència, per tant) i de la sedimentació sin-cinemàtica en el creixement i evolució de l’anticlinal. Basat en dades de camp i interpretacions de perfils sísmics s’ha portat a terme una reconstrucció i restitució geomecànica tridimensional de l’anticlinal del Pico del Águila. D’aquestes se’n deriva la coexistència de de múltiples mecanismes de plegament actuant simultàniament en diferents unitats i dominis estructurals, la qual implica al seu torn un patró i distribució de la deformació que no poden ser avaluats mitjançant aproximacions o tècniques de modelització cinemàtiques o/i bidimensionals. S’integren també els resultats obtinguts a partir de les esmentades tècniques de modelització amb les dades i coneixements previs de la regió, es discuteixen els beneficis, desavantatges i limitacions de cadascuna d’aquestes tècniques de modelització, i es presenta un model integrat d’evolució estructural del Pico del Águila. Aquest anàlisi crític dels resultats i aquest esforç d’integració porten sense dubte cap a una millor comprensió de l’estructura i dels processos que menaren l’evolució dels plecs de desenganxament N-S de les Sierras Exteriores Aragonesas dels Pirineus Meridionals.

Low-lying coastal plains comprised of unconsolidated infill are internally complex hydrogeological settings, due to the high level of heterogeneity in the infill material. In order to resolve the hydrogeological processes active in these complex settings, an integrated multi-disciplinary, geoscientific approach is required. This research determines quantitatively, the effects of sedimentary aquifer heterogeneity on groundwater flowpaths and groundwater processes within a heavily laterised, coastal plain setting. The study site is the Bells Creek catchment in southeast Queensland, Australia. The methodology developed in this study provides a new approach to enable the determination of groundwater flowpaths and groundwater processes at macroscale resolution within other shallow alluvial and coastal plain aquifers. The multi-disciplinary approach utilises sedimentological, geophysical, chronological and hydrogeological techniques (including hydrochemistry and groundwater flow modelling) to develop a high-resolution aquifer framework, and to determine accurately, both groundwater flowpaths and relative flow rates. Sedimentary framework is confirmed to be the principal factor controlling the distribution of aquifer permeability pathways in any given setting, and is therefore, the dominant control over groundwater flow and processes. For the Bells Creek catchment, interpretation of stratigraphic and sedimentary data allowed the compilation of a detailed sedimentary framework. This interpretation demonstrated that weathering of the low-lying arkose sandstone bedrock has developed thick lateritic profiles. Within the weathering profiles, cemented, iron-rich horizons have resisted erosion and developed raised and elongated ridges in the modern landscape, while other clay-rich weathered layers have submitted to erosion and downgraded around those iron-rich ridges. Consequently, alluvial deposition throughout the Late Quaternary has been restricted to narrow, and relatively deep valleys containing sandrich channels, and thin floodplains at shallow depth. From a hydrogeological perspective, there is significant macroscopic aquifer heterogeneity between fine-grained lateritic mixed clay layers, floodplain clays, ironcemented ferricrete horizons, and permeable sand-rich alluvial aquifers. This variability of aquifer material has created a complex subsurface arrangement of permeability pathways. Application of Ground Penetrating Radar (GPR) in this setting enabled accurate definition of alluvial channel boundaries and the high degree of connectedness within the channels themselves. Interpretation of a comprehensive GPR dataset (that covered the entire catchment) allowed refinement of the sedimentary framework previously established to develop a detailed threedimensional aquifer framework. Finite-difference groundwater modelling and particle tracking analysis (using MODFLOW and MODPATH) has clearly demonstrated that the macroscopic heterogeneity within the various aquifer materials of the plain has marked impacts on groundwater pathways, and especially groundwater travel times. The variability between a maximum residence time of 18 months for groundwater within the alluvium, compared to hundreds of years for groundwater within the mixed clay layers of the laterite, clearly demonstrates the importance of accurately defining the spatial distribution of the various aquifer materials in a groundwater flow investigation. In this setting, the interconnection of the narrow alluvial channels feeding into a deeper alluvial delta has provided an effective conduit for shallow groundwater flow. The role of the alluvial delta in discharging the bulk of fresh groundwater from the central plain into the coastal and estuarine aquifers to the east, is certainly critical in preventing saline intrusion from encroaching further west. Hydrochemical and isotopic indicators have identified the dominant recharge processes and groundwater flowpaths within the plain, and indicated that the processes are strongly related to sub-surface permeability distributions determined in the aquifer framework (and groundwater modelling), as well as seasonal fluctuations in rainfall. In the northwest of the plain, sandstone hills provide a delayed and slightly mineralized component of groundwater recharge into adjacent highly permeable, unconfined alluvial aquifers; these aquifers also recharge directly via precipitation. Aluminosilicate weathering in the bedrock hills and eastern peripheries of the laterised bedrock are a source of excess Na, SiO2, and HCO3 to the alluvial groundwater. As this groundwater flows down-gradient to the east, however, its chemical composition evolves by sulfate reduction, silica equilibrium and ion exchange processes into a more mature Na-Cl type. Within the shallow coastal aquifers proximal to the eastern shoreline, sulfate enrichment is occurring (associated with increases in Ca, HCO3, Fe and Al) resulting in major deterioration in groundwater quality. The deterioration is produced by saline intrusion from the adjacent estuary coupled with oxidation of sulfide materials in shallow marine and estuarine clays. Reverses in salinity in those coastal aquifers have been correlated with surges in fresh recharge waters from unconfined coastal dunes and semi-confined landward alluvium, following significant rainfall events. The multi-disciplinary methodology developed, provides an effective approach for accurately defining the three-dimensional distribution of shallow aquifer material of varying permeability via detailed stratigraphic interpretation and GPR analysis. Utilising this aquifer framework, finite-difference groundwater modelling aided by hydrogeological data and hydrochemical analysis, allows accurate determination of groundwater flowpaths and groundwater processes. This research provides a new hydrogeological analogue for alluvial channel aquifers within a laterised coastal plain setting. Key Words: groundwater flow, aquifer heterogeneity, numerical modelling, hydrochemistry, recharge, ground penetrating radar, coastal plain aquifers, weathering, alluvial channels.

Geodynamic modelling, via computer simulations, offers an easily controllable method for investigating the behaviour of an Earth system and providing feedback to conceptual models of geological evolution. However, most available computer codes have been developed for engineering or hydrological applications, where strains are small and post-failure deformation is not studied. Such codes cannot simultaneously model large deformation and porous fluid flow. To remedy this situation in the face of tectonic modelling, a numerical approach was developed to incorporate porous fluid flow into an existing high-deformation code called Ellipsis. The resulting software, with these twin capabilities, simulates the evolution of highly deformed tectonic regimes where fluid flow is important, such as in mineral provinces. A realistic description of deformation depends on the accurate characterisation of material properties and the laws governing material behaviour. Aside from the development of appropriate physics, it can be a difficult task to find a set of model parameters, including material properties and initial geometries, that can reproduce some conceptual target. In this context, an interactive system for the rapid exploration of model parameter space, and for the evaluation of all model results, replaces the traditional but time-consuming approach of finding a result via trial and error. The visualisation of all solutions in such a search of parameter space, through simple graphical tools, adds a new degree of understanding to the effects of variations in the parameters, the importance of each parameter in controlling a solution, and the degree of coverage of the parameter space. Two final applications of the software code and interactive parameter search illustrate the power of numerical modelling within the feedback loop to field observations. In the first example, vertical rheological contrasts between the upper and lower crust, most easily related to thermal profiles and mineralogy, exert a greater control over the mode of crustal extension than any other parameters. A weak lower crust promotes large fault spacing with high displacements, often overriding initial close fault spacing, to lead eventually to metamorphic core complex formation. In the second case, specifically tied to the history of compressional orogenies in northern Nevada, exploration of model parameters shows that the natural reactivation of early normal faults in the Proterozoic basement, regardless of basement topography or rheological contrasts, would explain the subsequent elevation and gravitationally-induced thrusting of sedimentary layers over the Carlin gold trend, providing pathways and ponding sites for mineral-bearing fluids.