Identifying the Structure and Fate of Wastewater Derived Organic Micropollutants by High-resolution Mass Spectrometry

Human activities represent a significant burden on the global water cycle, with large and increasing demands placed on limited water resources by manufacturing, energy production and domestic water use. In addition to changing the quantity of available water resources, human activities lead to changes in water quality by introducing a large and often poorly-characterized array of chemical pollutants, which may negatively impact biodiversity in aquatic ecosystems, leading to impairment of valuable ecosystem functions and services. Domestic and industrial wastewaters represent a significant source of pollution to the aquatic environment due to inadequate or incomplete removal of chemicals introduced into waters by human activities. Currently, incomplete chemical characterization of treated wastewaters limits comprehensive risk assessment of this ubiquitous impact to water. In particular, a significant fraction of the organic chemical composition of treated industrial and domestic wastewaters remains uncharacterized at the molecular level. Efforts aimed at reducing the impacts of water pollution on aquatic ecosystems critically require knowledge of the composition of wastewaters to develop interventions capable of protecting our precious natural water resources.

The goal of this dissertation was to develop a robust, extensible and high-throughput framework for the comprehensive characterization of organic micropollutants in wastewaters by high-resolution accurate-mass mass spectrometry. High-resolution mass spectrometry provides the most powerful analytical technique available for assessing the occurrence and fate of organic pollutants in the water cycle. However, significant limitations in data processing, analysis and interpretation have limited this technique in achieving comprehensive characterization of organic pollutants occurring in natural and built environments. My work aimed to address these challenges by development of automated workflows for the structural characterization of organic pollutants in wastewater and wastewater impacted environments by high-resolution mass spectrometry, and to apply these methods in combination with novel data handling routines to conduct detailed fate studies of wastewater-derived organic micropollutants in the aquatic environment.

In Chapter 2, chemoinformatic tools were implemented along with novel non-targeted mass spectrometric analytical methods to characterize, map, and explore an environmentally-relevant “chemical space” in municipal wastewater. This was accomplished by characterizing the molecular composition of known wastewater-derived organic pollutants and substances that are prioritized as potential wastewater contaminants, using these databases to evaluate the pollutant-likeness of structures postulated for unknown organic compounds that I detected in wastewater extracts using high-resolution mass spectrometry approaches. Results showed that application of multiple computational mass spectrometric tools to structural elucidation of unknown organic pollutants arising in wastewaters improved the efficiency and veracity of screening approaches based on high-resolution mass spectrometry. Furthermore, structural similarity searching was essential for prioritizing substances sharing structural features with known organic pollutants or industrial and consumer chemicals that could enter the environment through use or disposal.

I then applied this comprehensive methodological and computational non-targeted analysis workflow to micropollutant fate analysis in domestic wastewaters (Chapter 3), surface waters impacted by water reuse activities (Chapter 4) and effluents of wastewater treatment facilities receiving wastewater from oil and gas extraction activities (Chapter 5). In Chapter 3, I showed that application of chemometric tools aided in the prioritization of non-targeted compounds arising at various stages of conventional wastewater treatment by partitioning high dimensional data into rational chemical categories based on knowledge of organic chemical fate processes, resulting in the classification of organic micropollutants based on their occurrence and/or removal during treatment. Similarly, in Chapter 4, high-resolution sampling and broad-spectrum targeted and non-targeted chemical analysis were applied to assess the occurrence and fate of organic micropollutants in a water reuse application, wherein reclaimed wastewater was applied for irrigation of turf grass. Results showed that organic micropollutant composition of surface waters receiving runoff from wastewater irrigated areas appeared to be minimally impacted by wastewater-derived organic micropollutants. Finally, Chapter 5 presents results of the comprehensive organic chemical composition of oil and gas wastewaters treated for surface water discharge. Concurrent analysis of effluent samples by complementary, broad-spectrum analytical techniques, revealed that low-levels of hydrophobic organic contaminants, but elevated concentrations of polymeric surfactants, which may effect the fate and analysis of contaminants of concern in oil and gas wastewaters.

