The subject of this dissertation is artisanal mining in Namibia and is based on two hypotheses. The first is that if the artisanal mining process is well understood in Namibia, then effective solutions to improving the artisanal mining cycle of poverty can be developed. The second hypothesis states that if effective solutions to improving the artisanal mining cycle of poverty can be implemented, then the miners' income can be increased.Research for the dissertation is broken down into four areas. The first area includes the development of an understanding of the artisanal mining process. A general model that shows the relationship of four basic elements of artisanal mining (investment, production, sales and consumption) as well as several components that make up the elements was created. The model also describes the various environments (physical, social/cultural, and political) that artisanal operate.The Artisanal Mining Process Model was then used to evaluate artisanal mining in Namibia. Artisanal mining in Namibia is somewhat unique because of the material mined (collector minerals and semi-precious gemstones), but the outcome is consistent with other locations with most artisanal miners only making a subsistence living. One of the key outcomes of this part of the study is the identification of low sales revenue exacerbated by the miners having to sell in a local market with few buyers.Because the local market and buyers are a critical component of the cycle of poverty, the third area of research is an understanding of the local market and the supply chain that ultimately leads to a much larger international market. The international markets lead to the fourth area of research, which is the Tucson Mineral Show, the largest of its kind in the world.The conclusions of the paper discuss the applicability of the Artisanal Mining Process Model in helping to understand issues facing the artisanal miners as well as methods that could be used to help the artisanal miners participate in the international market for collector minerals.

Large opencast coal mines require a complex infrastructure to fulfill production demand and quality values. The distinct specifications required by each customer are achieved by blending adjustments. There is limited control in variability. With only partial information available, operation controllers blend coal by empirical approximation, trying to keep quality between acceptable ranges in order to avoid penalizations, shipment rejections or even contract suspensions. When a decision support system (DSS) centralized in a control room is used for blending control, crew operators visualize enhanced displays of the different sources of information, obtaining a holistic perspective of operations. Using a simulator to reproduce the blending sequence, crew operators can experiment with diverse what-if scenarios and develop blending strategies for an entire working shift, in which they also incorporate their own expertise and the knowledge obtained after interpreting the simulation results. The research focuses on the empirical analysis of the effectiveness of the DSS by studying the performance of crew users in different operating scenarios produced with a simulator. The development of a methodology for measuring this effectiveness and its impact in the quantification of controlling the variability of blending represents a significant contribution in the area of quality improvement for coal production. The effectiveness of the DSS for controlling the blending and load out processes has been numerically measured after experimenting diverse simulated scenarios, proving that the difference between estimated and actual quality delivered is narrower when using a DSS, in comparison with the BTU variability obtained from historical data. The strategies that produced better results in terms of control of coal quality variability, maximization of infrastructure utilization, time spent in making decisions and the minimization of risk for penalizations and rejections, were scored proportionally to the benefits obtained.

The research presented in this dissertation consists of four topics. The first of these topics is an experimental study of rock fracturing due to rapid thermal cooling, and the other three topics are related to mine-to-optimization. This includes the development and testing of a site-specific model for blast fragmentation, the development of a technique for utilizing digital image processing and ground-based LIDAR for rock mass characterization, and an experimental study of the effects of ore blending on mineral recovery. All four topics are related through the subject of rock fracturing and rock fragmentation. The results from this research are important and can be used to improve engineering design associated with rock excavation and rock fragmentation. First of all, a successful set of laboratory experiments and 3D numerical modeling was conducted, looking at the effects of rapid thermal cooling on rock mechanical properties. The results gave the unexpected finding that depending on the rock type and the thermal conditions, rapid cooling can result in either overall crack growth or crack closing. Secondly, a site-specific model for predicting blast fragmentation was developed and tested at an open-pit copper mine in Arizona. The results provide a practical technique for developing a calibrated blasting model using digital images and digital image processing software to estimate in-situ block size, and a calibrated Schmidt hammer to estimate intact tensile strength. Thirdly, a new technique was developed to conduct cell mapping in open-pit mines using the new technologies of digital image processing and ground-based LIDAR. The results show that the use of these new technologies provide an increased accuracy and the ability for more sophisticated slope stability analyses with no increase in field time only a moderate increase in data processing time. Finally, a successful set of laboratory experiments was conducted looking at the effects of ore blending and grinding times on mineral recovery from a set of six ore from a copper mine in Arizona. The results gave the unexpected finding that for a fixed grinding time, the mineral recovery of the blended ores exceeded the average of the individual recoveries of the same ores unblended.

