Radiocarbon wiggle matching of Swedish lake varves reveals asynchronous climate changes around the 8.2 kyr cold event

Snowball, I., Muscheler, R., Zillén, L., Sandgren, P., Stanton, T. & Ljung, K. 2010: Radiocarbon wiggle matching of Swedish lake varves reveals asynchronous climate changes around the 8.2 kyr cold event. Boreas, Vol. 39, pp. 720–733. 10.1111/j.1502‐3885.2010.00167.x. ISSN 0300‐9483. A series of nine radiocarbon measurements obtained at increments of 50 years was used to wiggle match varved sediments in Lake Kälksjön (west central Sweden) to the calendar year time scale provided by the tree‐ring‐derived radiocarbon calibration curve. The anchor point for a series of 400 varves known to span the ‘8.2 kyr cold event’ was derived from a combination of correlation analysis, Monte Carlo statistics and Bayesian age modelling. The GRIP and NorthGRIP ice‐cores were matched to the same absolute time scale by comparing ¹⁰Be data and tree‐ring Δ¹⁴C. The radiocarbon‐based wiggle match, organic carbon measurements, mineral magnetic parameters and XRF data reconstruct a distinct period of enhanced erosion in Kälksjön's catchment between 8066 ± 25 and 7920 ± 25 cal. yr BP. Prior to human impact, the increased erosion in the boreal environment is assigned to increased winter precipitation in the form of snow, which caused spring meltwater discharge to intensify. Our results suggest that the abrupt onset of increased winter precipitation in west central Sweden started at least 50 years after the onset of the ‘8.2 kyr cold event’ as defined by oxygen isotope data from Greenland. The study highlights the value of synchronized time scales in the reconstruction of abrupt climate changes based on proxy data.     

Anticipating the effects of groundwater withdrawal on seawater intrusion and soil settlement in urban coastal areas

Intensive pumping in urban coastal areas is a common threat to water resource quality due to seawater intrusion. In those areas where subsurface water resources are not usually used for human consumption or irrigation, intensive pumping is associated with other activities like the lowering of the water table necessary to support underground structures and building foundations. This activity also increases the likelihood of soil settlement that affects building stability and the corrosion of concrete structures due to groundwater salinity. Under these circumstances, the awareness of a certain municipality (Calonge, NE Spain) of the potential effects of groundwater withdrawal upon foundations has led to an integrated approach to anticipate seawater intrusion related to urban development. Geological mapping and correlation of borehole logs, electrical resistivity tomography, and hydrochemical data provide comprehensive knowledge of the geology and hydrogeology of the area and act as screening tools necessary to discern the influence of hydrological processes in coastal areas. Developing Strack's analytical solution, new comprehensive, dimensionless expressions are herein derived to determine the critical pumping rate necessary to prevent seawater intrusion, as well as to reproduce the evolution of the wedge toe and the water table stagnation point under different withdrawal rates. Furthermore, the Dupuit–Forchheimer well discharge formula allows the estimation of the effects of the water table lowering due to such critical pumping in the surrounding building foundations. Field data from the Calonge coastal plain illustrate this approach and provide assessment criteria for future urban development and planning. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.     

Non linear waves in two dimensional keplerian flows

Abstract

Accretion discs in astrophysics are fundamental for converting gravitational binding energy into observed electromagnetic radiation. We study the behavior of waves in a two dimensional supersonic Keplerian flow inside a given gravitational potential. We present the effects of shearing and rotation on short waves, and the numerical study of the dynamical stability of such flows with respect to various perturbations. We show that a large class of dynamical effects, due to pressure and associated to short time scales, may be excited.     

TITAN2D simulations of pyroclastic flows at Cerro Machín Volcano,Colombia: Hazard implications

Cerro Machín is a dacitic tuff ring located in the central part of the Colombian Andes. It lies at the southern end of the Cerro Bravo–Cerro Machín volcanic belt. This volcano has experienced at least six major explosive eruptions during the last 5000 years. These eruptions have generated pyroclastic flows associated with Plinian activity that have traveled up to 8 km from the crater, and pyroclastic flows associated with Vulcanian activity with shorter runouts of 5 km from the source. Today, some 21,000 people live within a 8 km radius of Cerro Machín. The volcano is active with fumaroles and has shown increasing seismic activity since 2004, and therefore represents a potentially increasing threat to the local population. To evaluate the possible effects of future eruptions that may generate pyroclastic density currents controlled by granular flow dynamics we performed flow simulations with the TITAN2D code. These simulations were run in all directions around the volcano, using the input parameters of the largest eruption reported. The results show that an eruption of 0.3 km³ of pyroclastic flows from a collapsing Plinian column would travel up to 9 km from the vent, emplacing a deposit thicker than 60 m within the Toche River valley. Deposits >45 m thick can be expected in the valleys of San Juan, Santa Marta, and Azufral creeks, while 30 m thick deposits could accumulate within the drainages of the Tochecito, Bermellón, and Coello Rivers. A minimum area of 56 km² could be affected directly by this kind of eruption. In comparison, Vulcanian column-collapse pyroclastic flows of 0.1 km³ would travel up to 6 km from the vent depositing >45 m thick debris inside the Toche River valley and more than 30 m inside the valleys of San Juan, Santa Marta, and Azufral creeks. The minimum area that could be affected directly by this kind of eruption is 33 km². The distribution and thickness of the deposits obtained by these simulations are consistent with the hazard map presented by INGEOMINAS (Geological Survey of Colombia) in 2002. The composite map of the simulated flow deposits suggests that after major explosive events such as these, the generation of lahars is probable.     

