Rates of tectonic and magmatic processes in the North Cascades continental magmatic arc

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, 2004.Includes bibliographical references.Continental magmatic arcs are among the most dynamic. geologic systems, and documentation of the magmatic, thermal, and tectonic evolution of arcs is essential for understanding the processes of magma generation, ascent and crustal growth. The primary goal of this research is to determine rates of tectonic and magmatic processes in the mid to deep crustal levels of the crystalline core of the Cretaceous North Cascades arc. This region was selected for study because it preserves a -10-40 km depth-section through the arc, which allows an assessment of magmatic and structural processes over a range of crustal levels. The relatively young age of the arc (ca. 100-45 Ma) and the inherent high-precision of U-Pb zircon dates permit absolute uncertainties of <100 ky. Meta-supracrustal rocks of the Cascades core record some of the highest pressures obtained in the North American Cordillera. The timing of deposition and metamorphism of the 9-12 kbar Swakane Gneiss constrain tectonic burial models and the timescales of large crustal displacements within an arc setting. These models involve rapid burial (-7 mm/yr) of a fore- or back-arc basin from ca. 73-68 Ma. Nd isotopic signatures of all meta-clastic terranes of the Cascades core reflect mixing of arc- and craton-derived sediment, and the Swakane Gneiss has the most isotopically-evolved signature of these terranes. Nd isotopic signatures of plutons that intrude the core lack evidence of melting of this isotopically-evolved unit. 40Ar/39Ar and U-Pb thermochronologic data define regional cooling patterns that suggest mid- to Late Cretaceous exhumation coincident with contraction and crustal thickening at the deepest levels of the core, followed by Early Tertiary extension.(cont.) High-precision U-Pb geochronology reveals internal complexities inherent in the construction of an intrusive magmatic system. The Mount Stuart batholith was constructed over a ca. 5.6 Myr time period with four punctuated intervals of magma emplacement, whereas the Tenpeak intrusion was emplaced in a more continuous process over ca. 2.7 Myr time period. U-Pb zircon dates from two elongate intrusions, the Seven- Fingered Jack and Entiat suites, suggest that they were constructed from multiple magmatic sheets that were partially homogenized at the level of emplacement.by Jennifer E. Piontek Matzel.Ph.D

Pioneering Robotic Liver Surgery in Germany: First Experiences with Liver Malignancies

Background Minimally invasive liver surgery is growing worldwide with obvious benefits for the treated patients. These procedures maybe improved by robotic techniques, which add several innovative features. In Germany, we were the first surgical department implementing robotic assisted minimally invasive liver resections. Material and methods Between June 2013 and March 2015, we performed robotic based minimally invasive liver resections in nine patients with malignant liver disease. Five off these patients suffered from primary and four from secondary liver malignancies. We retrospectively analyzed the perioperative variables of these patients and the oncological follow up. Results Mean age of the patients was 63 years (range 45–71). One patient suffered from intrahepatic cholangiocellular, four from hepatocellular carcinoma, and four patients from colorectal liver metastases. In six patients, left lateral liver resection, in two cases single segment resection, and in one case minimally invasive guided liver ablation were performed. Five patients underwent previous abdominal surgery. Mean operation time was 312 min (range 115–458 min). Mean weight of the liver specimens was 182 g (range 62–260 g) and mean estimated blood loss was 251 ml (range 10–650 ml). The mean tumor size was 4.4 cm (range 3.5–5.5 cm). In all cases, R0 status was confirmed with a mean margin of 0.6 cm (range 0.1–1.5 cm). One patient developed small bowel fistula on postoperative day 5, which could be treated conservatively. No patient died. Mean hospital stay of the patients was 6 days (range 3–10 days). During a mean follow up of 12 months (range 1–21 months), two patients developed tumor recurrence. Conclusion Robotic-based liver surgery is feasible in patients with primary and secondary liver malignancies. To achieve perioperative parameters comparable to open settings, the learning curve must be passed. Minor liver resections are good candidates to start this technique. But the huge benefits of robotic-based liver resections should be expected in extended procedures beyond minor liver resections with the currently available technology

Partitioned Waveform Inversion Applied to Eurasia and Northern Africa

