Numerical simulation of the thermal maturation, oil generation and migration in the Songliao Basin, Northeastern China

The discovery of the giant Daqing oil field in the Songliao Basin led to the realisation of the significant petroleum potential of non-marine basins. In order to reconstruct the basin evolution and oil formation, an integrated organic geochemical-basin modelling study along a regional transect across the Songliao Basin was conducted. It provided a regional heat flow evolution model, and revealed post-orogenic or late syn-orogenic maturation in the Central Depression and pre-orogenic maturation in the Southeast Uplift Zone. Kinetic parameters of petroleum generation for the lacustrine source formations are the basis for the simulation of oil generation and migration in the Songliao Basin. Using the principle activation energy peaking at 54 kcal/mol and a pre-exponential factor of about 4.2·10²⁷ Ma⁻¹, the simulation obtained a relatively good match with the measured transformation ratios. The Qingshankou Formation in the West and East Central Depressions constituted the major source in the basin. Major oil generation, migration and accumulation occurred during the Early Tertiary. In the West Central Depression, the generated oils migrated upwards into the Yaojia Formation followed by the updip migration into the Daqing Anticline and towards the local structural high along the West Slope. In contrast, the oil migration in the East Central Depression was dominated by the downward movement from the lower member of the Qingshankou Formation followed by the updip migration towards the Caoyang Anticline. The simulated oil accumulations are in good agreement with discovered oil fields, implying a potential application of the model for prediction and evaluation of new exploration targets in the basin.     

Modelling of a medium-term dynamics in a shallow tidal sea, based on combined physical and neural network methods

The paper presents an approach towards a medium-term (decades) modelling of water levels and currents in a shallow tidal sea by means of combined hydrodynamic and neural network models. The two-dimensional version of the hydrodynamic model Delft3D, forced with realistic water level and wind fields, is used to produce a two-year-database of water levels and currents in the study area. The linear principal component analysis (PCA) of the results is performed to reveal dominating spatial patterns in the analyzed dataset and to significantly reduce the dimensionality of the data. It is shown that only a few principal components (PCs) are necessary to reconstruct the data with high accuracy (over 95% of the original variance). Feed-forward neural networks are set up and trained to effectively simulate the leading PCs based on water level and wind speed and direction time series in a single, arbitrarily chosen point in the study area. Assuming that the spatial modes resulting from the PCA are ‘universally’ applicable to the data from time periods not modelled with Delft3D, the trained neural networks can be used to very effectively and reliably simulate temporal and spatial variability of water levels and currents in the study area. The approach is shown to be able to accurately reproduce statistical distribution of water levels and currents in various locations inside the study area and thus can be viewed as a reliable complementary tool e.g., for computationally expensive hydrodynamic modelling. Finally, a detailed analysis of the leading PCs is performed to estimate the role of tidal forcing and wind (including its seasonal and annual variability) in shaping the water level and current climate in the study area.     

Origin and emplacement of the impact formations at Chicxulub,Mexico, as revealed by the ICDP deep drilling at Yaxcopoil‐1 and by numerical modeling

We present and interpret results of petrographic, mineralogical, and chemical analyses of the 1511 m deep ICDP Yaxcopoil‐1 (Yax‐1) drill core, with special emphasis on the impactite units. Using numerical model calculations of the formation, excavation, and dynamic modification of the Chicxulub crater, constrained by laboratory data, a model of the origin and emplacement of the impact formations of Yax‐1 and of the impact structure as a whole is derived. The lower part of Yax‐1 is formed by displaced Cretaceous target rocks (610 m thick), while the upper part comprises six suevite‐type allochthonous breccia units (100 m thick). From the texture and composition of these lithological units and from numerical model calculations, we were able to link the seven distinct impact‐induced units of Yax‐1 to the corresponding successive phases of the crater formation and modification, which are as follows: 1) transient cavity formation including displacement and deposition of Cretaceous “megablocks;” 2) ground surging and mixing of impact melt and lithic clasts at the base of the ejecta curtain and deposition of the lower suevite right after the formation of the transient cavity; 3) deposition of a thin veneer of melt on top of the lower suevite and lateral transport and brecciation of this melt toward the end of the collapse of the transient cavity (brecciated impact melt rock); 4) collapse of the ejecta plume and deposition of fall‐back material from the lower part of the ejecta plume to form the middle suevite near the end of the dynamic crater modification; 5) continued collapse of the ejecta plume and deposition of the upper suevite; 6) late phase of the collapse and deposition of the lower sorted suevite after interaction with the inward flowing atmosphere; 7) final phase of fall‐back from the highest part of the ejecta plume and settling of melt and solid particles through the reestablished atmosphere to form the upper sorted suevite; and 8) return of the ocean into the crater after some time and minor reworking of the uppermost suevite under aquatic conditions. Our results are compatible with: a) 180 km and 100 km for the diameters of the final crater and the transient cavity of Chicxulub, respectively, as previously proposed by several authors, and b) the interpretation of Chicxulub as a peak‐ring impact basin that is at the transition to a multi‐ring basin.     