Geophysical imaging methods significantly enhance our knowledge of subsurface characteristics and their use has become prevalent over a range of subsurface investigations. These methods facilitate the detection and characterization of both metallic and nonmetallic subsurface targets, and can provide spatially extensive information on subsurface structure and characteristics that is often impractical to obtain using standard drilling and sampling procedures alone. Electrical imaging methods such as electrical resistivity tomography (ERT) have proven to be particularly useful in hydrogeologic and geotechnical investigations because of the strong dependence of the electrical properties of soils to water saturation, soil texture, and solute concentration. Given the available geophysical tools as well as their applications, the selection of the appropriate geophysical survey design is an essential part of every subsurface geophysical investigation. Where investigations are located in an area with subsurface information already available, this information may be used as a guide for the design of a geophysical survey. In some instances, no subsurface information is available and a survey must be designed to cover a range of possible circumstances. Yet, in other instances, there may be significant subsurface information available, but because of subsurface complexities, a geophysical survey must still be designed to cover a broad range of possibilities. Demonstrating the application and limitations of ERT in a specific field application, the first investigation presented in this document provides guidance for developing methods to improve the design and implementation of ERT surveys in a complex subsurface environment. The two investigations that follow present the development of a relatively simple optimization approach based on limited forward modeling of the geophysical response for both static and mobile surveys. This process is demonstrated through examples of selecting a limited number of ERT surveys to identify and discriminate subsurface target tunnels (with a simple cylindrical geometry). These examples provide insights into the practical application of the optimization process for improved ERT survey design for subsurface target detection. Because of their relative simplicity, the optimization procedures developed here may be used to rapidly identify optimal array configurations without the need for computationally expensive inversion techniques.

No description available.

Approved for public release, distribution is unlimited / Performance of high resolution mesoscale models has been in a continuous state of refinement since their inception. Mesoscale models have become quite skillful in forecasting synoptic scale events such as mid-latitude cyclones. However, atmospheric forcing becomes a much more complicated process when faced with the challenge of forecasting near topography along the coastline. Phenomena such as gap flows, blocked flow winds and low level stratification become important to predictability at these scales. The problem is further complicated by the dynamics of a frontal passage event. The skill of mesoscale models in predicting these winds is not as well developed. This study examines several forecasts by the Coupled Ocean Atmospheric Mesoscale Prediction System (COAMPS) during frontal passage events for the Winter of 2003-2004. An attempt is made to characterize the predictability of the wind speed and direction both before and after frontal passage along the California coast. Synoptic forcing during this time is strong due to the effects of the mid-latitude cyclones propagate across the Pacific. The study's results indicate that the wind field predictability is subject to several consistent errors associated with the passage of fronts over topography. These errors arise due to difficulty in the model capturing weak thermal advection events and topographic wind funneling. The deficiencies in model representation of topography contributes to these errors. / Lieutenant, United States Navy

The characterization of rock properties is a vital task in the challenge for hard rock mining operation. A simplified and straightforward characterization of rock properties provides information about the safety of ground structure (e.g. slope, tunnel, etc.), and the strategy to improve productivity in terms of rock breakage process. The penetration-rate of drilling has been proposed to quantify the comminution characteristics of rock by virtue of real-time logging of drilling performance otherwise the data is obtained from a time- and cost-consuming laboratory test; this is called measurement while drilling. In the mining industry, this technique can be a useful tool that has allowed for the meticulous and routine data collection of geological information from blasthole drilling operations. In this study, the mechanical performance of drill and its interaction with the rock properties is investigated in laboratory scale. The rock properties include tensile strength, hardness, and grindability, which is considered as the influential parameters of the required energy consumption for the comminution processes. For sandstone samples, the penetration-rate data shows a good correlation with tensile strength, hardness, and Bond work index; this implies that penetration-rate data can be a good indicator to estimate comminution characteristics. Additionally we carried out the same test with limestone samples. Second, field study is conducted to investigate the interaction between current blast design and rock fragmentation. Fabricating the blast design and fragmentation through the blast operation might enable to construct proper strategy to reduce the energy cost of downstream processes including crushing and grinding by using the rock characteristics measured from the blasthole drilling. The concept of this process is a part of Mine-to-Mill optimization. The thesis proposed the blueprint of Mine-to-Mill optimization, providing a guideline for further in-situ research.

