Geochemical reconstruction of late Holocene drainage and mixing in Kluane Lake,Yukon Territory

The level of Kluane Lake in southwest Yukon Territory, Canada, has fluctuated tens of metres during the late Holocene. Contributions of sediment from different watersheds in the basin over the past 5,000 years were inferred from the elemental geochemistry of Kluane Lake sediment cores. Elements associated with organic material and oxyhydroxides were used to reconstruct redox fluctuations in the hypolimnion of the lake. The data reveal complex relationships between climate and river discharge during the late Holocene. A period of influx of Duke River sediment coincides with a relatively warm climate around 1,300 years BP. Discharge of Slims River into Kluane Lake occurred when Kaskawulsh Glacier advanced to the present drainage divide separating flow to the Pacific Ocean via Kaskawulsh and Alsek rivers from flow to Bering Sea via tributaries of Yukon River. During periods when neither Duke nor Slims river discharged into Kluane Lake, the level of the lake was low and stable thermal stratification developed, with anoxic and eventually euxinic conditions in the hypolimnion.     

Is the ocean floor a fractal?

The topographic structure of the ocean bottom is investigated at different scales of resolution to answer the question: Can the seafloor be described as a fractal process? Methods from geostatistics, the theory of regionalized variables, are used to analyze the spatial structure of the ocean floor at different scales of resolution. The key to the analysis is the variogram criterion: Self-similarity of a stochastic process implies self-similarity of its variogram. The criterion is derived and proved here: it also is valid for special cases of self-affinity (in a sense adequate for topography). It has been proposed that seafloor topography can be simulated as a fractal (an object of Hausdorff dimension strictly larger than its topological dimension), having scaling properties (self-similarity or self-affinity). The objective of this study is to compare the implications of these concepts with observations of the seafloor. The analyses are based on SEABEAM bathymetric data from the East Pacific Rise at 13°N/104°W and at 9°N/104°W and use tracks that run both across the ridge crest and along the ridge flank. In the geostatistical evaluation, the data are considered as a stochastic process. The spatial continuity of this process is described by variograms that are calculated for different scales and directions. Applications of the variogram criterion to scale-dependent variogram models yields the following results: Although the seafloor may be a fractal in the sense of the definition involving the Hausdorff dimension, it is not self-similar, nor self-affine (in the given sense). Mathematical models of scale-dependent spatial structures are presented, and their relationship to geologic processes such as ridge evolution, crust formation, and sedimentation is discussed.     

Oxidation of CO on surface hematite in high CO2 atmospheres

We propose a mechanism for the oxidation of gaseous CO into CO₂ occurring on the surface mineral hematite (Fe₂O₃(s)) in hot, CO₂-rich planetary atmospheres, such as Venus. This mechanism is likely to constitute an important source of tropospheric CO₂ on Venus and could at least partly address the CO₂ stability problem in Venus’ stratosphere, since our results suggest that atmospheric CO₂ is produced from CO oxidation via surface hematite at a rate of 0.4 petagrammes (Pg) CO₂ per (Earth) year on Venus which is about 45% of the mass loss of CO₂ via photolysis in the Venusian stratosphere. We also investigated CO oxidation via the hematite mechanism for a range of planetary scenarios and found that modern Earth and Mars are probably too cold for the mechanism to be important because the rate-limiting step, involving CO(g) reacting onto the hematite surface, proceeds much slower at lower temperatures. The mechanism may feature on extrasolar planets such as Gliese 581c or CoRoT-7b assuming they can maintain solid surface hematite which, e.g. starts to melt above about 1200 K. The mechanism may also be important for hot Hadean-type environments and for the emerging class of hot Super-Earths with planetary surface temperatures between about 600 and 900 K.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

A COMPARISON OF HYDROLOGY AND VEGETATION BETWEEN A CHANNELIZED STREAM AND A NONCHANNELIZED STREAM IN WESTERN TENNESSEE

The purpose of this study was to provide baseline data on floodplain forest structure, composition, and function that would be needed to predict and monitor the consequences of a proposed stream restoration project. This project would involve the “dechannelization” of Stokes Creek, a stream in western Tennessee that was channelized and leveed in the first half of the 1900s. To this end, we collected data on surface hydrology, soil redox potential (Eh), and the structure and composition of the floodplain vegetation of Stokes Creek. To place our findings into a regional context, we also collected comparable vegetation data from plots located along a nonchannelized stream reach of the Wolf River near Moscow, Tennessee. While hydrologic fluctuations of floodplain sites were synchronous with river dynamics for the Wolf River, the hydrology of floodplain sites at Stokes Creek was constrained by the influence of beaver dams, backflooding, and ponding of overland flow behind levees. Consequently, composition of the forest overstory, understory, and herbaceous strata was significantly different between the two sites. For example, Stokes Creek had a noticeable lack of cypress and tupelo sites, and a greater abundance of red maple. Analyses of size-class structure and woody debris quantity reinforced the existing differences between the more natural and human-impacted systems. While the current hydrology apparently has a negative affect on bottomland hardwoods, scattered regeneration stems and soil redox measurements indicate that a dechannelization effort that yielded lower water tables in the Stokes Creek floodplain potentially could increase bottomland hardwood establishment. [Key words: channelization, western Tennessee, bottomland hardwoods, hydrology.]