Tracking Holocene environmental changes in an alpine lake sediment core: application of regional diatom calibration, geochemistry, and pollen

Diatom abundances in the surface sediment samples of 41 mountain lakes in the central Austrian Alps (Niedere Tauern) were related to environmental variables using multi-variate techniques. Canonical correspondence analysis (CCA) revealed that the pH, date of autumn mixing (A _mix), mean August water temperature (T _Aug), dissolved organic carbon (DOC), and relative water depth (Z _rel) made significant contributions to explain the diatom assemblage variation in the lakes of the training set. A weighted averaging partial least square regression and calibration model was used to establish Di-pH (R ² _boot= 0.72, RMSEP_boot= 0.131), and a thermistor measurements-based PLS model for A _mix (R ² _boot= 0.71, RMSEP_boot= 0.006 log₁₀ Julian days). The latter showed a better prediction than T _Aug, and was used in terms of climate change. These transfer functions, together with analyses of loss on ignition (LOI), the total carbon/nitrogen (C/N)-ratios, and selected pollen, were applied to an early to mid-Holocene (11.5–4 cal. ky BP) sediment core section from an Austrian Alpine treeline lake on crystalline bedrock. Additionally, passive sample scores in the CCA of the diatom training set were used to show trends in the variables DOC and Z _rel. During the early Holocene, diatoms indicative of increased pH, extended warm summers, and low water levels dominated. Between 10.2 and 7.6 cal. ky BP it was followed by diatom assemblages that indicated an increase in lake water depth and an earlier A _mix. The multi-proxy data suggest that the A _mix decline is the result of a series of snow-rich summer cool and wet climate fluctuations, which were divided by climate warming at ∼9 cal. ky BP. Increased A _mix, LOI and DOC, and the correspondent decline in the C/N-ratios, show subsequent climate warming between 7.3 and 6 cal. ky BP. The long-term trend in Di-pH indicates the impact of catchment-related processes during the early-Holocene, that were superimposed by climate.     

Coupling of early diagenetic processes and sedimentary dynamics in tropical shelf environments: the Gulf of Papua deltaic complex

Tropical mobile mud belts represent a major class of biogeochemical and diagenetic systems characterized by extensive and frequent physical reworking of fine-grained, organic-rich deposits underlying oxygenated waters. Large regions of the Gulf of Papua, Papua New Guinea deltaic complex are dominated by such conditions. A reworked mud belt lies within the inner shelf between 10 and 20 m depth on a sedimentary clinoform derived from coalescing deltas. Deposits across the topset are typically suboxic, nonsulfidic over the upper 0.5–1 m, and have low to moderate maximum pore water concentrations of dissolved Fe(II) and Mn(II) (100–200, but up to 800 μM). Sediments are reactive, with surficial ΣCO₂ production 0.1–0.3 mM d⁻¹ and benthic O₂ fluxes 23±15 mmol m⁻² d⁻¹ (upper 20 cm). The highest rates occur within inner topset deposits (10–20 m) and near the high accumulation rollover region of the topset–foreset beds (40–50 m). Lower rates are found inshore along intertidal channels—mangrove fringe and within scoured or exposed consolidated deposits of the middle topset region. Remineralization rate patterns are independent of relative dominance by terrestrial or marine carbon in sediments. Dissolved O₂ usually penetrates 2–5 mm into surface sediments when macrofaunal burrows are absent. More than 75% of the highly reactive sedimentary Fe(III) pool (350–400 μmol g⁻¹) is typically diagenetically reduced in the upper 0.5 m. Pore water can be measureably depleted at depths >0.5 m, but dissolved H₂S generally remains below detection over the upper 1–2 m. As in other deltaic topset regions, concentration gradients often indicate that compared to many marine deposits of similar sediment accumulation rates, relatively refractory C_org is supplied to the SO₄ reducing zone. Sedimentary C/S ratios are 4–6 within the suboxic topset regions but decrease to <3 in offshore foreset beds where sulfidic diagenesis dominates. Only 15–20% of the diagenetically reduced Fe(II) is pyritic and a maximum of 10–25% is carbonate, implying that most Fe(II) is associated with authigenic or lithogenic silicates or oxides. The dominance of suboxic, nonsulfidic diagenetic processes reflect coupling between delivery of oxide-rich terrestrial debris, remobilization and reoxidation of deposits, and repetitive entrainment/remineralization of both labile and refractory organics. Distinct sedimentary indicators of reactive, suboxic mobile mud belts within tropical climatic zones are: abundant total highly reactive Fe (ΣFeR )>300 μmol g⁻¹; most reactive Fe is diagenetically reduced (ΣFe(II)/ΣFeR0.7–0.8); the proportion of diagenetically reduced Fe present as pyrite is low (Py–Fe(II)<0.2); C/S 4–8; and C_org/particle surface area <0.4 (mg C m⁻²). These depositional environments must be most common in tropical climates during high sea stand.     

