Pollution history from 256 BC to AD 2005 inferred from the accumulation of elements in a varve record of Lake Korttajärvi in Finland     

Jarmo J. Meriläinen  Virve Kustula  Allan Witick  Eeva Haltia-Hovi  Timo Saarinen 《Journal of Paleolimnology》2010,44(2):531-545

In order to assess the importance of different pollution sources for the increase in element concentration and accumulation, historical changes in selected elements were studied in the annually laminated sediment of Lake Korttajärvi in Central Finland (62°20′N; 25°41′E). The sediment chronology based on varve counting (256 BC to AD 2005) provided a unique opportunity to explore and date signals of metal emissions, including the ancient metallurgical activities of the Roman Empire at the beginning of the Current Era. Records of this kind are mostly lacking in Finland and northernmost Europe. The stratigraphic sequence of element concentrations did not reflect any major changes in the lake, but changes in element accumulation rates provided distinct pollution signals caused by airborne fallout, catchment erosion, and to some extent municipal wastewater loading. The maximum bulk sedimentation recorded in the twentieth century was 11-fold and organic sedimentation 4-fold higher than the mean background sedimentation rate (256 BC to AD 1019). The increase in the accumulation rates of the majority of the elements, such as Cd, Sn, Pb, Si, Ni, B, Cu, Zn, Sr, Na, K, Sb, Ca, Cr, U and Mg, in descending order, was at least equal to that of bulk sedimentation or much greater, especially for Cd, Sn, and Pb. Changes in the accumulation of Co, Fe, Mn, Mo and As were small and mainly followed those of organic sedimentation. The earliest pollution signals were those of Pb recorded in AD 1055–1141. A weak signal of Pb pollution from the Roman Era was detected in metal concentrations, but this could not be confirmed by the accumulation rate data for Pb.  相似文献   

Geophysics of sea ice in the Baltic Sea: A review     

Timo Vihma  Jari Haapala   《Progress in Oceanography》2009,80(3-4):129-148

With improved observation methods, increased winter navigation, and increased awareness of the climate and environmental changes, research on the Baltic Sea ice conditions has become increasingly active. Sea ice has been recognized as a sensitive indicator for changes in climate. Although the inter-annual variability in the ice conditions is large, a change towards milder ice winters has been detected from the time series of the maximum annual extent of sea ice and the length of the ice season. On the basis of the ice extent, the shift towards a warmer climate took place in the latter half of the 19th century. On the other hand, data on the ice thickness, which are mostly limited to the land-fast ice zone, basically do not show clear trends during the 20th century, except that during the last 20 years the thickness of land-fast ice has decreased. Due to difficulties in measuring the pack-ice thickness, the total mass of sea ice in the Baltic Sea is, however, still poorly known. The ice extent and length of the ice season depend on the indices of the Arctic Oscillation and North Atlantic Oscillation. Sea ice dynamics, thermodynamics, structure, and properties strongly interact with each other, as well as with the atmosphere and the sea. The surface conditions over the ice-covered Baltic Sea show high spatial variability, which cannot be described by two surface types (such as ice and open water) only. The variability is strongly reflected to the radiative and turbulent surface fluxes. The Baltic Sea has served as a testbed for several developments in the theory of sea ice dynamics. Experiences with advanced models have increased our understanding on sea ice dynamics, which depends on the ice thickness distribution, and in turn redistributes the ice thickness. During the latest decade, advance has been made in studies on sea ice structure, surface albedo, penetration of solar radiation, sub-surface melting, and formation of superimposed ice and snow ice. A high vertical resolution has been found as a prerequisite to successfully model thermodynamic processes during the spring melt period. A few observations have demonstrated how the river discharge and ice melt affect the stratification of the oceanic boundary layer below the ice and the oceanic heat flux to the ice bottom. In general, process studies on ice–ocean interaction have been rare. In the future, increasingly multidisciplinary studies are needed with close links between sea ice physics, geochemistry and biology.  相似文献   

Autoregressive filtering versus linear detrending in estimation of fluxes by the eddy covariance method   总被引：10，自引：7，他引：3  

Üllar Rannik  Timo Vesala 《Boundary-Layer Meteorology》1999,91(2):259-280

