The Evolution of Massive Young Stellar Objects in the Large Magellanic Cloud

This thesis presents an analysis of the largest catalog to date of infrared spectra of massive young stellar objects in the Large Magellanic Cloud. Evidenced by their very different spectral features, the luminous objects span a range of evolutionary states from those most embedded in their natal molecular material to those that have dissipated and ionized their surroundings to form compact HII regions and photodissociation regions. We quantify the contributions of the various spectral features using the statistical method of principal component analysis. Using this analysis, we classify the YSO spectra into several distinct groups based upon their dominant spectral features: silicate absorption (S Group), silicate absorption and fine-structure line emission (SE), polycyclic aromatic hydrocarbon (PAH) emission (P Group), PAH and fine-structure line emission (PE), and only fine-structure line emission (E). Based upon the relative numbers of sources in each category, we are able to estimate the amount of time massive YSOs spend in each evolutionary stage. We find that approximately 50% of the sources have ionic fine-structure lines, indicating that a compact HII region forms about half-way through the YSO lifetime probed in our study. Of the 277 YSOs we collected spectra for, 41 have ice absorption features, indicating they are surrounded by cold ice-bearing dust particles. We have decomposed the shape of the ice features to probe the composition and thermal history of the ice. We find that most the CO2 ice is embedded a polar ice matrix that has been thermally processed by the embedded YSO. The amount of thermal processing may be correlated with the luminosity of the YSO. Using the Australia Telescope Compact Array, we imaged the dense gas around a subsample of our sources in the HII complexes N44, N105, N113, and N159 using HCO+ and HCN as dense gas tracers. We find that the molecular material in star forming environments is highly clumpy, with clumps that range from subparsec to ~2 parsecs in size and with masses between 10^2 to 10^4 solar masses. We find that there are varying levels of star formation in the clumps, with the lower-mass clumps tending to be without massive YSOs. These YSO-less clumps could either represent an earlier stage of clump to the more massive YSO-bearing ones or clumps that will never form a massive star. Clumps with massive YSOs at their centers have masses larger than those with massive YSOs at their edges, and we suggest that the difference is evolutionary: edge YSO clumps are more advanced than those with YSOs at their centers. Clumps with YSOs at their edges may have had a significant fraction of their mass disrupted or destroyed by the forming massive star. We find that the strength of the silicate absorption seen in YSO IR spectra feature is well-correlated with the on-source HCO+ and HCN flux densities, such that the strength of the feature is indicative of the embeddedness of the YSO. We estimate that ~40% of the entire spectral sample has strong silicate absorption features, implying that the YSOs are embedded in circumstellar material for about 40% of the time probed in our study.

Thermal and Barometric Constraints on the Intrusive and Unroofing History of the Black Mountains: Implications for Timing, Initial Dip, and Kinematics of Detachment Faulting in the Death-Valley Region, California

Unroofing of the Black Mountains, Death Valley, California, has resulted in the exposure of 1.7 Ga crystalline basement, late Precambrian amphibolite facies metasedimentary rocks, and a Tertiary magmatic complex. The Ar-40/Ar-39 cooling ages, obtained from samples collected across the entire length of the range (>55 km), combined with geobarometric results from synextensional intrusions, provide time-depth constraints on the Miocene intrusive history and extensional unroofing of the Black Mountains. Data from the southeastern Black Mountains and adjacent Greenwater Range suggest unroofing from shallow depths between 9 and 10 Ma. To the northwest in the crystalline core of the range, biotite plateau ages from approximately 13 to 6.8 Ma from rocks making up the Death Valley turtlebacks indicate a midcrustal residence (with temperatures >300-degrees-C) prior to extensional unroofing. Biotite Ar-40/Ar-39 ages from both Precambrian basement and Tertiary plutons reveal a diachronous cooling pattern of decreasing ages toward the northwest, subparallel to the regional extension direction. Diachronous cooling was accompanied by dike intrusion which also decreases in age toward the northwest. The cooling age pattern and geobarometric constraints in crystalline rocks of the Black Mountains suggest denudation of 10-15 km along a northwest directed detachment system, consistent with regional reconstructions of Tertiary extension and with unroofing of a northwest deepening crustal section. Mica cooling ages that deviate from the northwest younging trend are consistent with northwestward transport of rocks initially at shallower crustal levels onto deeper levels along splays of the detachment. The well-known Amargosa chaos and perhaps the Badwater turtleback are examples of this "splaying" process. Considering the current distance of the structurally deepest samples away from moderately to steeply east tilted Tertiary strata in the southeastern Black Mountains, these data indicate an average initial dip of the detachment system of the order of 20-degrees, similar to that determined for detachment faults in west central Arizona and southeastern California. Beginning with an initially listric geometry, a pattern of footwall unroofing accompanied by dike intrusion progress northwestward. This pattern may be explained by a model where migration of footwall flexures occur below a scoop-shaped banging wall block. One consequence of this model is that gently dipping ductile fabrics developed in the middle crust steepen in the upper crust during unloading. This process resolves the low initial dips obtained here with mapping which suggests transport of the upper plate on moderately to steeply dipping surfaces in the middle and upper crust.

