Travel time tomographic studies of seismic structures around subducted lithospheric slabs

The nature of the subducted lithospheric slab is investigated seismologically by tomographic inversions of ISC residual travel times. The slab, in which nearly all deep earthquakes occur, is fast in the seismic images because it is much cooler than the ambient mantle. High resolution three-dimensional P and S wave models in the NW Pacific are obtained using regional data, while inversion for the SW Pacific slabs includes teleseismic arrivals. Resolution and noise estimations show the models are generally well-resolved.

The slab anomalies in these models, as inferred from the seismicity, are generally coherent in the upper mantle and become contorted and decrease in amplitude with depth. Fast slabs are surrounded by slow regions shallower than 350 km depth. Slab fingering, including segmentation and spreading, is indicated near the bottom of the upper mantle. The fast anomalies associated with the Japan, Izu-Bonin, Mariana and Kermadec subduction zones tend to flatten to sub-horizontal at depth, while downward spreading may occur under parts of the Mariana and Kuril arcs. The Tonga slab appears to end around 550 km depth, but is underlain by a fast band at 750-1000 km depths.

The NW Pacific model combined with the Clayton-Comer mantle model predicts many observed residual sphere patterns. The predictions indicate that the near-source anomalies affect the residual spheres less than the teleseismic contributions. The teleseismic contributions may be removed either by using a mantle model, or using teleseismic station averages of residuals from only regional events. The slab-like fast bands in the corrected residual spheres are are consistent with seismicity trends under the Mariana Tzu-Bonin and Japan trenches, but are inconsistent for the Kuril events.

The comparison of the tomographic models with earthquake focal mechanisms shows that deep compression axes and fast velocity slab anomalies are in consistent alignment, even when the slab is contorted or flattened. Abnormal stress patterns are seen at major junctions of the arcs. The depth boundary between tension and compression in the central parts of these arcs appears to depend on the dip and topology of the slab.

Dynamic rupture simulation integrated with earthquake observations

Dynamic rupture simulations are unique in their contributions to the study of earthquake physics. The current rapid development of dynamic rupture simulations poses several new questions: Do the simulations reflect the real world? Do the simulations have predictive power? Which one should we believe when the simulations disagree? This thesis illustrates how integration with observations can help address these questions and reduce the effects of non-uniqueness of both dynamic rupture simulations and kinematic inversion problems. Dynamic rupture simulations with observational constraints can effectively identify non-physical features inferred from observations. Moreover, the integrative technique can also provide more physical insights into the mechanisms of earthquakes. This thesis demonstrates two examples of such kinds of integration: dynamic rupture simulations of the M_w 9.0 2011 Tohoku-Oki earthquake and of earthquake ruptures in damaged fault zones:

(1) We develop simulations of the Tohoku-Oki earthquake based on a variety of observations and minimum assumptions of model parameters. The simulations provide realistic estimations of stress drop and fracture energy of the region and explain the physical mechanisms of high-frequency radiation in the deep region. We also find that the overridding subduction wedge contributes significantly to the up-dip rupture propagation and large final slip in the shallow region. Such findings are also applicable to other megathrust earthquakes.

(2) Damaged fault zones are usually found around natural faults, but their effects on earthquake ruptures have been largely unknown. We simulate earthquake ruptures in damaged fault zones with material properties constrained by seismic and geological observations. We show that reflected waves in fault zones are effective at generating pulse-like ruptures and head waves tend to accelerate and decelerate rupture speeds. These mechanisms are robust in natural fault zones with large attenuation and off-fault plasticity. Moreover, earthquakes in damaged fault zones can propagate at super-Rayleigh speeds that are unstable in homogeneous media. Supershear transitions in fault zones do not require large fault stresses. In the end, we present observations in the Big Bear region, where variability of rupture speeds of small earthquakes correlates with the laterally variable materials in a damaged fault zone.

Complexity of rupture propagation in large earthquakes in relation to tectonic environment

Complexity in the earthquake rupture process can result from many factors. This study investigates the origin of such complexity by examining several recent, large earthquakes in detail. In each case the local tectonic environment plays an important role in understanding the source of the complexity.