A statistical retrieval algorithm for root zone soil moisture

An algorithm for the estimation of root zone soil moisture is presented. Global fields of the soil moisture within the uppermost metre of soil are derived with a temporal resolution of 10 days. For calibration, long-term soil moisture observations from the former Soviet Union are used. The variance of the measurements is largely dominated by the spatial variability of the long-term mean soil moisture, while the temporal variability gives comparatively small contribution. Consequently, the algorithm is organised into two steps. The first step concentrates on the retrieval of the spatial variance of the long-term means, which comprises more than 85% of the total soil moisture variability. A major part of the spatial variance can be explained by four easily available fields: the climatological precipitation, land use, soil texture, and terrain slope. The second step of the algorithm is dedicated to the local temporal variability. This part of variability is recovered by using passive microwave data from scanning multichannel microwave radiometre (SMMR) supported by monthly averaged fields of air temperature and precipitation. The 6-GHz channel of SMMR is shown to be severely disturbed by radio frequency interference, so that information from the 10-GHz channel is used instead. The algorithm provides reasonable soil moisture fields which is confirmed by a comparison with independent measurements from Illinois.     

The Decay of Supersonic and Super-Alfvénic Turbulence in Star-Forming Clouds

Supersonic and super-Alfvénic turbulent motions generate an intricate pattern of shock waves through which the turbulence decays. We here illustrate the spatial structure during shock formation and decay. This revised version was published online in July 2006 with corrections to the Cover Date.     

Selected topics in arctic atmosphere and climate

This paper summarizes the main elements of four IPY projects that examine the Arctic Atmosphere. All four projects focus on present conditions with a view to anticipating possible climate change. All four investigate the Arctic atmosphere, ocean, ice, and land interfacial surfaces. One project uses computer models to simulate the dynamics of the Arctic atmosphere, storms, and their interactions with the ocean and ice interface. Another project uses statistical methods to infer transports of pollutants as simulated in large-scale global atmospheric and oceanic models verifying results with available observations. A third project focuses on measurements of pollutants at the ice-ocean?Catmosphere interface, with reference to model estimates. The fourth project is concerned with multiple, high accuracy measurements at Eureka in the Canadian Archipelago. While these projects are distinctly different, led by different teams and interdisciplinary collaborators, with different technical approaches and methodologies, and differing objectives, they all strive to understand the processes of the Arctic atmosphere and climate, and to lay the basis for projections of future changes. Key findings include: ? Decreased sea ice leads to more intense storms, higher winds, reduced surface albedo, increased surface air temperature, and enhanced vertical mixing in the upper ocean. ? Arctic warming may affect toxic chemicals by remobilizing persistent organic pollutants and augmenting mercury deposition/retention in the environment. ? Changes in sea ice can dramatically change processes in and at the ice surface related to ozone, mercury and bromine oxide and related chemical/physical properties. ? Structure and properties of the Arctic atmospheric??troposphere to stratosphere??and tracking of transport of pollution and smoke plumes from mid-latitudes to the poles.     

Identification of extreme storm surges with high-impact potential along the German North Sea coastline

Planning and design of coastal protection rely on information about the probabilities of very severe storm tides and the possible changes that may occur in the course of climate change. So far, this information is mostly provided in the form of high percentiles obtained from frequency distributions or return values. More detailed information and assessments of events that may cause extreme damages or have extreme consequences at the coast are so far still unavailable. We describe and compare two different approaches that may be used to identify highly unlikely but still physically possible and plausible events from model simulations. Firstly, in the case when consistent wind and tide-surge data are available, different metrics such as the height of the storm surge can be derived directly from the simulated water levels. Secondly, in cases where only atmospheric data are available, the so called effective wind may be used. The latter is the projection of the horizontal wind vector on that direction which is most effective in producing surges at the coast. Comparison of events identified by both methods show that they can identify extreme events but that knowledge of the effective wind alone does not provide sufficient information to identify the highest storm surges. Tracks of the low-pressure systems over the North Sea need to be investigated to find those cases, where the duration of the high wind is too short to induce extreme storm tides. On the other hand, factors such as external surges or variability in mean sea level may enhance surge heights and are not accounted for in estimates based on effective winds only. Results from the analysis of an extended data set suggest that unprecedented storm surges at the German North Sea coast are possible even without taking effects from rising mean sea level into account. The work presented is part of the ongoing project “Extreme North Sea Storm Surges and Their Consequences” (EXTREMENESS) and represents the first step towards an impact assessment for very severe storm surges which will serve as a basis for development of adaptation options and evaluation criteria.     

