Gas hydrate systems at Hydrate Ridge offshore Oregon inferred from molecular and isotopic properties of hydrate-bound and void gases

We report and discuss molecular and isotopic properties of hydrate-bound gases from 55 samples and void gases from 494 samples collected during Ocean Drilling Program (ODP) Leg 204 at Hydrate Ridge offshore Oregon. Gas hydrates appear to crystallize in sediments from two end-member gas sources (deep allochthonous and in situ) as mixtures of different proportions. In an area of high gas flux at the Southern Summit of the ridge (Sites 1248-1250), shallow (0-40 m below the seafloor [mbsf]) gas hydrates are composed of mainly allochthonous mixed microbial and thermogenic methane and a small portion of thermogenic C₂₊ gases, which migrated vertically and laterally from as deep as 2- to 2.5-km depths. In contrast, deep (50-105 mbsf) gas hydrates at the Southern Summit (Sites 1248 and 1250) and on the flanks of the ridge (Sites 1244-1247) crystallize mainly from microbial methane and ethane generated dominantly in situ. A small contribution of allochthonous gas may also be present at sites where geologic and tectonic settings favor focused vertical gas migration from greater depth (e.g., Sites 1244 and 1245). Non-hydrocarbon gases such as CO₂ and H₂S are not abundant in sampled hydrates. The new gas geochemical data are inconsistent with earlier models suggesting that seafloor gas hydrates at Hydrate Ridge formed from gas derived from decomposition of deeper and older gas hydrates. Gas hydrate formation at the Southern Summit is explained by a model in which gas migrated from deep sediments, and perhaps was trapped by a gas hydrate seal at the base of the gas hydrate stability zone (GHSZ). Free gas migrated into the GHSZ when the overpressure in gas column exceeded sealing capacity of overlaying sediments, and precipitated as gas hydrate mainly within shallow sediments. The mushroom-like 3D shape of gas hydrate accumulation at the summit is possibly defined by the gas diffusion aureole surrounding the main migration conduit, the decrease of gas solubility in shallow sediment, and refocusing of gas by carbonate and gas hydrate seals near the seafloor to the crest of the local anticline structure.     

Contents to Volume 8

Volume Contents

Contents to Volume 8     

On-line estimation of nonlinear physical systems

Recursive algorithms for estimating states of nonlinear physical systems are presented. Orthogonality properties are rediscovered and the associated polynomials are used to linearize state and observation models of the underlying random processes. This requires some key hypotheses regarding the structure of these processes, which may then take account of a wide range of applications. The latter include streamflow forecasting, flood estimation, environmental protection, earthquake engineering, and mine planning. The proposed estimation algorithm may be compared favorably to Taylor series-type filters, nonlinear filters which approximate the probability density by Edgeworth or Gram-Charlier series, as well as to conventional statistical linearization-type estimators. Moreover, the method has several advantages over nonrecursive estimators like disjunctive kriging. To link theory with practice, some numerical results for a simulated system are presented, in which responses from the proposed and extended Kalman algorithms are compared.     

Vertical vibration of a rigid strip footing on viscoelastic half-space

This paper presents an analytical method for modeling the dynamic response of a rigid strip footing subjected to vertical-only loads. The footing is assumed to rest on the surface of a viscoelastic half-space; therefore, effects of hysteretic soil damping on the impedance of the foundation and the generated ground vibrations are considered in the solution. To solve the mixed boundary value problem, we use the Fourier transform to cast a pair of dual integral equations providing contact stresses, which are solved by means of Jacobi orthogonal polynomials. The resulting soil and footing displacements and stresses are obtained by means of the Fourier inverse transform. The solution provides more realistic estimates of footing impedance, compared to existing solutions for elastic soil, as well as of the attenuation of ground vibrations with distance of the footing. The latter is important for the estimation of machine vibration effects on nearby structures and installations.     

Comparative urbanization in Ghana and Kenya in time and space