Taken together, my work represents significant progress in the characterization of polar organic chemical pollutants associated with wastewater-impacted environments by high-resolution mass spectrometry. Application of these comprehensive methods to examine micropollutant fate processes in wastewater treatment systems, water reuse environments, and water applications in oil/gas exploration yielded new insights into the factors that influence transport, transformation, and persistence of organic micropollutants in these systems across an unprecedented breadth of chemical space.

Spectral sensitivity, spatial resolution and temporal resolution and their implications for conspecific signalling in cleaner shrimp.

Cleaner shrimp (Decapoda) regularly interact with conspecifics and client reef fish, both of which appear colourful and finely patterned to human observers. However, whether cleaner shrimp can perceive the colour patterns of conspecifics and clients is unknown, because cleaner shrimp visual capabilities are unstudied. We quantified spectral sensitivity and temporal resolution using electroretinography (ERG), and spatial resolution using both morphological (inter-ommatidial angle) and behavioural (optomotor) methods in three cleaner shrimp species: Lysmata amboinensis, Ancylomenes pedersoni and Urocaridella antonbruunii. In all three species, we found strong evidence for only a single spectral sensitivity peak of (mean ± s.e.m.) 518 ± 5, 518 ± 2 and 533 ± 3 nm, respectively. Temporal resolution in dark-adapted eyes was 39 ± 1.3, 36 ± 0.6 and 34 ± 1.3 Hz. Spatial resolution was 9.9 ± 0.3, 8.3 ± 0.1 and 11 ± 0.5 deg, respectively, which is low compared with other compound eyes of similar size. Assuming monochromacy, we present approximations of cleaner shrimp perception of both conspecifics and clients, and show that cleaner shrimp visual capabilities are sufficient to detect the outlines of large stimuli, but not to detect the colour patterns of conspecifics or clients, even over short distances. Thus, conspecific viewers have probably not played a role in the evolution of cleaner shrimp appearance; rather, further studies should investigate whether cleaner shrimp colour patterns have evolved to be viewed by client reef fish, many of which possess tri- and tetra-chromatic colour vision and relatively high spatial acuity.

Low-linewidth and tunable single frequency 1x2 multimode-interferometer-Fabry-Perot laser diode

In this paper, we present a novel 1x2 multi-mode-interferometer-Fabry-Perot (MMI-FP) laser diode, which demonstrated tunable single frequency operation with more than 30dB side mode suppression ratio (SMSR) and a tuning range of 25nm in the C and L bands, as well as a 750 kHz linewidth. These lasers do not require material regrowth and high resolution gratings; resulting in a simpler process that can significantly increase the yield and reduce the cost.

Elemental Analysis of Biological Matrices by Laser Ablation High Resolution Inductively Coupled Plasma Mass Spectrometry (LA-HR-ICP-MS) and High Resolution Inductively Coupled Plasma Mass Spectrometry (HR-ICP-MS)