Thermal and collisional history of Tishomingo iron meteorite: More evidence for early disruption of differentiated planetesimals

Tishomingo is a chemically and structurally unique iron with 32.5 wt.% Ni that contains 20% residual taenite and 80% martensite plates, which formed on cooling to between -75 and -200 °C, probably the lowest temperature recorded by any meteorite. Our studies using transmission (TEM) and scanning electron microscopy (SEM), X-ray microanalysis (AEM) and electron backscatter diffraction (EBSD) show that martensite plates in Tishomingo formed in a single crystal of taenite and decomposed during reheating forming 10-100 nm taenite particles with ∼50 wt.% Ni, kamacite with ∼4 wt.%Ni, along with martensite or taenite with 32 wt.% Ni. EBSD data and experimental constraints show that Tishomingo was reheated to 320-400 °C for about a year transforming some martensite to kamacite and to taenite particles and some martensite directly to taenite without composition change. Fizzy-textured intergrowths of troilite, kamacite with 2.7 wt.% Ni and 2.6 wt.% Co, and taenite with 56 wt.% Ni and 0.15 wt.% Co formed by localized shock melting. A single impact probably melted the sub-mm sulfides, formed stishovite, and reheated and decomposed the martensite plates. Tishomingo and its near-twin Willow Grove, which has 28 wt.% Ni, differ from IAB-related irons like Santa Catharina and San Cristobal that contain 25-36 wt.% Ni, as they are highly depleted in moderately volatile siderophiles and enriched in Ir and other refractory elements. Tishomingo and Willow Grove therefore resemble IVB irons but are chemically distinct. The absence of cloudy taenite in these two irons shows that they cooled through 250 °C abnormally fast at >0.01 °C/yr. Thus this grouplet, like the IVA and IVB irons, suffered an early impact that disrupted their parent body when it was still hot. Our noble gas data show that Tishomingo was excavated from its parent body about 100 to 200 Myr ago and exposed to cosmic rays as a meteoroid with a radius of ∼50-85 cm.

Thermal profile of a near-adiabatic compression process in a cylindrical tube and establishment of critical control elements for repeatable process control

The compressed gas industry and government agencies worldwide utilize "adiabatic compression" testing for qualifying high-pressure valves, regulators, and other related flow control equipment for gaseous oxygen service. This test methodology is known by various terms including adiabatic compression testing, gaseous fluid impact testing, pneumatic impact testing, and BAM testing as the most common terms. The test methodology will be described in greater detail throughout this document but in summary it consists of pressurizing a test article (valve, regulator, etc.) with gaseous oxygen within 15 to 20 milliseconds (ms). Because the driven gas1 and the driving gas2 are rapidly compressed to the final test pressure at the inlet of the test article, they are rapidly heated by the sudden increase in pressure to sufficient temperatures (thermal energies) to sometimes result in ignition of the nonmetallic materials (seals and seats) used within the test article. In general, the more rapid the compression process the more "adiabatic" the pressure surge is presumed to be and the more like an isentropic process the pressure surge has been argued to simulate. Generally speaking, adiabatic compression is widely considered the most efficient ignition mechanism for directly kindling a nonmetallic material in gaseous oxygen and has been implicated in many fire investigations. Because of the ease of ignition of many nonmetallic materials by this heating mechanism, many industry standards prescribe this testing. However, the results between various laboratories conducting the testing have not always been consistent. Research into the test method indicated that the thermal profile achieved (i.e., temperature/time history of the gas) during adiabatic compression testing as required by the prevailing industry standards has not been fully modeled or empirically verified, although attempts have been made. This research evaluated the following questions: 1) Can the rapid compression process required by the industry standards be thermodynamically and fluid dynamically modeled so that predictions of the thermal profiles be made, 2) Can the thermal profiles produced by the rapid compression process be measured in order to validate the thermodynamic and fluid dynamic models; and, estimate the severity of the test, and, 3) Can controlling parameters be recommended so that new guidelines may be established for the industry standards to resolve inconsistencies between various test laboratories conducting tests according to the present standards?

