Evolution of cool skin and direct air-sea gas transfer coefficient during daytime

Absorption of solar radiation within the thermal molecular sublayer of the ocean can modify the temperature difference across the cool skin as well as the air-sea gas transfer. Our model of renewal type is based on the assumption that the thermal and diffusive molecular sublayers below the ocean surface undergo cyclic growth and destruction, the heat and gas transfer between the successive burst events are performed by molecular diffusion. The model has been upgraded to include heating due to solar radiation. The renewal time is parameterized as a function of the surface Richardson number and the Keulegan number. A Rayleigh number criterion characterizes the convective instability of the cool skin under solar heating. Under low wind speed conditions, the solar heating can damp the convective instability, strongly increasing the renewal time and correspondingly decreasing the interfacial gas exchange. In the ocean, an additional convective instability caused by salinity flux due to evaporation becomes of importance in such cases. The new parameterization is compared with the cool skin data obtained in the western equatorial Pacific during the Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment in February 1993. In combination with a model of the diurnal thermocline it describes main features of the field data both in nighttime and daytime. Under low wind speed conditions (< 5 m s^-1) diurnal variations of the sea surface temperature due to the formation of a diurnal thermocline were substantially larger than those across the cool skin. Under wind speeds > 5 m s^-1, diurnal variations of the surface temperature due to the variations of the thermal molecular sublayer become more important.     

Gravity field determination using boundary element methods

The Boundary Element Method (BEM), a numerical technique for solving boundary integral equations, is introduced to determine the earth's gravity field. After a short survey on its main principles, we apply this method to the fixed gravimetric boundary value problem (BVP), i.e. the determination of the earth's gravitational potential from measurements of the intensity of the gravity field in points on the earth's surface. We show how to linearize this nonlinear BVP using an implicit function theorem and how to transform the linearized BVP into a boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation into a linear system of equations. We discuss the major problems of this approach for setting up and solving the linear system. The BVP is numerically solved for a bounded part of the earth's surface using a high resolution reference gravity model, measured gravity values of high density, and a 50 50 m² digital terrain model to describe the earth's surface. We obtain a gravity field resolution of 1 1 km² with an accuracy of the order 10^–3 to 10^–4 in about 1 CPU-hour on a Siemens/Fujitsu SIMD vector pipeline machine using highly sophisticated numerical integration techniques and fast equation solvers. We conclude that BEM is a powerful numerical tool for solving boundary value problems and may be an alternative to classical geodetic techniques.     

Intercomparison of measurements of stratospheric hydrogen fluoride

Observations of the vertical profile of hydrogen fluoride (HF) vapor in the stratosphere and of the vertical column amounts of HF above certain altitudes were made using a variety of spectroscopic instruments in the 1982 and 1983 Balloon Intercomparison Campaigns. Both emission instruments working in the far infrared spectral region and absorption instruments using solar occultation in the 2.5m region were employed. No systematic differences were seen in results from the two spectral regions. A mean profile from 20–45 km is presented, with uncertainties ranging from 20% to 50%. Total columns measured from ground and from 12 km are consistent with the profile if the mixing ratio for HF is small in the tropophere and low stratosphere.     

Oil-generation kinetics for organic facies with Type-II and -IIS kerogen in the Menilite Shales of the Polish Carpathians

The Menilite Shales (Oligocene) of the Polish Carpathians are the source of low-sulfur oils in the thrust belt and some high-sulfur oils in the Carpathian Foredeep. These oil occurrences indicate that the high-sulfur oils in the Foredeep were generated and expelled before major thrusting and the low-sulfur oils in the thrust belt were generated and expelled during or after major thrusting. Two distinct organic facies have been observed in the Menilite Shales. One organic facies has a high clastic sediment input and contains Type-II kerogen. The other organic facies has a lower clastic sediment input and contains Type-IIS kerogen. Representative samples of both organic facies were used to determine kinetic parameters for immiscible oil generation by isothermal hydrous pyrolysis and S₂ generation by non-isothermal open-system pyrolysis. The derived kinetic parameters showed that timing of S₂ generation was not as different between the Type-IIS and -II kerogen based on open-system pyrolysis as compared with immiscible oil generation based on hydrous pyrolysis. Applying these kinetic parameters to a burial history in the Skole unit showed that some expelled oil would have been generated from the organic facies with Type-IIS kerogen before major thrusting with the hydrous-pyrolysis kinetic parameters but not with the open-system pyrolysis kinetic parameters. The inability of open-system pyrolysis to determine earlier petroleum generation from Type-IIS kerogen is attributed to the large polar-rich bitumen component in S₂ generation, rapid loss of sulfur free-radical initiators in the open system, and diminished radical selectivity and rate constant differences at higher temperatures. Hydrous-pyrolysis kinetic parameters are determined in the presence of water at lower temperatures in a closed system, which allows differentiation of bitumen and oil generation, interaction of free-radical initiators, greater radical selectivity, and more distinguishable rate constants as would occur during natural maturation. Kinetic parameters derived from hydrous pyrolysis show good correlations with one another (compensation effect) and kerogen organic-sulfur contents. These correlations allow for indirect determination of hydrous-pyrolysis kinetic parameters on the basis of the organic-sulfur mole fraction of an immature Type-II or -IIS kerogen.     

