Lake Koronia, Greece: Shift from autotrophy to heterotrophy with cultural eutrophication and progressive water-level reduction

Lake Koronia, a Ramsar site, is shallow, polymictic, hypertrophic and until recently was aerially the fourth largest lake in Greece. Although exceeding 5 m in the past, lake depth has declined progressively from 3.8 m in 1980 to < 1 m in 1997, reducing surface area and water volume by 50% and 80%, respectively. Specific conductivityincreased from 1300 μS cm⁻¹ in 1977 to >6000 μS cm⁻¹ in 1991. Increased phosphate concentrations from the late 1970's (8–45 μg L⁻¹) to the late 1990's (100–1000 μg L⁻¹) document that the previously eutrophic system with a limited littoral zone switched to hypertrophy dominated by massive cyanobacteria blooms. Oxygen saturation of the water column increased progressively from about 80% in 1983 to full saturation about 1993, after which it decreased progressively to only 20% saturation in 1997. In spite of cyanobacteria dominance, community metabolism of the lake switched from progressively increasing autotrophy to rapidly advancing heterotrophy associated with progressive water-level reduction leading to fish extirpation in the lake.     

Torngat ultramafic lamprophyres and their relation to the North Atlantic Alkaline Province

Geological mapping and diamond exploration in northern Quebec and Labrador has revealed an undeformed ultramafic dyke swarm in the northern Torngat Mountains. The dyke rocks are dominated by an olivine-phlogopite mineralogy and contain varying amounts of primary carbonate. Their mineralogy, mineral compositional trends and the presence of typomorphic minerals (e.g. kimzeyitic garnet), indicate that these dykes comprise an ultramafic lamprophyre suite grading into carbonatite. Recognized rock varieties are aillikite, mela-aillikite and subordinate carbonatite. Carbonatite and aillikite have in common high carbonate content and a lack of clinopyroxene. In contrast, mela-aillikites are richer in mafic silicate minerals, in particular clinopyroxene and amphibole, and contain only small amounts of primary carbonate. The modal mineralogy and textures of the dyke varieties are gradational, indicating that they represent end-members in a compositional continuum.

The Torngat ultramafic lamprophyres are characterized by high but variable MgO (10–25 wt.%), CaO (5–20 wt.%), TiO₂ (3–10 wt.%) and K₂O (1–4 wt.%), but low SiO₂ (22–37 wt.%) and Al₂O₃ (2–6 wt.%). Higher SiO₂, Al₂O₃, Na₂O and lower CO₂ content distinguish the mela-aillikites from the aillikites. Whereas the bulk rock major and trace element concentrations of the aillikites and mela-aillikites overlap, there is no fractional crystallization relation between them. The major and trace element characteristics imply related parental magmas, with minor olivine and Cr-spinel fractionation accounting for intra-group variation.

The Torngat ultramafic lamprophyres have a Neoproterozoic age and are spatially and compositionally closely related with the Neoproterozoic ultramafic lamprophyres from central West Greenland. Ultramafic potassic-to-carbonatitic magmatism occurred in both eastern Laurentia and western Baltica during the Late Neoproterozoic. It can be inferred from the emplacement ages of the alkaline complexes and timing of Late Proterozoic processes in the North Atlantic region that this volatile-rich, deep-seated igneous activity was a distal effect of the breakup of Rodinia. This occurred during and/or after the rift-to-drift transition that led to the opening of the Iapetus Ocean.     

GIS Spatial-Temporal Modeling of Water Systems in Greater Toronto Area, Canada

