No description available.

Yes / Field measurements and numerical simulations demonstrate the existence of an upper saline plume in tidally dominated beaches. The effect of tides on the saltwater-freshwater mixing occurring at both the upper saline plume and lower salt wedge is well understood. However, it is poorly understood whether the tidal driven force acts equally on the mixing behaviours of above two regions and what factors control the mixing fluctuation features. In this study, variable-density, saturated-unsaturated, transient groundwater flow and solute transport numerical models are proposed and performed for saltwater-freshwater mixing subject to tidal forcing on a sloping beach. A range of tidal amplitude, fresh groundwater flux, hydraulic conductivity, beach slope and dispersivity anisotropy are simulated. Based on time sequential salinity data, the gross mixing features are quantified by computing the spatial moments in three different aspects, namely, the centre point, length and width, and the volume (or area in a two-dimensional case). Simulated salinity distribution varies significantly at saltwater-freshwater interfaces. Mixing characteristics of the upper saline plume greatly differ from those in the salt wedge for both the transient and quasi-steady state. The mixing of the upper saline plume largely inherits the fluctuation characteristics of the sea tide in both the transverse and longitudinal directions when the quasi-steady state is reached. On the other hand, the mixing in the salt wedge is relatively steady and shows little fluctuation. The normalized mixing width and length, mixing volume and the fluctuation amplitude of the mass centre in the upper saline plume are, in general, one-magnitude-order larger than those in the salt wedge region. In the longitudinal direction, tidal amplitude, fresh groundwater flux, hydraulic conductivity and beach slope are significant control factors of fluctuation amplitude. In the transverse direction, tidal amplitude and beach slope are the main control parameters. Very small dispersivity anisotropy (e.g., α_L⁄α_T <5) could greatly suppress mixing fluctuation in the longitudinal direction. This work underlines the close connection between the sea tides and the upper saline plume in the aspect of mixing, thereby enhancing understanding of the interplay between tidal oscillations and mixing mechanisms in tidally dominated sloping beach systems. / Shenzhen Key Laboratory for Coastal Ocean Dynamics and Environment (No. ZDSY20130402163735964), National High Technology Research & Development Program of China (No. 2012AA09A409).

When properly designed, hypolimnetic aeration and oxygenation systems can replenish dissolved oxygen in water bodies while preserving stratification. A comprehensive literature review of design methods for the three primary devices was completed. Using fundamental principles, a discrete-bubble model was first developed to predict plume dynamics and gas transfer for a circular bubble-plume diffuser. This approach has subsequently been validated in a large vertical tank and applied successfully at full-scale to an airlift aerator as well as to both circular and linear bubble-plume diffusers. The unified suite of models, all based on simple discrete-bubble dynamics, represents the current state-of-the-art for designing systems to add oxygen to stratified lakes and reservoirs. An existing linear bubble plume model was improved, and data collected from a full-scale diffuser installed in Spring Hollow Reservoir, Virginia (U.S.A.) were used to validate the model. The depth of maximum plume rise was simulated well for two of the three diffuser tests. Temperature predictions deviated from measured profiles near the maximum plume rise height, but predicted dissolved oxygen profiles compared very well to observations. Oxygen transfer within the hypolimnion was independent of all parameters except initial bubble radius. The results of this work suggest that plume dynamics and oxygen transfer can successfully be predicted for linear bubble plumes using the discrete-bubble approach. To model the complex interaction between a bubble plume used for hypolimnetic oxygenation and the ambient water body, a model for a linear bubble plume was coupled to two reservoir models, CE-QUAL-W2 (W2) and Si3D. In simulations with a rectangular basin, predicted oxygen addition was directly proportional to the update frequency of the plume model. W2 calculated less oxygen input to the basin than Si3D and significantly less mixing within the hypolimnion. The coupled models were then applied to a simplified test of a full-scale linear diffuser. Both the W2 and Si3D coupled models predicted bulk hypolimnetic DO concentrations well. Warming within the hypolimnion was overestimated by both models, but more so by W2. The lower vertical resolution of the reservoir grid in W2 caused the plume rise height to be over-predicted, enhancing erosion of the thermocline. / Ph. D.