Late Holocene sea-level changes and isostasy in western Denmark

Cores and exposed cliff sections in salt marshes around Ho Bugt, a tidal embayment in the northernmost part of the Danish Wadden Sea, were subjected to ¹⁴C dating and litho- and biostratigraphical analyses to reconstruct paleoenvironmental changes and to establish a late Holocene relative sea-level history. Four stages in the late Holocene development of Ho Bugt can be identified: (1) groundwater-table rise and growth of basal peat (from at least 2300 BC to AD 0); (2) salt-marsh formation (0 to AD 250); (3) a freshening phase (AD 250 to AD 1600?), culminating in the drying out of the marshes and producing a distinct black horizon followed by an aeolian phase with sand deposition; and (4) renewed salt-marsh deposition (AD 1600? to present). From 16 calibrated AMS radiocarbon ages on fossil plant fragments and 4 calibrated conventional radiocarbon ages on peat, we reconstructed a local relative sea-level history that shows a steady sea-level rise of 4 m since 4000 cal yr BP. Contrary to suggestions made in the literature, the relative sea-level record of Ho Bugt does not contain a late Holocene highstand. Relative sea-level changes at Ho Bugt are controlled by glacio-isostatic subsidence and can be duplicated by a glacial isostatic adjustment model in which no water is added to the world's oceans after ca. 5000 cal yr BP.     

T -distributed Random Fields: A Parametric Model for Heavy-tailedWell-log Data1

Histograms of observations from spatial phenomena are often found to be more heavy-tailed than Gaussian distributions, which makes the Gaussian random field model unsuited. A T-distributed random field model with heavy-tailed marginal probability density functions is defined. The model is a generalization of the familiar Student-T distribution, and it may be given a Bayesian interpretation. The increased variability appears cross-realizations, contrary to in-realizations, since all realizations are Gaussian-like with varying variance between realizations. The T-distributed random field model is analytically tractable and the conditional model is developed, which provides algorithms for conditional simulation and prediction, so-called T-kriging. The model compares favourably with most previously defined random field models. The Gaussian random field model appears as a special, limiting case of the T-distributed random field model. The model is particularly useful whenever multiple, sparsely sampled realizations of the random field are available, and is clearly favourable to the Gaussian model in this case. The properties of the T-distributed random field model is demonstrated on well log observations from the Gullfaks field in the North Sea. The predictions correspond to traditional kriging predictions, while the associated prediction variances are more representative, as they are layer specific and include uncertainty caused by using variance estimates.     

The Record of an Extreme Flood in the Sediments of Montane Lillooet Lake, British Columbia: Implications for Paleoenvironmental Assessment

Severe rainfall in mid October, 2003 produced the largest floods in almost a century of record on rivers in the Cordillera of southwestern British Columbia. Sediment deposited in Lillooet Lake as a result of this event is clearly distinguished by stratigraphy, colour, texture, magnetic properties, and organic content. Each of these physical properties is related to the lacustrine processes, especially turbid underflow, that distributed the sediment through the lake. The flood, which lasted less than a week, delivered 8–12 times the amount of sediment that accumulates in most entire years in the deepest, central parts of the lake. Recognition of events of this type in the stratigraphic record offers a means of assessing the changing nature of extreme hydroclimatic events, and their relation to more ubiquitous, lower-energy processes.     

Lake Sediment Core Records of Sulphur Accumulation and Sulphur Isotopic Composition in Central Ontario, Canada Lakes

Stable isotopic compositions and concentrations of total sedimentary sulphur (S) were determined in cores from 6 lakes in the acid-sensitive Muskoka-Haliburton region of south-central Ontario. The isotopic composition of S in deep sediment (> ~ 20 cm) was approximately constant in all lakes, and indicated a pre-industrial δ ³⁴S value between +4.0 and +5.3‰, which is similar to current bulk deposition. Similarly, total S concentrations in deep sediment were relatively low (1.9–5 mg S g⁻¹ dwt) and approximately constant with depth within cores. All lakes exhibited up-core increases in total S and decreases in δ ³⁴S at a depth corresponding to the beginning of industrialization in the Great Lakes region ( ~ 1900), resulting in a generally reciprocal depth pattern between total S concentration and δ ³⁴S ratios. While initial shifts in total S and δ ³⁴S were likely due to enhanced SO₄ reduction of newly available anthropogenic SO₄, both the magnitude and pattern of up-core S enrichment and shifts in δ ³⁴S varied greatly among lakes, and did not match changes in S deposition post 1900. Differences between lakes in total S and δ ³⁴S were not related to any single hydrologic (e.g., residence time) or physical (e.g., catchment-area-to-lake area ratio) lake characteristic. This work indicates that sediment cores do not provide consistent records of changes in post-industrial S deposition in this region, likely due to redox-related mobility of S in upper sediment.     

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

The   M _w γ 9.0  2004 December 26 Sumatra-Andaman and   M _w = 8.7  2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well.