Palaeolimnology of Lake Piene-Kuuppalanlampi (Kurkijoki, Karelian Republic, Russia): isolation history, lake ecosystem development and long-term agricultural impact

The post-glacial history of Lake Pieni-Kuuppalanlampi, western Lake Ladoga region, was studied by means of stratigraphic pollen and diatom analyses. Diatoms were analysed to track the isolation history of the basin and the limnological effects of the early land-use phases indicated by pollen analysis. Chrysophycean stomatocysts and Isoëtes spores were also employed in the limnological reconstructions. Sediment dating was provided by six conventional radiocarbon dates.The lower part of the 370-cm long sediment sequence represents early Holocene, large lake conditions: the (freshwater) Yoldia and Ancylus stages of the Baltic basin, with a short-term lagoonal or isolation phase at the end of Yoldia. The basin was isolated due to Ancylus lake regression at 9785 cal B.P. For the small-lake sequence (0–250 cm) we used diatom inferences (WA-method) for hindcasting water chemistry. The post-isolation limnology of Pieni-Kuuppalanlampi reflects the development of vegetation on its small hill-top catchment. In its early development, the lake was mesotrophic, but became more acidic (pH about 6) and oligotrophic after the decline of temperate hardwood trees and the spread of spruce to the area after around 5000 B.P. The lake ecosystem appears to have reacted sensitively to agricultural land-use in the catchment from AD 400–800 onwards: inferred total phosphorus levels and pH both rise during these periods.     

The boundary integral calculations of the forward problem for D.C. sounding and MMR methods for a 3-D body near a vertical contact

Summary The paper presents comprehensive theory based on the boundary integral method for calculations of the electric potential, electric field and corresponding magnetic field due to a pair of D.C. source electrodes near a vertical resistivity contact in the halfspace, indlucing a 3-D disturbing body in the vicinity of the contact. Special attention is paid to the case when the disturbing body touches the vertical contact. Results of numerical calculations are presented in the form of sounding curves and a set of isoline graphs for potential, components of the electric and magnetic field (total and anomalous) on the surface of the Earth. It is shown that the presence of the disturbing body at the contact is most pronounced in the electrical characteristics. Anomalies in the magnetic field are small in comparison to the field due to the electric current in the electrode cable and primary currents flowing from the electrodes.     

Severe weather affecting European transport systems: the identification,classification and frequencies of events

Severe weather can have serious repercussions in the transport sector as a whole by increasing the number of accidents, injuries and other damage, as well as leading to highly increased travel times. This study, a component of the EU FP7 Project EWENT, delineates a Europe-wide climatology of adverse and extreme weather events that can be expected to affect the transport network. We first define and classify the relevant severe weather events by investigating the effects of hazardous conditions on different transportation modes and the infrastructure. Consideration is given to individual phenomena such as snowfall, heavy precipitation, heat waves, cold spells, wind gusts; a combined phenomenon, the blizzard, is also considered. The frequency of severe weather events, together with the changes in their spatial extension and intensity, is analyzed based on the E-OBS dataset (1971–2000) and the ERA-Interim reanalysis dataset (1989–2010). Northern Europe and the Alpine region are the areas most impacted by winter extremes, such as snowfall, cold spells and winter storms, the frequency of heavy snowfall. The frequency of hot days is highest in Southern Europe. Severe winds and blizzards are the most common over the Atlantic and along its shores. Although heavy rainfall may affect the whole continent on an annual basis, extreme precipitation events are relative sparse, affecting particularly the Alps and the Atlantic coastline. A European regionalization covering similar impacts on the transport network is performed.     

Estimation of pH optima and tolerances of diatoms in lake sediments by the methods of weighted averaging,least squares and maximum likelihood,and their use for the prediction of lake acidity

Ecological optima and tolerances with respect to autumn pH were estimated for 63 diatom taxa in 47 Finnish lakes. The methods used were weighted averaging (WA), least squares (LS) and maximum likelihood (ML), the two latter methods assuming the Gaussian response model.WA produces optimum estimates which are necessarily within the observed lake pH range, whereas there is no such restriction in ML and LS. When the most extreme estimates of ML and LS were excluded, a reasonably close agreement among the results of different estimation methods was observed. When the species with unrealistic optima were excluded, the tolerance estimates were also rather similar, although the ML estimates were systematically greater.The parameter estimates were used to predict the autumn pH of 34 other lakes by weighted averaging. The ML and LS estimates including the extreme optima produced inferior predictions. A good prediction was obtained, however, when prediction with these estimates was additionally scaled with inverse squared tolerances, or when the extreme values were removed (censored). Tolerance downweighting was perhaps more efficient, and when it was used, no additional improvement was gained by censoring. The WA estimates produced good predictions without any manipulations, but these predictions tended to be biased towards the centroid of the observed range of pH values.At best, the average bias in prediction, as measured by mean difference between predicted and observed pH, was 0.082 pH units and the standard deviation of the differences, measuring the average random prediction error, was 0.256 pH units.     

Subarctic catchment water storage and carbon cycling – Leading the way for future studies using integrated datasets at Pallas,Finland

