Modeling ground water flow in alluvial mountainous catchments on a watershed scale

In large mountainous catchments, shallow unconfined alluvial aquifers play an important role in conveying subsurface runoff to the foreland. Their relatively small extent poses a serious problem for ground water flow models on the river basin scale. River basin scale models describing the entire water cycle are necessary in integrated water resources management and to study the impact of global climate change on ground water resources. Integrated regional-scale models must use a coarse, fixed discretization to keep computational demands low and to facilitate model coupling. This can lead to discrepancies between model discretization and the geometrical properties of natural systems. Here, an approach to overcome this discrepancy is discussed using the example of the German-Austrian Upper Danube catchment, where a coarse ground water flow model was developed using MODFLOW. The method developed uses a modified concept from a hydrological catchment drainage analysis in order to adapt the aquifer geometry such that it respects the numerical requirements of the chosen discretization, that is, the width and the thickness of cells as well as gradients and connectivity of the catchment. In order to show the efficiency of the developed method, it was tested and compared to a finely discretized ground water model of the Ammer subcatchment. The results of the analysis prove the applicability of the new approach and contribute to the idea of using physically based ground water models in large catchments.     

Influence of root characteristics and soil variables on the uprooting mechanics of Avena sativa and Medicago sativa seedlings

The success of seedlings and rejuvenated woody debris growing on river bedforms depends on the resistance to uprooting by flow provided by their simple root architecture. Avena sativa and Medicago sativa seedlings were used in flume experiments as prototypes for juvenile riparian plants. Very little is known about the magnitude of root anchoring forces and the role of secondary roots of such simple root systems. We performed 1550 vertical uprooting experiments on Avena sativa and Medicago sativa seedlings grown in quartz sand. Seedlings were pulled up by direct traction using a wheel driven by a computer‐controlled motor and the force was recorded. Roots were scanned and architectural parameters (root length and number of roots) determined. Uprooting force and work (the integral of the applied force times the distance over which it is applied) were then related to root architecture and soil variables. Resistance to uprooting increased with decreasing sediment size and sediment moisture content. The initial response of the root–soil system to uprooting showed linear elastic behaviour with modulus increasing with plant age. While the maximum uprooting force was found to increase linearly with total root length and be mainly dependent on the length of the main root, uprooting work followed a power law and has to be related to the whole root system. Thus, for the young plants we considered, secondary roots are responsible for the ability to withstand environmental disturbances in terms of duration rather than magnitude. This distinction between primary and secondary roots can be of crucial importance for seedlings of riparian species germinating on river bars and islands where inundation is a main cause of mortality. Beyond clarifying the biomechanical role of soil and root variables, the uprooting statistics obtained are useful in interpreting and designing ecomorphodynamic flume experiments. Copyright © 2014 John Wiley & Sons, Ltd.     

On the developments of spectral wave models: numerics and parameterizations for the coastal ocean

The development of numerical wave models for coastal applications, including coupling with ocean circulation models, has spurred an ongoing effort on theoretical foundations, numerical techniques, and physical parameterizations. Some important aspects of this effort are reviewed here, and results are shown in the case of the French Atlantic and Channel coast using version 4.18 of the WAVEWATCH III ^R model. Compared to previous results, the model errors have been strongly reduced thanks to, among other things, the introduction of currents, coastal reflection, and bottom sediment types. This last item is described here for the first time, allowing unprecedented accuracy at some sites along the French Atlantic Coast. The adequate resolution, necessary to represent strong gradients in tidal currents, was made possible by the efficiency brought by unstructured grids. A further increase in resolution, necessary to resolve surf zones and still cover vast regions,will require further developments in numerical methods.     

Fingerprinting feldspar phenocrysts using crystal isotopic composition stratigraphy: implications for crystal transfer and magma mingling in S-type granites

 Microsampling of cm-scale feldspar crystals within an S-type granite from the Lachlan Fold Belt of southeastern Australia has revealed complex internal Sr and Nd isotopic variations. The observed isotopic zonations are in part interpreted as recording feldspar crystallisation in a dynamically mixing magma system, the isotopic composition of which was varying in response to the influx of more mafic and isotopically more mantle-like magmas, the latter stages of which are now represented in modified form by microgranular enclaves. Similar core to rim isotopic variations in feldspar megacrysts from a microgranular enclave and the adjacent host granite strongly suggest megacrysts in the enclave were transferred from the granitic magma during crystallisation. Feldspar rims have higher ⁸⁷Sr/⁸⁶Sr_i and lower ɛ_Nd(i) than adjacent whole rock analyses, but match those of mineral separates from the surrounding enclave matrix. This suggests that the final stages of megacryst growth occurred in the presence of a component that had previously interacted with a high ⁸⁷Sr/⁸⁶Sr, low ɛ_Nd(i) component such as metasedimentary wall rocks. Isotopic heterogeneities are also presererved within different mineral phases in the enclave matrix, suggesting that differing phases grew at differing stages of equilibration between the enclave magma and its host granitic magma. Our results reveal major isotopic heterogeneities on a single crystal and also inter-mineral scale in a pluton which shows well constrained evidence for magma mingling. These results indicate the suitability of feldspars as recorders of isotopic change in magmatic systems, even those which have cooled slowly in the plutonic environment and suggest that much heterogeneity in plutonic systems may be overlooked on a whole rock scale. Received: 28 September 1998 / Accepted: 29 December 1999     