The Collisional Orogeny in the Scandinavian Caledonides (COSC) project aims to provide a deeper understanding of mountain belt dynamics through scientific deep drilling in the central parts of the mountain belt of western Sweden. The main targets include a subduction related allochthon, the basal orogenic detachment and the underlying partially subducted Precambrian basement. Research covered by this thesis, focusing primarily on reflection seismic data, was done within the framework of the COSC project. The 55 km long composite COSC Seismic Profile (CSP) images the upper crust in high resolution and established the basis for the selection of the optimum location for the two 2.5 km deep COSC boreholes. Together with potential field and magnetotelluric data, these profiles allowed the construction of a constrained regional interpretation of the major tectonic units. Non-conventional pseudo 3D processing techniques were applied to the 2D data prior to the drilling of the first borehole, COSC-1, to provide predictions about the 3D geometry of subsurface structures and potential zones of interest for the sampling programs. COSC-1 was drilled in 2014 and reached the targeted depth with nearly complete core recovery. A continuous geological section and a wealth of information from on-site and off-site scientific investigations were obtained. A major post-drilling seismic survey was conducted in and around the borehole and included a 3D reflection seismic experiment. The structurally and lithologically complex Lower Seve Nappe proved difficult to image in detail using standard processing techniques, but its basal mylonite zone and underlying structures are well resolved. The 3D data, from the surface down to the total drilled depth, show good correlation with the initial mapping of the COSC-1 core as well as with preliminary results from on-core and downhole logging.   Good correlation is also observed between the 2D and 3D reflection seismic datasets. These will provide a strong link between the two boreholes and a means to extrapolate the results from the cores and boreholes into the surrounding rock. Ultimately, they will contribute to the deeper understanding of the tectonic evolution of the region, the Scandinavian Caledonides and the formation of major orogens. / Collisional Orogeny in the Scandinavian Caledonides