Subarctic ecohydrological processes are changing rapidly, but detailed and integrated ecohydrological investigations are not as widespread as necessary. We introduce an integrated research catchment site (Pallas) for atmosphere, ecosystems, and ecohydrology studies in subarctic conditions in Finland that can be used for a new set of comparative catchment investigations. The Pallas site provides unique observational data and high-intensity field measurement datasets over long periods. The infrastructure for atmosphere- to landscape-scale research in ecosystem processes in a subarctic landscape has recently been complemented with detailed ecohydrological measurements. We identify three dominant processes in subarctic ecohydrology: (a) strong seasonality drives ecohydrological regimes, (b) limited dynamic storage causes rapid stream response to water inputs (snowmelt and intensive storms), and (c) hydrological state of the system regulates catchment-scale dissolved carbon dynamics and greenhouse (GHG) fluxes. Surface water and groundwater interactions play an important role in regulating catchment-scale carbon balances and ecosystem respiration within subarctic peatlands, particularly their spatial variability in the landscape. Based on our observations from Pallas, we highlight key research gaps in subarctic ecohydrology and propose several ways forward. We also demonstrate that the Pallas catchment meets the need for sustaining and pushing the boundaries of critical long-term integrated ecohydrological research in high-latitude environments.     

Micrometeorological Observations of a Microburst in Southern Finland

On the afternoon of 3 July 2004 in Hyytiälä (Juupajoki, Finland), convective cells produced a strong downburst causing forest damage. The SMEAR II field station, situated near the damage site, enabled a unique micrometeorological analysis of a microburst with differences above and inside the canopy. At the time of the event, a squall line associated with a cold front was crossing Hyytiälä with a reflectivity maximum in the middle of the squall line. A bow echo, rear-inflow notch, and probable mesovortex were observed in radar data. The bow echo moved west-north-west, and its apex travelled just north of Hyytiälä. The turbulence data were analysed at two locations above the forest canopy and at one location at sub-canopy. At 1412 EET (Eastern European Time, UTC+2), the horizontal and vertical wind speed increased and the wind veered, reflecting the arrival of a gust front. At the same time, the carbon dioxide concentration increased due to turbulent mixing, the temperature decreased due to cold air flow from aloft and aerosol particle concentration decreased due to rain scavenging. An increase in the number concentration of ultra-fine particles (< 10 nm) was detected, supporting the new particle formation either from cloud outflow or due to rain. Five minutes after the gust front (1417 EET), strong horizontal and downward vertical wind speed gusts occurred with maxima of 22 and 15 m s^?1, respectively, reflecting the microburst. The turbulence spectra before, during and after the event were consistent with traditional turbulence spectral theory.     

Testing a spatially distributed tracer‐aided runoff model in a snow‐influenced catchment: Effects of multicriteria calibration on streamwater ages

Integrating stable isotope tracers into rainfall‐runoff models allows investigation of water partitioning and direct estimation of travel times and water ages. Tracer data have valuable information content that can be used to constrain models and, in integration with hydrometric observations, test the conceptualization of catchment processes in model structure and parameterization. There is great potential in using tracer‐aided modelling in snow‐influenced catchments to improve understanding of these catchments' dynamics and sensitivity to environmental change. We used the spatially distributed tracer‐aided rainfall‐runoff (STARR) model to simulate the interactions between water storage, flux, and isotope dynamics in a snow‐influenced, long‐term monitored catchment in Ontario, Canada. Multiple realizations of the model were achieved using a combination of single and multiple objectives as calibration targets. Although good simulations of hydrometric targets such as discharge and snow water equivalent could be achieved by local calibration alone, adequate capture of the stream isotope dynamics was predicated on the inclusion of isotope data in the calibration. Parameter sensitivity was highest, and most local, for single calibration targets. With multiple calibration targets, key sensitive parameters were still identifiable in snow and runoff generation routines. Water ages derived from flux tracking subroutines in the model indicated a catchment where runoff is dominated by younger waters, particularly during spring snowmelt. However, resulting water ages were most sensitive to the partitioning of runoff sources from soil and groundwater sources, which was most realistically achieved when isotopes were included in the calibration. Given the paucity of studies where hydrological models explicitly incorporate tracers in snow‐influenced regions, this study using STARR is an important contribution to satisfactorily simulating snowpack dynamics and runoff generation processes, while simultaneously capturing stable isotope variability in snow‐influenced catchments.     

Two-dimensional non-linear inversion of VLF-R data using simulated annealing

The VLF-R (very low frequency-resistivity) data, i.e. the apparent resistivity ( ρ _a ) and phase ( φ ) data, were inverted individually and jointly using the VFSA (very fast simulated annealing) global inversion approach. Global inversion results for synthetic data without and with various amounts of random and normally distributed Gaussian noise reveal that the inversion of neither the ρ _a nor φ data alone yields the true parameters of the structures. However, the joint inversion of the ρ _a and φ data yields very good estimates of the model parameters. Five models, representing typical subsurface structures in the shield areas, are studied here. Various models achieved after 10 VFSA runs were used to compute the mean model and the corresponding covariance and correlation matrices, which were used to estimate the uncertainties in the mean model parameters and correlations between the model parameters. We observe that these correlations follow the physics associated with the problem. VLF-R field data due to a nearly vertical contact structure and a very thick dyke-like structure were also inverted to demonstrate the efficacy of the approach in the delineation of the parameters of 2-D structures.