The application of autoregressive running mean filtering (RMF) and linear detrending (LDT) in the estimation of turbulent fluxes by the eddy covariance method is analysed. The systematic, as well as the random, errors of the fluxes arising from filtering and/or limited observation time effects are described. To observe negligible systematic errors in fluxes, the RMF has to be applied with moderately long time constants. However, the obtained flux values are subject to increased random errors during periods of non-stationarity and the method leads to systematic overestimation of variances. These shortcomings are not inherent in the LDT approach, which is recommended for use. But the systematic errors of fluxes due to LDT are not negligible under certain experimental conditions and have to be accounted for. The corrections are important because the relatively small errors in short-period fluxes can translate to significant errors in long-period averages. The corrections depend on the turbulence time scales, which should be preferably estimated as ensemble mean variables for a particular site.  相似文献   

Locating regional seismic events with global optimization based on interval arithmetic   总被引：4，自引：0，他引：4  

Matti Tarvainen  Timo Tiira  & Eystein S. Husebye 《Geophysical Journal International》1999,138(3):879-885

More than 1000 seismic events in northern Europe at distances of up to 400  km from the detecting network are located using an optimization method in which the global minimum of the traveltime function residuals is searched for using an Interval Arithmetic (IA) method. Epicentres are determined using P waves detected by the Finnish national seismic network: up to 15 stations were used in the analysis. The IA results coincide with locations provided by the University of Helsinki bulletins with a median location bias of 7.6  km.
  A second data set of 59 explosions in the Siilinjärvi mine in central Finland was examined in detail, because the locations of the explosions were known exactly. In this case, the median difference of IA locations was 3.8  km from the average location of mine explosions, while all 59 events were located within 9  km of the 'true' epicentres. The corresponding median error of the University of Helsinki locations was smaller (3.2  km), but some Helsinki locations were well over 10  km from the mine. The convergence towards the global optimum using interval arithmetic was fast when compared with the conventional least-squares approaches for epicentre determinations.  相似文献   

Distribution of nitrate over sea-salt and soil derived particles — Implications from a field study   总被引：1，自引：0，他引：1  

Tuomo A. Pakkanen  Veli-Matti Kerminen  Risto E. Hillamo  Minna Màkinen  Timo Màkelà  Aki Virkkula 《Journal of Atmospheric Chemistry》1996,24(2):189-205

A method for determining the distribution of supermicrometer nitrate between size-segregated sea-salt and soil derived particles is presented. The analysis is based on field data from six measurements at a coastal site in southern Finland, and on a theoretical treatment taking into account the transfer of gaseous species onto particle surfaces and their subsequent reaction. Significant amounts of nitrate were found in both the particle types, with the fraction of nitrate associated with soil particles varying from 20–50% in the 1–2 m size to near 90% in particles larger than 10 m. Overall, the nitrate accumulation followed closely the relative abundances of these two particle types. Two overlapping modes in supermicron nitrate mass size distributions could be identified. The lower mode, associated with sea-salt, was located between the surface-area and volume distribution of sodium peaking at about 2–3 m of EAD. The upper mode peaked at 3–5 m and followed more closely the surface-area distribution of calcium in all samples. At our site, the accumulation of nitrate into both particle types was shown to be limited by an effective surface reaction rate rather than by gas-phase diffusion. This rate was estimated to be considerably larger for sea-salt particles. Strong evidence in support of the saturation of nitrate in sea-salt particles were obtained.  相似文献   

Rate-Dependent Embedded Discontinuity Approach Incorporating Heterogeneity for Numerical Modeling of Rock Fracture     

Timo Saksala 《Rock Mechanics and Rock Engineering》2015,48(4):1605-1622

  相似文献   

Damage–viscoplastic consistency model with a parabolic cap for rocks with brittle and ductile behavior under low‐velocity impact loading     

Timo Saksala 《国际地质力学数值与分析法杂志》2010,34(13):1362-1386

This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Mineralogical and strontium isotopic record of hydrothermal processes in the lower ocean crust at and near the East Pacific Rise     

Timo?M.?Kirchner  Kathryn?M.?GillisEmail author 《Contributions to Mineralogy and Petrology》2012,164(1):123-141