The need for elemental analysis of biological matrices such as bone, teeth, and plant matter for sourcing purposes has emerged within the forensic and geochemical laboratories. Trace elemental analyses for the comparison of aterials such as glass by inductively coupled plasma mass spectrometry (ICP-MS) and laser ablation ICP-MS has been shown to offer a high degree of discrimination between different manufacturing sources. Unit resolution ICP-MS instruments may suffer from some polyatomic interferences including 40Ar16O+, 40Ar16O1H+, and 40Ca16O+ that affect iron measurement at trace levels. Iron is an important element in the analysis of glass and also of interest for the analysis of several biological matrices. A comparison of the nalytical performance of two different ICP-MS systems for iron analysis in glass for determining the method detection limits (MDLs), accuracy, and precision of the measurement is presented. Acid digestion and laser ablation methods are also compared. Iron polyatomic interferences were reduced or resolved by using dynamic reaction cell and high resolution ICP-MS. MDLs as low as 0.03 ìg g-1 and 0.14 ìg g-1 for laser ablation and solution based analyses respectively were achieved. The use of helium as a carrier gas demonstrated improvement in the detection limits of both iron isotopes (56Fe and 57Fe) in medium resolution for the HR-ICP-MS and with a dynamic reaction cell (DRC) coupled to a quadrupole ICP-MS system. The development and application of robust analytical methods for the quantification of trace elements in biological matrices has lead to a better understanding of the potential utility of these measurements in forensic chemical analyses. Standard reference materials (SRMs) were used in the development of an analytical method using HR-ICP-MS and LA-HR-ICP-MS that was subsequently applied on the analysis of real samples. Bone, teeth and ashed marijuana samples were analyzed with the developed method. Elemental analysis of bone samples from 12 different individuals provided discrimination between individuals, when femur and humerus bones were considered separately. Discrimination of 14 teeth samples based on elemental composition was achieved with the exception of one case where samples from the same individual were not associated with each other. The discrimination of 49 different ashed plant (cannabis)samples was achieved using the developed method.

(Table 1) Age control and temporal resolution of paleoceanographic records of the North Atlantic

Short-term changes in sea surface conditions controlling the thermohaline circulation in the northern North Atlantic are expected to be especially efficient in perturbing global climate stability. Here we assess past variability of sea surface temperature (SST) in the northeast Atlantic and Norwegian Sea during Marine Isotope Stage (MIS) 2 and, in particular, during the Last Glacial Maximum (LGM). Five high-resolution SST records were established on a meridional transect (53°N-72°N) to trace centennial-scale oscillations in SST and sea-ice cover. We used three independent computational techniques (SIMMAX modern analogue technique, Artificial Neural Networks (ANN), and Revised Analog Method (RAM)) to reconstruct SST from planktonic foraminifer census counts. SIMMAX and ANN reproduced short-term SST oscillations of similar magnitude and absolute levels, while RAM, owing to a restrictive analog selection, appears less suitable for reconstructing "cold end" SST. The SIMMAX and ANN SST reconstructions support the existence of a weak paleo-Norwegian Current during Dansgaard-Oeschger (DO) interstadials number 4, 3, 2, and 1. During the LGM, two warm incursions of 7°C water to occurred in the northern North Atlantic but ended north of the Iceland Faroe Ridge. A rough numerical estimate shows that the near-surface poleward heat transfer from 53° across the Iceland-Faroe Ridge up to to 72° N dropped to less than 60% of the modern value during DO interstadials and to almost zero during DO stadials. Summer sea ice was generally confined to the area north of 70°N and only rarely expanded southward along the margins of continental ice sheets. Internal LGM variability of North Atlantic (>40°N) SST in the GLAMAP 2000 compilation (Sarnthein et al., 2003, doi:10.1029/2002PA000771; Pflaumann et al., 2003, doi:10.1029/2002PA000774) indicates maximum instability in the glacial subpolar gyre and at the Iberian Margin, while in the Nordic Seas, SST was continuously low.

LOTUS: A low cost, ultraviolet spectrograph

We describe the design, construction and commissioning of LOTUS; a simple, low-cost long-slit spectrograph for the Liverpool Telescope. The design is optimized for near-UV and visible wavelengths and uses all transmitting optics. It exploits the instrument focal plane field curvature to partially correct axial chromatic aberration. A stepped slit provides narrow (2.5x95 arcsec) and wide (5x25 arcsec) options that are optimized for spectral resolution and flux calibration respectively. On sky testing shows a wavelength range of 3200-6300 Angstroms with a peak system throughput (including detector quantum efficiency) of 15 per cent and wavelength dependant spectral resolution of R=225-430. By repeated observations of the symbiotic emission line star AG Peg we demonstrate the wavelength stability of the system is less than 2 Angstroms rms and is limited by the positioning of the object in the slit. The spectrograph is now in routine operation monitoring the activity of comet 67P/Churyumov-Gerasimenko during its current post-perihelion apparition.