A thesis submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg, in fulfilment of the requirements for the degree of Doctor of Philosophy Johannesburg, August 2015 School of Geosciences, University of the Witwatersrand / The possibility of extensive shale gas resources in the main Karoo Basin has resulted in a renewed focus on the basin, and particularly the Whitehill Formation. The main Karoo Basin has been the subject of geological studies since before the 1920s, but geophysical data provides an opportunity to shed new light on the basin architecture and formation. In this thesis, I use regional gravity, magnetic and borehole data over the basin, as well as vintage seismic data in the southern part of the basin. Modern computational capacity allows for more information to be extracted from these seismic data, and for these data to be better integrated with potential field data. The integration of datasets in a three-dimensional model (3D) has allowed for a better understanding of the shape of the basin and its internal structure, in turn shedding light on basin formation. A new depth map of the basin constructed using this extensive database confirms that the basin deepens from on- to off-craton. The basin is deepest along the northern boundary of the Cape Fold Belt (CFB), with a depth of ~4000 m in the southwestern Karoo and ~5000 m in the southeastern part of the basin. Sediment thickness ranges from ~5500 to 6000 m. The Whitehill Formation along this boundary reaches a depth of ~ 3000 m in the southwest and ~4000 m in the southeast. Despite limited boreholes in this region, the basin appears to broadly deepen to the southeast. These seismic and borehole data also allow for mapping of the Cape Supergroup pinch-out below the Karoo basin (32.6°S for the Bokkeveld and 32.4°S for the Table Mountain Group), with the basin reaching a thickness of around 4 km just north of the CFB. The gravity effect of these sediments in the south is not sufficient to account for the low of the Cape Isostatic Anomaly near Willowmore and Steytlerville. This ~45 mGal Bouguer gravity low dominates the central region of the southern Karoo at the northern border of the CFB. The seismic data for the first time show uplift of lower-density shales of the Ecca Group (1800 – 2650 kg/m3) in this region, and structural and seismic data suggest that these lower density sediments continue to depth of 11 to 12 km along normal and thrust faults in this region. Two-dimensional density models show that these shallow crustal features, as well as deeper lower crust compared to surrounding regions, account for the anomaly. These seismic and borehole data also allow for constraints to be placed on the distribution and geometry of the dolerite intrusions that intruded the basin after its formation, and in some cases impacted on the shale layer, to be constrained. The highest concentrations of dolerites are found in the northwest and east of the basin, pointing towards two magma sources. The region of lowest concentration is in the south-central part of the basin. Here the intrusions are confined to the Beaufort Group, ~1000 m shallower than the shale reservoir, suggesting it should be the focus of exploration efforts. These dolerite sills are shown to be between 5 and 30 km wide and are saucer-shaped with ~ 800 m vertical extent, and dips of between 2° and 8° on the edges. The sheets in the south of the basin extend for over 150 km, dipping at between 3° and 13°, and are imaged down to ~ 5 km. This change in dip of the sheets is linked to deformation within the Cape Fold Belt, with greater dips closer to the belt, although these sheets do not appear to intrude strata dipping at more than 15 to 20°. In order to understand the shape of the Karoo basin and construct a 3D model of the basin, an understanding is needed of the underlying basement rocks. The Beattie Magnetic Anomaly (BMA) that stretches across the entire southern part of the basin forms part of the basement Namaqua-Natal Belt. Filtered magnetic data confirm that the Namaqua and Natal Belts are two separate regions with different magnetic characteristics, which is taken into account during modelling. The BMA is shown to be part of a group of linear magnetic anomalies making up the Natal Belt. The anomaly itself will therefore not have an individual effect on basin formation, and the effect of the Natal Belt as a whole will have to be investigated. An in-depth study of outcrops associated with one of these linear magnetic anomalies on the east coast of South Africa suggest the BMA can be attributed to regions of highly magnetic (10 to 100 x 10-3 SI) supracrustal rocks in Proterozoic shear zones. Along two-dimensional magnetic models in the southwestern Karoo constrained by seismic data, these magnetic zones are modelled as dipping slabs with horizontal extents of ~20-60 km and vertical extents of ~10-15 km. Body densities range from 2800- 2940 kg/m3 and magnetic susceptibilities from 10 to 100 x 10-3 SI. These, as well as other geophysical and geological constraints, are used to construct a 3D model of the basin down to 300 km. Relatively well-constrained crustal structure allows for inversion modelling of lithospheric mantle densities using GOCE satellite gravity data, with results in-line with xenolith data. These results confirm the existence of lower density mantle below the craton (~3270 kg/m3) that could contribute to the buoyancy of the craton, and an almost 50 kg/m3 density increase in the lithospheric mantle below the surrounding Proterozoic belts. It is this change in lithospheric density along with changes in Moho depths that isostatically compensate a large portion of South Africa’s high topography (<1200 m). The topography higher than 1200 m along the edge of the plateau, along the Great Escarpment, are shown to be accommodated by an asthenospheric buoyancy anomaly with a density contrast of around 40 kg/m3, while still mimicking the Bouguer gravity field. These findings are in line with recent tomographic studies below Africa suggesting an “African Superplume” or “Large Low Velocity Seismic Province” in the deep mantle. The basin sediment thickness maps were further used to investigate the formation of the main Karoo Basin. This was accomplished by studying the past flexure of the Whitehill Formation using north-south two-dimensional (2D) profiles. Deepening of the formation from ~3000 m in the southwest to ~4000 m in the southeast is explained using the concept of isostasy, i.e., an infinite elastic beam that is subjected to an increasing load size across the Cape Fold Belt. Load height values increase from 4 km in the southwest to 8 km in the southeast. This larger load is attributed here to “locking” along a subduction zone further to the south. The effective elastic thickness (Te) of the beam also increases from around 50 km over the Namaqua and Natal Belts in the southwest to 80 km over the Kaapvaal Craton and Natal Belt in the southeast. The changes in Te values do not correlate with changes in terrane, i.e., a north to south change, as previously though. The large extent and shape of the Karoo basin can therefore, in general, be explained as a flexural basin, with the strength of the basement increasing towards the southeast. Therefore, while factors such as mantle flow could have contributed towards basin formation, reducing the load size needed, it is no longer necessary in order to account for the large extent of the basin. This flexure model breaks down further to the southeast, most likely due to a very high Te value. This could be the reason for later plate break in this region during Gondwana breakup. It is inferred that this increase in Te is linked to the buoyancy anomaly in the asthenospheric mantle.