This report summarizes the data analysis achieved during Heather Bedle's eleven-week Technical Scholar internship at Lawrence Livermore National Labs during the early summer 2006. The work completed during this internship resulted in constraints on the crustal and upper mantle S-velocity structure in Northern Africa, the Mediterranean, the Middle East, and Europe, through the fitting of regional waveform data. This data extends current raypath coverage and will be included in a joint inversion along with data from surface wave group velocity measurements, S and P teleseismic arrival time data, and receiver function data to create an improved velocity model of the upper mantle in this region. The tectonic structure of the North African/Mediterranean/Europe/Middle Eastern study region is extremely heterogeneous. This region consists of, among others, stable cratons and platforms such as the West Africa Craton, and Baltica in Northern Europe; oceanic subduction zones throughout the Mediterranean Sea where the African and Eurasian plate collide; regions of continental collision as the Arabian Plate moves northward into the Turkish Plate; and rifting in the Red Sea, separating the Arabian and Nubian shields. With such diverse tectonic structures, many of the waveforms were difficult to fit. This is not unexpected as the waveforms are fit using an averaged structure. In many cases the raypaths encounter several tectonic features, complicating the waveform, and making it hard for the software to converge on a 1D average structure. Overall, the quality of the waveform data was average, with roughly 30% of the waveforms being discarded due to excessive noise that interfered with the frequency ranges of interest. An inversion for the 3D S-velocity structure of this region was also performed following the methodology of Partitioned Waveform Inversion (Nolet, 1990; Van der Lee and Nolet, 1997). The addition of the newly fit waveforms drastically extends the range of the model. The model now extends as far east in Africa to cover Chad and Niger, and reaches south to cover Zambia. The model is also stretched eastward to cover the eastern half of India, and northward to cover the southern portion of Scandinavia

MODEL-BASED HYDROACOUSTIC BLOCKAGE ASSESSMENT AND DEVELOPMENT OF AN EXPLOSIVE SOURCE DATABASE

We are continuing the development of the Hydroacoustic Blockage Assessment Tool (HABAT) which is designed for use by analysts to predict which hydroacoustic monitoring stations can be used in discrimination analysis for any particular event. The research involves two approaches (1) model-based assessment of blockage, and (2) ground-truth data-based assessment of blockage. The tool presents the analyst with a map of the world, and plots raypath blockages from stations to sources. The analyst inputs source locations and blockage criteria, and the tool returns a list of blockage status from all source locations to all hydroacoustic stations. We are currently using the tool in an assessment of blockage criteria for simple direct-path arrivals. Hydroacoustic data, predominantly from earthquake sources, are read in and assessed for blockage at all available stations. Several measures are taken. First, can the event be observed at a station above background noise? Second, can we establish backazimuth from the station to the source. Third, how large is the decibel drop at one station relative to other stations. These observational results are then compared with model estimates to identify the best set of blockage criteria and used to create a set of blockage maps for each station. The model-based estimates are currently limited by the coarse bathymetry of existing databases and by the limitations inherent in the raytrace method. In collaboration with BBN Inc., the Hydroacoustic Coverage Assessment Model (HydroCAM) that generates the blockage files that serve as input to HABAT, is being extended to include high-resolution bathymetry databases in key areas that increase model-based blockage assessment reliability. An important aspect of this capability is to eventually include reflected T-phases where they reliably occur and to identify the associated reflectors. To assess how well any given hydroacoustic discriminant works in separating earthquake and in-water explosion populations it is necessary to have both a database of reference earthquake events and of reference in-water explosive events. Although reference earthquake events are readily available, explosive reference events are not. Consequently, building an in-water explosion reference database requires the compilation of events from many sources spanning a long period of time. We have developed a database of small implosive and explosive reference events from the 2003 Indian Ocean Cruise data. These events were recorded at some or all of the IMS Indian Ocean hydroacoustic stations: Diego Garcia, Cape Leeuwin, and Crozet Island. We have also reviewed many historical large in-water explosions and identified five that have adequate source information and can be positively associated to the hydrophone recordings. The five events are: Cannekin, Longshot, CHASE-3, CHASE-5, and IITRI-1. Of these, the first two are nuclear tests on land but near water. The latter three are in-water conventional explosive events with yields from ten to hundreds of tons TNT equivalent. The objective of this research is to enhance discrimination capabilities for events located in the world's oceans. Two research and development efforts are needed to achieve this: (1) improvement in discrimination algorithms and their joint statistical application to events, and (2) development of an automated and accurate blockage prediction capability that will identify all stations and phases (direct and reflected) from a given event that will have adequate signal to be used in a discrimination analysis. The strategy for improving blockage prediction in the world's oceans is to improve model-based prediction of blockage and to develop a ground-truth database of reference events to assess blockage. Currently, research is focused on the development of a blockage assessment software tool. The tool is envisioned to develop into a sophisticated and unifying package that optimally and automatically assesses both model and data based blockage predictions in all ocean basins, for all NDC stations, and accounting for reflected phases (Pulli et al., 2000). Currently, we have focused our efforts on the Diego Garcia, Cape Leeuwin and Crozet Island hydroacoustic stations in the Indian Ocean

