Comments on “Evidence for global runoff increase related to climate warming” by Labat et al.

We have examined the evidence presented by Labat et al. and found that (1) their claims for a 4% increase in global runoff arising from a 1 °C increase in air temperature and (2) that their article provides the “first experimental data-based evidence demonstrating the link between the global warming and the intensification of the global hydrological cycle” are not supported by the data presented. Our conclusions are based on the facts that (1) their discharge records exhibit non-climatic influences and trends, (2) their work cannot refute previous studies finding no relation between air temperature and runoff, (3) their conclusions cannot explain relations before 1925, and (4) the statistical significance of their results hinges on a single data point that exerts undue influence on the slope of the regression line. We argue that Labat et al. have not provided sufficient evidence to support their claim for having detected increases in global runoff resulting from climate warming.     

Image processing of 2D resistivity data for imaging faults

A methodology to locate automatically limits or boundaries between different geological bodies in 2D electrical tomography is proposed, using a crest line extraction process in gradient images. This method is applied on several synthetic models and on field data set acquired on three experimental sites during the European project PALEOSIS where trenches were dug. The results presented in this work are valid for electrical tomographies data collected with a Wenner-alpha array and computed with an l₁ norm (blocky inversion) as optimization method. For the synthetic cases, three geometric contexts are modelled: a vertical and a dipping fault juxtaposing two different geological formations and a step-like structure. A superficial layer can cover each geological structure. In these three situations, the method locates the synthetic faults and layer boundaries, and determines fault displacement but with several limitations. The estimated fault positions correlate exactly with the synthetic ones if a conductive (or no superficial) layer overlies the studied structure. When a resistive layer with a thickness of 6 m covers the model, faults are positioned with a maximum error of 1 m. Moreover, when a resistive and/or a thick top layer is present, the resolution significantly decreases for the fault displacement estimation (error up to 150%). The tests with the synthetic models for surveys using the Wenner-alpha array indicate that the proposed methodology is best suited to vertical and horizontal contacts. Application of the methodology to real data sets shows that a lateral resistivity contrast of 1:5–1:10 leads to exact faults location. A fault contact with a resistivity contrast of 1:0.75 and overlaid by a resistive layer with a thickness of 1 m gives an error location ranging from 1 to 3 m. Moreover, no result is obtained for a contact with very low contrasts (1:0.85) overlaid by a resistive soil. The method shows poor results when vertical gradients are greater than horizontal ones. This kind of image processing technique should be systematically used for improving the objectiveness of tomography interpretation when looking for limits between geological objects.     

Short‐term intertidal bar mobility on a ridge‐and‐runnel beach,Merlimont, northern France

Digital elevation models and topographic pro?les of a beach with intertidal bar and trough (ridge‐and‐runnel) morphology in Merlimont, northern France, were analysed in order to assess patterns of cross‐shore and longshore intertidal bar mobility. The beach exhibited a pronounced dual bar–trough system that showed cross‐shore stationarity. The bars and troughs were, however, characterized by signi?cant longshore advection of sand under the in?uence of suspension by waves and transport by strong tide‐ and wind‐driven longshore currents. Pro?le changes were due in part to the longshore migration of medium‐sized bedforms. The potential for cross‐shore bar migration appears to be mitigated by the large size of the two bars relative to incident wave energy, which is modulated by high vertical tidal excursion rates on this beach due to the large tidal range (mean spring tidal range = 8·3 m). Cross‐shore bar migration is also probably hindered by the well‐entrenched troughs which are maintained by channelled high‐energy intertidal ?ows generated by swash bores and by tidal discharge and drainage. The longshore migration of intertidal bars affecting Merlimont beach is embedded in a regional coastal sand transport pathway involving tidal and wind‐forced northward residual ?ows affecting the rectilinear northern French coast in the eastern English Channel. Copyright © 2004 John Wiley & Sons, Ltd.     

Dike intrusion associated with the 2000 eruption of Miyakejima Volcano, Japan

