Gas hydrate associated with mud diapirs in southern Okinawa Trough

Many mud diapirs have been recognized in southern Okinawa Trough by a multi-channel seismic surveying on R/V KEXUE I in 2001. Gas hydrates have been identified, by the seismic reflection characteristics, the velocity analysis and the impedance inversion. Geothermal heat flow around the central of the mud diapir has been determined theoretically by the Bottom Simulating Reflectors (BSRs). Comparing the BSR derived and the measured heat flow values, we infer that the BSR immediately at the top of the mud diapirs indicate the base of the saturated gas hydrate formation zone (BSGHFZ), but not, as we ordinarily know, the base of the gas hydrate stability zone (BGHSZ), which could be explained by the abnormal regional background heat flow and free gas flux associated with mud diapirs. As a result, it helps us to better understand the generation mechanism of the gas hydrates associated with mud diapirs and to predict the gas hydrate potential in the southern Okinawa Trough. (C) 2008 Elsevier Ltd. All rights reserved.

The characteristics of fluid potential in mud diapirs associated with gas hydrates in the Okinawa Trough

Many mud diapirs have been identified in the southern Okinawa Trough from a seismic survey using R/V KEXUE I in 2001. The movement and accumulation of free gas related to mud diapirs are discussed in detail by an analysis of fluid potential which is based upon velocity data. It can be found that free gas moves from the higher fluid potential strata to the lower ones and the gas hydrate comes into being during free gas movement meeting the proper criteria of temperature and pressure. In fact, gas hydrates have been found in the upper layers above the mud diapirs and in host rocks exhibiting other geophysical characteristics. As the result of the formation of the gas hydrate, the free gas bearing strata are enclosed by the gas hydrate bearing strata. Due to the high pressure anomalies of the free gas bearing strata the fluid potential increases noticeably. It can then be concluded that the high fluid potential anomaly on the low fluid potential background may be caused by the presence of the free gas below the gas hydrate bearing strata.

Dating and characterizing late holocene earthquakes using paleomagnetics

In this thesis I apply paleomagnetic techniques to paleoseismological problems. I investigate the use of secular-variation magnetostratigraphy to date prehistoric earthquakes; I identify liquefaction remanent magnetization (LRM), and I quantify coseismic deformation within a fault zone by measuring the rotation of paleomagnetic vectors.

In Chapter 2 I construct a secular-variation reference curve for southern California. For this curve I measure three new well-constrained paleomagnetic directions: two from the Pallett Creek paleoseismological site at A.D. 1397-1480 and A.D. 1465-1495, and one from Panum Crater at A.D. 1325-1365. To these three directions I add the best nine data points from the Sternberg secular-variation curve, five data points from Champion, and one point from the A.D. 1480 eruption of Mt. St. Helens. I derive the error due to the non-dipole field that is added to these data by the geographical correction to southern California. Combining these yields a secular variation curve for southern California covering the period A.D. 670 to 1910, with the best coverage in the range A.D. 1064 to 1505.

In Chapter 3 I apply this curve to a problem in southern California. Two paleoseismological sites in the Salton trough of southern California have sediments deposited by prehistoric Lake Cahuilla. At the Salt Creek site I sampled sediments from three different lakes, and at the Indio site I sampled sediments from four different lakes. Based upon the coinciding paleomagnetic directions I correlate the oldest lake sampled at Salt Creek with the oldest lake sampled at Indio. Furthermore, the penultimate lake at Indio does not appear to be present at Salt Creek. Using the secular variation curve I can assign the lakes at Salt Creek to broad age ranges of A.D. 800 to 1100, A.D. 1100 to 1300, and A.D. 1300 to 1500. This example demonstrates the large uncertainties in the secular variation curve and the need to construct curves from a limited geographical area.