Several large shallow earthquakes (M_s > 7.0) along the Middle American Trench have similarities and differences between them that may lead to a better understanding of fracture and subduction processes. They are predominantly thrust events consistent with the known subduction of the Cocos plate beneath N. America. Two events occurring along this subduction zone close to triple junctions show considerable complexity. This may be attributable to a more heterogeneous stress environment in these regions and as such has implications for other subduction zone boundaries.

An event which looks complex but is actually rather simple is the 1978 Bermuda earthquake (M_s ~ 6). It is located predominantly in the mantle. Its mechanism is one of pure thrust faulting with a strike N 20°W and dip 42°NE. Its apparent complexity is caused by local crustal structure. This is an important event in terms of understanding and estimating seismic hazard on the eastern seaboard of N. America.

A study of several large strike-slip continental earthquakes identifies characteristics which are common to them and may be useful in determining what to expect from the next great earthquake on the San Andreas fault. The events are the 1976 Guatemala earthquake on the Motagua fault and two events on the Anatolian fault in Turkey (the 1967, Mudurnu Valley and 1976, E. Turkey events). An attempt to model the complex P-waveforms of these events results in good synthetic fits for the Guatemala and Mudurnu Valley events. However, the E. Turkey event proves to be too complex as it may have associated thrust or normal faulting. Several individual sources occurring at intervals of between 5 and 20 seconds characterize the Guatemala and Mudurnu Valley events. The maximum size of an individual source appears to be bounded at about 5 x 10²⁶ dyne-cm. A detailed source study including directivity is performed on the Guatemala event. The source time history of the Mudurnu Valley event illustrates its significance in modeling strong ground motion in the near field. The complex source time series of the 1967 event produces amplitudes greater by a factor of 2.5 than a uniform model scaled to the same size for a station 20 km from the fault.

Three large and important earthquakes demonstrate an important type of complexity --- multiple-fault complexity. The first, the 1976 Philippine earthquake, an oblique thrust event, represents the first seismological evidence for a northeast dipping subduction zone beneath the island of Mindanao. A large event, following the mainshock by 12 hours, occurred outside the aftershock area and apparently resulted from motion on a subsidiary fault since the event had a strike-slip mechanism.

An aftershock of the great 1960 Chilean earthquake on June 6, 1960, proved to be an interesting discovery. It appears to be a large strike-slip event at the main rupture's southern boundary. It most likely occurred on the landward extension of the Chile Rise transform fault, in the subducting plate. The results for this event suggest that a small event triggered a series of slow events; the duration of the whole sequence being longer than 1 hour. This is indeed a "slow earthquake".

Perhaps one of the most complex of events is the recent Tangshan, China event. It began as a large strike-slip event. Within several seconds of the mainshock it may have triggered thrust faulting to the south of the epicenter. There is no doubt, however, that it triggered a large oblique normal event to the northeast, 15 hours after the mainshock. This event certainly contributed to the great loss of life-sustained as a result of the Tangshan earthquake sequence.

What has been learned from these studies has been applied to predict what one might expect from the next great earthquake on the San Andreas. The expectation from this study is that such an event would be a large complex event, not unlike, but perhaps larger than, the Guatemala or Mudurnu Valley events. That is to say, it will most likely consist of a series of individual events in sequence. It is also quite possible that the event could trigger associated faulting on neighboring fault systems such as those occurring in the Transverse Ranges. This has important bearing on the earthquake hazard estimation for the region.

Imaging the earth with ambient noise and earthquakes

In this thesis, I develop the velocity and structure models for the Los Angeles Basin and Southern Peru. The ultimate goal is to better understand the geological processes involved in the basin and subduction zone dynamics. The results are obtained from seismic interferometry using ambient noise and receiver functions using earthquake- generated waves. Some unusual signals specific to the local structures are also studied. The main findings are summarized as follows:

(1) Los Angeles Basin

The shear wave velocities range from 0.5 to 3.0 km/s in the sediments, with lateral gradients at the Newport-Inglewood, Compton-Los Alamitos, and Whittier Faults. The basin is a maximum of 8 km deep along the profile, and the Moho rises to a depth of 17 km under the basin. The basin has a stretch factor of 2.6 in the center decreasing to 1.3 at the edges, and is in approximate isostatic equilibrium. This "high-density" (~1 km spacing) "short-duration" (~1.5 month) experiment may serve as a prototype experiment that will allow basins to be covered by this type of low-cost survey.