Oxygen Isotope Measurements of a Rare Murchison Type A CAI and Its Rim

Ca-, Al-rich inclusions (CAIs) from CV chondrites commonly show oxygen isotope heterogeneity among different mineral phases within individual inclusions reflecting the complex history of CAIs in both the solar nebula and/or parent bodies. The degree of isotopic exchange is typically mineral-specific, yielding O-16-rich spinel, hibonite and pyroxene and O-16-depleted melilite and anorthite. Recent work demonstrated large and systematic variations in oxygen isotope composition within the margin and Wark-Lovering rim of an Allende Type A CAI. These variations suggest that some CV CAIs formed from several oxygen reservoirs and may reflect transport between distinct regions of the solar nebula or varying gas composition near the proto-Sun. Oxygen isotope compositions of CAIs from other, less-altered chondrites show less intra-CAI variability and 16O-rich compositions. The record of intra-CAI oxygen isotope variability in CM chondrites, which commonly show evidence for low-temperature aqueous alteration, is less clear, in part because the most common CAIs found in CM chondrites are mineralogically simple (hibonite +/- spinel or spinel +/- pyroxene) and are composed of minerals less susceptible to O-isotopic exchange. No measurements of the oxygen isotope compositions of rims on CAIs in CM chondrites have been reported. Here, we present oxygen isotope data from a rare, Type A CAI from the Murchison meteorite, MUM-1. The data were collected from melilite, hibonite, perovskite and spinel in a traverse into the interior of the CAI and from pyroxene, melilite, anorthite, and spinel in the Wark-Lovering rim. Our objectives were to (1) document any evidence for intra-CAI oxygen isotope variability; (2) determine the isotopic composition of the rim minerals and compare their composition(s) to the CAI interior; and (3) compare the MUM-1 data to oxygen isotope zoning profiles measured from CAIs in other chondrites

Source and Propagation Characteristics of Explosive and Other Seismic Sources

Understanding of the source and propagation characteristics of seismic events of different types including earthquakes, explosions and mining-induced events is essential for successful discrimination of nuclear explosions. We are compiling a data set of mining related seismic events in east Eurasia. Natural earthquake data in the same region are also collected for comparison study between mining related events and earthquakes. The ground-truth data set will provide a unique and valuable resource for monitoring research. We will utilize the data set to investigate the source and propagation characteristics of seismic sources of different types including mine blasts, tremors, collapses and earthquakes. We will use various seismological techniques including spectral analysis, and waveform modeling to conduct the investigation. The research will improve our understanding of the S-wave excitation and propagation characteristics of chemical explosions and other source types

Coordinated Oxygen Isotopic and Petrologic Studies of CAIS Record Varying Composition of Protosolar

Ca-, Al-rich inclusions (CAIs) record the O-isotope composition of Solar nebular gas from which they grew [1]. High spatial resolution O-isotope measurements afforded by ion microprobe analysis across the rims and margin of CAIs reveal systematic variations in (Delta)O-17 and suggest formation from a diversity of nebular environments [2-4]. This heterogeneity has been explained by isotopic mixing between the O-16-rich Solar reservoir [6] and a second O-16-poor reservoir (probably nebular gas) with a "planetary-like" isotopic composition [e.g., 1, 6-7], but the mechanism and location(s) where these events occur within the protoplanetary disk remain uncertain. The orientation of large and systematic variations in (Delta)O-17 reported by [3] for a compact Type A CAI from the Efremovka reduced CV3 chondrite differs dramatically from reports by [4] of a similar CAI, A37 from the Allende oxidized CV3 chondrite. Both studies conclude that CAIs were exposed to distinct, nebular O-isotope reservoirs, implying the transfer of CAIs among different settings within the protoplanetary disk [4]. To test this hypothesis further and the extent of intra-CAI O-isotopic variation, a pristine compact Type A CAI, Ef-1 from Efremovka, and a Type B2 CAI, TS4 from Allende were studied. Our new results are equally intriguing because, collectively, O-isotopic zoning patterns in the CAIs indicate a progressive and cyclic record. The results imply that CAIs were commonly exposed to multiple environments of distinct gas during their formation. Numerical models help constrain conditions and duration of these events

Heterogenous Oxygen Isotopic Composition of a Complex Wark-Lovering Rim and the Margin of a Refractory Inclusion from Leoville