Chapter 4 demonstrates that seismically induced liquefaction can cause resetting of detrital remanent magnetization and acquisition of a liquefaction remanent magnetization (LRM). I sampled three different liquefaction features, a sandbody formed in the Elsinore fault zone, diapirs from sediments of Mono Lake, and a sandblow in these same sediments. In every case the liquefaction features showed stable magnetization despite substantial physical disruption. In addition, in the case of the sandblow and the sandbody, the intensity of the natural remanent magnetization increased by up to an order of magnitude.

In Chapter 5 I apply paleomagnetics to measuring the tectonic rotations in a 52 meter long transect across the San Andreas fault zone at the Pallett Creek paleoseismological site. This site has presented a significant problem because the brittle long-term average slip-rate across the fault is significantly less than the slip-rate from other nearby sites. I find sections adjacent to the fault with tectonic rotations of up to 30°. If interpreted as block rotations, the non-brittle offset was 14.0+2.8, -2.1 meters in the last three earthquakes and 8.5+1.0, -0.9 meters in the last two. Combined with the brittle offset in these events, the last three events all had about 6 meters of total fault offset, even though the intervals between them were markedly different.

In Appendix 1 I present a detailed description of my standard sampling and demagnetization procedure.

In Appendix 2 I present a detailed discussion of the study at Panum Crater that yielded the well-constrained paleomagnetic direction for use in developing secular variation curve in Chapter 2. In addition, from sampling two distinctly different clast types in a block-and-ash flow deposit from Panum Crater, I find that this flow had a complex emplacement and cooling history. Angular, glassy "lithic" blocks were emplaced at temperatures above 600° C. Some of these had cooled nearly completely, whereas others had cooled only to 450° C, when settling in the flow rotated the blocks slightly. The partially cooled blocks then finished cooling without further settling. Highly vesicular, breadcrusted pumiceous clasts had not yet cooled to 600° C at the time of these rotations, because they show a stable, well clustered, unidirectional magnetic vector.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

Influência do Complexo Vulcânico de Abrolhos na evolução tectono-sedimentar da porção centro-sul da bacia do Espírito Santo

A busca por novas acumulações de hidrocarbonetos necessita de esforços exploratórios contínuos, gerando novas possibilidades e modelos geológicos e diminuindo os riscos associados à atividade exploratória. O interesse no entendimento da formação de armadilhas, migração e reservatórios de hidrocarbonetos, associado à halocinese motivou a realização deste trabalho. Apresenta-se como principal objetivo deste trabalho a caracterização e a descrição da evolução halocinética na porção centro-sul da bacia do Espírito Santo. Dados de poços, sísmicos e gravimetria foram utilizados com o intuito de gerar uma interpretação geológica integrada, possibilitando entender à influência do Complexo Vulcânico de Abrolhos (CVA) na evolução tectonossedimentar da área, por meio da técnica de restauração de seção geológica. Na área estudada, ocorreu uma intensa atividade halocinética, já a partir do Albiano, em resposta a distensão causada pela subsidência da bacia e a abertura do Atlântico Sul. Durante o Neocretáceo, cunhas clásticas do Rio Doce adentraram na bacia provocando um novo pulso halocinético, resultando num aumento da taxa de sedimentação nas mini-bacias. Em outras regiões esta progradação causou a migração da camada-mãe de sal para porções distais da bacia, acarretando uma deficiência no suprimento de sal. Isto ocasionou o colapso de alguns diápiros associados a uma quiescência tectônica na área. A principal fase tectônica na área ocorreu no Eoterciário, época em que ocorre a implantação do CVA, formando estilos estruturais característicos de terrenos compressivos, com falhas de empurrão, popups, dobras e gotas de sal. Esta nova configuração tectônica na área mudou os eixos dos principais depocentros, que passaram a ser controlados pelos altos estruturais gerados pela tectônica compressiva, e pelos seus baixos relativos, que passaram a receber os sedimentos sin-tectônicos. (As associações destas características de remobilização tectonossedimentar formou uma nova compartimentação, a saber: a) Zona de translação; b) Zona dobrada e c) Zona de Cavalgamento com falhas de empurrão . Esta nova configuração tectônica tem sua formação diretamente relacionada à implantação do CVA.