Thermal decomposition of alkali metal perchlorates

The effect of past mechanical history on the subsequent thermal decomposition kinetics of sodium, potassium, rubidium and caesium perchlorates, has been investigated. At low temperatures the decomposition of all these salts is significantly sensitized by pre-compression. At high temperatures, however, prior compression results in a lowered decomposition rate in the case of potassium, rubidium and caesium perchlorates and in an increase in the thermal reactivity of sodium perchlorate. The high temperature behaviour is shown to be an indirect consequence of the low temperature behaviour. The difference in behaviour between sodium perchlorate and the other alkali metal perchlorates is explained on the basis of the stability of the respective chlorates, formed during the low temperature decomposition. This is substantiated by experiments which show that the addition of sodium chlorate to sodium perchlorate brings about a sensitization while potassium perchlorate admixed with potassium chlorate results in a desensitization at high temperatures.

Effect of process variables on die-billet temperature history in a slow speed hot coining type process

A computer code is developed as a part of an ongoing project on computer aided process modelling of forging operation, to simulate heat transfer in a die-billet system. The code developed on a stage-by-stage technique is based on an Alternating Direction Implicit scheme. The experimentally validated code is used to study the effect of process specifics such as preheat die temperature, machine ascent time, rate of deformation, and dwell time on the thermal characteristics in a batch coining operation where deformation is restricted to surface level only.

Critical stress intensity factors for cracked hollow pipes under transient thermal loads

Finite element analyses of a long hollow cylinder having an axisymmetric circumferential internal edge crack, subjected to convective cooling on the inner surface are performed. The transient thermal stress intensity factor is estimated using a domain version of the J-integral method. The effect of the thickness of the cylinder, crack length, and heat transfer coefficient on the stress intensity factor history are studied. The variations of critical normalized stress intensity factor with crack length-to-thickness ratio for different parameters are presented. The results show that if a small inner surface crack begins to grow, its stress intensity factor will increase with increase in crack length, reach a maximum, and then begin to drop. Based on the results, a fracture-based design methodology for cracked hollow pipes under transient thermal loads is discussed.

On the thermal budget of Pahoehoe lava flows

In this thesis I investigate some aspects of the thermal budget of pahoehoe lava flows. This is done with a combination of general field observations, quantitative modeling, and specific field experiments. The results of this work apply to pahoehoe flows in general, even though the vast bulk of the work has been conducted on the lavas formed by the Pu'u 'O'o - Kupaianaha eruption of Kilauea Volcano on Hawai'i. The field observations rely heavily on discussions with the staff of the United States Geological Survey's Hawaiian Volcano Observatory (HVO), under whom I labored repeatedly in 1991-1993 for a period totaling about 10 months.

The quantitative models I have constructed are based on the physical processes observed by others and myself to be active on pahoehoe lava flows. By building up these models from the basic physical principles involved, this work avoids many of the pitfalls of earlier attempts to fit field observations with "intuitively appropriate" mathematical expressions. Unlike many earlier works, my model results can be analyzed in terms of the interactions between the different physical processes. I constructed models to: (1) describe the initial cooling of small pahoehoe flow lobes and (2) understand the thermal budget of lava tubes.

The field experiments were designed either to validate model results or to constrain key input parameters. In support of the cooling model for pahoehoe flow lobes, attempts were made to measure: (1) the cooling within the flow lobes, (2) the amount of heat transported away from the lava by wind, and (3) the growth of the crust on the lobes. Field data collected by Jones [1992], Hon et al. [1994b], and Denlinger [Keszthelyi and Denlinger, in prep.] were also particularly useful in constraining my cooling model for flow lobes. Most of the field observations I have used to constrain the thermal budget of lava tubes were collected by HVO (geological and geophysical monitoring) and the Jet Propulsion Laboratory (airborne infrared imagery [Realmuto et al., 1992]). I was able to assist HVO for part of their lava tube monitoring program and also to collect helicopterborne and ground-based IR video in collaboration with JPL [Keszthelyi et al., 1993].

The most significant results of this work are (1) the quantitative demonstration that the emplacement of pahoehoe and 'a'a flows are the fundamentally different, (2) confirmation that even the longest lava flows observed in our Solar System could have formed as low effusion rate, tube-fed pahoehoe flows, and (3) the recognition that the atmosphere plays a very important role throughout the cooling of history of pahoehoe lava flows. In addition to answering specific questions about the thermal budget of tube-fed pahoehoe lava flows, this thesis has led to some additional, more general, insights into the emplacement of these lava flows. This general understanding of the tube-fed pahoehoe lava flow as a system has suggested foci for future research in this part of physical volcanology.