Modeling landscape with high-resolution digital elevation model (DEM) in a geographic information system can provide essential morphological and structural information for modeling surface processes such as geomorphologic process and water systems. This paper introduces several DEM-based spatial analysis processes applied to characterize spatial distribution and their interactions of ground and surface water systems in the Greater Toronto Area (GTA), Canada. The stream networks and drainage basin systems were derived from the DEM with 30 m resolution and the regularities of the surface stream and drainage patterns were modeled from a statistical/multifractal point of view. Together with the elevation and slope of topography, other attributes defmed from modeling the stream system, and drainage networks were used to associate geological, hydrological and topographical features to water flow in river systems and the spatial locations of artesian aquifers in the study area. Stream flow data derived from daily flow measurements recorded at river gauging stations for multi-year period were decomposed into “drainage-area dependent“ and “drainage-area independent“ flow components by two-step “frequency“ and “spatial“ analysis processes. The latter component was further demonstrated to relate most likely to the ground water discharge. An independent analysis was conducted to model the distribution of aquifers with information derived from the records of water wells. The focus was given on quantification of the likelihood of ground water discharge to river and ponds through flowing wells, springs and seepages. It has been shown that the Oak Ridges Moraine (ORM) is a unique glacial deposit that serves as a recharge layer and that the aquifers in the ORM underlain by Hilton Tills and later deposits exposed near the steep slope zone of the ridges of ORM provide significant discharge to the surface water systems (river flow and ponds) through flowing wells, springs and seepages. Various statistics (cross- and auto-correlation coefficients, fractal R/S exponent) were used in conjunction with GIS to demonstrate the influence of land types, topography and geometry of drainage basins on short- and long-term persistence of river flows as well as responding time to precipitation events. The current study has provided not only insight in understanding the interaction of water systems in the GTA, but also a base for further establishment of an on-line GIS system for predicting spatial-temporal changes of river flow and groundwater level in the GTA.     

Seismic studies of the Brasília fold belt at the western border of the São Francisco Craton, Central Brazil, using receiver function, surface-wave dispersion and teleseismic tomography

The Tocantins Province in Central Brazil is composed of a series of SSW–NNE trending terranes of mainly Proterozoic ages, which stabilized in the Neoproterozoic in the final collision between the Amazon and São Francisco cratons. No previous information on crustal seismic properties was available for this region. Several broadband stations were used to study the regional patterns of crustal and upper mantle structure, extending the results of a recent E–W seismic refraction profile. Receiver functions and surface wave dispersion showed a thin crust (33–37 km) in the Neoproterozoic Magmatic Arc terrane. High average crustal Vp/Vs ratios (1.74–1.76) were consistently observed in this unit. The foreland domain of the Brasília foldbelt, on the other hand, is characterized by thicker crust (42–43 km). Low Vp/Vs ratios (1.70–1.72) were observed in the low-grade foreland fold and thrust zone of the Brasília belt adjacent to the São Francisco craton. Teleseismic P-wave tomography shows that the lithospheric upper mantle has lower velocities beneath the Magmatic Arc and Goiás Massif compared with the foreland zone of the belt and São Francisco craton. The variations in crustal thickness and upper mantle velocities observed with the broadband stations correlate well with the measurements along the seismic refraction profile. The integration of all seismic observations and gravity data indicates a strong lithospheric contrast between the Goiás Massif and the foreland domain of the Brasília belt, whereas little variation was found across the foldbelt/craton surface boundary. These results support the hypothesis that the Brasília foreland domain and the São Francisco craton were part of a larger São Francisco-Congo continental plate in the final collision with the Amazon plate.     

Pliocene and Quaternary regional uplift in western Turkey: the Gediz River terrace staircase and the volcanism at Kula