(2) Peruvian subduction zone

Two prominent mid-crust structures are revealed in the 70 km thick crust under the Central Andes: a low-velocity zone interpreted as partially molten rocks beneath the Western Cordillera – Altiplano Plateau, and the underthrusting Brazilian Shield beneath the Eastern Cordillera. The low-velocity zone is oblique to the present trench, and possibly indicates the location of the volcanic arcs formed during the steepening of the Oligocene flat slab beneath the Altiplano Plateau.

The Nazca slab changes from normal dipping (~25 degrees) subduction in the southeast to flat subduction in the northwest of the study area. In the flat subduction regime, the slab subducts to ~100 km depth and then remains flat for ~300 km distance before it resumes a normal dipping geometry. The flat part closely follows the topography of the continental Moho above, indicating a strong suction force between the slab and the overriding plate. A high-velocity mantle wedge exists above the western half of the flat slab, which indicates the lack of melting and thus explains the cessation of the volcanism above. The velocity turns to normal values before the slab steepens again, indicating possible resumption of dehydration and ecologitization.

(3) Some unusual signals

Strong higher-mode Rayleigh waves due to the basin structure are observed in the periods less than 5 s. The particle motions provide a good test for distinguishing between the fundamental and higher mode. The precursor and coda waves relative to the interstation Rayleigh waves are observed, and modeled with a strong scatterer located in the active volcanic area in Southern Peru. In contrast with the usual receiver function analysis, multiples are extensively involved in this thesis. In the LA Basin, a good image is only from PpPs multiples, while in Peru, PpPp multiples contribute significantly to the final results.

Studies of noble gases in meteorites and in the earth

The isotopic and elemental abundances of noble gases in the solar system are investigated, using simple mixing models and mass-spectrometric measurements of the noble gases in meteorites and terrestrial rocks and minerals.

Primordial neon is modeled by two isotopically distinct components from the interstellar gas and dust. Neon from the gas dominates solar neon, which contains about ten times more ²⁰Ne than ²²Ne. Neon from the dust is represented in meteorites by neon-E, with ²⁰Ne/²²Ne less than 0.6. Isotopic variations in meteorites require neon from both dust and gas to be present. Mixing dust and gas without neon loss generates linear correlation lines on three-isotope and composition-concentration diagrams. A model for solar wind implantation predicts small deviations from linear mixing, due to preferential sputtering of the lighter neon isotopes.

Neon in meteorites consists of galactic cosmic ray spallation neon and at least two primordial components, neon-E and neon-S. Neon was measured in several meteorites to investigate these end- members. Cosmogenic neon produced from sodium is found to be strongly enriched in ²²Ne. Neon measurements on sodium-rich samples must be interpreted with care so not to confuse this source of ²²Ne with neon-E, which is also rich in ²²Ne.

Neon data for the carbonaceous chondrite Mokoia show that the end member composition of neon-Si in meteorites is ²⁰Ne/²²Ne = 13.7, the same as the present solar wind. The solar wind composition evidently has remained constant since before the compaction of Mokoia.

Ca, Al-rich inclusions from the Allende meteorite were examined for correlation between neon-E and oxygen or magnesium isotopic anomalies. ²²Ne and ³⁶Ar enrichments found in some inclusions are attributed to cosmic- ray-induced reactions on Na and Cl, not to a primordial component. Neon-E is not detectably enriched in Allende.

Measurements were made to determine the noble gas contents of various terrestrial rocks and minerals, and to investigate the cycling of noble gases between different terrestrial reservoirs. Beryl crystals contain a characteristic suite of magmatic gases including nucleogenic ²¹Ne and ²²Ne from (α,n) reactions, radiogenic ⁴⁰Ar, and fissiogenic ^131-136Xe from the decay of K and U in the continental crust. Significant concentrations of atmospheric noble gases are also present in beryl.

Both juvenile and atmospheric noble gases are found in rocks from the Skaergaard intrusion. The ratio ⁴⁰Ar/³⁶Ar (corrected for in situ decay of ⁴⁰K) correlates with δ¹⁸O in plagioclase. Atmospheric argon has been introduced into samples that have experienced oxygen-isotope exchange with circulating meteoric hydrothermal fluids. Unexchanged samples contain juvenile argon with ⁴⁰Ar/³⁶Ar greater than 6000 that was trapped from the Skaergaard magma.