Methane hydrate growth velocity in seepage system of Qiongdongnan Basin, South China Sea

The occurrences of diapirs, gas-filled zones and gas plumes in seawater in Qiongdongnan Basin of South China Sea indicate that there may exist seepage system gas-hydrate reservoirs. Assuming there has a methane venting zone of 1500 m in diameter, and the methane flux is 1000 kmol/a, and the temperature of methane hydrate-bearing sediments ranges from 3 degrees C to 20 degrees C, then according to the hydrate film growth theory, by numerical simulation, this paper computes the temperatures and velocities in 0 mbsf, 100 mbsf, 200 mbsf, 425 mbsf over discrete length, and gives the change charts. The results show that the cementation velocity in sediments matrix of methane hydrate is about 0.2 nm/s, and the seepage system will evolve into diffusion system over probably 35000 years. Meanwhile, the methane hydrate growth velocity in leakage system is 20 similar to 40 times faster than in diffusion system.

The long-term stability of a possible aqueous ammonium sulfate ocean inside Titan

We model the thermal evolution of a subsurface ocean of aqueous ammonium sulfate inside Titan using a parameterized convection scheme. The cooling and crystallization of such an ocean depends on its heat flux balance, and is governed by the pressure-dependent melting temperatures at the top and bottom of the ocean. Using recent observations and previous experimental data, we present a nominal model which predicts the thickness of the ocean throughout the evolution of Titan; after 4.5 Ga we expect an aqueous ammonium sulfate ocean 56 km thick, overlain by a thick (176 km) heterogeneous crust of methane clathrate, ice I and ammonium sulfate. Underplating of the crust by ice I will give rise to compositional diapirs that are capable of rising through the crust and providing a mechanism for cryovolcanism at the surface. We have conducted a parameter space survey to account for possible variations in the nominal model, and find that for a wide range of plausible conditions, an ocean of aqueous ammonium sulfate can survive to the present day, which is consistent with the recent observations of Titan's spin state from Cassini radar data [Lorenz, R.D., Stiles, B.W., Kirk, R.L., Allison, M.D., del Marmo, P.P., Iess, L., Lunine, J.I., Ostro, S.J., Hensley, S., 2008. Science 319, 1649–1651].

Magnetic fabric and zircon U-Pb geochronology of the Itaoca pluton: Implications for the Brasiliano deformation of the southern Ribeira Belt (SE Brazil)

The Itaoca pluton consists of porphyritic monzogranite that intruded the upper crust into low-grade metasedimentary rocks of the Apiai Dornain (Ribeira Belt). Anisotropy of magnetic susceptibility and zircon U-Pb (Shrimp) geochronology were combined to determine pluton emplacement mechanisms and its chronology relative to the collision structures of the Paranapiacaba (Brasiliano II) orogenic system. Magnetic susceptibility ranges between 4 and 38 x 10(-3) SI, and thermomagnetic measurements indicate multidomain magnetite is the main carrier of anisotropy. The pluton shows an ""onion-skin"" structure roughly elongated to the northeast with its hinge zone including kilometer-wide roof-pendants. Magnetic lineations are variable in orientation in consistency with the dominant oblate symmetry of the magnetic fabric. A distinct NE-trending point-maxima, however, indicates the mean lineation is parallel to the stretching direction of the transpressive deformation that affected the regional host rocks. Prismatic zircon from the monzogranite, both in the core and in the finely-zoned margins, yielded an age of 623 +/- 10 Ma. These results suggest the magmatic fabric recorded the earlier strain increments of the regional shear deformation. It may correspond to the transition from continental arc to collision tectonics of the southern Ribeira Belt. (C) 2008 Elsevier Ltd. All rights reserved.

Sediment deformation in Miocene and post-Miocene strata, Northeastern Brazil: Evidence for paleoseismicity in a passive margin