Tectonic exposures of upper plutonics (>800 m) that are part of a contiguous section of young East Pacific Rise (EPR) crust at the Hess Deep Rift provide the first regional-scale constraints on hydrothermal processes in the upper plutonic crust formed at a fast-spreading ridge. Submersible-collected samples recovered over a 4-km-wide region show that the sheeted dike complex is largely underlain by a 150- to 200-m-thick gabbro unit, followed by a more primitive gabbronorite unit. Gabbroic samples are variably altered by pervasive fluid flow along fracture networks to amphibole-dominated assemblages. The gabbroic rocks are significantly less altered (average 11% hydrous phases) than the overlying sheeted dike complex (average 24%), and the percentage of hydrous alteration diminishes with depth. Incipient, pervasive fluid flow occurred at amphibolite facies conditions (average 720°C), with slightly higher temperatures in the lower 500 m of the section. The extent of subsequent lower-temperature alteration is generally low and regionally variable. The gabbroic samples are slightly elevated in ⁸⁷Sr/⁸⁶Sr relative to fresh rock values (0.7024) and less enriched than the overlying sheeted dike complex. ⁸⁷Sr/⁸⁶Sr for the pervasively altered gabbroic samples ranges from 0.70244 to 0.70273 (mean 0.70257), tonalites is 0.7038, and pyroxene hornfels ranges from 0.70259 to 0.70271. ⁸⁷Sr/⁸⁶Sr does not vary with depth, and there is a strong positive correlation with the percentage of hydrous phases. Strontium contents of igneous and hydrothermal minerals, combined with bulk rock ⁸⁷Sr/⁸⁶Sr, indicate that Sr-isotopic exchange is largely controlled by the uptake of fluid ⁸⁷Sr/⁸⁶Sr in hydrous minerals and does not require Sr gain or loss. The minimum, time-integrated fluid–rock ratio for the sheeted dike complex and upper plutonics is 0.55–0.66, and the fluid flux calculated by mass balance is ~2.1 to 2.5 × 10⁶ kg m⁻², 30–60% higher than fluid fluxes calculated in the same manner for sheeted dike complexes on their own at Hess and Pito Deeps, and Ocean Drilling Program Hole 504B. Alteration patterns within the upper plutonics evolved in response to axial magma chamber (AMC) dynamics at the EPR, such that magma replenishment led to assimilation and thermal metamorphism of the country rock, and the position of the hydrothermal root-zone tracked the vertical migration of the AMC. The freshness of the lowermost gabbroic rocks suggests that pervasive fluid flow does not lead to significant fluid and heat fluxes at and near fast-spreading ridges.  相似文献   

Numerical modelling of bit–rock fracture mechanisms in percussive drilling with a continuum approach     

Timo Saksala 《国际地质力学数值与分析法杂志》2011,35(13):1483-1505

This paper presents a numerical method for continuum modelling of the dynamic bit–rock interaction process in percussive drilling. The method includes a constitutive model based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. The interaction between the drill bit and rock is modelled using contact mechanics by treating the bit as a rigid body. As the bit–rock interaction in percussive drilling is a transient event, the method is implemented in explicit dynamics FEM. The rock strength heterogeneity is characterized at the mesoscopic level statistically using the Weibull distribution. The bit–rock interaction is simulated under axisymmetric conditions using cylindrical and hemispherical buttons. The choice of the quite complex constitutive model accounting, e.g. for plastic compaction, viscoplastic shear and tensile failure along with induced damage and rate dependency is justified by numerical simulations. Moreover, the quasi‐static and dynamic cases are compared in plane strain simulations. Finally, some results clarifying the discrepancy of opinions found in the literature concerning the side (lateral) crack formation are obtained. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

How to attribute market leakage to CDM projects     

Frank Vöhringer  Timo Kuosmanen  Rob Dellink 《Climate Policy》2013,13(5):503-516

Abstract

Economic studies suggest that market leakage rates of greenhouse gas abatement can reach the two-digit percentage range. Although the Marrakesh Accords require Clean Development Mechanism (CDM) projects to account for leakage, most projects neglect market leakage. Insufficient leakage accounting is facilitated by a lack of applicable methods regarding the quantification and attribution of project-related leakage effects. This article proposes a method for attributing CDM-related market leakage effects to individual projects. To this purpose, alternative attribution methods are analysed. We find that project-specific approaches fail to take account of market leakage effects. Consequently, we propose to estimate aggregate market leakage effects and attribute them proportionally to individual projects. We suggest that predetermined commodity-specific leakage factors are applied by project developers to any emission reductions that are associated with a project's leakage-relevant demand or supply changes. This approach is conservative, equitable, incentive-compatible and applicable at manageable costs.  相似文献