The lakes of Michoacan (Mexico): A brief history and alternative point of view

This article is intended to open a discussion about the historical development of lakes Zirahuen, Patzcuaro and Cuitzeo in the state of Michoacan, and the postulated relationships between lake ecology and evolution. Dr Fernando De Buen was the first man dedicated to limnology in Mexico who came to the country in the 1930s. He was adviser at the Estacion Limnologica de Patzcuaro and wrote outstanding papers dealing with Mexican lakes. The lakes of Michoacan probably formed in the late Pliocene or Holocene, and were part of a tributary to the Lerma River, which became isolated by successive volanic barriers to form lake basins. Lake Zirahuen is a warm monomictic waterbody with unique water dynamics amongst the Michoacan lakes. Because it is relatively deep (max depth 40m), seasonal patterns of alternating circulation and thermal stratification develop in the lake, a feature not shared by the other two polymictic shallow lakes, Patzcuaro and Cuitzeo.

Variability in thermal and UV-B energy fluxes through time and their influence on plant diversity and speciation

Throughout Earth's history there have been temporal and spatial variations in the amount of visible and ultraviolet radiation received by ecosystems. This paper examines if temporal changes in these forms of energy receipt could have influenced the tempo and mode of plant diversity and speciation, focusing in particular upon Cenozoic time-scales. Evidence for changing patterns of plant diversity and speciation apparent in various fossil records and molecular phylogenies are considered alongside calculated changes in thermal and solar ultraviolet energy (specifically UV-B) over the past 50 Myr. We suggest that changes in thermal energy influx (amount and variability) affected the tempo of evolution through its influence upon community dynamics (e.g. population size, diversity, turnover, extinctions). It was not only the amount of thermal energy but also variability in its flux that may have influenced these processes, and ultimately the rate of diversification. We suggest that variations in UV-B would have influenced the mode and tempo of speciation through changes to genome stability during intervals of elevated UV-B. We argue, therefore, that although variability in thermal energy and UV-B fluxes through time may lead to the same end-point (changing the rate of diversification), the processes responsible are very different and both need to be considered when linking evolutionary processes to energy flux.

Paragenetic history of the Ordovician Trenton Group carbonates, Southwestern Ontario

Geochemical examination of the rock matrix and cements from core material extracted from four oil wells within southwestern Ontario suggest various stages of diagenetic alteration and preservation of the Trenton Group carbonates. The geochemical compositions of Middle Ordovician (LMC) brachiopods reflect the physicochemical water conditions of the ambient depositional environment. The sediments appear to have been altered in the presence of mixed waters during burial in a relatively open diagenetic microenvironment. Conodont CAl determination suggests that the maturation levels of the Trenton Group carbonates are low and proceeded at temperatures of about 30 - 50°C within the shallow burial environment. The Trenton Group carbonates are characterized by two distinct stages of dolomitization which proceeded at elevated temperatures. Preexisting fracture patterns, and block faulting controlled the initial dolomitization of the precursor carbonate matrix. Dolomitization progressed In the presence of warm fluids (60 75°C) with physicochemical conditions characteristic of a progressively depleted basinal water. The matrix is mostly Idiotopic-S and Idiotopic-E dolomite, with Xenotopic-A dolomite dominating the matrix where fractures occur. The second stage of dolomitization involved hydrothermal basinal fluid(s) with temperatures of about 60 - 70°C. These are the postulated source for the saddle dolomite and blocky calcite cements occurring in pore space and fractures. Rock porosity was partly occluded by Idiotopic-E type dolomite. Late stage saddle dolomite, calcite, anhydrite, pyrite, marcasite and minor sphalerite and celestite cements effectively fill any remaining porosity within specific horizons. Based on cathode luminescence, precipitation of the different diagenetic phases probably proceeded in open diagenetic systems from chemically homogeneous fluids. Ultraviolet fluorescence of 11 the matrix and cements demonstrated that hydrocarbons were present during the earliest formation of saddle dolomite. Oxygen isotope values of -7.6 to -8.5 %0 (PDB), and carbon isotope values of - 0.5 and -3.0 %0 (PDB) from the latest stage dog-tooth calcite cement suggest that meteoric water was introduced into the system during their formation. This is estimated to have occurred at temperatures of about 25 - 40°C. Specific facies associations within the Trenton Group carbonates exhibit good hydrocarbon generating potential based on organic carbon preservation (1-3.5%). Thermal maturation and Lopatin burial-history evaluations suggest that hydrocarbons were generated within the Trenton Group carbonates some time after 300 Ma . Progressively depleted vanadium trends measured from hydrocarbon samples within southwestern Ontario suggests its potential use as a hydrocarbon migration indicator on local (within an oilfield) and on regional scales.