High resolution Thomson Parabola Spectrometer for full spectral capture of multi-species ion beams

We report on the experimental characterisation of laser-driven ion beams using a Thomson Parabola Spectrometer (TPS) equipped with trapezoidally shaped electric plates, proposed by Gwynne et al. [Rev. Sci. Instrum. 85, 033304 (2014)]. While a pair of extended (30 cm long) electric plates was able to produce a significant increase in the separation between neighbouring ion species at high energies, deploying a trapezoidal design circumvented the spectral clipping at the low energy end of the ion spectra. The shape of the electric plate was chosen carefully considering, for the given spectrometer configuration, the range of detectable ion energies and species. Analytical tracing of the ion parabolas matches closely with the experimental data, which suggests a minimal effect of fringe fields on the escaping ions close to the wedged edge of the electrode. The analytical formulae were derived considering the relativistic correction required for the high energy ions to be characterised using such spectrometer.

Convective Cold Pools over the Atlas Mountains and Their Influence on the Saharan Heat Low

The West African Monsoon (WAM) and its representation in numerical models are strongly influenced by the Saharan Heat Low (SHL), a low-pressure system driven by radiative heating over the central Sahara and ventilated by the cold and moist inflow from adjacent oceans. It has recently been shown that a significant part of the southerly moisture flux into the SHL originates from convective cold pools over the Sahel. These density currents driven by evaporation of rain are largely absent in models with parameterized convection. This crucial issue has been hypothesized to contribute to the inability of many climate models to reproduce the variability of the WAM. Here, the role of convective cold pools approaching the SHL from the Atlas Mountains, which are a strong orographic trigger for deep convection in Northwest Africa, is analyzed. Knowledge about the frequency of these events, as well as their impact on large-scale dynamics, is required to understand their contribution to the variability of the SHL and to known model uncertainties. The first aspect is addressed through the development of an objective and automated method for the generation of multi-year climatologies not available before. The algorithm combines freely available standard surface observations with satellite microwave data. Representativeness of stations and influence of their spatial density are addressed by comparison to a satellite-only climatology. Applying this algorithm to data from automated weather stations and manned synoptic stations in and south of the Atlas Mountains reveals the frequent occurrence. On the order of 6 events per month are detected from May to September when the SHL is in its northernmost position. The events tend to cluster into several-days long convectively active periods, often with strong events on consecutive days. This study is the first to diagnose dynamical impacts of such periods on the SHL, based on simulations of two example cases using the Weather Research and Forecast (WRF) model at convection-permitting resolution. Sensitivity experiments with artificially removed cold pools as well as different resolutions and parameterizations are conducted. Results indicate increases in surface pressure of more than 1 hPa and significant moisture transports into the desert over several days. This moisture affects radiative heating and thus the energy balance of the SHL. Even though cold pool events north of the SHL are less frequent when compared to their Sahelian counterparts, it is shown that they gain importance due to their temporal clustering on synoptic timescale. Together with studies focusing on the Sahel, this work emphasizes the need for improved parameterization schemes for deep convection in order to produce more reliable climate projections for the WAM.

Comparison of 2.3 & 5 mega pixel (MP) resolution monitors when detecting mammography image blurring