There are few inter-African country urban analyses because of the continent’s enormous size and socioeconomic diversity, language barriers, and wide variations in national and regional urban research capacity. Nevertheless, comparative urban studies are critical in understanding contemporary African urbanization. In this comparative spatial and temporal analysis of Ghana and Kenya’s urbanization, we find that both countries are urbanizing rapidly and are faced with many common urban problems. Moreover, Ghana is more urbanized than Kenya and has a larger indigenous urban imprint and a more widely dispersed urban pattern. Besides their physiographic and population conditions, we trace these countries’ convergent and divergent urban trends to their shared but unique experiences of colonialism, nationalism and globalization.     

Maritime tsunami evacuation guidelines for the Pacific Northwest coast of Oregon

Recent tsunamis affecting the West Coast of the USA have resulted in significant damage to ports and harbors, as well as to recreational and commercial vessels attempting to escape the tsunami. With the completion of tsunami inundation simulations for a distant tsunami originating from the Aleutian Islands and a locally generated tsunami on the Cascadia subduction zone (CSZ), the State of Oregon is now able to provide guidance on the magnitudes and directions of the simulated currents for the Oregon coast and shelf region. Our analyses indicate that first wave arrivals for an Aleutian Island event would take place on the north coast,?~?3 h 40 min after the start of the earthquake,?~?20 min later on the southern Oregon coast. The simulations demonstrated significant along-coast variability in both the tsunamis water levels and currents, caused by localized bathymetric effects (e.g., submarine banks and reefs). A locally generated CSZ event would reach the open coast within 7–13 min; maximum inundation occurs at?~?30–40 min. As the tsunami current velocities increase, the potential for damage in ports and harbors correspondingly increases, while also affecting a vessels ability to maintain control out on the ocean. Scientific consensus suggests that tsunami currents?<?1.54 m/s are unlikely to impact maritime safety in ports and harbors. No such guidance is available for boats operating on the ocean, though studies undertaken in Japan suggest that velocities in the region of 1–2 m/s may be damaging to boats. In addition to the effects of currents, there is the added potential for wave amplification of locally generated wind waves interacting with opposing tsunami currents in the offshore. Our analyses explore potential wave amplification effects for a range of generic sea states, ultimately producing a nomogram of wave amplification for a range of wave and opposing current conditions. These data will be useful for US Coast Guard and Port authorities as they evaluate maritime tsunami evacuation options for the Oregon coast. Finally, we identify three regions of hazard (high, moderate, and low) across the Oregon shelf, which can be used to help guide final designation of tsunami maritime evacuation zones for the coast.     

Thermal metamorphism of the Arunachal Himalaya, India: Raman thermometry and thermochronological constraints on the tectono-thermal evolution

Determination of the peak thermal condition is vital in order to understand tectono-thermal evolution of the Himalayan belt. The Lesser Himalayan Sequence (LHS) in the Western Arunachal Pradesh, being rich in carbonaceous material (CM), facilitates the determination of peak metamorphic temperature based on Raman spectroscopy of carbonaceous material (RSCM). In this study, we have used RSCM method of Beyssac et al. (J Metamorph Geol 20:859–871, 2002a) and Rahl et al. (Earth Planet Sci Lett 240:339–354, 2005) to estimate the thermal history of LHS and Siwalik foreland from the western Arunachal Pradesh. The study indicates that the temperature of 700–800 °C in the Greater Himalayan Sequence (GHS) decreases to 650–700 °C in the main central thrust zone (MCTZ) and decreases further to <200 °C in the Mio-Pliocene sequence of Siwaliks. The work demonstrates greater reliability of Rahl et al.’s (Earth Planet Sci Lett 240:339–354, 2005) RSCM method for temperatures >600 and <340 °C. We show that the higher and lower zones of Bomdila Gneiss (BG) experienced temperature of ~600 °C and exhumed at different stages along the Bomdila Thrust (BT) and Upper Main Boundary Thrust (U.MBT). Pyrolysis analysis of the CM together with the Fission Track ages from upper Siwaliks corroborates the RSCM thermometry estimate of ~240 °C. The results indicate that the Permian sequence north of Lower MBT was deposited at greater depths (>12 km) than the upper Siwalik sediments to its south at depths <8 km before they were exhumed. The ⁴⁰Ar/³⁹Ar ages suggest that the upper zones of Se La evolved ~13–15 Ma. The middle zone exhumed at ~11 Ma and lower zone close to ~8 Ma indicating erosional unroofing of the MCT sheet. The footwall of MCTZ cooled between 6 and 8 Ma. Analyses of P–T path imply that LHS between MCT and U.MBT zone falls within the kyanite stability field with near isobaric condition. At higher structural level, the temperatures increase gradually with P–T conditions in the sillimanite stability field. The near isothermal (700–800 °C) condition in the GHS, isobaric condition in the MCTZ together with T–t path evidence of GHS that experienced relatively longer duration of near peak temperatures and rapid cooling towards MCTZ, compares the evolution of GHS and inverted metamorphic gradient closely to channel flow predictions.     