The model for the 2000 dike intrusion event between Kozushima and Miyakejima volcano, Japan, was reinvestigated. After the sudden earthquake swarm in Miyakejima volcano, a dike intrusion of large volume was detected by the nationwide GPS network (Geonet). The displacements detected with GPS stations over an area with a radius of about 200 km shows a distribution that is consistent with the dike source being located near Miyakejima volcano.The dike was intruded northwestwards between Miyakejima and the neighboring Kozushima volcano. We searched for the parameters in the models that reproduce the regional displacements due to dike intrusion between Miyakejima and Kozushiima islands. We tested three models, (1) the model with a single dike, (2) the model with a dike and a point dislocation source which represents a creep dislocation source and (3) the model with a dike and a deflation source which represents a magma reservoir. Though all three models can match the horizontal displacements near the source area, model 1 fails to reproduce the regional displacements in the central part of Japan. Both models 2 and 3 can reproduce the regional displacement for horizontal components. Model 3 produces slightly better results than model 2 for vertical components. The balance in the volume budget for models 2 and 3 is also consistent with the observations. These results show that we cannot distinguish between the two models using only GPS observation. As there is no direct evidence for such a large creep or ductile source (corresponds to M7 or more) as proposed in model 2 and the active seismic region migrated back and forth within the linear swarm region, the model with a dike and a deep magma source is preferable. For the deflation point source, we obtained a deflation volume of 1.5 km³ at the depth of 20 km below the dike. An additional ~0.95 km³ of volume loss through caldera collapse and edifice deflation took place at Miyakejima. We conclude that the magma that intruded the dike came in part from below Miyakejima and in part from below the sea floor between Miyakejima and Kozushima, perhaps from reservoirs at the Moho.Editorial responsibility: S Nakada, T Druitt     

Evidence from the Northeastern Atlantic basin for variability in the rate of the meridional overturning circulation through the last deglaciation

A first study from the subtropical western Atlantic, using ²³¹Pa/²³⁰Th ratios as a kinematic proxy for deep water circulation, provided compelling evidence for a strong link between climate and the rate of Meridional Overturning Circulation (MOC) over the last deglaciation. However, these results warrant confirmation from additional locations and water depths because the interpretation of the sedimentary ²³¹Pa/²³⁰Th ratio in terms of circulation vigor can be biased by variations in particle flux and composition. We have measured ²³¹Pa/²³⁰Th in a core from the Iberian margin, in the Northeastern Atlantic basin, and have compared these new results to the data from the western Atlantic basin. We find that the reduction in the circulation during H1 and YD and the subsequent increases first recognized in the sediment deposited on Bermuda Rise are also evident in the eastern basin, in a totally different sedimentary regime, confirming that sedimentary ²³¹Pa/²³⁰Th ratios record basin-wide changes in deep water circulation. However, some differences between the eastern and western records are also recognized, providing preliminary evidence to differentiate between renewal rates in the two North Atlantic basins and between shallower and deeper overturning. Our results suggest the possible existence of two sources of Glacial North Atlantic Intermediate Deep Water (GNAIW), one in the south Labrador Sea and another west of Rockall Plateau. Both sources contributed to the meridional overturning but the two had different sensitivity to meltwater from the Laurentide and the Fennoscandian ice sheets during the deglaciation. These results indicate that additional information on the geometry and strength of the ventilation of the deep Atlantic can be obtained by contrasting the evolution of sediment ²³¹Pa/²³⁰Th in different sections of the Atlantic Ocean.     

Assessment of regional air pollution distribution by point cumulative semivariogram method at Erzurum urban center, TURKEY