Along the upper reaches of the Gediz River in western Turkey, in the eastern part of the Aegean extensional province, the land surface has uplifted by 400 m since the Middle Pliocene. This uplift is revealed by progressive gorge incision, and its rate can be established because river terraces are capped by basalt flows that have been K–Ar and Ar–Ar dated. At present, the local uplift rate is 0.2 mm a⁻¹. Uplift at this rate began around the start of the Middle Pleistocene, following a span of time when the uplift was much slower. This was itself preceded by an earlier uplift phase, apparently in the late Late Pliocene and early Early Pleistocene, when the uplift rate was comparable to the present. The resulting regional uplift history resembles what is observed in other regions and is analogously interpreted as the isostatic response to changing rates of surface processes linked to global environmental change. We suggest that this present phase of surface uplift, amounting so far to 150 m, is being caused by the nonsteady-state thermal and isostatic response of the crust to erosion, following an increase in erosion rates in the late Early Pleistocene, most likely as a result of the first large northern-hemisphere glaciation during oxygen isotope stage 22 at 870 ka. We suggest that the earlier uplift phase, responsible for the initial 250 m of uplift, resulted from a similar increase in erosion rates caused by the deterioration in local climate at 3.1 Ma. This uplift thus has no direct relationship to the crustal extension occurring in western Turkey, the rate and sense of which are thought not to have changed significantly on this time scale. Our results thus suggest that the present, often deeply incised, landscape of western Turkey has largely developed from the Middle Pleistocene onwards, for reasons not directly related to the active normal faulting that is also occurring. The local isostatic consequences of this active faulting are instead superimposed onto this “background” of regional surface uplift. Modelling of this surface uplift indicates that the effective viscosity of the lower continental crust beneath this part of Turkey is of the order of 10¹⁹ Pa s, similar to a recent estimate for beneath central Greece. The lower uplift rates observed in western Turkey, compared with central Greece, result from the longer typical distances of fluvial sediment transport, which cause weaker coupling by lower-crustal flow between offshore depocentres and eroding onshore regions that provide the sediment source.     

Construction of detrital mineral populations: insights from mixing of U–Pb zircon ages in Himalayan rivers

Fission‐track, U–Pb and Pb–Pb analyses of detrital heavy mineral populations in depositional basins and modern river sediments are widely used to infer the exhumational history of mountain belts. However, relatively few studies address the underlying assumption that detrital mineral populations provide an accurate representation of their entire source region. Implicit in this assumption is the idea that all units have equal potential to contribute heavy minerals in proportion to their exposure area in the source region. In reality, the detrital mineral population may be biased by variable concentrations of minerals in bedrock and differential erosion rates within the source region. This study evaluates the relative importance of these two variables by using mixing of U–Pb zircon ages to trace zircon populations from source units, through the fluvial system, and into the foreland. The first part of the study focuses on the Marsyandi drainage in central Nepal, using tributaries that drain single formations to define the U–Pb age distributions of individual units and using trunk river samples to evaluate the relative contributions from each lithology. Observed mixing proportions are compared with proportions predicted by a simple model incorporating lithologic exposure area and zircon concentration. The relative erosion rates that account for the discrepancy between the observed and predicted mixing proportions are then modelled and compared with independent erosional proxies. The study also compares U–Pb age distributions from four adjacent drainages spanning ～250 km along the Himalayan front using the Kolmogorov–Smirnov statistic and statistical estimates of the proportion of zircon derived from each upstream lithology. Results show that, along this broad swath of rugged mountains, the U–Pb age distributions are remarkably similar, thereby allowing data from more localized sources to be extrapolated along strike.     

Continental rift systems and anorogenic magmatism

Precambrian Laurentia and Mesozoic Gondwana both rifted along geometric patterns that closely approximate truncated-icosahedral tessellations of the lithosphere. These large-scale, quasi-hexagonal rift patterns manifest a least-work configuration. For both Laurentia and Gondwana, continental rifting coincided with drift stagnation, and may have been driven by lithospheric extension above an insulated and thermally expanded mantle. Anorogenic magmatism, including flood basalts, dike swarms, anorthosite massifs and granite-rhyolite provinces, originated along the Laurentian and Gondwanan rift tessellations. Long-lived volcanic regions of the Atlantic and Indian Oceans, sometimes called hotspots, originated near triple junctions of the Gondwanan tessellation as the supercontinent broke apart. We suggest that some anorogenic magmatism results from decompression melting of asthenosphere beneath opening fractures, rather than from random impingement of hypothetical deep-mantle plumes.     

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.