An algorithm to classify and monitor seasonal variations in vegetation phenologies and their inter-annual change

The algorithm presented in this paper classifies vegetation from annual Normalized Difference Vegetation Index (NVDI) time series according to the shape of the temporal cycle. Shape is described using the Fourier components’ magnitude and phase. The degree of an NDVI cycle’s similarity to a predefined reference cycle is measured by the similarity in their amplitude ratios and in their phase differences. Tolerable deviations from the ideal ratios and differences can be set by the user depending on individual accuracy requirements. Tolerable vegetation coverage variation within a shape class is another user defined variable. The algorithm is invariant to cycle modifications including temporal shifts, vertical displacements, and intensity variations, modifications that may be caused by differences in climate, soil (-type, -water, -fertility), or topography, but are unrelated to the vegetation type. The output is a highly consistent clustering of NDVI cycles according to their shapes, which can be linked to distinct vegetation types or land use practices. Intra-class coverage variations in the form of continuous fields measured relative to the reference cycle provide additional information about vegetation covers. Based on the same principles, inter-annual vegetation changes can be monitored with the possibility to distinguish between coverage fluctuations and phenological variations/changes.The algorithms are independent from scene statistics and can be used to create spatially and temporally comparable classifications. Their potential is demonstrated using a 250 m MODIS NDVI time series (version 4) from the Middle East.     

A comparison of parameterized vs. measured transilient mixing coefficients for a convective mixed layer

The parameterization for transilient turbulence coefficients suggested by Stull and Driedonks (SD. 1987) is tested against the large-eddy simulations (LES) of Ebert et al. (ESS, 1989) for the special case of an idealized convective boundary layer. The SD parameterization is based on a nonlocal approximation to the turbulence kinetic energy (TKE) equation, and requires turbulent exchange (i.e., the matrix of transilient mixing coefficients is assumed to be symmetric) and dominance of the smaller eddies (i.e., elements closer to the main diagonal of the matrix are greater). Measurements from the LES model, however, show that the transilient matrix is asymmetric in convective situations, with larger eddies dominating.Mean-state conditions such as the deep convective mixing and mixed-layer growth are satisfactorily described by the parameterization, but the surface layer is too deep and the entrainment zone thickness is poorly defined. Turbulence properties such as skewed vertical velocity distributions are not possible within the constraint of a symmetric matrix, and partial convective overturning is also not possible because of the restriction that small eddies dominate. Future improved parametcrizations might continue to be based on the TKE equation, but should allow transilient matrix asymmetry.Research was performed while the author worked at the Institute of Atmospheric Physics, German Aerospace Research Establishment (DLR), D-8031 Weling. Germany.     

Stand structure,floristic composition and species diversity along altitudinal gradients of a Bornean mountain range 30 years after selective logging

The mountain ecosystem is a hotspot of biodiversity and provides a significant contribution to the humans' livelihoods. However, there is still a lack of information regarding floral diversity and how it has been affected due to a logging operation in Sarawak. A study was conducted to determine vegetation recovery at the Heart of Borneo conservation area at different elevations: 750–1200 m(S1), 1200–1500 m(S2) and 1500 m(S3). A total of 891 trees were enumerated. The trees comprised 55 families, 98 genera and 203 species. The lowest elevation site, S1, contained the highest number of species(117) and families(42), while the highest elevation site, S3, contained the lowest diversity(25 families and 59 species). A similar forest structure was recorded, whereby diameter at breast height(dbh), basal area and height decreased as elevation increased. The most speciose families in S1 were Dipterocarpaceae(11 species) and Lauraceae(10 species), while Lauraceae and Fagaceae were the dominant tree families at S2 and S3. Secondary tree species are still abundantly found in the logged plots, especially at S2(45.5%), indicating that after 30 years, these stands are still in the recovery stages. The logged plots retained typical secondary forest appearances with lower canopy diameters and heights and dbh than the unlogged forest. Although the species richness approached that of the unlogged plots, especially at S1(90.9%) and S3(94.9%), there were still major differences in the floristic compositions between logged and unlogged plots(35.4%, 21.5% and 44.7% for S1, S2 and S3, respectively). Although potential seed sources were in close proximity, succession at S2 has proceeded much slower than at S1 and S3. Lower concentrations of soil nutrients, especially N, may be the main reason for the slow recovery. In conclusion, the results suggested that the selectively logged forest at the Payeh Maga Highland required more than three decades to recover and achieve a climate forest stage.     

Disentangling late-Holocene climate and land use impacts on an Austrian alpine lake using seasonal temperature anomalies,ice-cover,sedimentology, and pollen tracers