Iterative selective spatial variance reduction of MYD11A2 LST data

The aim of this research effort is to develop a method that will allow to map and evaluate thermal anomalies in SW USA from the MYD11A2 night land surface temperature (LST) imagery being available for the year 2014, that present higher spatial (1 km) and temporal (46 images per year) resolution than the MYD11C3 LST data (12 images per year at 5.6 km spatial resolution). The fact that is MYD11A2 LST imagery is projected to a rectangular grid did not affect the X, Y and elevation (H) spatial decorrelation stretch. Principal component analysis and linear regression models isolated and removed the X, Y, H (spatial) dependent variance included in the data while metrics devised verified the selective spatial variance reduction. The reconstructed 46 LST images represent the amount the LST deviates from the X, Y and H predicted for the year 2014. The thematic information content of the reconstructed LST images is verified by cluster analysis and mapped the spatial extend and the temporal variability of thermal anomalies within the study area. The positive thermal anomaly clusters are spatially arranged mainly west of Sierra Nevada in Great Basin Section where extensional tectonics create a series of titled elongated mountain blocks along the N to S direction in between basins bounded by normal faults, while the negative thermal anomaly clusters are spatially arranged along the coastal region, further north and in the western region far from the tilted mountain tectonic blocks of the Great Basin Section. The spatial maps that define regions with (positive or negative) thermal anomalies and distinct mean land response could assist landcover studies and support urban and rural planning in the context of emerging climatic change.     

Correlation between the structural pattern and the development of the hydrographic network in a portion of the Western Thessaly Basin (Greece)

In the context of the present study the structural pattern in the Western Thessaly Basin (Greece) has been examined, based on structural data collected from the entire study area that were further correlated to the hydrographic network. The structural pattern of the area was revealed from tectonic analysis. Additionally, the topography, stratigraphy and sedimentology of the study area were taken into account. GIS techniques were used to map the spatial distribution of the geological and tectonic features on the topographic relief of the area. The oldest structures are contractional in nature, deformed by normal faulting related to the extensional episodes initiated in Serravallian times. It is inferred that the orientation of the stress field in the area has changed several times: the N-S stress field which was dominant during Late Serravallian times changed to NW-SE (Messinian-Zanclian) up to E-W in Zanclian and subsequently to roughly NNW-SSE (in late Piacencian). The NE-SW stress that was dominant in Pleistocene became N-S in later times. In addition, some changes in orientation are also indicated for the transitional periods of the pre-mentioned extensional episodes, possibly related to local events, or as a block-related deformation. The development of the 7^th order streams is probably related to the N-S extensional faulting initiated in Pliocene times, while the dominant direction of the 6^th, 5^th and 4^th order streams is possibly connected with the presence of the NNE-SSW and NW-SE extensional faults. Finally, the lower order streams are probably related to the most recent E-W striking normal faults.