Source rock potential of the Middle Jurassic to Middle Pliocene,onshore Nile Delta Basin,Egypt

Organic geochemical characterization of cutting samples from the Abu Hammad-1 and Matariya-1 wells elucidates the depositional environment and source rock potential of the Jurassic and Lower Cretaceous successions and the Middle Miocene to Pleistocene section in the southern and eastern Nile Delta Basin. The burial and thermal histories of the Mesozoic and Miocene sections were modeled using 1D basin modeling based on input data from the two wells. This study reveals fair to good gas-prone source rocks within the Upper Jurassic and Lower Cretaceous sections with total organic carbon (TOC) averaging 2.7% and hydrogen index (HI) up to 130 mg HC/g TOC. The pristane/n-C₁₇ versus phytane/n-C₁₈ correlation suggests mixed marine and terrestrial organic matter with predominant marine input. Burial and thermal history modeling reveals low thermal maturity due to low heat flow and thin overburden. These source rocks can generate gas in the western and northern parts of the basin where they are situated at deeper settings. In contrast, the thick Middle Miocene shows fair source rock quality (TOC averaging at 1.4%; HI maximizing at 183 mg HC/g TOC). The quality decreases towards the younger section where terrestrial organic matter is abundant. This section is similar to previously studied intervals in the eastern Nile Delta Basin but differs from equivalents in the central parts where the quality is better. Based on 1D modeling, the thick Middle Miocene source rocks just reached the oil generation stage, but microbial gas, however, is possible.     

Collision of the Kronotskiy arc at the NE Eurasia margin and structural evolution of the Kamchatka–Aleutian junction

Structural evolution of the Kamchatka–Aleutian junction area in late Mesozoic and Tertiary was generally controlled by (1) the processes of subduction in Kronotskiy and Proto-Kamchatka subduction zones and (2) collision of the Kronotskiy arc against NE Eurasia margin. Two structural zones of the pre-Pliocene age and six structural assemblages are recognized in studied region. 1: Eastern ranges zone comprises SE-vergent thrust folded belt, which evolved in accretionary and collisional setting. Two structural assemblages (ER1 and ER2), developed there, document shortening in the NW–SE direction and in the N–S direction, respectively. 2: Eastern Peninsulas zone generally corresponds to Kronotskiy arc terrane. Four structural assemblages are recognized in this zone. They characterize (1) precollisional deformations in the accretionary wedge (EP1) and in the fore-arc basin and volcanic belt (EP2), and (2) syn-collisional deformation of the entire Kronotskiy terrane in plunging folds (EP3) and deformations in the foreland basin (EP4). Analysis of paleomagnetic declinations versus present day structural strike in the Kronotskiy arc terrane shows that originally the arc was trending from west to east. Relative position of the accretionary wedge, fore-arc basin and volcanic belt, as well as northward dipping thrusts in accretionary wedge indicate, that a northward dipping subduction zone was located south of the arc. The accretionary wedge developed from the Late Cretaceous through the Eocene, and it implies that the subduction zone maintained its direction and position during this time. It implies that Kronotskiy arc was neither a part of the Pacific nor Kula plates and was located on an individual smaller plate, which included the arc and Vetlovka back-arc basin. Motion of the Kronotskiy arc towards Eurasia was connected only with NW-directed subduction at Kamchatka margin since Middle Eocene (42–44 Ma). Emplacement of the Kronotskiy arc at the Kamchatka margin occurred between Late Eocene and Early Miocene. This is based on the age of syn-collisional plunging folds in Kronotskiy terrane, and provenance data for the Upper Eocene to Middle Miocene Tyushevka basin, which indicate in situ evolution of the basin with respect to Kamchatka. Collision was controlled by the common motion of the Kronotskiy arc with Pacific plate towards the northwest, and by the motion of the Eurasian margin towards the south. The latter motion was responsible for the southward deflection of the western part of the Kronotskiy arc (EP3 structures), and for oblique transpressional structures in the collisional belt (ER2 structures).