Juvenile and atmospheric gases have been measured in the glassy rims of mid-ocean ridge (MOR) pillow basalts. Evidence is presented that three samples contain excess radiogenic ¹²⁹Xe and fission xenon, in addition to the excess radiogenic ⁴⁰Ar found in all samples. These juvenile gases are being outgassed from the upper-mantle source region of the MOR magma. No isotopic evidence has been found here for juvenile primordial noble gases accompanying the juvenile radiogenic gases in the MOR glasses. Large argon isotopic variations in a single specimen provide a clear indication of the late-stage addition of atmospheric argon, probably from seawater.

The Skaergaard data demonstrate that atmospheric noble gases dissolved in ground water can be transferred into crustal rocks. Subduction of oceanic crust altered by seawater can transport atmospheric noble gases into the upper mantle. A substantial portion of the noble gases in mantle derived rocks may represent subducted gases, not a primordial component as is often assumed.

Dynamical models for the Earth's geoid

The Earth's largest geoid anomalies occur at the lowest spherical harmonic degrees, or longest wavelengths, and are primarily the result of mantle convection. Thermal density contrasts due to convection are partially compensated by boundary deformations due to viscous flow whose effects must be included in order to obtain a dynamically consistent model for the geoid. These deformations occur rapidly with respect to the timescale for convection, and we have analytically calculated geoid response kernels for steady-state, viscous, incompressible, self-gravitating, layered Earth models which include the deformation of boundaries due to internal loads. Both the sign and magnitude of geoid anomalies depend strongly upon the viscosity structure of the mantle as well as the possible presence of chemical layering.

Correlations of various global geophysical data sets with the observed geoid can be used to construct theoretical geoid models which constrain the dynamics of mantle convection. Surface features such as topography and plate velocities are not obviously related to the low-degree geoid, with the exception of subduction zones which are characterized by geoid highs (degrees 4-9). Recent models for seismic heterogeneity in the mantle provide additional constraints, and much of the low-degree (2-3) geoid can be attributed to seismically inferred density anomalies in the lower mantle. The Earth's largest geoid highs are underlain by low density material in the lower mantle, thus requiring compensating deformations of the Earth's surface. A dynamical model for whole mantle convection with a low viscosity upper mantle can explain these observations and successfully predicts more than 80% of the observed geoid variance.

Temperature variations associated with density anomalies in the man tie cause lateral viscosity variations whose effects are not included in the analytical models. However, perturbation theory and numerical tests show that broad-scale lateral viscosity variations are much less important than radial variations; in this respect, geoid models, which depend upon steady-state surface deformations, may provide more reliable constraints on mantle structure than inferences from transient phenomena such as postglacial rebound. Stronger, smaller-scale viscosity variations associated with mantle plumes and subducting slabs may be more important. On the basis of numerical modelling of low viscosity plumes, we conclude that the global association of geoid highs (after slab effects are removed) with hotspots and, perhaps, mantle plumes, is the result of hot, upwelling material in the lower mantle; this conclusion does not depend strongly upon plume rheology. The global distribution of hotspots and the dominant, low-degree geoid highs may correspond to a dominant mode of convection stabilized by the ancient Pangean continental assemblage.

Petrogênese das rochas máficas alcalinas do litoral entre São Sebastião (SP) e Parati (RJ)