Coastal cliff retreat rates at Beit‐Yannay,Israel, in the 20th century

Research indicates that the aeolianite (Kurkar) cliffs along the Israeli Mediterranean coastline have continuously retreated eastward during the last few decades. There seems to be no dispute among Earth scientists regarding the general trend of cliff retreat. However the majority of papers displaying cliff retreat rates are based upon comparison of aerial photographs. Their lack of advanced geometric measurement methods causes a high margin of error. Public attention is focused upon the Beit‐Yannay coastal cliff since private homes are located along the southern section of the cliff crest. The current research compares the historic location of the cliff crest edge at Beit‐Yannay as observed in a series of aerial photographs taken during the period 1918–2000. Quantitative measurement methods included applications of satellite geodesy and digital photogrammetry and mapping. Research results offer quantitative, consecutive and highly accurate data regarding retreat rates over a relatively long period of 82 years. It is concluded that: 1. Annual average cliff retreat rates of the cliff crest is 20 cm/year. 2. Categorization of the study time span reveals periods displaying varying retreat rates such as 27 cm/year during 1918–1946, 21 cm/year during 1946–1973 and 10 cm/year during 1973–2000. 3. Maximum retreat distances of the cliff crest, over the study period were found to be approximately 25 m along the northern, lowest section of the cliff. Minimum distances of 11 m were identi?ed at the highest, southern section of the cliff. 4. The eolianite (Kurkar) cliffs along the Israeli Mediterranean coast throughout the 20th century have been an important source of sediment, contributing approximately 24 × 10⁶ m³ of sediments to the sediment balance of Israeli beaches. Copyright © 2004 John Wiley & Sons, Ltd.     

Impacts of Spatial Partitioning in Hydroecological Models: Predicting Grassland Productivity with RHESSys

Environmental models constructed with a spatial domain require choices about the representation of space. Decisions in the adaptation of a spatial data model can have significant consequences on the ability to predict environmental function as a result of changes to levels of aggregation of input parameters and scaling issues in the processes being modelled. In some cases, it is possible to construct a systematic framework to evaluate the uncertainty in predictions using different spatial models; in other cases, the realm of possibilities plus the complexity of the environmental model in question may inhibit numeric uncertainty estimates. We demonstrate a range of potential spatial data models to parameterize a landscape‐level hydroecological model (RHESSys). The effects of data model choice are illustrated, both in terms of input parameter distributions and resulting ecophysiological predictions. Predicted productivity varied widely, as a function of both the number of modelling units, and of arbitrary decisions such as the origin of a raster grid. It is therefore important to use as much information about the modelled environment as possible. Combinations of adaptive methods to evaluate distributions of input data, plus knowledge of dominant controls of ecosystem processes, can help evaluate potential representations. In this case, variance‐based delineation of vegetation patches is shown to improve the ability to intelligently choose a patch distribution that minimizes the number of patches, while maintaining a degree of aggregation that does not overly bias the predictions.     

The geology and Re–Os geochronology of the Palaeoproterozoic Vaikijaur Cu–Au–(Mo) porphyry-style deposit in the Jokkmokk granitoid, northern Sweden