Major and trace elements, minerals and grain-size were analysed from a sediment core covering the last 4,000 years of an Austrian Alpine lake (Oberer Landschitzsee, ObLAN, 2,076 m a.s.l.). These analyses were combined with autumn and spring temperature anomalies and ice-cover estimated from a diatom and chrysophyte cyst, thermistor-based regional calibration dataset and selected pollen markers published previously. Diatom-inferred pH (Di-pH) and DOC (Di-DOC) completed the multi-proxy approach, which, together with ordination techniques (PCA and RDA), helped (1) to confirm proposed climatic patterns and hypotheses, and (2) to disentangle the complex interactions between climatic and anthropogenic impacts. Shifts in chemical and physical weathering, erosion, production, lake stratification, redox potentials and air pollution were the major processes that have affected Oberer Landschitzsee in relation to climate, lake/catchment interactions, and human impact. Geochemistry supported the finding that the four waves of high-altitude land use (Early to Middle Bronze, ca. 1,800–1,300 B.C.; Late Bronze to Hallstatt, ca. 1,000–500 B.C.; Celtic to Roman, ca. 300 B.C.–400 A.D.; and Medieval, ca. 1,000–1,600 A.D.) were coupled mainly with warm periods. Increased production, onset of lake stratification, and the deterioration of hypolimnetic oxygen conditions were the major changes governing in-lake processes during climate warming. They resulted in specific element assemblages associated with organic matter accumulation (Br), oxygen depletion (As, Ga, Cu, S), and changes in redox (Fe/Mn). The Se/Di-DOC ratio was introduced to track shifts between in-lake production and allochthonous sources. Nutrient loading from pastures, coupled with climate warming, could explain that in-lake production was highest during Roman and Medieval land use. Lithogenic elements mainly originated from chemical weathering of silicate bedrock and they increased when intense land-use was coupled with climate deterioration and/or increased humidity. These perturbations were highest during a High Medieval climate fluctuation around 1,000 A.D. The association of sand with LOI and C/N and the decoupling of sand from quartz and feldspar separates erosion from physical weathering. Di-DOC, S, and C/N showed hybrid characters in relation to climate and human impact. The mineral proportions indicated gradients in relation to weathering, snow-cover and running waters, as well as vegetation. Air pollution by metallurgic industries, starting during the High Medieval and culminating during the Late Medieval, caused lead accumulation and could have contributed to As and S enrichment. Corresponding features in the distribution of selected elements and pollen tracers, as well as changes in mineral proportions, supported the hypotheses of shifts in seasonal climate and an overall trend towards more continental climate conditions since Medieval times.     

A lost Tethyan evaporitic basin: Evidence from a Cretaceous hemipelagic meta‐selenite – red chert association in the Eastern Mediterranean realm

Ancient evaporite deposits are geological archives of depositional environments characterized by a long‐term negative precipitation balance and bear evidence for global ocean element mass balance calculations. Here, Cretaceous selenite pseudomorphs from western Anatolia (‘Rosetta Marble’) — characterized by their exceptional morphological preservation — and their ‘marine’ geochemical signatures are described and interpreted in a process‐oriented context. These rocks recorded Late Cretaceous high‐pressure/low‐temperature, subduction‐related metamorphism with peak conditions of 1·0 to 1·2 GPa and 300 to 400°C. Metre‐scale, rock‐forming radiating rods, now present as fibrous calcite marble, clearly point to selenitic gypsum as the precursor mineral. Stratigraphic successions are recorded along a reconstructed proximal to distal transect. The cyclical alternation of selenite beds and radiolarian ribbon‐bedded cherts in the distal portions are interpreted as a two type of seawater system. During arid intervals, shallow marine brines cascaded downward into basinal settings and induced precipitation. During more humid times, upwelling‐induced radiolarian blooms caused the deposition of radiolarite facies. Interestingly, there is no comparable depositional setting known from the Cenozoic world. Meta‐selenite geochemical data (δ¹³C, δ¹⁸O and ⁸⁷Sr/⁸⁶Sr) plot within the range of reconstructed middle Cretaceous seawater signatures. Possible sources for the ¹³C‐enriched (mean 2·2‰) values include methanogenesis, gas hydrates and cold seep fluid exhalation. Spatially resolved component‐specific analysis of a rock slab displays isotopic variances between meta‐selenite crystals (mean δ¹³C 2·2‰) and host matrix (mean δ¹³C 1·3‰). The Cretaceous evaporite‐pseudomorphs of Anatolia represent a basin wide event coeval with the Aptian evaporites of the Proto‐Atlantic and the pseudomorphs share many attributes, including lateral distribution of 600 km and stratigraphic thickness of 1·5 to 2·0 km, with the evaporites formed during the younger Messinian salinity crisis. The Rosetta Marble of Anatolia may represent the best‐preserved selenite pseudomorphs worldwide and have a clear potential to act as a template for the study of meta‐selenite in deep time.