Lamprófiros e diabásios alcalinos afloram no litoral dos estados de São Paulo e Rio de Janeiro e integram o Enxame de Diques da Serra do Mar (EDSM). Essas rochas ocorrem sob a forma de diques e intrudem o Orógeno Ribeira, de idade Neoproterozóica/Cambro-Ordoviciana, inserindo-se no contexto geodinâmico de abertura do Oceano Atlântico Sul durante o Cretáceo Superior. Essas intrusões são subverticais e orientam-se preferencialmente a NE-SW, seguindo a estruturação das rochas encaixantes. Os lamprófiros são classificados como monchiquitos e camptonitos e exibem, respectivamente, textura hipocristalina e holocristalina. Apresentam também textura panidiomórfica, fenocristais de clinopiroxênio e olivina, imersos em matriz formada essencialmente por esses mesmos minerais, além de biotita, kaersutita e minerais opacos. O camptonito apresenta ainda plagioclásio na matriz. Os diabásios alcalinos são hipocristalinos a holocristalinos, equigranulares a inequigranulares, com fenocristais de olivina e/ou clinopiroxênio e/ou plagioclásio, em uma matriz composta essencialmente por esses minerais. As rochas estudadas caracterizam séries alcalinas miaskíticas, com os lamprófiros sendo tanto sódicos, potássicos e ultrapotássicos e os diabásios alcalinos como predominantemente sódicos. Modelagens petrogenéticas envolvendo possíveis processos evolutivos mostram que é improvável que os lamprófiros sejam cogenéticos por processos evolutivos envolvendo tanto cristalização fracionada, com ou sem assimilação concomitante, quanto hibridização. O mesmo ocorre para os diabásios alcalinos. A discriminação de fontes mantélicas foi feita com base nos teores de elementos traços de amostras representativas de líquidos parentais e indica que esse magmatismo alcalino está relacionado a fontes lherzolíticas com fusão parcial na zona de estabilidade do espinélio, isto é, a poucas profundidades. Os dados litogeoquímicos e isotópicos do sistema Sr-Nd das rochas estudadas sugerem tanto o envolvimento de fontes férteis, associadas ao manto sublitosférico, quanto de fontes enriquecidas, relacionadas ao manto litosférico subcontinental. Modelagens de mistura binária revelam que a petrogênese dos lamprófiros e diabásios alcalinos envolveu uma grande participação de um componente fértil misturado com contribuições menores de um componente enriquecido. Idades TDM (760-557 Ma) obtidas sugerem remobilização do manto litosférico no Neoproterozóico, talvez relacionadas à subducção da Placa São Francisco preteritamente à colisão do Orógeno Ribeira. Altas razões CaO/Al2O3 para os líquidos lamprofíricos menos evoluídos, altos teores de Zr, correlações negativas Zr/Hf e Ti/Eu e associação com carbonatitos indicam condições metassomáticas de alto CO2/H2O. Em escala local, modelos geodinâmicos baseados na astenosfera não isotérmica parecem mais aplicáveis. No entanto, modelos geodinâmicos baseados na astenosfera isotérmica (com o envolvimento de plumas) parecem mais indicados num contexto regional, considerando-se outras províncias alcalinas contemporâneas e correlatas.

Petrogênese e geocronologia das intrusões alcalinas de Morro Redondo, Mendanha e Morro de São João: caracterização do magmatismo alcalino no Estado do Rio de Janeiro e implicações geodinâmicas