Our understanding of how magma-poor rifts accommodate strain remains limited largely due to sparse geophysical observations from these rift systems. To better understand magma-poor rifting processes, chapter 1 of this dissertation is focused on investigating the lithosphere-asthenosphere interactions beneath the Malawi Rift, a segment of the magma-poor Western Branch of the East African Rift (EAR). Chapter 2 and 3 are focused on investigating the sources of melt beneath the Rungwe Volcanic Province (RVP), an anomalous volcanic center located at the northern tip of the Malawi Rift. In chapter 1, we use the lithospheric structure of the Malawi Rift derived from the World Gravity Model 2012 to constrain three-dimensional (3D) numerical models of lithosphere-asthenosphere interactions, which indicate ~3 cm/yr asthenospheric upwelling beneath the thin lithosphere (115-125 km) of the northern Malawi Rift and the RVP from lithospheric modulated convection (LMC) that is decoupling from surface motions. We suggest that the asthenospheric upwelling may generate decompression melts which weakens the lithosphere thereby enabling extension. The source of asthenospheric melt for the RVP is still contentious. Some studies suggest the asthenospheric melt beneath the RVP arises from thermal perturbations in the upper mantle associated with plume head materials, while others propose decompression melting from upwelling asthenosphere due to LMC where the lithosphere is thin. Chapter 2 of this dissertation is focused on testing the hypothesis that asthenospheric melt feeding the RVP can be generated from LMC using realistic constraints on the mantle potential temperature (Tp). We develop a 3D thermomechanical model of LMC beneath the RVP and the entire Malawi Rift that incorporates melt generation. We find decompression melt associated with LMC upwelling (~3 cm/yr) occurs at a maximum depth of ~150 km localized beneath the RVP. Studies of volcanic rock samples from the RVP indicate plume signatures which are enigmatic since the RVP is highly localized, unlike the large igneous provinces in the Eastern Branch of the EAR. In chapter 3, we test the hypothesis that the melt beneath the RVP is generated from plume materials. We investigate melt generation from plume-lithosphere interactions (PLI) beneath the RVP by developing a 3D seismic tomography-based convection (TBC) model beneath the RVP. The seismic constraints indicate excess temperatures of ~250 K in the sublithospheric mantle beneath the RVP suggesting the presence of a plume. We find a relatively fast upwelling (~10 cm/yr) beneath the RVP which we interpret as a rising plume. The TBC upwelling generates decompression melt (~0.25 %) at a maximum depth of ~200 km beneath the RVP where the lithosphere is thinnest (~100 km). Our results demonstrate that an excess heat source from may be plume materials is necessary for melt generation in the sublithospheric mantle beneath the RVP because passive asthenospheric upwelling of ambient mantle will require a higher than normal Tp to generate melt. Studies of volcanic rock samples from the RVP indicate plume signatures which are enigmatic since the RVP is highly localized, unlike the large igneous provinces in the Eastern Branch of the EAR. In chapter 3, we test the hypothesis that the melt beneath the RVP is generated from plume materials. We investigate melt generation from plume-lithosphere interactions (PLI) beneath the RVP by developing a 3D seismic tomography-based convection (TBC) model beneath the RVP. The seismic constraints indicate excess temperatures of ≈ 250K in the sublithospheric mantle beneath the RVP suggesting the presence of a plume. We find a relatively fast upwelling (≈10 cm/yr) beneath the RVP which we interpret as a rising plume. The TBC upwelling generates decompression melt (≈0.25 %) at a maximum depth of ≈200 km beneath the RVP where the lithosphere is thinnest (≈100 km). Our results demonstrate that an excess heat source from may be plume materials is necessary for melt generation in the sublithospheric mantle beneath the RVP because passive asthenospheric upwelling of ambient mantle will require a higher than normal Tp to generate melt. / Doctor of Philosophy / Studies suggest the presence of hot, melted rock deep in the continents makes them weaker and easier to break apart, however, our understanding of how continents with less melted rock break apart remains limited largely due to sparse geophysical observations from these dry areas. To better understand how continents with less melted rock break apart, chapter 1 of this dissertation is focused on investigating the interactions between the rigid part of the Earth, called lithosphere, and the underlying lower viscosity rock layer called asthenosphere beneath the Malawi Rift, a segment of the magma-poor Western Branch of the East African Rift (EAR). Chapter 2 and 3 are focused on investigating the sources of melt beneath the Rungwe Volcanic Province (RVP), an anomalous volcanic center located at the northern tip of the Malawi Rift. In chapter 1, we use the lithospheric structure of the Malawi Rift derived from gravity data to constrain three-dimensional (3-D) numerical models of lithosphere-asthenosphere interactions, which indicate ~3 cm/yr asthenospheric upwelling beneath the thin lithosphere (115-125 km) of the northern Malawi Rift and the RVP that does not seem to drive movements at the surface. We suggest that the asthenospheric upwelling may generate melted rock which weakens the lithosphere thereby enabling extension. However, the source of asthenospheric melt for the RVP is still contentious. Some studies suggest the asthenospheric melt beneath the RVP arises from thermal perturbations in the upper mantle associated with rising mantle rocks or plume head materials, while others propose melting occurs from upwelling asthenosphere due to lithospheric modulated convection (LMC) where the lithosphere is thin. Chapter 2 of this dissertation is focused on testing the hypothesis that asthenospheric melt feeding the RVP can be generated from LMC. We develop a 3D thermomechanical model of LMC beneath the RVP and the entire Malawi Rift that incorporates melt generation. We find decompression melt associated with LMC upwelling (~3 cm/yr) occurs at a maximum depth of ~150 km localized beneath the RVP. Studies of volcanic rock samples from the RVP indicate plume signatures which are enigmatic since the RVP is highly localized, unlike the large igneous provinces in the Eastern Branch of the EAR. In chapter 3, we investigate melt generation from plume-lithosphere interactions (PLI) beneath the RVP. We develop a 3D model of convection using information from seismology we call tomography-based convection (TBC) beneath the RVP. The seismic data indicate excess temperatures of ~250 K beneath the RVP suggesting the presence of a plume. We find a relatively fast upwelling (~10 cm/yr) beneath the RVP which we interpret as a rising plume. The TBC upwelling generates decompression melt at a maximum depth of ~200 km beneath the RVP. Our results demonstrate that an excess heat source from may be plume materials is necessary for melt generation in the sublithospheric mantle beneath the RVP because passive asthenospheric upwelling of ambient mantle will require a higher than normal mantle potential temperatures to generate melt.