Wark-Lovering (WL) rims [1] surrounding many refractory inclusions represent marker events in the early evolution of the Solar System in which many inclusions were exposed to changes in pressure [2], temperature [3], and isotopic reservoirs [4-7]. The effects of these events can be complex, not only producing mineralogical variability of WL rims [2], but also leading to mineralogical [8-10] and isotopic [7, 11, 12] changes within inclusion interiors. Extreme oxygen isotopic heterogeneity measured in CAIs has been explained by mixing between distinct oxygen gas reservoirs in the nebula [13]. Some WL rims contain relatively simple mineral layering and/or are isotopically homogeneous [14, 15]. As part of a larger effort to document and understand the modifications observed in some CAIs, an inclusion (L6) with a complex WL rim from Leoville, a member of the reduced CV3 subgroup was studied. Initial study of the textures and mineral chemistry was presented by [16]. Here we present NanoSIMS oxygen isotopic measurements to complement these petrologic observations

Does Oxygen Isotopic Heterogeneity in Refractory Inclusions and Their Wark-Lovering Rims Record Nebular Repressing?

Large systematic variations in O-isotopic compositions found within individual mineral layers of rims surrounding Ca-, Al-rich inclusions (CAIs) and at the margins of some CAIs imply formation from distinct environments [e.g., 1-3]. The O-isotope compositions of many CAIs preserve a record of the Solar nebula gas believed to initially be O-16-rich (delta O-17 less than or equal to -25%0) [4-5]. Data from a recent study of the compact Type A Allende CAI, A37, preserve a diffusion profile in the outermost ~70 micrometers of the inclusion and show greater than 25%0 variations in delta O-17 within its ~100 micrometer-thick Wark-Lovering rim (WL-rim) [3]. This and comparable heterogeneity measured in several other CAIs have been explained by isotopic mixing between the O-16-rich Solar reservoir and a second O-16-poor reservoir (probably nebular gas) with a planetary-like isotopic composition, e.g., [1,2,3,6]. However, there is mineralogical and isotopic evidence from the interiors of CAIs, in particular those from Allende, for parent body alteration. At issue is how to distinguish the record of secondary reprocessing in the nebula from that which occurred on the parent body. We have undertaken the task to study a range of CAI types with varying mineralogies, in part, to address this problem

REGIONAL SEISMIC CHEMICAL AND NUCLEAR EXPLOSION DISCRIMINATION: WESTERN U.S. EXAMPLES

We continue exploring methodologies to improve regional explosion discrimination using the western U.S. as a natural laboratory. The western U.S. has abundant natural seismicity, historic nuclear explosion data, and widespread mine blasts, making it a good testing ground to study the performance of regional explosion discrimination techniques. We have assembled and measured a large set of these events to systematically explore how to best optimize discrimination performance. Nuclear explosions can be discriminated from a background of earthquakes using regional phase (Pn, Pg, Sn, Lg) amplitude measures such as high frequency P/S ratios. The discrimination performance is improved if the amplitudes can be corrected for source size and path length effects. We show good results are achieved using earthquakes alone to calibrate for these effects with the MDAC technique (Walter and Taylor, 2001). We show significant further improvement is then possible by combining multiple MDAC amplitude ratios using an optimized weighting technique such as Linear Discriminant Analysis (LDA). However this requires data or models for both earthquakes and explosions. In many areas of the world regional distance nuclear explosion data is lacking, but mine blast data is available. Mine explosions are often designed to fracture and/or move rock, giving them different frequency and amplitude behavior than contained chemical shots, which seismically look like nuclear tests. Here we explore discrimination performance differences between explosion types, the possible disparity in the optimization parameters that would be chosen if only chemical explosions were available and the corresponding effect of that disparity on nuclear explosion discrimination. Even after correcting for average path and site effects, regional phase ratios contain a large amount of scatter. This scatter appears to be due to variations in source properties such as depth, focal mechanism, stress drop, in the near source material properties (including emplacement conditions in the case of explosions) and in variations from the average path and site correction. Here we look at several kinds of averaging as a means to try and reduce variance in earthquake and explosion populations and better understand the factors going into a minimum variance level as a function of epicenter (see Anderson ee et al. this volume). We focus on the performance of P/S ratios over the frequency range from 1 to 16 Hz finding some improvements in discrimination as frequency increases. We also explore averaging and optimally combining P/S ratios in multiple frequency bands as a means to reduce variance. Similarly we explore the effects of azimuthally averaging both regional amplitudes and amplitude ratios over multiple stations to reduce variance. Finally we look at optimal performance as a function of magnitude and path length, as these put limits the availability of good high frequency discrimination measures