Background - Image blurring in Full Field Digital Mammography (FFDM) is reported to be a problem within many UK breast screening units resulting in significant proportion of technical repeats/recalls. Our study investigates monitors of differing pixel resolution, and whether there is a difference in blurring detection between a 2.3 MP technical review monitor and a 5MP standard reporting monitor. Methods - Simulation software was created to induce different magnitudes of blur on 20 artifact free FFDM screening images. 120 blurred and non-blurred images were randomized and displayed on the 2.3 and 5MP monitors; they were reviewed by 28 trained observers. Monitors were calibrated to the DICOM Grayscale Standard Display Function. T-test was used to determine whether significant differences exist in blurring detection between the monitors. Results - The blurring detection rate on the 2.3MP monitor for 0.2, 0.4, 0.6, 0.8 and 1 mm blur was 46, 59, 66, 77and 78% respectively; and on the 5MP monitor 44, 70, 83 , 96 and 98%. All the non-motion images were identified correctly. A statistical difference (p <0.01) in the blurring detection rate between the two monitors was demonstrated. Conclusions - Given the results of this study and knowing that monitors as low as 1 MP are used in clinical practice, we speculate that technical recall/repeat rates because of blurring could be reduced if higher resolution monitors are used for technical review at the time of imaging. Further work is needed to determine monitor minimum specification for visual blurring detection.

Atomic resolution studies of oxide superlattices and ultrathin films

The ability to grow ultrathin films layer-by-layer with well-defined epitaxial relationships has allowed research groups worldwide to grow a range of artificial films and superlattices, first for semiconductors, and now with oxides. In the oxides thin film research community, there have been concerted efforts recently to develop a number of epitaxial oxide systems grown on single crystal oxide substrates that display a wide variety of novel interfacial functionality, such as enhanced ferromagnetic ordering, increased charge carrier density, increased optical absorption, etc, at interfaces. The magnitude of these novel properties is dependent upon the structure of thin films, especially interface sharpness, intermixing, defects, and strain, layering sequence in the case of superlattices and the density of interfaces relative to the film thicknesses. To understand the relationship between the interfacial thin film oxide atomic structure and its properties, atomic scale characterization is required. Transmission electron microscopy (TEM) offers the ability to study interfaces of films at high resolution. Scanning transmission electron microscopy (STEM) allows for real space imaging of materials with directly interpretable atomic number contrast. Electron energy loss spectroscopy (EELS), together with STEM, can probe the local chemical composition as well as local electronic states of transition metals and oxygen. Both techniques have been significantly improved by aberration correctors, which reduce the probe size to 1 Å, or less. Aberration correctors have thus made it possible to resolve individual atomic columns, and possibly probe the electronic structure at atomic scales. Separately, using electron probe forming lenses, structural information such as the crystal structure, strain, lattice mismatches, and superlattice ordering can be measured by nanoarea electron diffraction (NED). The combination of STEM, EELS, and NED techniques allows us to gain a fundamental understanding of the properties of oxide superlattices and ultrathin films and their relationship with the corresponding atomic and electronic structure. In this dissertation, I use the aforementioned electron microscopy techniques to investigate several oxide superlattice and ultrathin film systems. The major findings are summarized below. These results were obtained with stringent specimen preparation methods that I developed for high resolution studies, which are described in Chapter 2. The essential materials background and description of electron microscopy techniques are given in Chapter 1 and 2. In a LaMnO3-SrMnO3 superlattice, we demonstrate the interface of LaMnO3-SrMnO3 is sharper than the SrMnO3-LaMnO3 interface. Extra spectral weights in EELS are confined to the sharp interface, whereas at the rougher interface, the extra states are either not present or are not confined to the interface. Both the structural and electronic asymmetries correspond to asymmetric magnetic ordering at low temperature. In a short period LaMnO3-SrTiO3 superlattice for optical applications, we discovered a modified band structure in SrTiO3 ultrathin films relative to thick films and a SrTiO3 substrate, due to charge leakage from LaMnO3 in SrTiO3. This was measured by chemical shifts of the Ti L and O K edges using atomic scale EELS. The interfacial sharpness of LaAlO3 films grown on SrTiO3 was investigated by the STEM/EELS technique together with electron diffraction. This interface, when prepared under specific conditions, is conductive with high carrier mobility. Several suggestions for the conductive interface have been proposed, including a polar catastrophe model, where a large built-in electric field in LaAlO3 films results in electron charge transfer into the SrTiO3 substrate. Other suggested possibilities include oxygen vacancies at the interface and/or oxygen vacancies in the substrate. The abruptness of the interface as well as extent of intermixing has not been thoroughly investigated at high resolution, even though this can strongly influence the electrical transport properties. We found clear evidence for cation intermixing through the LaAlO3-SrTiO3 interface with high spatial resolution EELS and STEM, which contributes to the conduction at the interface. We also found structural defects, such as misfit dislocations, which leads to increased intermixing over coherent interfaces.