Bubble plume oxygenation systems are commonly used to mitigate anoxia and its deleterious effects on water quality in thermally stratified lakes and reservoirs. Following installation, increases in sediment oxygen flux (JO2) are typically observed during oxygenation and are positively correlated with the bubble plume gas flow rate. Studies show that JO2 is controlled by the thickness of the diffusive boundary layer (DBL) at the sediment-water interface (SWI), which is in turn controlled by turbulence. As a result, JO2 can be quite spatially and temporally variable. Accurately predicting oxygenation-induced JO2 is vitally important for ensuring successful oxygenation system design and operation. Yet despite the current understanding of physical and chemical controls on JO2, methods for predicting oxygenation-induced JO2 are still based on empirical correlations and factors of safety. As hypolimnetic oxygenation becomes more widely used as a lake management tool for improving and maintaining water quality, there is a need to move from the current empirically based approach to a mechanistic approach and improve the ability to predict induced JO2. This work details field campaigns to investigate and identify appropriate models of oxygen supply to the SWI and oxygen demand exerted from the sediment, with the intent to use these models to predict oxygenation-induced JO2. Oxygen microprofiles across the SWI and near-sediment velocity measurements were collected in situ during three field campaigns on two oxygenated lakes, providing simultaneous measurements of JO2 and turbulence. Field observations show that oxygenation can increase JO2 by increasing bulk hypolimnetic oxygen concentrations, which increases the concentration gradient across the SWI. Oxygenation can also enhance turbulence, which decreases the DBL thickness and increases JO2. Existing models of interfacial flux were compared to field measurements to determine which model best predicted the observed JO2. Models based on the Batchelor scale, friction velocity, and film-renewal theory all agree reasonably well with field observations in both lakes. Additionally, the oxygen microprofiles were used to fit a transient model of oxygen kinetics in lake sediment and determine the appropriate kinetic model. Oxygen microprofiles in both lakes can be described using zero-order kinetics, rather than first-order kinetics. The interfacial flux and sediment kinetic models are incorporated into a coupled bubble plume and 3-D hydrodynamic lake model, allowing for spatial and temporal variation in simulated JO2. This comprehensive model was calibrated and validated to field data from two separate field campaigns on Carvin's Cove Reservoir, Virginia. Simulated temperature profiles agreed quite well with field observations, while simulated oxygen profiles differed from observed profiles, particularly in the bottom 1 m of the water column. The model overestimates oxygen concentrations near the sediment, which results in higher simulated JO2 than was observed during the field campaigns. These discrepancies are attributed to oxygen-consuming chemical processes, such as oxidation of soluble metals, which are not accounted for in the hydrodynamic model. Despite this, the model is still able to capture the impact of bubble plume operation on JO2, as simulated JO2 is higher when the diffusers are operating. With some additional improvements to the water quality modeling aspects of the model, as well as further calibration and validation, the model should be able to reproduce observed JO2 provided oxygen concentrations near the SWI are accurately reproduced as well. The current work is an attempt to push toward a comprehensive lake oxygenation model. A comprehensive model such as this should improve the ability to predict oxygenation-induced JO2 and lead to improvements in the design and operation of hypolimnetic oxygenation systems. / Ph. D.