The Vaikijaur Cu–Au–(Mo) deposit is located in the ca. 1.88 Ga calc-alkaline Jokkmokk granitoid near the Archaean–Proterozoic palaeoboundary within the Fennoscandian shield of northern Sweden. The Skellefte VMS district lies immediately to the south and the northern Norrbotten Fe-oxide–Cu–Au deposits to the north. The Vaikijaur deposit occupies an area of 2×3 km within the Jokkmokk granitoid and includes stockwork quartz-sulphide veinlets and disseminated chalcopyrite, pyrite, gold, molybdenite, magnetite, and pyrrhotite. Porphyritic mafic dykes were emplaced along fractures in a ring dyke pattern. The Jokkmokk granitoid, dykes, and the mineralized area are foliated, indicating that mineralization predated the main regional deformation. The mineralized area is characterized by strong potassic alteration. Phyllic and propylitic alteration zones are also present. A pyrite-rich inner core is surrounded by a concentric zone with pyrite, chalcopyrite, and gold. Molybdenite is distributed irregularly throughout the chalcopyrite zone. Geophysical data indicate a strongly conductive central zone in the mineralized area bordered by conductive and high magnetic zones. Five high precision Re–Os age determinations for three molybdenite occurrences from outcrop and drill core samples constrain the age of porphyry-style Cu–Au–(Mo) mineralization to between 1889±10 and 1868±6 Ma. A younger molybdenite is associated with a much later metamorphic event at about 1750 Ma. These data suggest that primary porphyry-style mineralization was associated with calc-alkaline magmatism within the Archaean–Proterozoic boundary zone at ca. 1.89–1.87 Ga.     

Monitoring changes in shoreline position adjacent to the Hadera power station, Israel

The paper focuses on the shoreline changes following the construction of a power station cooling harbor in Hadera, Israel, and uses a statistical analysis method capable of mitigating the effect of short-term changes on the measurement of long-term changes. The changes in shoreline position detected 25 years after the harbor's construction were compared with forecasts of the Central Lab for Hydraulics of France (LCHF) in 1976, which correctly predicted the effects of the construction on the nearby coast and with studies conducted a few years after the construction. A comparison between later studies and the study described herein shows a spatial effect of the harbor on shoreline position of 750 m north of the northern breakwater. No effect was found south of the harbor. The harbor's temporal effects extend over a few years from the beginning of the construction until the coast reaches a new sedimentary equilibrium. The statistical analysis method employed was found useful for the detection of long-term changes in shoreline position such as those caused by the construction of coastal structures; it will not be suitable for the detection of long-term changes resulting from a rise in sea-level.     

Fluid inclusion and stable isotope (O, H, C, and S) constraints on the genesis of the Serrinha gold deposit, Gurupi Belt, northern Brazil