Holocene environmental changes in the lower Thames Valley, London, UK: implications for our understanding of the history of Taxus woodland

A radiocarbon-dated multiproxy palaeoenvironmental record from the Lower Thames Valley at Hornchurch Marshes has provided a reconstruction of the timing and nature of vegetation succession against a background of Holocene climate change, relative sea level movement and human activities. The investigation recorded widespread peat formation between c. 6300 and 3900 cal. yr BP (marine ‘regression’), succeeded by evidence for marine incursion. The multiproxy analyses of these sediments, comprising pollen, Coleoptera, diatoms, and plant and wood macrofossils, have indicated significant changes in both the wetland and dryland environment, including the establishment of Alnus (Alder) carr woodland, and the decline of both Ulmus (Elm; c. 5740 cal. yr BP) and Tilia (Lime; c. 5600 cal. yr BP, and 4160–3710 cal. yr BP). The beetle faunas from the peat also suggest a thermal climate similar to that of the present day. At c. 4900 cal. yr BP, Taxus (L.; Yew) woodland colonised the peatland forming a plant community that has no known modern analogue in the UK. The precise reason, or reasons, for this event remain unclear, although changes in peatland hydrology seem most likely. The growth of Taxus on peatland not only has considerable importance for our knowledge of the vegetation history of southeast England, and NW Europe generally, but also has wider implications for the interpretation of Holocene palaeobotanical records. At c. 3900 cal. yr BP, Taxus declined on the peatland surface during a period of major hydrological change (marine incursion), an event also strongly associated with the decline of dryland woodland taxa, including Tilia and Quercus, and the appearance of anthropogenic indicators.

Quantitative genetic divergence and standing genetic (co)variance in thermal reaction norms along latitude

Although the potential to adapt to warmer climate is constrained by genetic trade-offs, our understanding of how selection and mutation shape genetic (co)variances in thermal reaction norms is poor. Using 71 isofemale lines of the fly Sepsis punctum, originating from northern, central, and southern European climates, we tested for divergence in juvenile development rate across latitude at five experimental temperatures. To investigate effects of evolutionary history in different climates on standing genetic variation in reaction norms, we further compared genetic (co)variances between regions. Flies were reared on either high or low food resources to explore the role of energy acquisition in determining genetic trade-offs between different temperatures. Although the latter had only weak effects on the strength and sign of genetic correlations, genetic architecture differed significantly between climatic regions, implying that evolution of reaction norms proceeds via different trajectories at high latitude versus low latitude in this system. Accordingly, regional genetic architecture was correlated to region-specific differentiation. Moreover, hot development temperatures were associated with low genetic variance and stronger genetic correlations compared to cooler temperatures. We discuss the evolutionary potential of thermal reaction norms in light of their underlying genetic architectures, evolutionary histories, and the materialization of trade-offs in natural environments.

Experimental evolution for generalists and specialists reveals multivariate genetic constraints on thermal reaction norms

Theory predicts the emergence of generalists in variable environments and antagonistic pleiotropy to favour specialists in constant environments, but empirical data seldom support such generalist–specialist trade-offs. We selected for generalists and specialists in the dung fly Sepsis punctum (Diptera: Sepsidae) under conditions that we predicted would reveal antagonistic pleiotropy and multivariate trade-offs underlying thermal reaction norms for juvenile development. We performed replicated laboratory evolution using four treatments: adaptation at a hot (31 °C) or a cold (15 °C) temperature, or under regimes fluctuating between these temperatures, either within or between generations. After 20 generations, we assessed parental effects and genetic responses of thermal reaction norms for three correlated life-history traits: size at maturity, juvenile growth rate and juvenile survival. We find evidence for antagonistic pleiotropy for performance at hot and cold temperatures, and a temperature-mediated trade-off between juvenile survival and size at maturity, suggesting that trade-offs associated with environmental tolerance can arise via intensified evolutionary compromises between genetically correlated traits. However, despite this antagonistic pleiotropy, we found no support for the evolution of increased thermal tolerance breadth at the expense of reduced maximal performance, suggesting low genetic variance in the generalist–specialist dimension.