Two series of experiments were performed to determine the flow characteristics and to quantify the heat transfer components from a propane diffusion flame impinging onto a ceiling. A 0.3 m square sand burner with propane as the fuel type provided a steady-state fire. In the first series of experiments, measurements of gas temperature and velocity were made at 76 mm vertical intervals above the burner up to the ceiling. Fire heat release rates (HRRs) of 50 kW and 90 kW with free flame length to ceiling height ratios, Lf/H, of 2, 1.5, 1, 0.8, 0.85 were used to determine their effects on the measured parameters. Gas temperatures within the continuous flaming region were relatively constant, and measured to be independent of ceiling height and HRR, while velocities increased with elevation and were independent of ceiling height yet weakly dependent on HRR. Within the intermittent region, gas temperature was weakly affected by the presence of the ceiling at various heights, while the effect on velocity was more pronounced. HRR had an effect on both temperature and velocity within the intermittent region of the fire plume. Comparisons with existing fire plume correlations showed that the unbounded correlations can be used to provide a good approximation of the gas temperature for the ceiling bounded case; while the correlations for the velocity can only be used for elevations up to approximately 60% of the ceiling height. Elevations above this cutoff were significantly affected by the presence of the ceiling. The second series of experiments investigated HRRs of 50 kW and 90 kW with free flame length to ceiling height ratios, Lf/H, of 2, 1.5, and 1. Heat flux and gas temperature at the stagnation point of the ceiling were measured using hybrid heat flux gauges and an aspirated Type K thermocouple. Four methods of calculating the convective heat transfer coefficient, h, were developed and adapted; two reference methods and two slope methods. The components of heat transfer at the impingement point were separated using these calculated h values. The reference method 2, and both slope methods only required the use of the non-cooled hybrid gauge measurements and were in overall good agreement with one another. The reference method 1 differed significantly, being up to 15.8 times lower than the others. The trends in the two groups were contradictory, with the h calculated using the reference method 1 increasing with ceiling height while the others showed no strong trend with ceiling height. The disagreements between the methods greatly affected the components of heat transfer, particularly at the lowest ceiling heights. Convection calculated using the h from reference method 1 contributed only 2-5% of the total exposure heat flux at the lowest ceiling heights, whereas with the other methods convection contributed 20-50% of the total exposure heat flux. The limitations of each method are discussed. Further investigation is required for all methods to determine their applicability within the flaming region of a fire. / Master of Science