Low temperature magneto-photoluminescence characterization of high purity gallium arsenide and indium phosphide

Low-temperature magneto-photoluminescence is a very powerful technique to characterize high purity GaAs and InP grown by various epitaxial techniques. These III-V compound semiconductor materials are used in a wide variety of electronic, optoelectronic and microwave devices. The large binding energy differences of acceptors in GaAs and InP make possible the identification of those impurities by low-temperature photoluminescence without the use of any magnetic field. However, the sensitivity and resolution provided by this technique rema1ns inadequate to resolve the minute binding energy differences of donors in GaAs and InP. To achieve higher sensitivity and resolution needed for the identification of donors, a magneto-photoluminescence system 1s installed along with a tunable dye laser, which provides resonant excitation. Donors 1n high purity GaAs are identified from the magnetic splittings of "two-electron" satellites of donor bound exciton transitions 1n a high magnetic field and at liquid helium temperature. This technique 1s successfully used to identify donors 1n n-type GaAs as well as 1n p-type GaAs in which donors cannot be identified by any other technique. The technique is also employed to identify donors in high purity InP. The amphoteric incorporation of Si and Ge impurities as donors and acceptors in (100), (311)A and (3ll)B GaAs grown by molecular beam epitaxy is studied spectroscopically. The hydrogen passivation of C acceptors in high purity GaAs grown by molecular beam epitaxy (MBE) and metalorganic chemical vapor deposition (MOCVD) 1s investigated using photoluminescence. Si acceptors ~n MBE GaAs are also found to be passivated by hydrogenation. The instabilities in the passivation of acceptor impurities are observed for the exposure of those samples to light. Very high purity MOCVD InP samples with extremely high mobility are characterized by both electrical and optical techniques. It is determined that C is not typically incorporated as a residual acceptor ~n high purity MOCVD InP. Finally, GaAs on Si, single quantum well, and multiple quantum well heterostructures, which are fabricated from III-V semiconductors, are also measured by low-temperature photoluminescence.

Effects of efficient fronto-temporal circuitry on lexical ambiguity resolution: converging evidence from cross-age comparisons in eye-tracking and ERP data

Eye-tracking was used to examine how younger and older adults use syntactic and semantic information to disambiguate noun/verb (NV) homographs (e.g., park). We find that young adults exhibit inflated first fixations to NV-homographs when only syntactic cues are available for disambiguation (i.e., in syntactic prose). This effect is eliminated with the addition of disambiguating semantic information. Older adults (60+) as a group fail to show the first fixation effect in syntactic prose; they instead reread NV homographs longer. This pattern mirrors that in prior event-related potential work (Lee & Federmeier, 2009, 2011), which reported a sustained frontal negativity to NV-homographs in syntactic prose for young adults, which was eliminated by semantic constraints. The frontal negativity was not observed in older adults as a group, although older adults with high verbal fluency showed the young-like pattern. Analyses of individual differences in eye-tracking patterns revealed a similar effect of verbal fluency in both young and older adults: high verbal fluency groups of both ages show larger first fixation effects, while low verbal fluency groups show larger downstream costs (rereading and/or refixating NV homographs). Jointly, the eye-tracking and ERP data suggest that effortful meaning selection recruits frontal brain areas important for suppressing contextually inappropriate meanings, which also slows eye movements. Efficacy of fronto-temporal circuitry, as captured by verbal fluency, predicts the success of engaging these mechanisms in both young and older adults. Failure to recruit these processes requires compensatory rereading or leads to comprehension failures (Lee & Federmeier, in press).

Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets

International audience

Proteomic analysis of Rhizobium freirei PRF 81 responses to low pH.

Rhizobium freirei PRF 81 is employed in common bean commercial inoculants in Brazil, due to its outstanding efficiency in fixing nitrogen, competitiveness and tolerance to abiotic stresses. Among the environmental conditions faced by rhizobia in soils, acidity is perhaps the encountered most, especially in Brazil. So, we used proteomics based approaches to study the responses of PRF 81 to a low pH condition. R. freirei PRF 81 was grown in TY medium until exponential phase in two treatments: pH 6,8 and pH 4,8. Whole-cell proteins were extracted and separated by two-dimensional gel electrophoresis, using IPG-strips with pH range 4-7 and 12% polyacrilamide gels. The experiment was performed in triplicate. Protein spots were detected in the high-resolution digitized gel images and analyzed by Image Master 2D Platinum v 5.0 software. Relative volumes (%vol) of compared between the two conditions tested and were statistically evaluated (p ≤ 0.05). Even knowing that R. freirei PRF 81 can still grow in more acid conditions, pH 4.8 was chosen because didn´t affect significantly the bacterial growth kinetics, a factor that could compromise the analysis. Using a narrow pH range, the gel profiles displayed a better resolution and reprodutibility than using broader pH range. Spots were mostly concentrated between pH 5-7 and molecular masses between 17-95 kDa. From the six hundred well-defined spots analyzed, one hundred and sixty-three spots presented a significant change in % vol, indicating that the pH led to expressive changes in the proteome of R. freirei PRF 81. Of these, sixty-one were up-regulated and one hundred two was downregulated in pH 4.8 condition. Also, fourteen spots were only identified in the acid condition, while seven spots was exclusively detected in pH 6.8. Ninety-five differentially expressed spots and two exclusively detected in pH 4,8 were selected for Maldi-Tof identification. Together with the genome sequencing and the proteome analysis of heat stress, we will search for molecular determinants of PRF 81 related to capacity to adapt to stressful tropical conditions.

Investigating the mesh dependency and upscaling of 3D grain-based models for the simulation of brittle fracture processes in low-porosity crystalline rock

The application of 3D grain-based modelling techniques is investigated in both small and large scale 3DEC models, in order to simulate brittle fracture processes in low-porosity crystalline rock. Mesh dependency in 3D grain-based models (GBMs) is examined through a number of cases to compare Voronoi and tetrahedral grain assemblages. Various methods are used in the generation of tessellations, each with a number of issues and advantages. A number of comparative UCS test simulations capture the distinct failure mechanisms, strength profiles, and progressive damage development using various Voronoi and tetrahedral GBMs. Relative calibration requirements are outlined to generate similar macro-strength and damage profiles for all the models. The results confirmed a number of inherent model behaviors that arise due to mesh dependency. In Voronoi models, inherent tensile failure mechanisms are produced by internal wedging and rotation of Voronoi grains. This results in a combined dependence on frictional and cohesive strength. In tetrahedral models, increased kinematic freedom of grains and an abundance of straight, connected failure pathways causes a preference for shear failure. This results in an inability to develop significant normal stresses causing cohesional strength dependence. In general, Voronoi models require high relative contact tensile strength values, with lower contact stiffness and contact cohesional strength compared to tetrahedral tessellations. Upscaling of 3D GBMs is investigated for both Voronoi and tetrahedral tessellations using a case study from the AECL’s Mine-by-Experiment at the Underground Research Laboratory. An upscaled tetrahedral model was able to reasonably simulate damage development in the roof forming a notch geometry by adjusting the cohesive strength. An upscaled Voronoi model underestimated the damage development in the roof and floor, and overestimated the damage in the side-walls. This was attributed to the discretization resolution limitations.