Os modelos para a formação de plútons alcalinos da Província Alcalina do Sudeste Brasileiro ou Alinhamento Poços de Caldas-Cabo Frio associam a gênese destas rochas a grandes reativações ou a passagem de uma pluma mantélica, registrada pelo traço de um hot spot. O objetivo desta tese é, apresentar novos dados e interpretações para contribuir com a melhor elucidação e discussão destes modelos. Os estudos incluem mapeamento, petrografia, litogeoquímica, geoquímica isotópica de Sr, Nd e Pb e datação 40Ar/39Ar. As intrusões selecionadas correspondem ao Morro Redondo, Mendanha e Morro de São João, no Rio de Janeiro, localizados em posições distintas no alinhamento Poços de Caldas-Cabo Frio. A intrusão alcalina do Morro Redondo é composta majoritariamente de nefelina sienitos e sienitos com nefelina, com rara ocorrência de rochas máficas e é caracterizada por uma suíte alcalina sódica insaturada em sílica, de caráter metaluminosa a peralcalina. Esta intrusão foi datada em aproximadamente 74 Ma (idade-platô 40Ar/39Ar). A intrusão alcalina do Mendanha é composta por diversos tipos de rochas sieníticas, além de brechas e estruturas subvulcânicas, como rochas piroclásticas e diques e caracteriza-se por ser uma suíte alcalina sódica saturada em sílica, de caráter metaluminosa, diferente do que ocorre no Marapicu, este subsaturado em sílica. Esta intrusão apresentou duas idades-platô 40Ar/39Ar distintas de magmatismo: 64 Ma para as rochas do Mendanha e 54 Ma em dique de lamprófiro, registrando magmatismo policíclico. O Morro do Marapicu foi datado em aproximadamente 80 Ma. Já a intrusão alcalina do Morro de São João possui uma ampla variedade de litotipos saturados a subsaturados em sílica, tais como sienitos, álcali-sienitos e monzossienitos (alguns portadores de pseudoleucita), com variedades melanocráticas, tais como malignitos e fergustios. Estas rochas definem suas distintas suítes alcalinas subsaturadas em sílica: Uma de composição sódica e outra potássica. Há também uma suíte alcalina saturada em sílica, definida por gabros alcalinos e shonkinitos. A petrogênese destas intrusões corresponde ao modelo de cristalização fracionada, com assimilação de rochas encaixantes (AFC) como indicado pela alta variabilidade de razões isotópicas de estrôncio. No Morro de São João é sugerido o modelo de mistura magmática. Estas intrusões foram geradas a partir de magmas mantélicos enriquecidos, possivelmente associados à antiga zona de subducção relacionada ao orógeno Ribeira. Em razão das novas idades obtidas, o modelo de hot spot proposto fica prejudicado, visto que o Marapicu é de idade mais antiga das intrusões analisadas, o que era esperado para o Morro Redondo. Alguns modelos projetam plumas mantélicas com aproximadamente 1000 km de diâmetro, o que poderia explicar o Mendanha ser contemporâneo ao Morro de São João. As assinaturas isotópicas obtidas para as intrusões não se associam à assinatura isotópica de Trindade e, caso o modelo de plumas mantélicas seja o correto, a pluma que teria maior semelhança de assinatura isotópica é a pluma de Tristão da Cunha.

Proveniência sedimentar do grupo São Fidélis, terreno oriental da Faixa Ribeira com base em dados U-Pb

O Domínio Costeiro integra o Terreno Oriental, no segmento central da Faixa Ribeira e abriga rochas ortoderivadas com afinidade de arcos magmáticos (Complexo Rio Negro, ca. 790-605 Ma). Os ortognaisses deste complexo possuem clara assinatura para ambientes de zonas de subducção, encaixados em rochas metassedimentares de alto grau, integrantes do Grupo São Fidélis. O conjunto acima descrito é ainda intrudido por várias de rochas granitóides sin a tardi- colisionais, relacionadas às várias etapas de desenvolvimento da Orogenia Brasiliana neste setor do orógeno (ca. 605-480 Ma). Idades U-Pb (LA-ICP-MS) em zircões detríticos de rochas quartzíticas do Grupo São Fidélis indicam um amplo espectro com modas significativas no Mesoproterozóico e Paleoproterozoico, além de zircões do Neoproterozóico e do Arqueano. Sinteticamente os resultados obtidos foram: a) Idades concordantes Arqueanas com ca. 2,85, 2,84 e 2,70 Ga; b) zircões Paleoproterozóicos (ca. 2,3 a 1,7 Ga), com máxima concentração em torno de ca. 2,2 Ga, representando a segunda maior moda; c) Idades Mesoproterozóicas (ca. 1,3 -1,1 Ga) com idades de espectro dominantes, com moda em ca. 1,5 Ga; d) Zircões Neoproterozóicos com idades de ca. 0,95-90 Ga e 0,86-0,61 Ga. Em vários grãos detríticos observou-se sobrecrescimento metamórfico em ca. 602-570 Ma. Dados U-Pb (LA-ICP-MS) obtidos para zircões para Ortognaisse Rio Grande e o Biotita Ortognaisse, intrudidos na unidade basal do Grupo São Fidélis, apresentam idades em ca. 620 Ma e são equivalentes ao período pré-colisional de geração de rochas do arco magmático Rio Negro. Combinando estas idades com os núcleos de zircões detríticos mais jovens, com assinatura do Arco Rio Negro em ca. 613 Ma, pode-se definir o intervalo máximo de sedimentação da unidade superior do Grupo São Fidélis no Neoproterozóico. Cristais de monazitas selecionadas para análise U-Pb (ID-TIMS) apresentam relações com os principais episódios tectono-metamórficos da Faixa Ribeira. Dois cristais de uma amostra quartzítica e dois do ortognaisse Rio Grande alinham-se em uma discórdia que gerou idade de 603 Ma, referente ao metamorfismo progressivo descrito na literatura, durante a Orogenia Brasiliana. Enquanto a idade concordante obtida em 535 Ma, adquirida em uma amostra quartzítica, é correspontente ao último metamorfismo colisional da Faixa Ribeira.