THE HYDROGEOLOGY OF GHANA

SUMMARY

Many hydrological systems exhibit long periods of quiescent behaviour, broken only occasionally by short periods of comparatively intense activity. Clear examples arise in stream run-off measurements and in ground water table observations. The number of measurements recorded can then be considerably reduced by reducing the frequency of recording over the quiescent periods. The paper describes an algorithm whereby this reduction is attained in an optimum manner. When recordings can be accepted directly by a computer, the algorithm simply provides a translation programme Whereby the quality of data may be very considerably reduced. Other applications are, however, discussed.

The problems of Optimum representations is first presented, and choices of norm indicated. Using a particular norm, it is shown that the optimum representation is characterised, under certain limiting conditions, by the property that the second moments of the second derivatives are constant over all optimum intervals. This characterisation is used to construct an iterative difference scheme, whereby a certain truncation error is accepted. Sufficient conditions for stability of the scheme are determined and these incorporated into a computer programme. The programme has been tested on a variety of data with satisfactory results. In conclusion, a further development of the method is indicated.     

Formation of intracratonic basins by lithospheric shortening and phase changes: a case study from the ultra‐deep East Barents Sea basin

The origin of large subsidence in intracratonic basins is still under debate. We propose a new and self‐consistent model for the formation of those basins, where lithospheric shortening/buckling triggers metamorphism and densification of crustal mafic heterogeneities. We use a forward thermo‐mechanical finite element technique to evaluate this mechanism for the typical example of the East Barents Sea basin (EBB) where a very large and compensated subsidence, accommodating an up to 20‐km‐thick sediment succession, is observed. The lower crust in the dynamic model is modelled with petrologic‐consistent densities for a wet mafic gabbroic composition that depend on pressure and temperature taking into account dehydration at high PT conditions. The model successfully explains the main characteristics of the EBB, notably the large anomalous and fast subsidence during the Late Permian–Early Triassic, its present‐day geometry and the absence of a significant gravity anomaly.     

Oxygen and organic carbon fluxes in sediments of the Bay of Biscay

The relationship between particulate organic carbon (POC) concentrations measured in modern sediment and fluxes of exported POC to the sediment surface needs to be understood in order to use POC content as a proxy of paleo-environmental conditions. The objective of our study was to compare POC concentrations, POC mineralization rates calculated from O₂ consumption and POC burial rates. Benthic O₂ distributions were determined in 58 fine-grained sediment cores collected at different periods at 14 stations in the southeastern part of the Bay of Biscay with depths ranging from 140 to 2800 m. Depth-dependent volume-specific oxygen consumption rates were used to assess rates of aerobic oxidation of organic matter (OM), assuming that O₂ consumption solely was related to heterotrophic activity at the sediment–water interface. Heterogeneity of benthic O₂ fluxes denoted changes in time and space of fresh organic material sedimentation. The most labile fraction of exported POC engendered a steep decrease in concentration in the upper 5 mm of vertical O₂ profiles. The rupture in the gradient of O₂ microprofile may be related to the bioturbation-induced mixing depth of fast-decaying carbon. Average diffusive O₂ fluxes showed that this fast-decaying OM flux was much higher than buried POC, although diffusive O₂ fluxes underestimated the total sediment oxygen demand, and thus the fast-decaying OM flux to the sediment surface. Sedimentary POC burial was calculated from sediment mass accumulation rate and the organic carbon content measured at the top of the sediment. The proportion of buried POC relative to total exported POC ranged at the most between 50% and 10%, depending on station location. Therefore, for a narrow geographic area like the Bay of Biscay, burial efficiency of POC was variable. A fraction of buried POC consisted of slow-decaying OM that was mineralized within the upper decimetres of sediment through oxic and anoxic processes. This fraction was deduced from the decrease with depth in POC concentration. At sites located below 500 m water depth, where the fast-decaying carbon did not reach the anoxic sediment, the slow-decaying pool may control the O₂ penetration depth. Only refractory organic material was fossilized in sedimentary records at locations where labile OM did not reach the anoxic portion of the sediment.     

Integration of remote sensing data into hydrological models for reservoir management

Abstract

The purpose of this paper is to present the methodology set up to derive catchment soil moisture from Earth Observation (EO) data using microwave spaceborne Synthetic Aperture Radar (SAR) images from ERS satellites and to study the improvements brought about by an assimilation of this information into hydrological models. The methodology used to derive EO data is based on the appropriate selection of land cover types for which the radar signal is mainly sensitive to soil moisture variations. Then a hydrological model is chosen, which can take advantage of the new information brought by remote sensing. The assimilation of soil moisture deduced from EO data into hydrological models is based principally on model parameter updating. The main assumption of this method is that the better the model simulates the current hydrological system, the better the following forecast will be. Another methodology used is a sequential one based on Kalman filtering. These methods have been put forward for use in the European AIMWATER project on the Seine catchment upstream of Paris (France) where dams are operated to alleviate floods in the Paris area.