Stratigraphic intervals characterized by varied and complex styles of soft-sediment deformation structures are well preserved in Miocene and Late Pleistocene to Holocene deposits of a sedimentary basin located in Northeastern Brazil. The Miocene strata, represented by the Barreiras Formation, record only brittle structures, including numerous faults and fractures with straight and high angle-dipping planes that are often filled with sands derived from overlying beds. Folds consisting of broad anticlines and synclines are also present in this unit. The late Pleistocene to Holocene deposits, named Post-Barreiras Sediments, contain an indurated sandy package with a large variety of ductile and brittle deformation structures (i.e., massive sandstones with isolated sand fragments and breccias, undulatory strata, sand dykes and diapirs, sinks and bowls, pebbly pockets, plunged sediment mixtures, fitted sand masses, cone-shaped cracks, fault grading and sedimentary enclaves). These features, confined to sharp-based stratigraphic horizons that progressively grade downward into undisturbed deposits, are related to seismic shocks of high surface-wave magnitude (i.e., Ms>5 or 6). Amalgamated seismites suggest that previously formed seismites were affected by subsequent seismic-wave propagation. Seismic waves caused by activity along one, or most likely, several tectonic structures would have propagated throughout the depositional environment, producing laterally extensive seismites. The close proximity to earthquake epicenters would have promoted pervasive re-sedimentation due to pore overpressure, resulting high volumes of massive sandstones and breccia. The similarity between deposits with correlatable strata from many other areas along the Brazilian coast allows raise the hypothesis that the seismic episodes might have affected sedimentation patterns in a large (i.e., extension of several hundreds of kilometers) geographic area. Thus, the modern seismicity recorded along Northeastern Brazil was recurrent during the Quaternary and, perhaps, also in the Pliocene. The estimated high magnitude of the seismic events and the great regional extent of the affected area demonstrate that the Brazilian coast experienced tectonic stress through the last geological episodes of its evolution, which would have favored sediment accumulation and penecontemporaneous re-sedimentation. This geological context is unexpected in a passive margin, inducing to revisit the debate on how active is a passive margin. (C) 2010 Elsevier B.V. All rights reserved.

Caracterização tectonoestratigráfica da porção adjacente ao epicentro do sismo São Vicente na bacia de Santos

Pós-graduação em Geologia Regional - IGCE

Evolução halocinética da região centro-norte da Bacia de Santos, Brasil

Pós-graduação em Geociências e Meio Ambiente - IGCE

2D and 3D numerical modelling of multilayer detachment folding and salt tectonics

Numerical modelling was performed to study the dynamics of multilayer detachment folding and salt tectonics. In the case of multilayer detachment folding, analytically derived diagrams show several folding modes, half of which are applicable to crustal scale folding. 3D numerical simulations are in agreement with 2D predictions, yet fold interactions result in complex fold patterns. Pre-existing salt diapirs change folding patterns as they localize the initial deformation. If diapir spacing is much smaller than the dominant folding wavelength, diapirs appear in fold synclines or limbs.rnNumerical models of 3D down-building diapirism show that sedimentation rate controls whether diapirs will form and influences the overall patterns of diapirism. Numerical codes were used to retrodeform modelled salt diapirs. Reverse modelling can retrieve the initial geometries of a 2D Rayleigh-Taylor instability with non-linear rheologies. Although intermediate geometries of down-built diapirs are retrieved, forward and reverse modelling solutions deviate. rnFinally, the dynamics of fold-and-thrusts belts formed over a tilted viscous detachment is studied and it is demonstrated that mechanical stratigraphy has an impact on the deformation style, switching from thrust- to folding-dominated. The basal angle of the detachment controls the deformation sequence of the fold-and-thrust belt and results are consistent with critical wedge theory.rn

Complex geomorphologic assemblage of terrains in association with the banded terrain in Hellas basin, Mars