Variações dos tempos de percurso de fases P e PKP foram usados para tomografia sísmica do manto superior sob o sudeste e centro-oeste do Brasil. Os principais objetivos foram: Melhorar a resolução obtida pelos estudos anteriores (VanDecar et al,. 1995; Escalante, 2002; Schimmel et al., 2003, feitos de 1992 a 2001) com a inclusão de novos dados e estações, mapear as áreas que não tinham sido cobertas pelas estações anteriores, tentar mostrar a consistência da base de dados, principalmente relacionada às novas estações e verificar a robustez das anomalias encontradas. Foi estudada a influência das estações localizadas nas regiões anômalas através de inversões secundárias (retirando estações). Nesta base foram incluídos registros recentes do ano 2002 e também em registros de 2000 e 2001 para eventos utilizados nos trabalhos anteriores. Os resultados confirmaram as estruturas observadas nos trabalhos anteriores e também revelaram novas regiões anômalas, particularmente no sul do estado de Mato Grosso. Das anomalias observadas nos trabalhos anteriores, confirmamos em nossos resultados o Cráton do São Francisco como uma anomalia de alta velocidade, com suas raízes chegando a 250 km de profundidade, a anomalia de baixa velocidade interpretada como resto fóssil da pluma de Tristan da Cunha VanDecar et al. (1995), a boa correlação das anomalias rasas de baixa velocidade com as intrusões alcalinas do Cretáceo Superior Schimmel et al., (2003), a anomalia de baixa velocidade inferida por Escalante (2002) na região de Iporá em Goiás (possível região do impacto inicial da pluma de Trindade), uma anomalia de alta velocidade sob a região da Bacia do Paraná (profundidades rasas) interpretada inicialmente por Schimmel et al. (2003) como possível núcleo cratônico da Bacia do Paraná e também, nesta mesma região (maiores profundidades), uma anomalia de alta velocidade interpretada como a litosfera subduzida da Placa de Nazca (Schimmel et al., 2003 e Escalante, 2002). A nova base permitiu a expansão da área de estudo e o imageamento de anomalias de baixa velocidade na província ígnea de Poxoréu em Mato Grosso, as quais são consistentes com o afinamento litosferico proposto no modelo de Thompson et al. (1998). / Variations of P and PKP travel times were used for seismic tomography of the upper mantle beneath southeast and central Brazil. Our principal objectives were: To improve the resolution obtained by the previous studies (carried out by VanDecar et al., 1995; Escalante, 2002; Schimmel et al., 2003 between 1992 and 2001) with inclusion of new data and stations, to map areas not covered by previous stations, to show the consistency of the database, mainly of the related to the new stations and to verify the robustness of the detected anomalies. The influence of the stations located in the anomalous areas was studied through secondary inversions (removing stations). The new data set includes recent records of 2002 and also new records from 2000 and 2001 for events used in the previous works. Our results confirm the structures observed in the previous works and also revealed new anomalous regions, particularly in the south of the Mato Grosso state. We confirmed the anomalies observed in the previous works: The São Francisco cráton has as a high-velocity anomaly, with roots down to 250 km depth, the low velocity anomaly interpreted as a fossil remnant of the Tristão da Cunha plume (VanDecar et al. 1995), the good correlation of the shallow low velocity anomalies with the alkaline intrusions of the Late-Cretaceous (Schimmel et al. 2003), the low velocity anomaly inferred by Escalante (2002) in the Iporá igneous province in Goiás (possible area of the initial impact of the plume of Trindade), a high-velocity anomaly under the Paraná Basin (shallow depths) interpreted initially by Schimmel et al. (2003) as possible cratonic nucleus of the Paraná Basin and also, in this same area (larger depths), a high-velocity anomaly interpreted as the slab of the Nazca Plate (Schimmel et al., 2003 and Escalante, 2002). The new data set allowed the expansion of the study area and the imaging of low velocity anomalies in the igneous province of Poxoréu in Mato Grosso, which are consistent with the model of litospheric thinning proposed by Thompson et al. (1998).