The Serrinha gold deposit of the Gurupi Belt, northern Brazil, belongs to the class of orogenic gold deposits. The deposit is hosted in highly strained graphitic schist belonging to a Paleoproterozoic (∼2,160 Ma) metavolcano-sedimentary sequence. The ore-zones are up to 11 m thick, parallel to the regional NW–SE schistosity, and characterized by quartz-carbonate-sulfide veinlets and minor disseminations. Textural and structural data indicate that mineralization was syn- to late-tectonic and postmetamorphic. Fluid inclusion studies identified early CO₂ (CH₄-N₂) and CO₂ (CH₄-N₂)-H₂O-NaCl inclusions that show highly variable phase ratios, CO₂ homogenization, and total homogenization temperatures both to liquid and vapor, interpreted as the product of fluid immiscibility under fluctuating pressure conditions, more or less associated with postentrapment modifications. The ore-bearing fluid typically has 18–33mol% of CO₂, up to 4mol% of N₂, and less than 2mol% of CH₄ and displays moderate to high densities with salinity around 4.5wt% NaCl equiv. Mineralization occurred around 310 to 335°C and 1.3 to 3.0 kbar, based on fluid inclusion homogenization temperatures and oxygen isotope thermometry with estimated oxygen fugacity indicating relatively reduced conditions. Stable isotope data on quartz, carbonate, and fluid inclusions suggest that veins formed from fluids with δ¹⁸O_H2O and δD_H2O (310–335°C) values of +6.2 to +8.4‰ and −19 to −80‰, respectively, which might be metamorphic and/or magmatic and/or mantle-derived. The carbon isotope composition (δ¹³C) varies from −14.2 to −15.7‰ in carbonates; it is −17.6‰ in fluid inclusion CO₂ and −23.6‰ in graphite from the host rock. The δ³⁴S values of pyrite are −2.6 to −7.9‰. The strongly to moderately negative carbon isotope composition of the carbonates and inclusion fluid CO₂ reflects variable contribution of organic carbon to an originally heavier fluid (magmatic, metamorphic, or mantle-derived) at the site of deposition and sulfur isotopes indicate some oxidation of the originally reduced fluid. The deposition of gold is interpreted to have occurred mainly in response to phase separation and fluid-rock interactions such as CO₂ removal and desulfidation reactions that provoked variations in the fluid pH and redox conditions.     

Regression model for generating time series of daily precipitation amounts for climate change impact studies

The precipitation amounts on wet days at De Bilt (the Netherlands) are linked to temperature and surface air pressure through advanced regression techniques. Temperature is chosen as a covariate to use the model for generating synthetic time series of daily precipitation in a CO₂ induced warmer climate. The precipitation-temperature dependence can partly be ascribed to the phenomenon that warmer air can contain more moisture. Spline functions are introduced to reproduce the non-monotonous change of the mean daily precipitation amount with temperature. Because the model is non-linear and the variance of the errors depends on the expected response, an iteratively reweighted least-squares technique is needed to estimate the regression coefficients. A representative rainfall sequence for the situation of a systematic temperature rise is obtained by multiplying the precipitation amounts in the observed record with a temperature dependent factor based on a fitted regression model. For a temperature change of 3°C (reasonable guess for a doubled CO₂ climate according to the present-day general circulation models) this results in an increase in the annual average amount of 9% (20% in winter and 4% in summer). An extended model with both temperature and surface air pressure is presented which makes it possible to study the additional effects of a potential systematic change in surface air pressure on precipitation.     

Inverse theory in the earth sciences—an introductory overview with emphasis on Gandin's method of optimum interpolation

In 1963, Gandin published a monograph on “optimum interpolation for the objective analysis of meteorological fields, ? a method that is similar mathematically to geodetical least-squares prediction and collocation, simple kriging, and spectral interpolation. The common problem is the interpolation or extrapolation or estimation of a continuous spatial property from finitely many observations. Gandin 's method is presented in an inverse-theoretical context with focus on a methodological comparison with related methods. Underlying mathematical assumptions as well as geological implications are discussed. An introductory overview of inverse methods in the earth sciences is given, with emphasis on methods with a structure analysis step.     

Simple Temperature Scenario for a Gulf Stream Induced Climate Change

Consequences of a Gulf Stream induced ocean surface cooling for the temperature climate of Western Europe were studied by means of a conditional perturbation of the observed daily temperature time series of the Netherlands. On days with advection of airmasses of maritime origin, the observed temperatures in the series were lowered with a fixed value, representing the influence of a cooler Atlantic Ocean. On the other days, the observed temperatures were left unchanged. The perturbation results in a decrease in the mean temperature that is almost constant over the year, and in a change in the standard deviation of the daily temperatures that is seasonally dependent. Due to preferential cooling of warm winter days, the standard deviation decreases in the winter, whereas in the other seasons the standard deviation increases as a result of preferential cooling of days with low temperatures. Although this ocean cooling scenario indicates an increase of the relative frequency of cold winters and cool summers, it is neither characterized by the occurrence of winters with unprecedented low temperatures nor by the disappearance of summer heatwaves.