Hellas basin acts as a major sink for the southern highlands of Mars and is likely to have recorded several episodes of sedimentation and erosion. The north-western part of the basin displays a potentially unique Amazonian landscape domain in the deepest part of Hellas, called “banded terrain”, which is a deposit characterized by an alternation of narrow band shapes and inter-bands displaying a sinuous and relatively smooth surface texture suggesting a viscous flow origin. Here we use high-resolution (HiRISE and CTX) images to assess the geomorphological interaction of the banded terrain with the surrounding geomorphologic domains in the NW interior of Hellas to gain a better understanding of the geological evolution of the region as a whole. Our analysis reveals that the banded terrain is associated with six geomorphologic domains: a central plateau named Alpheus Colles, plain deposits (P1 and P2), reticulate (RT1 and RT2) and honeycomb terrains. Based on the analysis of the geomorphology of these domains and their cross-cutting relationships, we show that no widespread deposition post-dates the formation of the banded terrain, which implies that this domain is the youngest and latest deposit of the interior of Hellas. Therefore, the level of geologic activity in the NW Hellas during the Amazonian appears to have been relatively low and restricted to modification of the landscape through mechanical weathering, aeolian and periglacial processes. Thermophysical data and cross-cutting relationships support hypotheses of modification of the honeycomb terrain via vertical rise of diapirs such as ice diapirism, and the formation of the plain deposits through deposition and remobilization of an ice-rich mantle deposit. Finally, the observed gradual transition between honeycomb and banded terrain suggests that the banded terrain may have covered a larger area of the NW interior of Hellas in the past than previously thought. This has implications on the understanding of the evolution of the deepest part of Hellas.

Chemistry of orthopyroxenes, clinopyroxenes, and spinels of ODP Hole 103-637A peridotites

Mineral compositions of the plagioclase-bearing ultramafic tectonites dredged and cored seaward of the continental slope of the Galicia margin (Leg 103, Site 637) were compared to mineral compositions from onshore low-pressure ultramafic bodies (southeastern Ronda, western Pyrenees, and Lizard Point), on the basis of standardized (30-s counting time) probe analyses. The comparison was extended to some plagioclase-free harzburgites related to ophiolites (Santa Elena in Costa Rica, north Oman, and the Humboldt body in New Caledonia) on the basis of new analytical data and data from the literature. The behavior of Cr, Na, Al, Mg, Fe, Ni, and Ti in olivine, pyroxenes, and spinel was examined in order to distinguish between the effects of partial melting and mineral facies change, from the spinel to plagioclase stability fields. The peridotite from the Galicia margin appears slightly depleted in major incompatible elements and experienced a minor partial melting. However, it experienced large scale but heterogeneous recrystallization in the plagioclase field. These features are very similar to those observed in Ronda, whereas in the western Pyrenees the minerals exemplify a very minor partial-melting event (or none at all) and have retained compositions corresponding to those of the relatively high-pressure Seiland sub facies. The minerals from the Lizard Point peridotite have characteristics (low Mg/(Mg + Fe) ratio; high Cr/(Cr + Al) ratio in spinel) more related to cumulate from a differentiated tholeiitic melt than related to ophiolitic tectonite. Diffusion profiles of Al and Cr across pyroxenes and spinel show that recrystallization features occurred at different speeds or temperatures in the different bodies. The pyroxenes from Ronda would have experienced recrystallization about 14 times faster than the peridotite from the Galicia margin. The western Pyrenean lherzolites also experienced rapid recrystallization; nevertheless, because they are of a different mineral facies, the data are not directly comparable to that from Ronda and Galicia. The harzburgite at Santa Elena as well as a xenolith from alkali basalt exemplify rapid cooling characterized by very weak re-equilibration. Recrystallization speed is related to emplacement speed in the present geological environment. The slow-rising Galicia margin peridotite was emplaced by thinning of the lithospheric subcontinental mantle near an incipient mid-oceanic ridge. The fast-rising peridotites from Ronda and the western Pyrenees were hot diapirs emplaced from the asthenosphere along transcurrent faults, possibly related to the opening of the Atlantic Ocean.