The Late Archean Magmatic Event (LAME), ca. 2.7 Ga, was the greatest crustal addition event in Earth history. My focus is the Superior Province of Canada, where LAME occurred ca. 2.75 – 2.65 Ga. Mantle plumes impinged on the Abitibi subprovince, where ~ 16 km regional thickness of tonalite-trondhjemite-granodiorite (TTG) melt was produced. Granites (sensu stricto) were the last magmatic phase of LAME, with a Superior-wide regional thickness of ~ 1 – 3 km. Assuming a crustal source for both TTG and granites, I use 1-D thermal models to investigate the origin of TTG in the Abitibi subprovince and that of late granites in the Superior Province. Melting curves appropriate to the source of TTG and granites are used to determine the thickness of melt produced in the models. I show that the incorporation of upward melt transfer into a standard model of lower crustal melting may increase the amount of predicted melt by ~ 1/(1-f), where f denotes the fraction of melt that is on average being extracted from the source rocks. Partitioning of heat producing elements between melt and restite reduces the amount of melt produced, but the effect is secondary compared to the increase in melt production through upward melt transfer. For the Abitibi subprovince, I show that the emplacement of a single plume coupled with the emplacement of a 12-km-thick greenstone cover can generate a maximum of ~ 9-km-thickness of TTG melt. However, the emplacement of a series of plumes, each coupled with the emplacement of a 3-km-thick greenstone cover and a 10-km-thick sill results in ~ 20-km-thickness of TTG melt. My model incorporates delamination of restitic eclogite. Finally, I show that late granites in the Superior Province may have resulted from thickening of a crust that had been “pre-heated” during earlier arc activity and that prolonged granitic magmatism observed in some areas of the Superior Province may be explained by late underthrusting of fertile source rocks into deeper and hotter regions of the crust.

The Late Archean Magmatic Event (LAME), ca. 2.7 Ga, was the greatest crustal addition event in Earth history. My focus is the Superior Province of Canada, where LAME occurred ca. 2.75 – 2.65 Ga. Mantle plumes impinged on the Abitibi subprovince, where ~ 16 km regional thickness of tonalite-trondhjemite-granodiorite (TTG) melt was produced. Granites (sensu stricto) were the last magmatic phase of LAME, with a Superior-wide regional thickness of ~ 1 – 3 km. Assuming a crustal source for both TTG and granites, I use 1-D thermal models to investigate the origin of TTG in the Abitibi subprovince and that of late granites in the Superior Province. Melting curves appropriate to the source of TTG and granites are used to determine the thickness of melt produced in the models. I show that the incorporation of upward melt transfer into a standard model of lower crustal melting may increase the amount of predicted melt by ~ 1/(1-f), where f denotes the fraction of melt that is on average being extracted from the source rocks. Partitioning of heat producing elements between melt and restite reduces the amount of melt produced, but the effect is secondary compared to the increase in melt production through upward melt transfer. For the Abitibi subprovince, I show that the emplacement of a single plume coupled with the emplacement of a 12-km-thick greenstone cover can generate a maximum of ~ 9-km-thickness of TTG melt. However, the emplacement of a series of plumes, each coupled with the emplacement of a 3-km-thick greenstone cover and a 10-km-thick sill results in ~ 20-km-thickness of TTG melt. My model incorporates delamination of restitic eclogite. Finally, I show that late granites in the Superior Province may have resulted from thickening of a crust that had been “pre-heated” during earlier arc activity and that prolonged granitic magmatism observed in some areas of the Superior Province may be explained by late underthrusting of fertile source rocks into deeper and hotter regions of the crust.

Oil well blow outs are investigated to determine methods to quickly and accurately respond to an emergency situation. Flow rate is needed to guide containment and dispersal operations. The Stratified Integral Multiphase Plume, SIMP, model was used to investigate the range of initial conditions available to integral modeling. Sensitivity to initial conditions is modest, but without experimental data at the appropriate scale the most accurate condition is unable to be determined. Flow rates are difficult to directly measure in blow out situations, so another method must be determined; therefore, sensitivity of several parameters to flow rate was also evaluated. Methane concentration in the first intrusion can be used in conjunction with velocity and trap height measurements to determine flow rate using an integral model. Plume width and temperature were determined to have little sensitivity. Separately, a containment dome was tested in the laboratory to determine if a full scale dome can be used to contain an oil leak in the field. The dome was found to have satisfactory entrapment in the designed position.