Change in Pacific Northwest coastal ecosystems. Proceedings of the Pacific Northwest Coastal Ecosystems Regional Study Workshop, August 73-14, 1996, Troutdale, Oregon

Over the past one hundred and fifty years, the landscape and ecosystems of the Pacific Northwest coastal region, already subject to many variable natural forces, have been profoundly affected by human activities. In virtually every coastal watershed from the Strait of Juan de Fuca to Cape Mendocino, settlement, exploitation and development of resou?-ces have altered natural ecosystems. Vast, complex forests that once covered the region have been largely replaced by tree plantations or converted to non-forest conditions. Narrow coastal valleys, once filled with wetlands and braided streams that tempered storm runoff and provided salmon habitat, were drained, filled, or have otherwise been altered to create land for agriculture and other uses. Tideflats and saltmarshes in both large and small estuaries were filled for industrial, commercial, and other urban uses. Many estuaries, including that of the Columbia River, have been channeled, deepened, and jettied to provide for safe, reliable navigation. The prodigious rainfall in the region, once buffered by dense vegetation and complex river and stream habitat, now surges down sirfiplified stream channels laden with increased burdens of sediment and debris. Although these and many other changes have occurred incrementally over time and in widely separated areas, their sum can now be seen to have significantly affected the natural productivity of the region and, as a consequence, changed the economic structure of its human communities. This activity has taken place in a region already shaped by many interacting and dynamic natural forces. Large-scale ocean circulation patterns, which vary over long time periods, determine the strength and location of currents along the coast, and thus affect conditions in the nearshore ocean and estuaries throughout the region. Periodic seasonal differences in the weather and ocean act on shorter time scales; winters are typically wet with storms from the southwest while summers tend to be dry with winds from the northwest. Some phenomena are episodic, such as El Nifio events, which alter weather, marine habitats, and the distribution and survival of marine organisms. Other oceanic and atmospheric changes operate more slowly; over time scales of decades, centuries, and longer. Episodic geologic events also punctuate the region, such as volcanic eruptions that discharge widespread blankets of ash, frequent minor earthquakes, and major subduction zone earthquakes each 300 to 500 years that release accumulated tectonic strain, dropping stretches of ocean shoreline, inundating estuaries and coastal valleys, and triggering landslides that reshape stream profiles. While these many natural processes have altered, sometimes dramatically, the Pacific Northwest coastal region, these same processes have formed productive marine and coastal ecosystems, and many of the species in these systems have adapted to the variable environmental conditions of the region to ensure their long-term survival.

Sri Lanka: is it a mid-plate platelet?

Two observations suggest the possibility that Sri Lanka is acting as a small-mid-plate platelet moving very slowly within and relative to the larger Indian plate. First, sediments of the Bengal Deep-Sea Fan off the SSE continental margin are folded and uplifted in a manner similar to the deformation from front of accretionary prisms where thick sediment columns are passing into subduction zones. And second, subsidence rates in the area of presumed spreading or continuing stretching of continental crust, the Cauvery-Palk Strait-Gulf of Mannar Basin, have not decreased during the Cenozoic as would be predicted by an aborted rift or aulacogen model, but instead appear to have accelerated during the Neogene. Information available on other phenomena which re predicted by the model is at the present time inadequate for evaluation.

Offshore geology of Pakistan and non-living resource prospects

Offshore geology of Pakistan is characterized by active and passive continental margins. These continental margins show very unique features such as an active Makran subduction zone in the west and the Indus delta and a submarine fan in the east. The geology of these features of Pakistan EEZ is inadequately known. This is a major obstacle in exploring mineral resources. Detailed study of coastal and shelf geology is needed for better understanding of the geology of the area and comprehensive evolution of its non-living resource potential. An understanding of the geological events including tectonic movements, sedimentation processes and geochemical history that comprise the geological history is very important to help in identification and estimation of resources. In Pakistan EEZ applying the current technology and undertaking research work to understand the seafloor features and mineral deposits associated with it will be very fruitful.