A dissertation submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg, in fulfilment of the requirements for the degree of Master of Science. Johannesburg, 2016 / The Mesoproterozoic Pilanesberg Complex, South Africa, is the world’s largest alkaline intrusive complex. Mapped geological field relationships suggest the Complex has circular inward dipping layers. However, it is unclear how the dipping layers extend at depth. As a result, the 3D geometry of the Pilanesberg Complex is unknown. Modelling of the Pilanesberg Complex uses 2D forward models as well as 3D forward and inversion, gravity and magnetic data models, to set limits on the 3D geometry of the Pilanesberg Complex. The 2D Bouguer gravity models and geology maps indicate that some of the Bushveld Complex Main Zone shifted to the west of the Pilanesberg Complex during emplacement. This, and a highly faulted country rock, accounts for a portion of how the host rock was able to accommodate the Pilanesberg Complex intrusion. The geometry of the Complex is explored with test gravity models where the model of outward dipping and vertically dipping cylinders are unable to match the Bouguer gravity signal over the Complex, but the inward dipping model matched the data to provide a possible solution for the geometry of the Complex. The Pilanesberg Complex geometry is modelled with 3D magnetic inversion, 3D forward gravity models and 2.5D gravity test profiles that were all constrained by the surface geology. The different models correlate so that best data fit for the Complex is represented by an overall inward dipping structure. Surface geological measurements indicate that the northern edge of the Complex dip out to the north. The 3D forward modelling was able to produce a positive solution that matched the gravity data with a northward dipping northern edge. The dipping northern edge is also observed on the University of British Columbia, UBC, 3D gravity inversion and the Euler deconvolution gravity profile solutions. The depth of the Pilanesberg Complex from 3D forward gravity modelling is estimated to be between 5 and 6 km. The Complex is suggested to have undergone block movement where the northern block and southern block are separated by the 30 km long Vlakfontein fault, which bisects the Complex from the north-east to the south-west. The image processing contact depth, Euler deconvolution solutions and the 3D Voxi inversion model suggest that the fresh bedrock is closer to surface in the north, while the southern block appears to be approximately 1km deeper than the northern block. The northern dip and block movement are explained by complicated structural events that include trap door graben settling which hinged on the northern edge as well as faulting and external block movement during a regional lateral extensional event. The Pilanesberg Complex intruded during a larger system of alkaline intrusions, known as the Pilanesberg Alkaline Province. The intrusions are associated with the Province due to their ages and chemical affinity. This Province includes two dyke swarms that radiate to the north-west and south of the Pilanesberg Complex, as well as smaller circular clinopyroxenite intrusions throughout the Bushveld Complex. The Pilanesberg dyke system and the circular clinopyroxenite intrusions are reversely magnetised with IGRF corrected values ranging between -150 to -320 nT compared to the normally magnetised 166 to 330 nT values of the Pilanesberg Complex. This suggests that a magnetic reversal occurred between the emplacement of the Pilanesberg Complex and the dyke System. The age data of the Complex and dyke Swarm suggest a magnetic reversal could have occurred between the emplacement of the Pilanesberg Complex and the Pilanesberg dyke System. The Complex is dated at 1602 ± 38 Ma and 1583 ± 10 Ma, from two white foyaite samples from the southern edge (using 40Ar/39Ar amphibole spectrum analysis). These ages are vastly different from previously reported ages, which ranged between 1200 Ma and 1450 Ma (Harmer R., 1992; Hansen et al., 2006). The error analysis has improved considerably from the published dates making the proposed dates plausible for the intrusion of the Pilanesberg Complex as the first and main intrusion of the Pilanesberg Alkaline Province. The Pilanesberg dyke System intruded much later between 1219 ± 6 Ma to 1268 ± 10 Ma for the red syenite dyke samples (using 40Ar/39Ar on feldspars spectrum analysis) and 1139 ± 18 Ma obtained for the grey syenite dyke (using 40Ar/39Ar on amphiboles inverse isochronal analysis). The dyke Swarm dates are significantly younger than the previously published ages for the dykes, which were between 1290 Ma and 1330 Ma (Van Niekerk, 1962; Emerman, 1991). / LG2017