In the assessment of air quality, regional distribution and dispersion with distance are important, together with the variations of pollutants in time. On this occasion, the point cumulative semi-variogram (PCSV) method is used in order to find simply regional distribution of pollutants of Erzurum urban centre. This method is based simply on the summation of square differences in air pollutant concentrations between different sites. Monthly regional variation maps of Erzurum are constructed by finding radius of influence (for SO₂, from 1000 m to 3500 m and, for TSP, 1000–2000 m) and PCSV scattering diagram data at different levels by using monthly average sulphur dioxide (SO₂) and total suspended particulate (TSP) matter concentrations in 2001–2002 winter season. Consequently, the air pollution distribution of Erzurum is assessed.     

Pseudo‐dynamic seismic response of reinforced concrete frames infilled with non‐structural brick masonry

This paper presents pseudo‐dynamic test results on the in‐plane seismic behaviour of infilled frames. Thirteen single‐storey, single‐bay, half‐size‐scale, reinforced concrete‐frame specimens, most of which infilled with non‐structural masonry made of perforated bricks and cement mortar are tested. The infills are in contact with frames, without any connector; openings are not covered. The frames are different in their strength and details, reinforcement grade, and aspect ratio. Seismic input is the 1976 Tolmezzo (Friuli, Italy) ground acceleration, to which specimens are subjected two times: virgin and damaged by the previous test. The global seismic response of initially virgin infilled specimens considerably differs from that of bare specimens. This follows a dramatic change of properties: compared to a bare frame, the initial stiffness increases by one order of magnitude, and the peak strength more than doubles. The peak drift lessens; however, the displacement ductility demand does not. The energy demand is greater. Nevertheless, the influence of infill decreases as damage proceeds. Displacement time histories of damaged specimens are quite similar. At the local level, infill causes asymmetry and concentration of the frame deformation. Copyright © 2005 John Wiley & Sons, Ltd.     

Evolutionary aseismic design and retrofit of structures with passive energy dissipation

A new computational framework is developed for the design and retrofit of building structures by considering aseismic design as a complex adaptive process. For the initial phase of the development within this framework, genetic algorithms are employed for the discrete optimization of passively damped structural systems. The passive elements may include metallic plate dampers, viscous fluid dampers and viscoelastic solid dampers. The primary objective is to determine robust designs, including both the non‐linearity of the structural system and the uncertainty of the seismic environment. Within the present paper, this computational design approach is applied to a series of model problems, involving sizing and placement of passive dampers for energy dissipation. In order to facilitate our investigations and provide a baseline for further study, we introduce several simplifications for these initial examples. In particular, we employ deterministic lumped parameter structural models, memoryless fitness function definitions and hypothetical seismic environments. Despite these restrictions, some interesting results are obtained from the simulations and we are able to gain an understanding of the potential for the proposed evolutionary aseismic design methodology. Copyright © 2005 John Wiley & Sons, Ltd.     

Sediment Distribution Around Glacially Abraded Bedrock Landforms (Whalebacks) at Lago Tranquilo, Chile

Whalebacks are convex landforms created by the smoothing of bedrock by glacial processes. Their formation is attributed to glacial abrasion either by bodies of subglacial sediment sliding over bedrock or by individual clasts contained within ice. This paper reports field measurements of sediment depth around two whaleback landforms in order to investigate the relationship between glacigenic deposits and whaleback formation. The study site, at Lago Tranquilo in Chilean Patagonia, is situated within the Last Glacial Maximum (LGM) ice limits. The two whalebacks are separated by intervening depressions in which sediment depths are generally 0.2 to 0.3 m. Two facies occur on and around the whalebacks. These facies are: (1) angular gravel found only on the surface of the whalebacks, interpreted as bedrock fracturing in response to unloading of the rock following pressure release after ice recession, and (2) sandy boulder‐gravel in the sediment‐filled depressions between the two whalebacks, interpreted as an ice‐marginal deposit, with a mixture of sediment types including basal glacial and glaciofluvial sediment. Since the whalebacks have heavily abraded and striated surfaces but are surrounded by only a patchy and discontinuous layer of sediment, the implication is that surface abrasion of the whalebacks was achieved primarily by clasts entrained in basal ice, not by subglacial till sliding.