Mechanical and hydraulic properties of Nankai accretionary prism sediments: Effect of stress path

We have conducted triaxial deformation experiments along different loading paths on prism sediments from the Nankai Trough. Different load paths of isotropic loading, uniaxial strain loading, triaxial compression (at constant confining pressure, Pc), undrained Pc reduction, drained Pc reduction, and triaxial unloading at constant Pc, were used to understand the evolution of mechanical and hydraulic properties under complicated stress states and loading histories in accretionary subduction zones. Five deformation experiments were conducted on three sediment core samples for the Nankai prism, specifically from older accreted sediments at the forearc basin, underthrust slope sediments beneath the megasplay fault, and overthrust Upper Shikoku Basin sediments along the frontal thrust. Yield envelopes for each sample were constructed based on the stress paths of Pc-reduction using the modified Cam-clay model, and in situ stress states of the prism were constrained using the results from the other load paths and accounting for horizontal stress. Results suggest that the sediments in the vicinity of the megasplay fault and frontal thrust are highly overconsolidated, and thus likely to deform brittle rather than ductile. The porosity of sediments decreases as the yield envelope expands, while the reduction in permeability mainly depends on the effective mean stress before yield, and the differential stress after yield. An improved understanding of sediment yield strength and hydromechanical properties along different load paths is necessary to treat accurately the coupling of deformation and fluid flow in accretionary subduction zones. © 2012 American Geophysical Union All Rights Reserved.

Discovery of Late Cretaceous garnet two-pyroxene granulite in the southern Lhasa terrane, Tibet and its tectonic significances

The Lhasa terrane, located between the Bangonghu-Nujiang suture zone and the Indus-Yalung Tsangpo suture zone in the southern Tibetan Plateau, was considered previously as a Precambrian continental block. Mesozoic and Cenozoic tectonic evolution of the Lhasa terrane is closely related to the subduction of the Tethys ocean and the collision between the Indian and European continents; so it is one of the keys to reveal the formation and evolution of the Tibetan plateau. The garnet two-pyroxene granulite which was found at the Nyingtri rock group of the southeastern Lhasa terrene consists of garnet, clinopyroxene, orthopyroxene, labradorite, Ti-rich amphibolite and biotite, with a chemical composition of mafic rock. The metamorphic conditions were estimated to be at T = 747 similar to 834 degrees C and P = 0.90 similar to 1.35GPa, suggesting a formation depth of 45km. The zircon U-Pb dating for the garnet amphibolite and marble associated with the granulite give a metamorphic age of 85 similar to 90Ma. This granulite-facies metamorphic event together with a contemporaneous magmatism demonstrated that the southern Lhasa terrane has undergone an Andean-type orogeny at Late Mesozoic time.

Trace element composition of peridotites from the southern Mariana forearc: Insights into the geochemical effects of serpentinization and/or seafloor weathering

Peridotites from the southern Mariana forearc were sampled on the landward trench slope of the Izu-Bonin-Mariana (IBM) subduction zone by dredging. These mantle wedge peridotites underwent hydration by fluid derived from a dehydrated descending slab, and later interacted with seawater after emplacement at or near the seafloor. This study investigates how these two different rock-fluid interaction processes influenced trace element distribution in the southern Mariana forearc peridotites. We measured trace element concentrations of peridotites from the southern Mariana forearc. The southern Mariana forearc peridotites are characterized by a distinct seawater-like REE pattern with an obvious negative Ce anomaly, and La shows good correlation with other REEs (except Ce). In addition, there is a great enrichment of U, Pb, Sr and Li elements, which show a distinct positive anomaly relative to adjacent elements in the multi-element diagram. For the seawater-like REE pattern, we infer that REEs are mainly influenced by seawater during peridotite-seawater interactions after their emplacement at or near the seafloor, by serpentinization or by marine weathering. Furthermore, the anomalous behavior of Ce, compared with other rare earth elements in these samples, may indicate that they have undergone reactions involving Ce (IV) when the peridotites interacted with seawater. Positive U, Pb, Sr and Li anomalies are inferred to be related to seawater and/or fluids released during dehydration of the subducting slab.