ALBERT

Description: Review Article: On the relation between the seismic activity and the Hurst exponent of the geomagnetic field at the time of the 2000 Izu swarm Natural Hazards and Earth System Science, 13, 2189-2194, 2013 Author(s): F. Masci and J. N. Thomas Many papers document the observation of earthquake-related precursory signatures in geomagnetic field data. However, the significance of these findings is ambiguous because the authors did not adequately take into account that these signals could have been generated by other sources, and the seismogenic origin of these signals have not been validated by comparison with independent datasets. Thus, they are not reliable examples of magnetic disturbances induced by the seismic activity. Hayakawa et al. (2004) claim that at the time of the 2000 Izu swarm the Hurst exponent of the Ultra-Low-Frequency (ULF: 0.001–10 Hz) band of the geomagnetic field varied in accord with the energy released by the seismicity. The present paper demonstrates that the behaviour of the Hurst exponent was insufficiently investigated and also misinterpreted by the authors. We clearly show that during the Izu swarm the changes of the Hurst exponent were strongly related to the level of global geomagnetic activity and not to the increase of the local seismic activity.

Description: Forest fire danger rating in complex topography – results from a case study in the Bavarian Alps in autumn 2011 Natural Hazards and Earth System Science, 13, 2157-2167, 2013 Author(s): C. Schunk, C. Wastl, M. Leuchner, C. Schuster, and A. Menzel Forest fire danger rating based on sparse meteorological stations is known to be potentially misleading when assigned to larger areas of complex topography. This case study examines several fire danger indices based on data from two meteorological stations at different elevations during a major drought period. This drought was caused by a persistent high pressure system, inducing a pronounced temperature inversion and its associated thermal belt with much warmer, dryer conditions in intermediate elevations. Thus, a massive drying of fuels, leading to higher fire danger levels, and multiple fire occurrences at mid-slope positions were contrasted by moderate fire danger especially in the valleys. The ability of fire danger indices to resolve this situation was studied based on a comparison with the actual fire danger as determined from expert observations, fire occurrences and fuel moisture measurements. The results revealed that, during temperature inversion, differences in daily cycles of meteorological parameters influence fire danger and that these are not resolved by standard meteorological stations and fire danger indices (calculated on a once-a-day basis). Additional stations in higher locations or high-resolution meteorological models combined with fire danger indices accepting at least hourly input data may allow reasonable fire danger calculations under these circumstances.

Description: Preface "Landslide hazard and risk assessment at different scales" Natural Hazards and Earth System Science, 13, 2169-2171, 2013 Author(s): P. Reichenbach, A. Günther, and T. Glade

Description: Linking ICT and society in early warning and adaptation to hydrological extremes in mountains Natural Hazards and Earth System Science, 13, 2253-2270, 2013 Author(s): C. de Jong The assessment of the societal impact of hydrological extremes is particularly important in mountain regions, since mountains can be considered both as the generators and victims of extreme events. ICT can provide a powerful tool for transmitting hydro-meteorological information to predict, prepare and adapt to such events. However, in remote regions, such as mountains, the poles, deserts and islands, preventive and adaptive measures are often restricted by data availability and lack and/or incoherence of data networks. This paper distinguishes between early warning of floods and droughts, emphasising the latter in particular in mountains and explores the possibilities of enhancing the role of society in data collection, the identification, activation and application of stakeholder knowledge and transferral of data from gauged to ungauged catchments.

Description: Similarity estimators for irregular and age uncertain time series Climate of the Past Discussions, 9, 5299-5346, 2013 Author(s): K. Rehfeld and J. Kurths Paleoclimate time series are often irregularly sampled and age uncertain, which is an important technical challenge to overcome for successful reconstruction of past climate variability and dynamics. Visual comparison and interpolation-based linear correlation approaches have been used to infer dependencies from such proxy time series. While the first is subjective, not measurable and not suitable for the comparison of many datasets at a time, the latter introduces interpolation bias, and both face difficulties if the underlying dependencies are nonlinear. In this paper we investigate similarity estimators that could be suitable for the quantitative investigation of dependencies in irregular and age uncertain time series. We compare the Gaussian-kernel based cross correlation (gXCF, Rehfeld et al., 2011) and mutual information (gMI, Rehfeld et al., 2013) against their interpolation-based counterparts and the new event synchronization function (ESF). We test the efficiency of the methods in estimating coupling strength and coupling lag numerically, using ensembles of synthetic stalagmites with short, autocorrelated, linear and nonlinearly coupled proxy time series, and in the application to real stalagmite time series. In the linear test case coupling strength increases are identified consistently for all estimators, while in the nonlinear test case the correlation-based approaches fail. The lag at which the time series are coupled is identified correctly as the maximum of the similarity functions in around 60–55% (in the linear case) to 53–42% (for the nonlinear processes) of the cases when the dating of the synthetic stalagmite is perfectly precise. If the age uncertainty increases beyond 5% of the time series length, however, the true coupling lag is not identified more often than the others for which the similarity function was estimated. Age uncertainty contributes up to half of the uncertainty in the similarity estimation process. Time series irregularity contributes less, particularly for the adapted Gaussian-kernel based estimators and the event synchronization function. The introduced link strength concept summarizes the hypothesis test results and balances the individual strengths of the estimators: while gXCF is particularly suitable for short and irregular time series, gMI and the ESF can identify nonlinear dependencies. ESF could, in particular, be suitable to study extreme event dynamics in paleoclimate records. Programs to analyze paleoclimatic time series for significant dependencies are included in a freely available software toolbox.

Description: Impact of precession on the climate, vegetation and fire activity in southern Africa during MIS4 Climate of the Past Discussions, 9, 5391-5438, 2013 Author(s): M.-N. Woillez, G. Levavasseur, A.-L. Daniau, M. Kageyama, D. H. Urrego, and M.-F. Sánchez-Goñi The relationships between climate, vegetation and fires are a major subject of investigation in the context of climate change. In southern Africa, fire is known to play a crucial role in the existence of grasslands and Mediterranean-like biomes. Microcharcoal-based reconstructions of past fire activity in that region have shown a tight correlation between grass-fueled fires and the precessional cycle, with maximum fire activity during maxima of the climatic precession index. These changes have been interpreted as the result of changes in fuel load in response to precipitation changes in eastern southern Africa. Here we use the general circulation model IPSL_CM5A and the dynamical vegetation model LPJ-LMfire to investigate the response of climate, vegetation and fire activity to precession changes in southern Africa during Marine Isotopic Stage 4. We perform two climatic simulations, for a maximum and minimum of the precession index, and use a statistical downscaling method to increase the spatial resolution of the IPSL_CM5A outputs over southern Africa and perform high-resolution simulations of the vegetation and fire activity. Our results show an anti-correlation between the North and South African monsoons in response to precession changes. A decrease of the precession climatic index leads to a precipitation decrease in the summer rainfall area of southern Africa. The drying of climate leads to a decrease of vegetation cover and fire activity. Our results are in qualitative agreement with data and confirm that fire activity in southern Africa is strongly dependent on the vegetation type.

Description: The environmental impact of the Puyehue–Cordon Caulle 2011 volcanic eruption on Buenos Aires Natural Hazards and Earth System Science, 13, 2319-2330, 2013 Author(s): G. B. Raga, D. Baumgardner, A. G. Ulke, M. Torres Brizuela, and B. Kucienska On 4 June 2011, the volcanic complex Puyehue–Cordon Caulle located in the Chilean Andes erupted, producing a plume of gases and particles that eventually circled the Southern Hemisphere, disrupting air travel and depositing ash in large quantities. On eight occasions, the plume passed over the city of Buenos Aires, Argentina, leading local authorities to close the two international airports. The eruption occurred during an on-going field campaign when measurements of the properties of atmospheric aerosol particles were being made in Buenos Aires as part of a year-long study of the concentration and optical properties of aerosol at one site in the city. The suite of instruments deployed in Buenos Aires were not tailored to measurements of volcanic ash, but were designed to characterize urban conditions. Nevertheless, these measurements were analysed for periods when vertical profiles of aerosol backscatter, made with a ceilometer, clearly showed the presence of the volcano plume over the research site and resulted in airport closure. Aerosol optical thickness derived from AERONET, MODIS and a ceilometer at our research site, all show enhanced values clearly indicating that the three platforms identified the volcanic plume simultaneously. However, a quantitative comparison of the different estimates proves difficult, suggesting large spatial and temporal variability of the plume. Our results indicate that the number concentration of condensation nuclei (CN), the mass concentration of particle-bound polycyclic aromatic hydrocarbons (PPAH) and the light absorption coefficient exceeded the average background values by more than one standard deviation during the events of volcanic plume. The anomalous concentrations of CN suggest new particle formation, presumably from the conversion of SO 2 , while the anomalous concentrations of PPAH may come from the uptake of PAHs on the plume particles or from chemical reactions on the surface of plume particles. The anomalous absorption coefficients indicate that plume particles may contain certain compounds that can absorb radiation at 550 nm. Another possible explanation consistent with the observations is the scavenging of black carbon from urban sources as the plume descends through the boundary layer to the surface. In addition, the volcanic plume influenced the local meteorology resulting in a decrease of the temperature when compared to the average temperature during days with no plume present.

Description: Contribution of land use changes to future flood damage along the river Meuse in the Walloon region Natural Hazards and Earth System Science, 13, 2301-2318, 2013 Author(s): A. Beckers, B. Dewals, S. Erpicum, S. Dujardin, S. Detrembleur, J. Teller, M. Pirotton, and P. Archambeau Managing flood risk in Europe is a critical issue because climate change is expected to increase flood hazard in many european countries. Beside climate change, land use evolution is also a key factor influencing future flood risk. The core contribution of this paper is a new methodology to model residential land use evolution. Based on two climate scenarios ("dry" and "wet"), the method is applied to study the evolution of flood damage by 2100 along the river Meuse. Nine urbanization scenarios were developed: three of them assume a "current trend" land use evolution, leading to a significant urban sprawl, while six others assume a dense urban development, characterized by a higher density and a higher diversity of urban functions in the urbanized areas. Using damage curves, the damage estimation was performed by combining inundation maps for the present and future 100 yr flood with present and future land use maps and specific prices. According to the dry scenario, the flood discharge is expected not to increase. In this case, land use changes increase flood damages by 1–40%, to €334–462 million in 2100. In the wet scenario, the relative increase in flood damage is 540–630%, corresponding to total damages of €2.1–2.4 billion. In this extreme scenario, the influence of climate on the overall damage is 3–8 times higher than the effect of land use change. However, for seven municipalities along the river Meuse, these two factors have a comparable influence. Consequently, in the "wet" scenario and at the level of the whole Meuse valley in the Walloon region, careful spatial planning would reduce the increase in flood damage by no more than 11–23%; but, at the level of several municipalities, more sustainable spatial planning would reduce future flood damage to a much greater degree.

Description: Estimating soil suction from electrical resistivity Natural Hazards and Earth System Science, 13, 2369-2379, 2013 Author(s): E. Piegari and R. Di Maio Soil suction and resistivity strongly depend on the degree of soil saturation and, therefore, both are used for estimating water content variations. The main difference between them is that soil suction is measured using tensiometers, which give point information, while resistivity is obtained by tomography surveys, which provide distributions of resistivity values in large volumes, although with less accuracy. In this paper, we have related soil suction to electrical resistivity with the aim of obtaining information about soil suction changes in large volumes, and not only for small areas around soil suction probes. We derived analytical relationships between soil matric suction and electrical resistivity by combining the empirical laws of van Genuchten and Archie. The obtained relationships were used to evaluate maps of soil suction values in different ashy layers originating in the explosive activity of the Mt Somma-Vesuvius volcano (southern Italy). Our findings provided a further example of the high potential of geophysical methods in contributing to more effective monitoring of soil stress conditions; this is of primary importance in areas where rainfall-induced landslides occur periodically.

Description: Sulfur dioxide emissions from Papandayan and Bromo, two Indonesian volcanoes Natural Hazards and Earth System Science, 13, 2399-2407, 2013 Author(s): P. Bani, Surono, M. Hendrasto, H. Gunawan, and S. Primulyana Indonesia hosts 79 active volcanoes, representing 14% of all active volcanoes worldwide. However, little is known about their SO 2 contribution into the atmosphere, due to isolation and access difficulties. Existing SO 2 emission budgets for the Indonesian archipelago are based on extrapolations and inferences as there is a considerable lack of field assessments of degassing. Here, we present the first SO 2 flux measurements using differential optical absorption spectroscopy (DOAS) for Papandayan and Bromo, two of the most active volcanoes in Indonesia. Results indicate mean SO 2 emission rates of 1.4 t d −1 from the fumarolic activity of Papandayan and more than 22–32 t d −1 of SO 2 released by Bromo during a declining eruptive phase. These DOAS results are very encouraging and pave the way for a better evaluation of Indonesian volcanic emissions.

Description: Ensemble meteorological reconstruction using circulation analogues of 1781–1785 Climate of the Past Discussions, 9, 5157-5182, 2013 Author(s): P. Yiou, M. Boichu, R. Vautard, M. Vrac, S. Jourdain, E. Garnier, F. Fluteau, and L. Menut This paper uses a method of atmospheric flow analogues to reconstruct an ensemble of atmospheric variables (namely sea-level pressure, surface temperature and wind speed) between 1781 and 1785. The properties of this ensemble are investigated and tested against observations of temperature. The goal of the paper is to assess whether the atmospheric circulation during the Laki volcanic eruption (in 1783) and the subsequent winter were similar to the conditions that prevailed in the winter 2009/2010 and during spring 2010. We find that the three months following the Laki eruption in June 1783 barely have analogues in 2010. The cold winter of 1783/1784 yields circulation analogues in 2009/2010. Our analysis suggests that it is unlikely that the Laki eruption was responsible for the cold winter of 1783/1784, of the relatively short memory of the atmospheric circulation.

Description: Seasonal changes in glacial polynya activity inferred from Weddell Sea varves Climate of the Past Discussions, 9, 5123-5156, 2013 Author(s): D. Sprenk, M. E. Weber, G. Kuhn, V. Wennrich, T. Hartmann, and K. Seelos The Weddell Sea and the associated Filchner-Rønne Ice Shelf constitute key regions for global bottom-water production today. However, little is known about bottom-water production under different climate and ice-sheet conditions. Therefore, we studied core PS1795, which consists primarily of fine-grained siliciclastic varves that were deposited on contourite ridges in the southeastern Weddell Sea during the Last Glacial Maximum (LGM). We conducted high-resolution X-ray fluorescence (XRF) analysis and grain-size measurements with the RADIUS tool (Seelos and Sirocko, 2005) using thin sections to characterize the two seasonal components of the varves at sub-mm resolution to distinguish the seasonal components of the varves. Bright layers contain coarser grains that can mainly be identified as quartz in the medium to coarse silt grain size. They also contain higher amounts of Si, Zr, Ca, and Sr, as well as more ice-rafted debris (IRD). Dark layers, on the other hand, contain finer particles such as mica and clay minerals from the chlorite and illite groups. In addition, chemical elements, Fe, Ti, Rb, and K are elevated as well. Based on these findings as well as on previous analyses on neighbouring cores, we propose a model of glacially enhanced thermohaline convection in front of a grounded ice sheet that is supported by seasonally variable coastal polynya activity. Accordingly, katabatic (i.e. offshore blowing) winds removed sea ice from the ice edge, leading to coastal polynya formation. We suggest that glacial processes were similar to today with stronger katabatic winds and enhanced coastal polynya activity during the winter season. If this is correct, silty layers are likely glacial winter deposits, when brine rejection was increased, leading to enhanced bottom water formation and increased sediment transport. Vice versa, finer-grained clayey layers were then deposited during summer, when coastal polynya activity was likely reduced.

Description: Interviewing insights regarding the fatalities inflicted by the 2011 Great East Japan Earthquake Natural Hazards and Earth System Science, 13, 2173-2187, 2013 Author(s): M. Ando, M. Ishida, Y. Hayashi, C. Mizuki, Y. Nishikawa, and Y. Tu One hundred fifty survivors of the 11 March 2011 Great East Japan Earthquake (Tohoku-oki earthquake) ( M w = 9.0) were interviewed to study the causes of deaths from the associated tsunami in coastal areas of Tohoku. The first official tsunami warning underestimated the height of the tsunami and 40% of the interviewees did not obtain this warning due to immediate blackouts and a lack of communication after the earthquake. Many chose to remain in dangerous locations based on the underestimated warning and their experiences with previous smaller tsunamis and/or due to misunderstanding the mitigating effects of nearby breakwaters in blocking incoming tsunamis. Some delayed their evacuation to perform family safety checks, and in many situations, the people affected misunderstood the risks involved in tsunamis. In this area, three large tsunamis have struck in the 115 yr preceding the 2011 tsunami. These tsunamis remained in the collective memory of communities, and numerous measures against future tsunami damage, such as breakwaters and tsunami evacuation drills, had been implemented. Despite these preparedness efforts, approximately 18 500 deaths and cases of missing persons occurred. The death rate with the age of 65 and above was particularly high, four times higher than that with other age groups. These interviews indicate that deaths resulted from a variety of reasons, but if residents had taken immediate action after the major ground motion stopped, most residents might have been saved. Education about the science behind earthquakes and tsunamis could help save more lives in the future.

Description: Is blue intensity ready to replace maximum latewood density as a strong temperature proxy? A tree-ring case study on Scots pine from northern Sweden Climate of the Past Discussions, 9, 5227-5261, 2013 Author(s): J. A. Björklund, B. E. Gunnarson, K. Seftigen, J. Esper, and H. W. Linderholm At high latitudes, where low temperatures mainly limit tree-growth, measurements of wood density (e.g. Maximum Latewood Density, MXD) using the X-Ray methodology provide a temperature proxy that is superior to that of TRW. Density measurements are however costly and time consuming and have lead to experimentation with optical flatbed scanners to produce Maximum Blue Intensity (BI max ). BI max is an excellent proxy for density on annual scale but very limited in skill on centennial scale. Discolouration between samples is limiting BI max where specific brightnesses can have different densities. To overcome this, the new un-exploited parameter Δ blue intensity (ΔBI) was constructed by using the brightness in the earlywood (BI EW ) as background, (BI max − BI EW = ΔBI). This parameter was tested on X-Ray material (MXD − earlywood density = ΔMXD) and showed great potential both as a quality control and as a booster of climate signals. Unfortunately since the relationship between grey scale and density is not linear, and between-sample brightness can differ tremendously for similar densities, ΔBI cannot fully match ΔMXD in skill as climate proxy on centennial scale. For ΔBI to stand alone, the range of brightness/density offset must be reduced. Further studies are needed to evaluate this possibility, and solutions might include heavier sample treatment (reflux with chemicals) or image-data treatment (digitally manipulating base-line levels of brightness).

Description: Review Article: Potential geomorphic consequences of a future great ( M w = 8.0+) Alpine Fault earthquake, South Island, New Zealand Natural Hazards and Earth System Science, 13, 2279-2299, 2013 Author(s): T. R. Robinson and T. R. H. Davies The Alpine Fault in New Zealand's South Island has not sustained a large magnitude earthquake since ca. AD 1717. The time since this rupture is close to the average inferred recurrence interval of the fault (~300 yr). The Alpine Fault is therefore expected to generate a large magnitude earthquake in the near future. Previous ruptures of this fault are inferred to have generated M w = 8.0 or greater earthquakes and to have resulted in, amongst other geomorphic hazards, large-scale landslides and landslide dams throughout the Southern Alps. There is currently 85% probability that the Alpine Fault will cause a M w = 8.0+ earthquake within the next 100 yr. While the seismic hazard is fairly well understood, that of the consequential geomorphic activity is less well studied, and these consequences are explored herein. They are expected to include landsliding, landslide damming, dam-break flooding, debris flows, river aggradation, liquefaction, and landslide-generated lake/fiord tsunami. Using evidence from previous events within New Zealand as well as analogous international examples, we develop first-order estimates of the likely magnitude and possible locations of the geomorphic effects associated with earthquakes. Landsliding is expected to affect an area 〉 30 000 km 2 and involve 〉 1billion m 3 of material. Some tens of landslide dams are expected to occur in narrow, steep-sided gorges in the affected region. Debris flows will be generated in the first long-duration rainfall after the earthquake and will continue to occur for several years as rainfall (re)mobilises landslide material. In total more than 1000 debris flows are likely to be generated at some time after the earthquake. Aggradation of up to 3 m will cover an area 〉 125 km 2 and is likely to occur on many West Coast alluvial fans and floodplains. The impact of these effects will be felt across the entire South Island and is likely to continue for several decades.

Description: Integrating spatial, temporal, and size probabilities for the annual landslide hazard maps in the Shihmen watershed, Taiwan Natural Hazards and Earth System Science, 13, 2353-2367, 2013 Author(s): C. Y. Wu and S. C. Chen Landslide spatial, temporal, and size probabilities were used to perform a landslide hazard assessment in this study. Eleven intrinsic geomorphological, and two extrinsic rainfall factors were evaluated as landslide susceptibility related factors as they related to the success rate curves, landslide ratio plots, frequency distributions of landslide and non-landslide groups, as well as probability–probability plots. Data on landslides caused by Typhoon Aere in the Shihmen watershed were selected to train the susceptibility model. The landslide area probability, based on the power law relationship between the landslide area and a noncumulative number, was analyzed using the Pearson type 5 probability density function. The exceedance probabilities of rainfall with various recurrence intervals, including 2, 5, 10, 20, 50, 100 and 200 yr, were used to determine the temporal probabilities of the events. The study was conducted in the Shihmen watershed, which has an area of 760 km 2 and is one of the main water sources for northern Taiwan. The validation result of Typhoon Krosa demonstrated that this landslide hazard model could be used to predict the landslide probabilities. The results suggested that integration of spatial, area, and exceedance probabilities to estimate the annual probability of each slope unit is feasible. The advantage of this annual landslide probability model lies in its ability to estimate the annual landslide risk, instead of a scenario-based risk.

Description: On the puzzling feature of the silence of precursory electromagnetic emissions Natural Hazards and Earth System Science, 13, 2381-2397, 2013 Author(s): K. Eftaxias, S. M. Potirakis, and T. Chelidze It has been suggested that fracture-induced MHz–kHz electromagnetic emissions (EME), which emerge from a few days up to a few hours before the main seismic shock occurrence permit a real-time monitoring of the damage process during the last stages of earthquake preparation, as it happens at the laboratory scale. Despite fairly abundant evidence, electromagnetic (EM) precursors have not been adequately accepted as credible physical phenomena. These negative views are enhanced by the fact that certain "puzzling features" are repetitively observed in candidate fracture-induced pre-seismic EME. More precisely, EM silence in all frequency bands appears before the main seismic shock occurrence, as well as during the aftershock period. Actually, the view that "acceptance of "precursive" EM signals without convincing co-seismic signals should not be expected" seems to be reasonable. In this work we focus on this point. We examine whether the aforementioned features of EM silence are really puzzling ones or, instead, reflect well-documented characteristic features of the fracture process, in terms of universal structural patterns of the fracture process, recent laboratory experiments, numerical and theoretical studies of fracture dynamics, critical phenomena, percolation theory, and micromechanics of granular materials. Our analysis shows that these features should not be considered puzzling.

Description: Migrating subtropical front and Agulhas Return Current affect the southwestern Indian Ocean during the late Quaternary Climate of the Past Discussions, 9, 5521-5551, 2013 Author(s): D. K. Naik, R. Saraswat, N. Khare, A. C. Pandey, and R. Nigam The position of sub-tropical front (STF), Agulhas Current (AC) and Agulhas Return Current (ARC) controls the hydrography of southwestern Indian Ocean. Although, equator-ward migration of STF and reduction in Agulhas leakage has been reported during the last glacial period, the fate of ARC during the last glacial–interglacial cycle is not clear. Therefore, in order to understand changes in the position and strength of ARC during the last glacial–interglacial cycle, here we reconstruct hydrographic changes in the southwestern Indian Ocean from temporal variation in planktic foraminiferal abundance, stable isotopic ratio (δ 18 O) and trace metal ratio (Mg / Ca) of planktic foraminifera Globigerina bulloides in a core collected from the Agulhas Retroflection Region (ARR) in the southwestern Indian Ocean. Increased abundance of G. bulloides suggests that the productivity in the southwestern Indian Ocean increased during glacial period which confirms previous reports of high glacial productivity in the Southern Ocean. The increased productivity was likely driven by a combination of equator-ward migration of subtropical front and westerlies. Increase in relative abundance of Neogloboquadrina pachyderma Dextral suggests warming of ARR leading to strong thermocline in the southwestern Indian Ocean during the last glacial period. We suggest that the warming of Agulhas Retroflection Region was driven by strengthened ARC which shifted to the east of its present location, thus bringing warmer and saltier water to the southwestern Indian Ocean. Therefore, it is inferred that over the last glacial–interglacial cycle, the hydrography of southwestern Indian Ocean was driven by an eastward shift of retroflection region as well as migrating subtropical front.

Description: The response of the Peruvian Upwelling Ecosystem to centennial-scale global change during the last two millennia Climate of the Past Discussions, 9, 5479-5519, 2013 Author(s): R. Salvatteci, D. Gutiérrez, D. Field, A. Sifeddine, L. Ortlieb, I. Bouloubassi, M. Boussafir, H. Boucher, and F. Cetin The Tropical Pacific ocean-atmosphere system influences global climate on interannual, decadal, as well as at longer timescales. Given the uncertainties in the response of the Tropical Pacific to the ongoing greenhouse effect, it is important to assess the natural range of the Tropical Pacific climate variability in response to global natural changes, and to understand the underlying mechanisms. The Peruvian Upwelling Ecosystem (PUE) represents an ideal area to reconstruct past changes in ocean-atmosphere systems because productivity and subsurface oxygenation are strongly linked to changes in the strength of the Walker circulation. Throughout the last 2000 yr, warmer (the Roman Warm Period [RWP], the Medieval Climate Anomaly [MCA] and the Current Warm Period [CWP]), and colder (the Dark Ages Cold Period [DACP] and Little Ice Age [LIA]) intervals occurred with considerable changes around the globe. In order to reconstruct the PUE response to these climatic periods and reveal the underlying mechanisms, we use a multi-proxy approach including organic and inorganic proxies in finely laminated sediments retrieved off Pisco (~ 14° S), Peru. Our results indicate that the PUE exhibited a La Niña-like mean state during the warm periods, characterized by an intense OMZ and high marine productivity. During cold periods the PUE exhibited an El Niño-like mean state, characterized by a weak OMZ and low marine productivity. Comparing our results with other relevant paleoclimatic reconstructions revealed that changes in the strength of the Walker circulation and the expansion/contraction of the South Pacific Sub-tropical High controlled productivity and subsurface oxygenation in the PUE during the last two millennia. This indicate that large scale circulation changes are the driving forces in maintaining productivity and subsurface oxygenation off Peru at centennial time scales during the past two millennia.

Description: Cumulated insolation: a simple explanation of Milankovitch's forcing on climate changes Climate of the Past Discussions, 9, 5553-5568, 2013 Author(s): F. Marra The occurrence of the sudden melting of the ice sheets during the glacial terminations is explained in this paper as the consequence of the combined role of the minima and the maxima of mean summer insolation on the Northern Hemisphere, providing a new contribution to understand the mechanisms ruling glacial forcing. Indeed, no satisfactory answer has been provided so far to the question why one specific maximum, after a series of consecutive maxima of insolation, has the potentiality to trigger a deglaciation. The explanation proposed in this paper accounts for a pre-conditioning factor, represented by "mild" (warmer) minimum, followed by a sufficiently warm maximum as the conditions that cause the end of a glacial cycle. These conditions are realized whenever the sum of the values of each consecutive minima and maxima ("cumulated insolation") on the curve of mean summer insolation at 65° N exceeds 742 Watt m −2 . The comparison of the succession of these cumulated insolation values with the astronomically tuned Oxygen isotopes record provides a satisfactory match with the occurrence of all the glacial terminations in the last 800 ka.

Description: Early warning of snow-caused disasters in pastoral areas on the Tibetan Plateau Natural Hazards and Earth System Science, 13, 1411-1425, 2013 Author(s): W. Wang, T. Liang, X. Huang, Q. Feng, H. Xie, X. Liu, M. Chen, and X. Wang This study develops a model for early warning of snow-caused livestock disasters on a county basis and proposes a method of qualitative risk assessment of snow disasters at 500 m resolution for pastoral areas on the Tibetan Plateau (TP). Data used for the model development include remote sensing data, statistical data of weather, livestock, and social economy, and 45 typical snow disaster cases from 2000 to 2010. The principal component analysis (PCA) approach is used to choose 7 crucial factors that contribute over 85% of information for early warning snow disasters on the TP. They are mean annual probability of snow disaster, number of snow-covered days, livestock stocking rate, continual days of mean daily temperature below −10 °C, grassland burial index, rate of snow-covered grassland, and per livestock gross domestic product. The chosen 411 cases from 2008 to 2010 are used to validate the prediction results from the developed early warning model, with an overall accuracy of 85.64% in predicting snow disasters and no disasters. This suggests that the early warning approach developed in the study has operational potential for predicting snow disasters on the TP.

Description: How much does participatory flood management contribute to stakeholders' social capacity building? Empirical findings based on a triangulation of three evaluation approaches Natural Hazards and Earth System Science, 13, 1427-1444, 2013 Author(s): M. Buchecker, S. Menzel, and R. Home Recent literature suggests that dialogic forms of risk communication are more effective to build stakeholders' hazard-related social capacities. In spite of the high theoretical expectations, there is a lack of univocal empirical evidence on the relevance of these effects. This is mainly due to the methodological limitations of the existing evaluation approaches. In our paper we aim at eliciting the contribution of participatory river revitalisation projects on stakeholders' social capacity building by triangulating the findings of three evaluation studies that were based on different approaches: a field-experimental, a qualitative long-term ex-post and a cross-sectional household survey approach. The results revealed that social learning and avoiding the loss of trust were more relevant benefits of participatory flood management than acceptance building. The results suggest that stakeholder involvements should be more explicitly designed as tools for long-term social learning.

Description: Stakeholders' issues for action during the warning process and the interpretation of forecasts' uncertainties Natural Hazards and Earth System Science, 13, 1469-1479, 2013 Author(s): L. Créton-Cazanave and C. Lutoff This article considers the socio-geographical approach carried out as part of the MedUp program. It presents a study of the ways that the French "actors" manage forecast uncertainties during a flash flood warning process. In order to better understand the role of forecasts' uncertainties in decision making, we focused on the actions people took and how what they say explains their actions. The practices of actors involved in warnings for the Vidourle watershed (Gard, France), in particular, are analyzed using a practice-based approach. A set of categories of the "actors" was developed based on their descriptions of the problems they faced during the flash flood warning, independent of their socio-professional status and position in the warning chain. Five actor profiles result from this: Translators, Managers, Committed, Navigators and Vulnerable. For each profile, specific action contexts are defined, determining how each deals with uncertainty.

Description: Augmenting Austrian flood management practices through geospatial predictive analytics: a study in Carinthia Natural Hazards and Earth System Science, 13, 1445-1455, 2013 Author(s): S. M. Ward and G. Paulus The Danube River basin has long been the location of significant flooding problems across central Europe. The last decade has seen a sharp increase in the frequency, duration and intensity of these flood events, unveiling a dire need for enhanced flood management policy and tools in the region. Located in the southern portion of Austria, the state of Carinthia has experienced a significant volume of intense flood impacts over the last decade. Although the Austrian government has acknowledged these issues, their remedial actions have been primarily structural to date. Continued focus on controlling the natural environment through infrastructure while disregarding the need to consider alternative forms of assessing flood exposure will only act as a provisional solution to this inescapable risk. In an attempt to remedy this flaw, this paper highlights the application of geospatial predictive analytics and spatial recovery index as a proxy for community resilience, as well as the cultural challenges associated with the application of foreign models within an Austrian environment.

Description: Chronology of Lake El'gygytgyn sediments Climate of the Past Discussions, 9, 3061-3102, 2013 Author(s): N. R. Nowaczyk, E. M. Haltia, D. Ulbricht, V. Wennrich, M. A. Sauerbrey, P. Rosén, H. Vogel, A. Francke, C. Meyer-Jacob, A. A. Andreev, and A. V. Lozhkin A 318 m long sedimentary profile drilled by the International Continental Scientific Drilling Program (ICDP) at Site 5011-1 in Lake El'gygytgyn, Far East Russian Arctic, has been analysed for its sedimentologic response to global climate modes by chrono-stratigraphic methods. The 12 km wide lake is sited in an 18 km large crater that was created by the impact of a meteorite 3.58 Ma ago. Since then sediments have been continuously deposited. For establishing their chronology, major reversals of the Earth's magnetic field provided initial tie points for the age model, confirming that the impact occurred in the earliest Gauss chron. Various stratigraphic parameters, reflecting redox conditions at the lake floor and climatic conditions in the catchment were tuned synchronously to Northern Hemisphere insolation variations and the marine oxygen isotope stack, respectively. Thus, a robust age model comprising more than 600 tie points could be defined. It could be shown that deposition of sediments in Lake El'gygytgyn occurred in concert with global climatic cycles. The upper ~160 m of sediments represent the past 3.3 Ma, equivalent to sedimentation rates of 4 to 5 cm ka −1 , whereas the lower 160 m represent just the first 0.3 Ma after the impact, equivalent to sedimentation rates in the order of 45 cm ka −1 .

Description: Late Cenozoic continuous aridification in the western Qaidam Basin: evidence from sporopollen records Climate of the Past Discussions, 9, 1485-1508, 2013 Author(s): Y. F. Miao, X. M. Fang, F. L. Wu, M. T. Cai, C. H. Song, Q. Q. Meng, and L. Xu Cenozoic climate changes in inner Asia provide a basis for understanding linkages between global cooling, the Tibetan Plateau uplift, and possibly the development of the East Asian monsoon. Based on the compiled palynological results from the western Qaidam Basin, this study reconstructed an 18 Ma record of changing vegetation and paleoclimates since the middle Miocene. Thermophilic taxa percentages were highest between 18 and 14 Ma and decreased after 14 Ma, corresponding closely with the Middle Miocene Climatic Optimum (MMCO) between 18 and 14 Ma and the following global climatic cooling. After 3.6 Ma, the thermophilic taxa percentages further decreased, showing the inevitable relations with the ice-sheets enlargement in the North Hemisphere. During the same period of time, the increase in xerophytic taxa percentages and decrease in conifers percentages imply aridification in both the basin and surrounding mountains since 18 Ma. These results indicate that global cooling mainly controlled the climate change from a relative warm-wet stage to a cold-dry stage during the late Cenozoic at the western Qaidam Basin, and that the Tibetan Plateau uplift also contributed in contrast to the East Asian summer monsoon.

Description: Post-Pliocene establishment of the present monsoonal climate in SW China: evidence from the late Pliocene Longmen megaflora Climate of the Past Discussions, 9, 1675-1701, 2013 Author(s): T. Su, F. M. B. Jacques, R. A. Spicer, Y.-S. Liu, Y.-J. Huang, Y.-W. Xing, and Z.-K. Zhou The paleoclimate of the late Pliocene Longmen flora from Yongping County located at the southeastern boundary of the Qinghai-Tibet Plateau was reconstructed using two leaf physiognomy based methods, i.e. Leaf Margin Analysis (LMA) and Climate Leaf Analysis Multivariate Program (CLAMP), to understand the paleoclimate condition and geographical pattern of monsoonal climate in southwestern China during the late Pliocene. The mean annual temperatures (MATs) estimated by LMA and CLAMP are 17.4 ± 3.3 °C and 17.4 ± 1.3 °C, respectively, compared with 15.9 °C at present. Meanwhile, the growing season precipitation (GSP) estimated by CLAMP is 1735.5 ± 217.7 mm in the Longmen flora, compared with 986.9 mm nowadays. The calculated monsoon index (MSI) of the Longmen flora is significantly lower than that of today. These results appear consistent with previous studies based on the coexistence approach (CA), and further suggest that there was a slightly warmer and much wetter climate during the late Pliocene than the present climate in western Yunnan. We conclude that the significant change of the monsoonal climate might have been resulted from the continuous uplift of mountains in western Yunnan, as well as the intensification of eastern Asian winter monsoon, both occurring concurrently in the post-Pliocene period.

Description: CDRD and PNPR satellite passive microwave precipitation retrieval algorithms: EuroTRMM/EURAINSAT origins and H-SAF operations Natural Hazards and Earth System Science, 13, 887-912, 2013 Author(s): A. Mugnai, E. A. Smith, G. J. Tripoli, B. Bizzarri, D. Casella, S. Dietrich, F. Di Paola, G. Panegrossi, and P. Sanò Satellite Application Facility on Support to Operational Hydrology and Water Management (H-SAF) is a EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites) program, designed to deliver satellite products of hydrological interest (precipitation, soil moisture and snow parameters) over the European and Mediterranean region to research and operations users worldwide. Six satellite precipitation algorithms and concomitant precipitation products are the responsibility of various agencies in Italy. Two of these algorithms have been designed for maximum accuracy by restricting their inputs to measurements from conical and cross-track scanning passive microwave (PMW) radiometers mounted on various low Earth orbiting satellites. They have been developed at the Italian National Research Council/Institute of Atmospheric Sciences and Climate in Rome (CNR/ISAC-Rome), and are providing operational retrievals of surface rain rate and its phase properties. Each of these algorithms is physically based, however, the first of these, referred to as the Cloud Dynamics and Radiation Database (CDRD) algorithm, uses a Bayesian-based solution solver, while the second, referred to as the PMW Neural-net Precipitation Retrieval (PNPR) algorithm, uses a neural network-based solution solver. Herein we first provide an overview of the two initial EU research and applications programs that motivated their initial development, EuroTRMM and EURAINSAT (European Satellite Rainfall Analysis and Monitoring at the Geostationary Scale), and the current H-SAF program that provides the framework for their operational use and continued development. We stress the relevance of the CDRD and PNPR algorithms and their precipitation products in helping secure the goals of H-SAF's scientific and operations agenda, the former helpful as a secondary calibration reference to other algorithms in H-SAF's complete mix of algorithms. Descriptions of the algorithms' designs are provided including a few examples of their performance. This aspect of the development of the two algorithms is placed in the context of what we refer to as the TRMM era, which is the era denoting the active and ongoing period of the Tropical Rainfall Measuring Mission (TRMM) that helped inspire their original development. In 2015, the ISAC-Rome precipitation algorithms will undergo a transformation beginning with the upcoming Global Precipitation Measurement (GPM) mission, particularly the GPM Core Satellite technologies. A few years afterward, the first pair of imaging and sounding Meteosat Third Generation (MTG) satellites will be launched, providing additional technological advances. Various of the opportunities presented by the GPM Core and MTG satellites for improving the current CDRD and PNPR precipitation retrieval algorithms, as well as extending their product capability, are discussed.

Description: Holocene vegetation and climate changes in central Mediterranean inferred from a high-resolution marine pollen record (Adriatic Sea) Climate of the Past Discussions, 9, 1969-2014, 2013 Author(s): N. Combourieu-Nebout, O. Peyron, V. Bout-Roumazeilles, S. Goring, I. Dormoy, S. Joannin, L. Sadori, G. Siani, and M. Magny To understand the effects of future climate change on the ecology of the central Mediterranean we can look to the impacts of long-term, millennial to centennial-scale climatic variability on vegetation in the basin. Pollen data from the Adriatic Marine core MD 90-917 allows us to reconstruct vegetation and regional climate changes over the south central Mediterranean during the Holocene. Clay mineral ratios from the same core reflect the relative contributions of riverine (illite and smectite) and eolian (kaolinite) contributions to the site, and thus act as an additional proxy with which to test precipitation changes in the Holocene. Vegetation reconstruction shows vegetation responses to the late-Glacial Preboreal oscillation, most likely driven by changes in seasonal precipitation. Pollen-inferred temperature declines during the early-mid Holocene, but increases during the mid-late Holocene, similar to southern-western Mediterranean climatic patterns during the Holocene. Several short climatic events appear in the record, indicating the sensitivity of vegetation in the region to millennial-scale variability. Reconstructed summer precipitation shows a regional maximum between 8000 and 7000 cal yr BP similar to the general pattern across southern Europe. Two important shifts in vegetation occur at 7700 and between 7500 and 7000 yr. These vegetation shifts are linked to changes in seasonal precipitation and are correlated to increased river inputs respectively from the north (7700 event) and from the central Adriatic borderlands (7500–7000 event). These results reinforce the strengths of multi-proxy analysis and provide a deeper understanding of the role of precipitation and particularly the seasonality of precipitation in mediating vegetation change in the central Mediterranean during the Holocene.

Description: Southern Hemisphere orbital forcing and its effects on CO 2 and tropical Pacific climate Climate of the Past Discussions, 9, 1869-1900, 2013 Author(s): K. Tachikawa, A. Timmermann, L. Vidal, C. Sonzogni, and O. E. Timm The western Pacific warm pool (WPWP) is an important heat source for the atmospheric circulation and influences climate conditions worldwide. Understanding its sensitivity to past radiative perturbations may help better contextualize the magnitudes and patterns of current and projected tropical climate change. Here we present a new Mg/Ca-based sea surface temperature (SST) reconstruction over the past 400 kyr from the Bismarck Sea, off Papua New Guinea, along with results from a transient earth system model simulation. Our results document the primary influence of CO 2 forcing on glacial/interglacial WPWP SSTs and secondary effects due to changes in wind-driven tropical boundary currents. In addition to the SST, deep ocean temperature reconstructions from this core are linked with Southern Ocean temperature and sea-ice variations on timescales of ~ 23 kyr. It is proposed that Southern Hemisphere insolation changes serve as pacemaker for sea-ice variations in the Southern Ocean, which in turn modulate windstress curl-driven upwelling of carbon-rich waters, hence controlling atmospheric CO 2 and tropical WPWP temperatures.

Description: A regional climate palaeosimulation for Europe in the period 1500–1990 – Part 1: Model validation Climate of the Past Discussions, 9, 1803-1839, 2013 Author(s): J. J. Gómez-Navarro, J. P. Montávez, S. Wagner, and E. Zorita We present and analyse a high-resolution regional climate palaeosimulation encompassing the European region for the period 1500–1990. We use the regional model MM5 coupled to the global model ECHO-G. Both models were driven by reconstructions of three external factors: greenhouse gas concentrations, Total Solar Irradiance and volcanic activity. The simulation has been assessed in a recent period by comparing the model results with the Climate Research Unit (CRU) database. The results show that although the regional model is tightly driven by the boundary conditions, it is able to improve the reliability of the simulations, narrowing the differences to the observations, especially in areas of complex topography. Additionally, the evolution of the spatial distributions of temperature and precipitation through the last five centuries has been analysed. The mean values of temperature reflects the influence of the external forcings but, contrary to the results obtained under climate change scenario conditions, we found that higher-order momenta of the probability distribution of seasonal temperature and precipitation are hardly affected by changes in the external forcings

Description: North–south palaeohydrological contrasts in the central Mediterranean during the Holocene: tentative synthesis and working hypotheses Climate of the Past Discussions, 9, 1901-1967, 2013 Author(s): M. Magny, N. Combourieu Nebout, J. L. de Beaulieu, V. Bout-Roumazeilles, D. Colombaroli, S. Desprat, A. Francke, S. Joannin, O. Peyron, M. Revel, L. Sadori, G. Siani, M. A. Sicre, S. Samartin, A. Simonneau, W. Tinner, B. Vannière, B. Wagner, G. Zanchetta, F. Anselmetti, E. Brugiapaglia, E. Chapron, M. Debret, M. Desmet, J. Didier, L. Essallami, D. Galop, A. Gilli, J. N. Haas, N. Kallel, L. Millet, A. Stock, J. L. Turon, and S. Wirth On the basis of a multi-proxy approach and a strategy combining lacustrine and marine records along a north–south transect, data collected in the Central Mediterranean within the framework of a collaborative project have led to reconstruction of high-resolution and well-dated palaeohydrological records and to assessment of their spatial and temporal coherency. Contrasting patterns of palaeohydrological changes have been evidenced in the Central Mediterranean: south (north) of around 40° N of latitude, the middle part of the Holocene was characterised by lake-level maxima (minima), during an interval dated to ca. 10 300–4500 cal BP to the south and 9000–4500 cal BP to the north. Available data suggest that these contrasting palaeohydrological patterns operated throughout the Holocene, both on millennial and centennial scales. Regarding precipitation seasonality, maximum humidity in the Central Mediterranean during the middle part of the Holocene was characterised by humid winters and dry summers north of ca. 40° N, and humid winters and summers south of ca. 40° N. This may explain an apparent conflict between palaeoclimatic records depending on the proxies used for reconstruction as well as the synchronous expansion of tree species taxa with contrasting climatic requirements. In addition, south of ca. 40° N, the first millennium of the Holocene was characterised by very dry climatic conditions not only in the Eastern, but also in the Central and the Western Mediterranean zones as reflected by low lake levels and delayed reforestation. These results suggest that, in addition to the influence of the Nile discharge reinforced by the African monsoon, the deposition of Sapropel 1 has been favoured (1) by an increase in winter precipitation in the northern Mediterranean borderlands, and (2) by an increase in winter and summer precipitation in the southern Mediterranean area. The climate reversal following the Holocene climate optimum appears to have been punctuated by two major climate changes around 7500 and 4500 cal BP. In the Central Mediterranean, the Holocene palaeohydrological changes developed in response to a combination of orbital, ice-sheet and solar forcing factors. The maximum humidity interval in the south-central Mediterranean started at ca. 10 300 cal BP, in correlation with the decline (1) of the possible blocking effects of the North Atlantic anticyclone linked to maximum insolation, and/or (2) of the influence of the remnant ice sheets and fresh water forcing in the North Atlantic Ocean. In the north-central Mediterranean, the lake-level minimum interval began only around 9000 cal BP when the Fennoscandian ice-sheet disappeared and a prevailing positive NAO-type circulation developed in the North Atlantic area. The major palaeohydrological oscillation around 4500–4000 cal BP may be a non-linear response to the gradual decrease, with additional key seasonal and interhemispherical changes, in insolation. On a centennial scale, the successive climatic events which punctuated the entire Holocene in the central Mediterranean coincided with cooling events associated with deglacial outbursts in the North Atlantic area and decreases in solar activity during the interval 11 700–7000 cal BP, and to a possible combination of NAO-type circulation and solar forcing since ca. 7000 cal BP onwards. Thus, regarding the centennial-scale climatic oscillations, the Mediterranean Basin appears to have been strongly linked to the North Atlantic area and affected by solar activity over the entire Holocene. In addition to model experiments, a better understanding of forcing factors and past atmospheric circulation patterns behind the Holocene palaeohydrological changes in the Mediterranean area will require further investigation to establish additional high-resolution and well-dated records in selected locations around the Mediterranean Basin and in adjacent regions. Special attention should be paid to greater precision in the reconstruction, on millennial and centennial time scales, of changes in the latitudinal location of the limit between the northern and southern palaeohydrological Mediterranean sectors, depending on (1) the intensity and/or characteristics of climatic periods/oscillations (e.g. Holocene thermal maximum versus Neoglacial, as well as, for instance, the 8.2 ka event versus the 4 ka event or the Little Ice Age), and (2) on varying geographical conditions from the western to the eastern Mediterranean areas (longitudinal gradients).

Description: The last 7 millennia of vegetation and climate changes at Lago di Pergusa (central Sicily, Italy) Climate of the Past Discussions, 9, 2059-2094, 2013 Author(s): L. Sadori, E. Ortu, O. Peyron, G. Zanchetta, B. Vannière, M. Desmet, and M. Magny The aim of this study is to investigate climate changes and human activities under the lens of palynology. Based on a new high-resolution pollen sequence (PG2) from Lago di Pergusa (667 m a.s.l., central Sicily, Italy) covering the last 6700 yr, we propose a reconstruction of climate and landscape changes over the recent past in central Sicily. Compared to former studies from Lago di Pergusa (Sadori and Narcisi, 2001), this work provides a reconstruction of the evolution of vegetation and climate over the last millennia in central Sicily, indeed completing previous results with new data which is particularly detailed on the last 3000 yr. Joint actions of increasing dryness, climate oscillations, and human impact shaped the landscape of this privileged site. Lago di Pergusa, in fact, besides being the main inland lake of Sicily, is very sensitive to climate change and its territory was inhabited and exploited continuously since the prehistory. The lake sediments turned out to be a good observatory for the natural phenomena occurred in the last thousands of years. Results of the pollen-based study are integrated with changes in magnetic susceptibility and a tephra layer characterization. The tephra layer was shown to be related to the Sicanians' event, radiocarbon dated at 3055 ± 75 yr BP (Sadori and Narcisi, 2001). We performed palaeoclimate reconstructions by MAT and WA-PLS. Palaeoclimate reconstructions based on the core show important climate fluctuations throughout the Holocene. Climate reconstruction points out four phases of cooling and enhanced wetness in the last three millennia (2600–2000, 1650–1100, 850–550, 400–200 cal BP). This appears to be the evidence of local responses to global climate oscillations during the recent past.

Description: Dust and associated trace element fluxes in a firn core from the coastal East Antarctica and its linkages with the Southern Hemisphere climate variability over the last ~ 50 yr Climate of the Past Discussions, 9, 1841-1867, 2013 Author(s): C. M. Laluraj, M. Thamban, and K. Satheesan High-resolution records of dust and trace element fluxes were studied in a firn core from the coastal Dronning Maud Land (cDML) in East Antarctica to identify the influence of climate variability on accumulation of these components over the past ~ 50 yr. A doubling of dust deposition was observed since 1985, coinciding with a shift in the Southern Annular Mode (SAM) index to positive values and associated increase in the wind speed. Back-trajectories showed that an increase in dust deposition is associated with the air parcels originating from north-west of the site, possibly indicating its origin from the Patagonian region. Our results suggest that while multiple processes could have influenced the increased dust formation, shift in SAM had a dominant influence on its transport. It is observed that since the 1985s the strength of easterlies increased significantly over the cDML region, which could sink air and dust material to the region that were brought by the westerlies through mass compensation. The correlation between the dust flux and δ 18 O records further suggest that enhanced dust flux in the firn core occurred during periods of colder atmospheric temperature, which reduced the moisture content and increased dust fall. Interestingly, the timing and amplitude of the insoluble dust peaks matched remarkably well with the fluxes of Ba, Cr, Cu, and Zn confirming that dust was the main carrier/source of atmospheric trace elements to East Antarctica during the recent past.

Description: An inter-laboratory investigation of the Arctic sea ice biomarker proxy IP 25 in marine sediments: key outcomes and recommendations Climate of the Past Discussions, 9, 5263-5298, 2013 Author(s): S. T. Belt, T. A. Brown, L. Ampel, P. Cabedo-Sanz, K. Fahl, J. J. Kocis, G. Massé, A. Navarro-Rodriguez, J. Ruan, and Y. Xu We describe the results of an inter-laboratory investigation into the identification and quantification of the Arctic sea ice biomarker proxy IP 25 in marine sediments. 7 laboratories took part in the study, which consisted of the analysis of IP 25 in a series of sediment samples from different regions of the Arctic, sub-Arctic and Antarctic, additional sediment extracts and purified standards. The results obtained allowed 4 key outcomes to be determined. First, IP 25 was identified by all laboratories in sediments from the Canadian Arctic with inter-laboratory variation in IP 25 concentration being substantially larger than within individual laboratories. This greater variation between laboratories was attributed to the difficulty in accurately determining instrumental response factors for IP 25 , despite provision of appropriate standards. Second, the identification of IP 25 by 3 laboratories in sediment from SW Iceland that was believed to represent a blank, was interpreted as representing a better limit of detection or quantification for such laboratories, contamination or mis-identification. These alternatives could not be distinguished conclusively with the data available, although it is noted that the precision of these data was significantly poorer compared with the other IP 25 concentration measurements. Third, 3 laboratories reported the occurrence of IP 25 in a sediment sample from the Antarctic Peninsula even though this biomarker is believed to be absent from the Southern Ocean. This anomaly is attributed to a combined chromatographic and mass spectrometric interference that results from the presence of a di-unsaturated highly branched isoprenoid (HBI) pseudo-homologue of IP 25 that occurs in Antarctic sediments. Finally, data are presented that suggest that extraction of IP 25 is consistent between Automated Solvent Extraction (ASE) and sonication methods and that IP 25 concentrations based on 7-hexylnonadecane as an internal standard are comparable using these methods. Recoveries of some more unsaturated HBIs and the internal standard 9-octylheptadecene, however, were lower with the ASE procedure, possibly due to partial degradation of these more reactive chemicals as a result of higher temperatures employed with this method. For future measurements, we recommend the use of reference sediment material with known concentration(s) of IP 25 for determining and routinely monitoring instrumental response factors. Given the significance placed on the presence (or otherwise) of IP 25 in marine sediments, some further recommendations pertaining to quality control are made that should also enable the two main anomalies identified here to be addressed.

Description: Expressions of climate perturbations in western Ugandan crater lake sediment records during the last 1000 yr Climate of the Past Discussions, 9, 5183-5226, 2013 Author(s): K. Mills, D. B. Ryves, N. J. Anderson, C. L. Bryant, and J. J. Tyler Equatorial East Africa has a complex, regional patchwork of climate regimes, with multiple interacting drivers. Recent studies have focussed on large lakes and reveal signals that are smoothed in both space and time, and, whilst useful at a continental scale, are of less relevance when understanding short-term, abrupt or immediate impacts of climate and environmental changes. Smaller-scale studies have highlighted spatial complexity and regional heterogeneity of tropical palaeoenvironments in terms of responses to climatic forcing (e.g. the Little Ice Age [LIA]) and questions remain over the spatial extent and synchroneity of climatic changes seen in East African records. Sediment cores from paired crater lakes in western Uganda were examined to assess ecosystem response to long-term climate and environmental change as well as testing responses to multiple drivers using redundancy analysis. These archives provide annual to sub-decadal records of environmental change. The records from the two lakes demonstrate an individualistic response to external (e.g. climatic) drivers, however, some of the broader patterns observed across East Africa suggest that the lakes are indeed sensitive to climatic perturbations such as a dry Mediaeval Climate Anomaly (MCA; 1000–1200 AD) and a relatively drier climate during the main phase of the LIA (1500–1800 AD); though lake levels in western Uganda do fluctuate. The relationship of Ugandan lakes to regional climate drivers breaks down c. 1800 AD, when major changes in the ecosystems appear to be a response to sediment and nutrient influxes as a result of increasing cultural impacts within the lake catchments. The data highlight the complexity of individual lake response to climate forcing, indicating shifting drivers through time. This research also highlights the importance of using multi-lake studies within a landscape to allow for rigorous testing of climate reconstructions, forcing and ecosystem response.

Description: Influence of management of variables, sampling zones and land units on LR analysis for landslide spatial prevision Natural Hazards and Earth System Science, 13, 2209-2221, 2013 Author(s): R. Greco and M. Sorriso-Valvo Several authors, according to different methodological approaches, have employed logistic Regression (LR), a multivariate statistical analysis adopted to assess the spatial probability of landslide, even though its fundamental principles have remained unaltered. This study aims at assessing the influence of some of these methodological approaches on the performance of LR, through a series of sensitivity analyses developed over a test area of about 300 km 2 in Calabria (southern Italy). In particular, four types of sampling (1 – the whole study area; 2 – transects running parallel to the general slope direction of the study area with a total surface of about 1/3 of the whole study area; 3 – buffers surrounding the phenomena with a 1/1 ratio between the stable and the unstable area; 4 – buffers surrounding the phenomena with a 1/2 ratio between the stable and the unstable area), two variable coding modes (1 – grouped variables; 2 – binary variables), and two types of elementary land (1 – cells units; 2 – slope units) units have been tested. The obtained results must be considered as statistically relevant in all cases (Aroc values 〉 70%), thus confirming the soundness of the LR analysis which maintains high predictive capacities notwithstanding the features of input data. As for the area under investigation, the best performing methodological choices are the following: (i) transects produced the best results (0 〈 P(y) ≤ 93.4%; Aroc = 79.5%); (ii) as for sampling modalities, binary variables (0 〈 P(y) ≤ 98.3%; Aroc = 80.7%) provide better performance than ordinated variables; (iii) as for the choice of elementary land units, slope units (0 〈 P(y) ≤ 100%; Aroc = 84.2%) have obtained better results than cells matrix.

Description: Risk evolution: how can changes in the built environment influence the potential loss of natural hazards? Natural Hazards and Earth System Science, 13, 2195-2207, 2013 Author(s): B. Schwendtner, M. Papathoma-Köhle, and T. Glade Alpine areas often suffer significant loss and damage due to a range of natural processes such as landslides, debris flows, snow avalanches or floods. Sealing of the soil surface, settling in endangered areas and enhanced human intervention in the natural settings, as well as socio-economic changes, increase the risk and susceptibility of built environments to natural hazards and the costs of the consequences in a spatio-temporal context. The present study examines the loss estimation of a particular debris flow event for different points in time. The event occurred in August 1987, affected the municipality Martell in South Tyrol, Italy, and resulted in a total cost of € 25 million. The approach presented in this paper focuses on the changes of the land use and settlement expansion in the area since 1954 and attempts to assess the monetary impact of a similar event, which could have happened before (1954, 1985) or following the actual event (1992, 1999, 2006). The method applied is based on the use of a vulnerability curve which was developed for the specific area, based on the documentation of the damage of the 1987 event. Based on this curve, a loss estimation was carried out in order to visualise the risk evolution in a period of 52 yr (1954 to 2006). The results show a significant increase in the extent of the built environment (number, size and value of buildings) which consequently reflect an increase of the potential overall loss through the years. The method can be used in order to assess the potential loss for future scenarios based on different spatial patterns of the built environment.

Description: Wind shear over the Nice Côte d'Azur airport: case studies Natural Hazards and Earth System Science, 13, 2223-2238, 2013 Author(s): A. Boilley and J.-F. Mahfouf The Nice Côte d'Azur international airport is subject to horizontal low-level wind shears. Detecting and predicting these hazards is a major concern for aircraft security. A measurement campaign took place over the Nice airport in 2009 including 4 anemometers, 1 wind lidar and 1 wind profiler. Two wind shear events were observed during this measurement campaign. Numerical simulations were carried out with Meso-NH in a configuration compatible with near-real time applications to determine the ability of the numerical model to predict these events and to study the meteorological situations generating an horizontal wind shear. A comparison between numerical simulation and the observation dataset is conducted in this paper.

Description: Explosive development of winter storm Xynthia over the subtropical North Atlantic Ocean Natural Hazards and Earth System Science, 13, 2239-2251, 2013 Author(s): M. L. R. Liberato, J. G. Pinto, R. M. Trigo, P. Ludwig, P. Ordóñez, D. Yuen, and I. F. Trigo In winter of 2009–2010 south-western Europe was hit by several destructive windstorms. The most important was Xynthia (26–28 February 2010), which caused 64 reported casualties and was classified as the 2nd most expensive natural hazard event for 2010 in terms of economic losses. In this work we assess the synoptic evolution, dynamical characteristics and the main impacts of storm Xynthia, whose genesis, development and path were very uncommon. Wind speed gusts observed at more than 500 stations across Europe are evaluated as well as the wind gust field obtained with a regional climate model simulation for the entire North Atlantic and European area. Storm Xynthia was first identified on 25 February around 30° N, 50° W over the subtropical North Atlantic Ocean. Its genesis occurred on a region characterized by warm and moist air under the influence of a strong upper level wave embedded in the westerlies. Xynthia followed an unusual SW–NE path towards Iberia, France and central Europe. The role of moist air masses on the explosive development of Xynthia is analysed by considering the evaporative sources. A lagrangian model is used to identify the moisture sources, sinks and moisture transport associated with the cyclone during its development phase. The main supply of moisture is located over an elongated region of the subtropical North Atlantic Ocean with anomalously high SST, confirming that the explosive development of storm Xynthia had a significant contribution from the subtropics.

Description: Seismic behavior of buried pipelines constructed by design criteria and construction specifications of both Korea and the US Natural Hazards and Earth System Science, 13, 2271-2278, 2013 Author(s): S.-S. Jeon Earthquake loss estimation systems in the US, for example HAZUS (Hazard in US), have been established based on sufficient damage records for the purpose of prevention and efficient response to earthquake hazards; however, in Korea, insufficient data sets of earthquakes and damage records are currently available. In this study, the earthquake damages to pipelines in Korea using the pipeline repair rate (RR) recommended in HAZUS was reevaluated with the degree of confidence when RR is used without modification for the damage estimation of pipelines in Korea. The numerical analyses using a commercial finite element model, ABAQUS, were carried out to compare stresses and strains mobilized in both brittle and ductile pipelines constructed by the design criteria and construction specifications of both Korea and the US. These pipelines were embedded in dense sand overlying three different in situ soils (clay, sand, and gravel) subjected to earthquake excitations with peak ground accelerations (PGAs) of 0.2 to 1.2 g and 1994 Northridge and 1999 Chi-Chi earthquake loadings. The numerical results show that differences in the stress and strain rates are less than 10%. This implies that RR in HAZUS can be used for earthquake damage estimation of pipelines with a 90% confidence level in Korea.

Description: Evaluation of seasonal climates of the Mediterranean and nothern Africa in the CMIP5 simulations Climate of the Past Discussions, 9, 5347-5389, 2013 Author(s): A. Perez-Sanz, G. Li, P. González-Sampériz, and S. P. Harrison We analyze the spatial expression of seasonal climates of the Mediterranean and northern Africa in pre-Industrial ( piControl ) and mid-Holocene ( midHolocene , 6 ka) simulations from the fifth phase of the Coupled Model Intercomparison Project (CMIP5). Modern observations show four distinct precipitation regimes characterized by differences in the seasonal distribution and total amount of precipitation: an equatorial band characterized by a double peak in rainfall, the monsoon zone characterized by summer rainfall, the desert characterized by low seasonality and total precipitation, and the Mediterranean zone characterized by summer drought. Most models correctly simulate the position of the Mediterranean and the equatorial climates in the piControl simulations, but over-estimate the extent of monsoon influence and underestimate the extent of desert. However, most models fail to reproduce the amount of precipitation in each zone. Model biases in the simulated magnitude of precipitation are unrelated to whether the models reproduce the correct spatial patterns of each regime. In the midHolocene , the models simulate a reduction in winter rainfall in the equatorial zone, and a northward expansion of the monsoon with a significant increase in summer and autumn rainfall. Precipitation is slightly increased in the desert, mainly in summer and autumn, with northward expansion of the monsoon. Changes in the Mediterranean are small, although there is an increase in spring precipitation consistent with palaeo-observations of increased growing-season rainfall. Comparison with reconstructions shows that most models under-estimate the mid-Holocene changes in annual precipitation, except in the equatorial zone. Biases in the piControl have only a limited influence on midHolocene anomalies in ocean-atmosphere models; carbon-cycle models show no relationship between piControl bias and midHolocene anomalies. Biases in the prediction of the midHolocene monsoon expansion are unrelated to how well the models simulate changes in Mediterranean climate.

Description: Orbital and millennial-scale environmental changes between 64 and 25 ka BP recorded in Black Sea sediments Climate of the Past Discussions, 9, 5439-5477, 2013 Author(s): L. S. Shumilovskikh, D. Fleitmann, N. R. Nowaczyk, H. Behling, F. Marret, A. Wegwerth, and H. W. Arz High-resolution pollen and dinoflagellate cyst records from marine sediment core 25-GC1 were used to reconstruct vegetation dynamics in Northern Anatolia and surface conditions of the Black Sea between 64 and 25 ka BP. During this period, the dominance of Artemisia in the pollen record indicates a steppe landscape and arid climate conditions. However, the presence of temperate and warm-temperate arboreal pollen suggests the existence of glacial refugia in Northern Anatolia. A general cooling trend towards 25 ka BP is evidenced by the decrease of Quercus and increase of Pinus . There is evidence of orbital-driven vegetation dynamics in Northern Anatolia during 64–25 ka BP with spread of steppe during precession minima (insolation maxima) and development of forests during precession maxima (insolation minima). Dansgaard–Oeschger (D–O) events are characterized by a marked increase in temperate tree pollen, indicating a spread of forests due to warm and wet conditions in Northern Anatolia. The dominance of Pyxidinopsis psilata and Spiniferites cruciformis in the dinocyst record indicates a rather brackish Black Sea during the last glacial period. The decrease of marine indicators (marine dinocysts, acritachs) at ~ 54 ka BP and increase of freshwater algae ( Pediastrum, Botryococcus ) from 32 to 25 ka BP reveals freshening of the Black Sea surface water, related to orbital-driven arid/humid phases in the region, influencing hydrology and level changes of the Black Sea. D–O interstadials are characterized by high dinocyst concentrations and calcium carbonate content, as a result of an increase in primary productivity in the Black Sea. Heinrich events show a similar impact on the environment in Northern Anatolia/Black Sea region as D–O stadials.

Description: Seismic zones for Azores based on statistical criteria Natural Hazards and Earth System Science, 13, 2337-2351, 2013 Author(s): M. C. M. Rodrigues and C. S. Oliveira The objective of this paper is to define seismic zones in the Azores based on statistical criteria. These seismic zones will likely be used in seismic simulations of occurrences in the Azores Archipelago. The data used in this work cover the time period from 1915 to 2011. The Azores region was divided into 1° × 1° area units, for which the seismicity and the maximum magnitudes of events were calculated. The seismicity, the largest earthquakes recorded and the geological characteristics of the region were used to group these area units because similar seismic zones must delineate areas with homogeneous seismic characteristics. We have identified seven seismic zones. To verify that the defined areas differ statistically, we considered the following dissimilarity measures (variables): time, size and seismic conditions – the number of seismic events with specific characteristics. Statistical tests, particularly goodness-of-fit tests, allowed us to conclude that, considering these three variables, the seven earthquake zones defined here are statistically distinct.

Description: Sub-ionospheric very low frequency perturbations associated with the 12 May 2008 M = 7.9 Wenchuan earthquake Natural Hazards and Earth System Science, 13, 2331-2336, 2013 Author(s): A. K. Maurya, R. Singh, B. Veenadhari, S. Kumar, and A. K. Singh The present study reports the VLF (very low frequency) sub-ionospheric perturbations observed on transmitter JJI (22.1 kHz), Japan, received at the Indian low-latitude station, Allahabad ( geographic lat. 25.41° N, long 81.93° E), due to Wenchuan earthquake (EQ) that occurred on 12 May 2008 with the magnitude 7.9 and at the depth of 19 km in Sichuan province of Southwest China, located at 31.0° N, 103.4° E. The nighttime amplitude fluctuation analysis gives a significant increase in fluctuation and dispersion two days before EQ, when it crosses 2σ criterion. However, there was no significant change observed in the amplitude trend. The diurnal amplitude variation shows a significant increase in the amplitude of JJI signal on 11 and 12 May 2008. The gravity wave channel and changes in the electric field associated with this EQ seem to be the potential factors of the observed nighttime amplitude fluctuation, dispersion, and significant increase in the signal strength.

Description: A catalog of high-impact windstorms in Switzerland since 1859 Natural Hazards and Earth System Science, 14, 2867-2882, 2014 Author(s): P. Stucki, S. Brönnimann, O. Martius, C. Welker, M. Imhof, N. von Wattenwyl, and N. Philipp In recent decades, extremely hazardous windstorms have caused enormous losses to buildings, infrastructure and forests in Switzerland. This has increased societal and scientific interest in the intensity and frequency of historical high-impact storms. However, high-resolution wind data and damage statistics mostly span recent decades only. For this study, we collected quantitative (e.g., volumes of windfall timber, losses relating to buildings) and descriptive (e.g., forestry or insurance reports) information on the impact of historical windstorms. To define windstorm severity, normalized and declustered quantitative data were processed by extreme value statistics. Descriptive information was classified using a conceptual guideline. Validation with independent damage information, as well as comparison with wind measurements and a reanalysis, indicates that the most hazardous winter storms are captured, while too few moderate windstorms are detected. Strong storms in the wind measurements and reanalysis are thus added to the catalog. The final catalog encompasses approximately 240 high-impact windstorms in Switzerland since 1859. It features three robust severity classes and contains eight extreme windstorms. Evidence of high winter storm activity in the early and late 20th century compared to the mid-20th century in both damage and wind data indicates a co-variability of hazard and related damage on decadal timescales.

Description: A GIS-based model to estimate flood consequences and the degree of accessibility and operability of strategic emergency response structures in urban areas Natural Hazards and Earth System Science, 14, 2847-2865, 2014 Author(s): R. Albano, A. Sole, J. Adamowski, and L. Mancusi Efficient decision-making regarding flood risk reduction has become a priority for authorities and stakeholders in many European countries. Risk analysis methods and techniques are a useful tool for evaluating costs and benefits of possible interventions. Within this context, a methodology to estimate flood consequences was developed in this paper that is based on GIS, and integrated with a model that estimates the degree of accessibility and operability of strategic emergency response structures in an urban area. The majority of the currently available approaches do not properly analyse road network connections and dependencies within systems, and as such a loss of roads could cause significant damages and problems to emergency services in cases of flooding. The proposed model is unique in that it provides a maximum-impact estimation of flood consequences on the basis of the operability of the strategic emergency structures in an urban area, their accessibility, and connection within the urban system of a city (i.e. connection between aid centres and buildings at risk), in the emergency phase. The results of a case study in the Puglia region in southern Italy are described to illustrate the practical applications of this newly proposed approach. The main advantage of the proposed approach is that it allows for defining a hierarchy between different infrastructure in the urban area through the identification of particular components whose operation and efficiency are critical for emergency management. This information can be used by decision-makers to prioritize risk reduction interventions in flood emergencies in urban areas, given limited financial resources.

Description: Wind waves in the Black Sea: results of a hindcast study Natural Hazards and Earth System Science, 14, 2883-2897, 2014 Author(s): V. S. Arkhipkin, F. N. Gippius, K. P. Koltermann, and G. V. Surkova In this study we describe the wind wave fields in the Black Sea. The general aims of the work were the estimation of statistical wave parameters and the assessment of interannual and seasonal wave parameter variability. The domain of this study was the entire Black Sea. Wave parameters were calculated by means of the SWAN wave model on a 5 × 5 km rectangular grid. Initial conditions (wind speed and direction) for the period between 1949 and 2010 were derived from the NCEP/NCAR reanalysis. According to our calculations the average significant wave height on the Black Sea does not exceed 0.7 m. Areas of most significant heavy sea are the southwestern and the northeastern parts of the sea as expressed in the spatial distribution of significant wave heights, wave lengths and periods. Besides, long-term annual variations of wave parameters were estimated. Thus, linear trends of the annual total duration of storms and of their quantity are nearly stable over the hindcast period. However, an intensification of storm activity is observed in the 1960s–1970s.

Description: Factors affecting flood insurance purchase in residential properties in Johor, Malaysia Natural Hazards and Earth System Science, 14, 3297-3310, 2014 Author(s): U. G. Aliagha, T. E. Jin, W. W. Choong, M. Nadzri Jaafar, and H. M. Ali High-impact floods have become a virtually annual experience in Malaysia, yet flood insurance has remained a grossly neglected part of comprehensive integrated flood risk management. Using discriminant analysis, this study seeks to identify the demand-side variables that best predict flood insurance purchase and risk aversion between two groups of residential homeowners in three districts of Johor State, Malaysia: those who purchased flood insurance and those who did not. Our results revealed an overall 34% purchase rate, with Kota Tinggi district having the highest (44%) and thus the highest degree of flood risk aversion. The Wilks' lambda F test for equality of group means, standardised discriminant function coefficients, structure correlation, and canonical correlation has clearly shown that there are strong significant attribute differences between the two groups of homeowners, based on the measures of objective flood risk exposure, subjective risk perception, and socio-economic cum demographic variables. However, the measures of subjective risk perception were found to be more predictive of flood insurance purchase and flood risk aversion.

Description: Kinematics of a mass movement constrained by sparse and inhomogeneous data Natural Hazards and Earth System Science, 11, 1609-1618, 2011 Author(s): M. Karbon, E. Brückl, E. Hegedüs, and A. Preh On 12 February 2008, a landslide occurred along a 50 m high bank of the Danube river near Dunaszekcsö, Hungary. The initial state is only incompletely documented and the geodetic data acquired after the mass movement are sparse. A generalized 3-D topographic model of the landslide and its surrounding area was assembled and a representative longitudinal profile extracted. The reconstruction of the original surface is based on an orthophoto as well as on morphological considerations. Recorded observations include the locations of the outcrops of basal sliding surfaces, displacements at the main scarp and in the lower part of the slide, and a value to describe the total mass transport. Such sparse and inhomogeneous data were insufficient to derive a comprehensive documentation of the landslide or obtain adequate constraints for an accurate numerical analysis. Therefore, slider block models were fitted to the field data, which have only a small number of free parameters. A general view on the morphology of the mass movement justifies its classification as a rotational slide. A double slider block model fits all observational parameters within their error margin and supplies valuable information on the geometry of the slide. Estimates of the residual friction angles were derived and the question of reactivation was addressed. Finite Difference (FD) modelling and the application of conventional stability analysis support the geometry of the slider blocks and the computed average residual friction angles. Generally, the results are assumed to represent preliminary information, which could only be attained by the combination of the thinly distributed geodetic data with qualitative morphological observations and the implementation of a model. This type of information can be gained quickly and may be valuable for preliminary hazard mitigation measures or the planning of a comprehensive exploration and monitoring program.

Description: Looking for evidence of climate change impacts in the eastern Irish Sea Natural Hazards and Earth System Science, 11, 1641-1656, 2011 Author(s): L. S. Esteves, J. J. Williams, and J. M. Brown Although storminess is often cited as a driver of long-term coastal erosion, a lack of suitable datasets has only allowed objective assessment of this claim in a handful of case studies. This reduces our ability to understand and predict how the coastline may respond to an increase in "storminess" as suggested by global and regional climate models. With focus on 16 km of the Sefton coastline bordering the eastern Irish Sea (UK), this paper analyses available measured datasets of water level, surge level, wave height, wind speed and barometric pressure with the objective of finding trends in metocean climate that are consistent with predictions. The paper then examines rates of change in shoreline position over the period 1894 to 2005 with the aim of establishing relationships with climatic variability using a range of measured and modelled metocean parameters (with time spans varying from two to eight decades). With the exception of the mean monthly wind speed, available metocean data do not indicate any statistically significant changes outside seasonal and decadal cycles. No clear relationship was found between changes in metocean conditions and rates of shoreline change along the Sefton coast. High interannual variability and the lack of long-term measurements make unambiguous correlations between climate change and shoreline evolution problematic. However, comparison between the North Atlantic Oscillation winter index (NAOw) and coastline changes suggest increased erosion at times of decreasing NAOw values and reduced erosion at times of increasing NAOw values. Erosion tends to be more pronounced when decreasing NAOw values lead to a strong negative NAO phase. At present, anthropogenic changes in the local sediment budget and the short-term impact of extreme events are still the largest threat likely to affect coastal flooding and erosion risk in the short- and medium-term. Nevertheless, the potential impacts of climate change in the long-term should not be ignored.

Description: Comparing historical-hydrogeomorphological reconstitution and hydrological-hydraulic modelling in the estimation of flood-prone areas – a case study in Central Portugal Natural Hazards and Earth System Science, 11, 1669-1681, 2011 Author(s): P. P. Santos, A. O. Tavares, and A. I. A. S. S. Andrade The Arunca River basin in Central Portugal has a historical record of hazardous events related to floods, causing widespread disturbance. This article describes the application of two approaches based on well-known methods for the estimation of flood-prone areas: (i) historical-hydrogeomorphological reconstitution, applied to the entire Arunca River basin, and (ii) hydrological-hydraulic modelling, applied to four sections selected from different (upper, middle and lower) sectors of the basin and including urban and rural areas along the Arunca River. The mapping of the flood-prone areas obtained by these two methods was compared in order to identify the main differences and similarities. Human interventions (river channel and floodplain morphological changes) were found to be the main factor explaining the differences and similarities between the results obtained by both methods. The application of hydrological-hydraulic modelling proved important in reinforcing the results of the historical-hydrogeomorphological method; it also helped in complementing the results produced by the latter method in urban areas and in areas with insufficient historical records. The application of the historical-hydrogeomorphological method, in turn, allowed for the size of the flood-prone areas to be determined where the primary data (e.g. geometry, roughness and flow) was not accurate enough for hydrological-hydraulic modelling. The methodological approach adopted demonstrates the strong complementary relationship between the different existing methods for estimating flood-prone areas, and may be reproduced for other drainage basins.

Description: The contribution of engineering surveys by means of GPS to the determination of crustal movements in Istanbul Natural Hazards and Earth System Science, 11, 1705-1713, 2011 Author(s): M. Özyaşar and M. T. Özlüdemir Global Navigation Satellite Systems (GNSS) are space based positioning techniques and widely used in geodetic applications. Geodetic networking accomplished by engineering surveys constitutes one of these tasks. Geodetic networks are used as the base of all kinds of geodetic implementations, Co from the cadastral plans to the relevant surveying processes during the realization of engineering applications. Geodetic networks consist of control points positioned in a defined reference frame. In fact, such positional information could be useful for other studies as well. One of such fields is geodynamic studies that use the changes of positions of control stations within a network in a certain time period to understand the characteristics of tectonic movements. In Turkey, which is located in tectonically active zones and struck by major earthquakes quite frequently, the positional information obtained in engineering surveys could be very useful for earthquake related studies. For this purpose, a GPS (Global Positioning System) network of 650 stations distributed over Istanbul (Istanbul GPS Triangulation Network; abbreviated IGNA) covering the northern part of the North Anatolian Fault Zone (NAFZ) was established in 1997 and measured in 1999. From 1998 to 2004, the IGNA network was extended to 1888 stations covering an area of about 6000 km 2 , the whole administration area of Istanbul. All 1888 stations within the IGNA network were remeasured in 2005. In these two campaigns there existed 452 common points, and between these two campaigns two major earthquakes took place, on 17 August and 12 November 1999 with a Richter scale magnitude of 7.4 and 7.2, respectively. Several studies conducted for estimating the horizontal and vertical displacements as a result of these earthquakes on NAFZ are discussed in this paper. In geodynamic projects carried out before the earthquakes in 1999, an annual average velocity of 2–2.5 cm for the stations along the NAFZ were estimated. Studies carried out using GPS observations in the same area after these earthquakes indicated that point displacements vary depending on their distance to the epicentres of the earthquakes. But the directions of point displacements are similar. The results obtained through the analysis of the IGNA network also show that there is a common trend in the directions of point displacements in the study area. In this paper, the past studies about the tectonics of Marmara region are summarised and the results of the displacement analysis on the IGNA network are discussed.

Description: Tropical climate and vegetation changes during Heinrich Event 1: comparing climate model output to pollen-based vegetation reconstructions with emphasis on the region around the tropical Atlantic Ocean Climate of the Past Discussions, 7, 1973-2019, 2011 Author(s): D. Handiani, A. Paul, and L. Dupont Abrupt climate changes associated with Heinrich Event 1 (HE1) about 18 to 15 thousand years before present (ka BP) strongly affected climate and vegetation patterns not only in the Northern Hemisphere, but also in tropical regions in the South Atlantic Ocean. We used the University of Victoria (UVic) Earth System-Climate Model (ESCM) with dynamical vegetation and land surface components to simulate four scenarios of climate-vegetation interaction: the pre-industrial era (PI), the Last Glacial Maximum (LGM), and a Heinrich-like event with two different climate backgrounds (interglacial and glacial). The HE1-like simulation with a glacial climate background produced sea surface temperature patterns and enhanced interhemispheric thermal gradients in accordance with the "bipolar seesaw" hypothesis. It allowed us to investigate the vegetation changes that result from a transition to a drier climate as predicted for northern tropical Africa due to a southward shift of the Intertropical Convergence Zone (ITCZ). We found that a cooling of the Northern Hemisphere caused a southward shift of those plant-functional types (PFTs) in Northern Tropical Africa that are indicative of an increased desertification, and a retreat of broadleaf forests in Western Africa and Northern South America. We used the PFTs generated by the model to calculate mega-biomes to allow for a direct comparison between paleodata and palynological vegetation reconstructions. Our calculated mega-biomes for the pre-industrial period and the LGM corresponded well to the modern and LGM sites of the BIOME6000 (v.4.2) reconstruction, except that our present-day simulation predicted the dominance of grassland in Southern Europe and our LGM simulation simulated more forest cover in tropical and sub-tropical South America. The mega-biomes from the HE1 simulation with glacial background climate were in agreement with paleovegetation data from land and ocean proxies in West, Central, and Northern Tropical Africa as well as Northeast South America. However, our model did not agree well with predicted biome distributions in Eastern South America.

Description: Quantitative estimation of landslide risk from rapid debris slides on natural slopes in the Nilgiri hills, India Natural Hazards and Earth System Science, 11, 1723-1743, 2011 Author(s): P. Jaiswal, C. J. van Westen, and V. Jetten A quantitative procedure for estimating landslide risk to life and property is presented and applied in a mountainous area in the Nilgiri hills of southern India. Risk is estimated for elements at risk located in both initiation zones and run-out paths of potential landslides. Loss of life is expressed as individual risk and as societal risk using F-N curves, whereas the direct loss of properties is expressed in monetary terms. An inventory of 1084 landslides was prepared from historical records available for the period between 1987 and 2009. A substantially complete inventory was obtained for landslides on cut slopes (1042 landslides), while for natural slopes information on only 42 landslides was available. Most landslides were shallow translational debris slides and debris flowslides triggered by rainfall. On natural slopes most landslides occurred as first-time failures. For landslide hazard assessment the following information was derived: (1) landslides on natural slopes grouped into three landslide magnitude classes, based on landslide volumes, (2) the number of future landslides on natural slopes, obtained by establishing a relationship between the number of landslides on natural slopes and cut slopes for different return periods using a Gumbel distribution model, (3) landslide susceptible zones, obtained using a logistic regression model, and (4) distribution of landslides in the susceptible zones, obtained from the model fitting performance (success rate curve). The run-out distance of landslides was assessed empirically using landslide volumes, and the vulnerability of elements at risk was subjectively assessed based on limited historic incidents. Direct specific risk was estimated individually for tea/coffee and horticulture plantations, transport infrastructures, buildings, and people both in initiation and run-out areas. Risks were calculated by considering the minimum, average, and maximum landslide volumes in each magnitude class and the corresponding minimum, average, and maximum run-out distances and vulnerability values, thus obtaining a range of risk values per return period. The results indicate that the total annual minimum, average, and maximum losses are about US$ 44 000, US$ 136 000 and US$ 268 000, respectively. The maximum risk to population varies from 2.1 × 10 −1 for one or more lives lost to 6.0 × 10 −2 yr −1 for 100 or more lives lost. The obtained results will provide a basis for planning risk reduction strategies in the Nilgiri area.

Description: Bridging the Faraoni and Selli oceanic anoxic events: short and repetitive dys- and anaerobic episodes during the late Hauterivian to early Aptian in the central Tethys Climate of the Past Discussions, 7, 2021-2059, 2011 Author(s): K. B. Föllmi, M. Bôle, N. Jammet, P. Froidevaux, A. Godet, S. Bodin, T. Adatte, V. Matera, D. Fleitmann, and J. E. Spangenberg A detailed stratigraphical and geochemical analysis was performed on the upper part of the Maiolica Formation outcropping in the Breggia (southern Switzerland) and Capriolo sections (northern Italy). In these localities, the Maiolica Formation consists of well-bedded, partly siliceous, pelagic, micritic carbonate, which lodges numerous thin, dark and organic-rich layers. Stable-isotope, phosphorus, organic-carbon and a suite of redox-sensitive trace-metal contents (RSTE: Mo, U, Co, V and As) were measured. Higher densities of organic-rich layers were identified in the uppermost Hauterivian, lower Barremian and the Barremian-Aptian boundary intervals, whereas the upper Barremian interval and the interval immediately following the Barremian-Aptian boundary interval are characterized by lower densities of organic-rich layers. TOC contents, RSTE pattern and C org :P tot ratios indicate that most layers were deposited under dysaerobic rather than anaerobic conditions and that latter conditions were likely restricted to short intervals in the latest Hauterivian, the early Barremian and the pre-Selli early Aptian. Correlations are possible with organic-rich intervals in central Italy (the Gorgo a Cerbara section) and the Boreal northwest German Basin, and with the facies and drowning pattern in the evolution of the Helvetic segment of the northern Tethyan carbonate platform. Our data and correlations suggest that the latest Hauterivian witnessed the progressive installation of dysaerobic conditions in the Tethys, which went along with the onset in sediment condensation, phosphogenesis and platform drowning on the northern Tethyan margin, and which culminated in the Faraoni anoxic episode. This brief episode is followed by further episodes of dysaerobic conditions in the Tethys and the northwest German Basin, which became more frequent and progressively stronger in the late early Barremian. Platform drowning persisted and did not halt before the latest early Barremian. The late Barremian witnessed diminishing frequencies and intensities in dysaerobic conditions, which went along with the progressive installation of the Urgonian carbonate platform. Near the Barremian-Aptian boundary, the increasing density in dysaerobic episodes in the Tethyan and northwest German Basins is paralleled by a change towards heterozoan carbonate production on the northern Tethyan shelf. The following return to more oxygenated conditions is correlated with the second phase of Urgonian platform growth and the period immediately preceding and corresponding to the Selli anoxic episode is characterized by renewed platform drowning and the change to heterozoan carbonate production. Changes towards more humid climate conditions were likely the cause for the repetitive installation of dys- to anaerobic conditions in the Tethyan and Boreal basins and the accompanying changes in the evolution of the carbonate platform towards heterozoan carbonate-producing ecosystems and platform drowning.

Description: Degree-day melt models for paleoclimate reconstruction from tropical glaciers: calibration from mass balance and meteorological data of the Zongo glacier (Bolivia, 16 ° S) Climate of the Past Discussions, 7, 2119-2158, 2011 Author(s): P.-H. Blard, P. Wagnon, J. Lavé, A. Soruco, J.-E. Sicart, and B. Francou This paper describes several simple positive degree-day models (hereafter referred as "PDD models") designed to provide past climatic reconstruction from tropical glacier paleo-equilibrium altitude lines (paleo-ELA). Several ablation laws were tested and calibrated using the monthly ablation and meteorological data recorded from 1997 to 2006 on the Zongo glacier (Cordillera Real, Bolivia, 16° S). The performed inversion analyses indicate that the model provides a better reconstruction of the mass balance if the ablation is modeled with different melting factors for snow and ice. The inclusion of short-wave solar radiations does not induce a substantial improvement. However, this type of model may be very useful to quantify the effects of local topographic (orientation, shading) and to take into account incoming solar radiation changes at geological timescale. The performed sensitivity test indicates that, in spite of the uncertainty in the calibrated snow-ice ablation factors, all models are able to provide paleotemperatures with ~1 °C uncertainty for a given paleoprecipitation. This error includes a 50 m uncertainty in the estimate of the paleoELA. Finally, the models are characterized by different precipitation-temperature sensitivities: if a similar warming is applied, model including different ablation factors for snow and ice requires a lower precipitation increase (by ∼15 %) than others to maintain the ELA.

Description: Trends in climatic variables and future reference evapotranspiration in Duero Valley (Spain) Natural Hazards and Earth System Science, 11, 1795-1805, 2011 Author(s): R. Moratiel, R. L. Snyder, J. M. Durán, and A. M. Tarquis The impact of climate change and its relation with evapotranspiration was evaluated in the Duero River Basin (Spain). The study shows possible future situations 50 yr from now from the reference evapotranspiration (ET o ). The maximum temperature ( T max ), minimum temperature ( T min ), dew point ( T d ), wind speed ( U ) and net radiation ( R n ) trends during the 1980–2009 period were obtained and extrapolated with the FAO-56 Penman-Montheith equation to estimate ET o . Changes in stomatal resistance in response to increases in CO 2 were also considered. Four scenarios were done, taking the concentration of CO 2 and the period analyzed (annual or monthly) into consideration. The scenarios studied showed the changes in ET o as a consequence of the annual and monthly trends in the variables T max , T min , T d , U and R n with current and future CO 2 concentrations (372 ppm and 550 ppm). The future ET o showed increases between 118 mm (11 %) and 55 mm (5 %) with respect to the current situation of the river basin at 1042 mm. The months most affected by climate change are May, June, July, August and September, which also coincide with the maximum water needs of the basin's crops.

Description: Possible interrelation between the lead time of precursory seismic electric signals (SES) and geodynamics in Aegean Sea Natural Hazards and Earth System Science, 11, 1599-1603, 2011 Author(s): E. Dologlou The seismicity of the last 15 years in the Aegean Sea revealed that earthquakes ( M w 〉 5) with epicentres falling within the Sporades basin and the confined area north of Samos island were preceded by electric seismic signals (SES) with a remarkably long lead time. A possible explanation of this behaviour by means of specific tectonics and geodynamics which characterise these two regions, such as a significant small crustal thickness and a high heat flow rate, has been attempted. New data seem to strengthen the above hypothesis.

Description: Deglaciation records of 17 O- excess in East Antarctica: reliable reconstruction of oceanic relative humidity from coastal sites Climate of the Past Discussions, 7, 1845-1886, 2011 Author(s): R. Winkler, A. Landais, H. Sodemann, L. Dümbgen, F. Prié, V. Masson-Delmotte, B. Stenni, and J. Jouzel We measured δ 17 O and δ 18 O in two Antarctic ice cores at EPICA Dome C (EDC) and TALDICE (TD), respectively and computed 17 O- excess with respect to VSMOW. The comparison of our 17 O- excess data with the previous record obtained at Vostok (Landais et al., 2008) revealed differences up to 35 ppm in 17 O- excess mean level and evolution for the three sites. Our data showed that the large increase depicted at Vostok (20 ppm) during the last deglaciation, is a regional and not a general pattern in the temporal distribution of 17 O- excess in East Antarctica. The EDC data display an increase of 13 ppm, whereas the TD data show no significant variation from the Last Glacial Maximum (LGM) to the Early Holocene (EH). Lagrangian moisture source diagnostic revealed very different source regions for Vostok and EDC compared to TD. These findings combined with the results of a sensitivity analysis, using a Rayleigh-type isotopic model, suggest that relative humidity (RH) at the oceanic source region (OSR) are a determining factor for the spatial differences of 17 O- excess in East Antarctica. However, 17 O- excess in remote sites of continental Antarctica (e.g. Vostok) may be highly sensitive to local effects. Hence, we consider 17 O- excess in coastal East Antarctic ice cores (TD) to be more reliable as a proxy for RH at the OSR.

Description: Impact of oceanic processes on the carbon cycle during the last termination Climate of the Past Discussions, 7, 1887-1934, 2011 Author(s): N. Bouttes, D. Paillard, D. M. Roche, C. Waelbroeck, M. Kageyama, A. Lourantou, E. Michel, and L. Bopp During the last termination (from ~18 000 yr ago to ~9000 yr ago) the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO 2 increased from ~190 ppm to ~260 ppm. Although this CO 2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a climate-carbon coupled model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependant diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO 2 , and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ 13 C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 yr ago, when the Antarctic ice sheet extent was at its maximum. Sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario it is possible to simulate both the amplitude and timing of the CO 2 increase during the last termination in agreement with data. The atmospheric δ 13 C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

Description: Productivity response of calcareous nannoplankton in the South Atlantic to the Eocene Thermal Maximum 2 (ETM2) Climate of the Past Discussions, 7, 2089-2118, 2011 Author(s): M. Dedert, H. M. Stoll, D. Kroon, N. Shimizu, and P. Ziveri The Early Eocene Thermal Maximum 2 (ETM2) at ~53.7 Ma is one of multiple hyperthermal events that followed the Paleocene-Eocene Thermal Maximum (PETM, ~55 Ma). In order to reconstruct the primary productivity response to the ETM2 in the South Atlantic, we have analyzed Sr/Ca ratios in various size fractions of bulk sediments and in picked monogeneric populations of calcareous nannofossils. The latter technique of measuring selected nannofossil populations using the ion probe circumvents possible contamination with secondary calcite. Avoiding such contamination is important for interpretation of the nannoplankton productivity record, since diagenetic processes can bias the productivity signal, as we demonstrate for Sr/Ca measurements in the fine (

Description: Daytime variations of foE connected to earthquakes Natural Hazards and Earth System Science, 11, 1807-1812, 2011 Author(s): E. V. Liperovskaya, V. V. Bogdanov, P.-F. Biagi, C.-V. Meister, and V. A. Liperovsky In the present work it is shown that, in accordance with the observations of the vertical sounding station "Tashkent", the critical foE-frequency of the daytime E-layer increases about one day before winter-earthquakes with magnitudes M 〉 5 and depths of the epicentre of h 〈 60 km, which appeared at distances of R 〈 2000 km from the station. The reliability of the result is larger than 99 %. The phenomenon is not observed for summer-earthquakes. It seems to be determined by the atmospheric wind system. Further, the variations of the foE-frequency are compared with possible simultaneous variations of the critical frequency foF2 of the F2-layer. First results show that only very large changes of the ionisation density in the E-layer influence the ionisation density in the F-region. Therefore, no synchronous growth of the foE- and foF2-frequencies 1–2 days before seismic shocks could be observed.

Description: Reconstruction of southeast Tibetan Plateau summer cloud cover over the past two centuries using tree ring δ 18 O Climate of the Past Discussions, 7, 1825-1844, 2011 Author(s): C. Shi, V. Daux, C. Risi, S.-G. Hou, M. Stievenard, M. Pierre, Z. Li, and V. Masson-Delmotte A tree-ring δ 18 O chronology of Linzhi spruce, spanning from AD 1781 to 2005, was developed in Bomi, Southeast Tibetan Plateau (TP). During the period with instrumental data (1961–2005), this record is strongly correlated with regional summer cloud cover, which is supported by a precipitation δ 18 O simulation conducted with the isotope-enabled atmospheric general circulation model LMDZiso. A 225-yr regional cloud cover reconstruction was therefore achieved. The observed cloud cover increased in the 1980s and this increase is not unprecedented in the entire reconstruction. The reconstructed cloud cover appears smaller and more stable in the 20th century than previously. A late 19th century decrease in our reconstructed cloud cover is consistent with a decrease in the TP glacier accumulation recorded in ice cores. Our data reveal a strong anomaly in the 1810s, which coincides with volcanic eruption in 1809 and the 1815 Tambora volcanic eruption.

Description: Shallow circulation groundwater – the main type of water containing hazardous radon concentration Natural Hazards and Earth System Science, 11, 1695-1703, 2011 Author(s): T. A. Przylibski The main factors affecting the value of 222 Rn activity concentration in groundwater are the emanation coefficient of reservoir rocks ( K em ), the content of parent 226 Ra in these rocks ( q ), changes in the volume and flow velocity as well as the mixing of various groundwater components in the circulation system. The highest values of 222 Rn activity concentration are recorded in groundwaters flowing towards an intake through strongly cracked reservoir rocks undergoing weathering processes. Because of these facts, waters with hazardous radon concentration levels, i.e. containing more than 100 Bq dm −3 222 Rn, could be characterised in the way that follows. They are classified as radon waters, high-radon waters and extreme-radon waters. They belong to shallow circulation systems (at less than a few dozen metres below ground level) and are contemporary infiltration waters, i.e. their underground flow time ranges from several fortnights to a few decades. Because of this, these are usually poorly mineralised waters (often below 0.2–0.5 g dm −3 ). Their resources are renewable, but also vulnerable to contamination. Waters of this type are usually drawn from private intakes, supplying water to one or at most a few households. Due to an increased risk of developing lung tumours, radon should be removed from such waters when still in the intake. To achieve this aim, appropriate legislation should be introduced in many countries.

Description: Precipitation change in Southern Italy linked to global scale oscillation indexes Natural Hazards and Earth System Science, 11, 1683-1694, 2011 Author(s): T. Caloiero, R. Coscarelli, E. Ferrari, and M. Mancini This study investigates precipitation variability in five regions of Southern Italy (Campania, Apulia, Basilicata, Calabria and Sicily) using a homogeneous database of about 70 rain gauges with more than 50 years of observation. First, a statistical analysis was performed through the Mann-Kendall non-parametric test in order to determine rainfall the trend on both yearly and seasonal scales. Then, the relationship between the rainfall and some teleconnection pattern indexes was investigated using Spearman's test. The results show remarkable statistically significant negative trends for annual and winter aggregations in most part of the series. Moreover, a strong correlation has emerged between the teleconnection patterns and precipitation in Southern Italy, particularly in winter and on the Tyrrhenian side of the study area.

Description: Reconstruction of Atlantic historical winter coastal storms in the Spanish coasts of the Gulf of Cadiz, 1929–2005 Natural Hazards and Earth System Science, 11, 1715-1722, 2011 Author(s): P. Ribera, D. Gallego, C. Pena-Ortiz, L. Del Rio, T. A. Plomaritis, and J. Benavente This paper presents the reconstruction of a climatological series of winter coastal storms on the northern coasts of the Gulf of Cadiz. This series has been put together using information extracted from regional and local Spanish newspapers. It includes all the storms coming from the Atlantic sector that have been detected during the winter season, from October to March, between 1929 and 2005. In order to validate this historical storm series, it has been compared with storms series identified from quasi-observational data and using different wave heights as thresholds to decide what is to be considered as a coastal storm. Nearly 2.6 reports per year about coastal storms are published in the press which correspond to waves of 3.6 m high or more and to prevailing winds from a direction ranging between SSW and WNW. A long- term positive trend has been detected for the complete storm series. If only the instrumental period is analysed, no significant trend is detected. It is suggested that this difference might be associated with the impact of the North Atlantic Oscillation over the occurrence of storms in this area.

Description: The neutral temperature in the ionospheric dynamo region and the ionospheric F region density during Wenchuan and Pingtung Doublet earthquakes Natural Hazards and Earth System Science, 11, 1759-1768, 2011 Author(s): Y. Y. Sun, K.-I. Oyama, J. Y. Liu, H. K. Jhuang, and C. Z. Cheng One of the possible candidates which modifies the ionosphere before large earthquake is electric field. We presume that the electric field associated with large earthquakes is generated in the ionosphere dynamo region (100–120 km). This paper tries to identify the evidence of the contribution of the neutral atmosphere in the dynamo region. The relationship between the critical frequency at the F2 peak (foF2) and the height profile of the neutral atmosphere temperature was studied for two large earthquakes: Wenchuan, 2008 and Pingtung Doublet, 2006. It is found that the wave amplitude of the vertical wavelength 20–30 km which is usually superposed on the height profile of the neutral atmosphere temperature enhances when the foF2 increases. The correlation between the wave amplitude and foF2 is found better along a longitudinal direction than along latitude direction.

Description: The relationship between the performance of soil conditions and damage following an earthquake: a case study in Istanbul, Turkey Natural Hazards and Earth System Science, 11, 1745-1758, 2011 Author(s): G. Ç. İnce Istanbul city has experienced many strong earthquakes throughout its history and suffered extensive damage. The old Istanbul (Fatih and Eminonu districts), one of the most densely populated locations of Istanbul city, is a commercial centre and has many significant historical buildings. In this study, the data pertaining to the damage sustained by historical artifacts and structures from past earthquakes, are examined along with the soil amplification of the region, the liquefaction and slope stability risk in relation to the seismic microzonation maps which were prepared using geographic information system techniques. The relationship between soil behaviour and the damage resulting from previous earthquakes. The structural damage observed in the region was in accordance with the microzonation maps. Since the area does not have a high risk in terms of slope stability, this does not have much impact on the level of damage. However, it was found that soil amplification and partial liquefaction contributed to the damage to historical artifacts and structures.

Description: Evaluation of future hydrological cycle under climate change scenarios in a mesoscale Alpine watershed of Italy Natural Hazards and Earth System Science, 11, 1769-1785, 2011 Author(s): B. Groppelli, A. Soncini, D. Bocchiola, and R. Rosso We investigate future (2045–2054) hydrological cycle of the snow fed Oglio (≈1800 km 2 ) Alpine watershed in Northern Italy. A Stochastic Space Random Cascade (SSRC) approach is used to downscale future precipitation from three general circulation models, GCMs (PCM, CCSM3, and HadCM3) available within the IPCC's data base and chosen for this purpose based upon previous studies. We then downscale temperature output from the GCMs to obtain temperature fields for the area. We also consider a projected scenario based upon trends locally observed in former studies, LOC scenario. Then, we feed the downscaled fields to a minimal hydrological model to build future hydrological scenarios. We provide projected flow duration curves and selected flow descriptors, giving indication of expected modified (against control run for 1990–1999) regime of low flows and droughts and flood hazard, and thus evaluate modified peak floods regime through indexed flood. We then assess the degree of uncertainty, or spread, of the projected water resources scenarios by feeding the hydrological model with ensembles projections consistent with our deterministic (GCMs + LOC) scenarios, and we evaluate the significance of the projected flow variables against those observed in the control run. The climate scenarios from the adopted GCMs differ greatly from one another with respect to projected precipitation amount and temperature regimes, and so do the projected hydrological scenarios. A relatively good agreement is found upon prospective shrinkage and shorter duration of the seasonal snow cover due to increased temperature patterns, and upon prospective increase of hydrological losses, i.e. evapotranspiration, for the same reason. However, precipitation patterns are less consistent, because HadCM3 and PCM models project noticeably increased precipitation for 2045–2054, whereas CCSM3 provides decreased precipitation patterns therein. The LOC scenario instead displays unchanged precipitation. The ensemble simulations indicate that several projected flow variables under the considered scenarios are significantly different from their control run counterparts, and also that snow cover seems to significantly decrease in duration and depth. The proposed hydrological scenarios eventually provide a what-if analysis, giving a broad view of the possible expected impacts of climate change within the Italian Alps, necessary to trigger the discussion about future adaptation strategies.

Description: Volcanic and ENSO effects in China in simulations and reconstructions: Tambora eruption 1815 Climate of the Past Discussions, 7, 2061-2088, 2011 Author(s): D. Zhang, R. Blender, and K. Fraedrich The co-operative effects of volcanic eruptions and ENSO (El Niño/Southern Oscillation) on the climate in China are analyzed in a millennium simulation for 800–2005 AD using the earth system model (ESM) ECHAM5/MPIOM/JSBACH subject to anthropogenic and natural forcings. The experiment includes two ensembles with weak (5 members) and strong (3 members) total solar irradiance variability. In the absence of El Niño and La Niña events, volcanoes, which are the dominant forcing in both ensembles, cause a dramatic cooling in West China (−2 °C) and a drought in East China during the year after the eruption. The recovery times for the volcano induced cooling vary globally between one and 12 yr; in China these values are mostly within 1–4 yr, but reach 10 yr in the Northeast. Without volcanoes, after El Niño events the summer precipitation is reduced in the North, while South China becomes wetter (indicated by the Standardized Precipitation Index, SPI, for summers, JJA); La Niña events cause opposite effects. El Niño events in the winters after eruptions compensate the cooling in most regions of China, while La Niña events intensify the cooling (up to −2.5 °C). The simulated impact of the eruption of the Tambora in 1815, which caused the "year without summer" 1816 in Europe and North America and coldness and famines for several years in the Chinese province Yunnan, depends crucially on the ENSO state of the coupled model. A comparison with reconstructed El Niño events shows a moderate cool climate with wet (in the South) and extreme dry anomalies (in the North) persisting for several years.

Description: A tool for assessing the quality of the Mediterranean cyclone forecast: a numerical index Natural Hazards and Earth System Science, 11, 1787-1794, 2011 Author(s): M. A. Picornell, A. Jansà, and A. Genovés Cyclones affecting the Mediterranean region, sometimes related to severe weather events, are often not well represented enough in numerical model predictions. Assessing the quality of the forecast of these cyclonic structures would be a significant advance in better knowing the goodness of the weather forecast in this region, and particularly the quality of predictions of high impact phenomena. In order to estimate the cyclone forecast uncertainty in operational models, in this work we compare two cyclone databases for the period 2006–2007: one from the operational analyses of the T799 ECMWF deterministic model; and the other from the forecasts provided by the same model in three ranges, H+12, H+24, and H+48. The skill of the model to detect cyclones and its accuracy in describing their features are assessed. An index is presented as an indicator of the quality of the prediction, derived from the frequency distribution of errors in the prediction of four characteristics of the cyclone: position, central pressure value, geostrophic circulation, and domain. Some sub-indexes are derived to verify each of the variables separately in order to analyse the most frequent sources of error. Other sub-indexes are also defined to indicate possible biases in the numerical prediction model.

Description: Occurrence analysis of daily rainfalls through non-homogeneous Poissonian processes Natural Hazards and Earth System Science, 11, 1657-1668, 2011 Author(s): B. Sirangelo, E. Ferrari, and D. L. De Luca A stochastic model based on a non-homogeneous Poisson process, characterised by a time-dependent intensity of rainfall occurrence, is employed to explain seasonal effects of daily rainfalls exceeding prefixed threshold values. The data modelling has been performed with a partition of observed daily rainfall data into a calibration period for parameter estimation and a validation period for checking on occurrence process changes. The model has been applied to a set of rain gauges located in different geographical areas of Southern Italy. The results show a good fit for time-varying intensity of rainfall occurrence process by 2-harmonic Fourier law and no statistically significant evidence of changes in the validation period for different threshold values.

Description: Regional frequency analysis of extreme storm surges along the French coast Natural Hazards and Earth System Science, 11, 1627-1639, 2011 Author(s): L. Bardet, C.-M. Duluc, V. Rebour, and J. L'Her A good knowledge of extreme storm surges is necessary to ensure protection against flood. In this paper we introduce a methodology to determine time series of skew surges in France as well as a statistical approach for estimating extreme storm surges. With the aim to cope with the outlier issue in surge series, a regional frequency analysis has been carried out for the surges along the Atlantic coast and the Channel coast. This methodology is not the current approach used to estimate extreme surges in France. First results showed that the extreme events identified as outliers in at-site analyses do not appear to be outliers any more in the regional empirical distribution. Indeed the regional distribution presents a curve to the top with these extreme events that a mixed exponential distribution seems to recreate. Thus, the regional approach appears to be more reliable for some sites than at-site analyses. A fast comparison at a given site showed surge estimates with the regional approach and a mixed exponential distribution are higher than surge estimates with an at-site fitting. In the case of Brest, the 1000-yr return surge is 167 cm in height with the regional approach instead of 126 cm with an at-site analysis.

Description: The Middle Miocene climate as modelled in an atmosphere-ocean-biosphere model Climate of the Past Discussions, 7, 1935-1972, 2011 Author(s): M. Krapp and J. H. Jungclaus We present simulations with a coupled ocean-atmosphere-biosphere model for the Middle Miocene 15 million years ago. The Middle Miocene topography, which alters both large-scale ocean and atmospheric circulations, causes a global warming of 0.7 K compared to present-day. Higher than present-day CO 2 levels of 480 and 720 ppm cause a global warming of 2.8 and 4.9 K, thereby matching proxy-based Middle Miocene global temperature estimates of 3–6 K warming. Higher CO 2 levels and the associated water vapour feedback enhance the greenhouse effect and lead to a polar amplification of the warming. Although oceanic and atmospheric poleward heat transport are individually altered by 10–30 % in the mid and high latitudes, changes of the total heat transport account only for 4–8 %, pointing toward a compensation between oceanic and atmospheric heat transport. Our model reproduces a denser vegetation in agreement with fossil records. These results suggest that higher than present-day CO 2 levels are essential to drive the warm Middle Miocene climate.

Description: Synchronicity of the East Asian Summer Monsoon variability and Northern Hemisphere climate change since the last deglaciation Climate of the Past Discussions, 7, 2159-2192, 2011 Author(s): T. Shinozaki, M. Uchida, K. Minoura, M. Kondo, S. F. Rella, and Y. Shibata Understanding of the mechanism of the East Asian Summer Monsoon (EASM) is required for the prediction of climate change in East Asia in a scenario of modern global warming. In this study, we present high-resolution climate records from peat sediments in Northeast Japan to reconstruct the EASM variability based on peat bulk cellulose δ 13 C since the last deglaciation. We used a 8.8 m long peat sediment core collected from the Tashiro Bog, Northeast Japan. Based on 42 14 C measurements, the core bottom reaches ~15.5 ka. δ 13 C, accumulation rate and accumulation flux time-series correlate well to Greenland ice core δ 18 O variability, suggesting that the climate record in Northeast Japan is linked to global climate changes. The δ 13 C record at Tashiro Bog and other paleo-EASM records at Northeast and Southern China consistently demonstrate that hydrological environments were spatially different in mid-high and mid-low latitude regions over the last 15.5 kyr. During global cooling (warming) periods, mid-high and mid-low latitude regions were characterized by wet (dry) and dry (wet) environments, respectively. We suggest that these climatic patterns are related to the migration of the EASM-related rain belt during global climate changes, as a consequence of variations in intensity and location of both the Intertropical Convergence Zone (ITCZ) and the Western Pacific Subtropical High (STH). The location of the rain belt largely influences the East Asian hydrological environment. Our δ 13 C time-series are characterized by a 1230 yr throughout the Holocene and a 680 yr periodicity during the early Holocene. The 1230 yr periodicity is in agreement with North Atlantic ice-rafted debris (IRD) events, suggesting a teleconnection between the Northeast Japan and the North Atlantic during the Holocene. In addition, it is the first evidence that the Bond events were recorded in terrestrial sediment in Japan. On the other hand, the 680 yr periodicity between 10.0 and 8.0 kyr is consistent with a prominent 649 yr solar activity cycle, suggesting that solar activity affected EASM precipitation during the Hypsithermal, when orbital-scale solar insolation was at a maximum in the Northern Hemisphere.

Description: Probabilistic assessments of climate change impacts on durum wheat in the Mediterranean region Natural Hazards and Earth System Science, 11, 1293-1302, 2011 Author(s): R. Ferrise, M. Moriondo, and M. Bindi Recently, the availability of multi-model ensemble prediction methods has permitted a shift from a scenario-based approach to a risk-based approach in assessing the effects of climate change. This provides more useful information to decision-makers who need probability estimates to assess the seriousness of the projected impacts. In this study, a probabilistic framework for evaluating the risk of durum wheat yield shortfall over the Mediterranean Basin has been exploited. An artificial neural network, trained to emulate the outputs of a process-based crop growth model, has been adopted to create yield response surfaces which are then overlaid with probabilistic projections of future temperature and precipitation changes in order to estimate probabilistic projections of future yields. The risk is calculated as the relative frequency of projected yields below a selected threshold. In contrast to previous studies, which suggest that the beneficial effects of elevated atmospheric CO 2 concentration over the next few decades would outweigh the detrimental effects of the early stages of climatic warming and drying, the results of this study are of greater concern.

Description: Debris-flow activity in abandoned channels of the Manival torrent reconstructed with LiDAR and tree-ring data Natural Hazards and Earth System Science, 11, 1247-1257, 2011 Author(s): J. Lopez Saez, C. Corona, M. Stoffel, A. Gotteland, F. Berger, and F. Liébault Hydrogeomorphic processes are a major threat in many parts of the Alps, where they periodically damage infrastructure, disrupt transportation corridors or even cause loss of life. Nonetheless, past torrential activity and the analysis of areas affected during particular events remain often imprecise. It was therefore the purpose of this study to reconstruct spatio-temporal patterns of past debris-flow activity in abandoned channels on the forested cone of the Manival torrent (Massif de la Chartreuse, French Prealps). A Light Detecting and Ranging (LiDAR) generated Digital Elevation Model (DEM) was used to identify five abandoned channels and related depositional forms (lobes, lateral levees) in the proximal alluvial fan of the torrent. A total of 156 Scots pine trees ( Pinus sylvestris L.) with clear signs of debris flow events was analyzed and growth disturbances (GD) assessed, such as callus tissue, the onset of compression wood or abrupt growth suppression. In total, 375 GD were identified in the tree-ring samples, pointing to 13 debris-flow events for the period 1931–2008. While debris flows appear to be very common at Manival, they have only rarely propagated outside the main channel over the past 80 years. Furthermore, analysis of the spatial distribution of disturbed trees contributed to the identification of four patterns of debris-flow routing and led to the determination of three preferential breakout locations. Finally, the results of this study demonstrate that the temporal distribution of debris flows did not exhibit significant variations since the beginning of the 20th century.

Description: Assessment of possible damaged areas due to landslide-induced waves at a constructed reservoir using empirical approaches: Kurtun (North Turkey) Dam reservoir area Natural Hazards and Earth System Science, 11, 1341-1350, 2011 Author(s): A. Akgün Landslide-induced wave in lakes or offshore is a natural hazard of significant concern throughout the world. In Turkey, several dam reservoir areas suffer from this problem, motivating for the exploration of this phenomenon. In this study, the potential for landslide-induced wave generation and the possible adverse effects of such an event were investigated for the Kurtun Dam reservoir area. This area is prone to translational debris-type landslides, and such slides represent a threat with respect to wave generation. To investigate this problem, a potential landslide was initially inspected with respect to key geometrical and index parameters, such as the internal friction angle and unit weight of the landslide material. After obtaining these data, the potential of sub-aerial landslide-induced wave characteristics such as wave height, wave run-up on the opposite slope and wave velocity were calculated using existing empirical relationships. Based on the obtained wave properties, a potential damage assessment was performed for vulnerable areas in recognition of the fact that the Kurtun Dam and Kurtun district are threatened by potential wave occurrence. According to the findings obtained from the potential damage assessment, it was determined that the Kurtun district and the Kurtun portion of the Kurtun-Gumushane highway are located within the run-up impact area. However, the Kurtun Dam was determined to be relatively safe due to the distance between the landslide area and the dam.

Description: The sensitivity of warm period precipitation forecasts to various modifications of the Kain-Fritsch Convective Parameterization scheme Natural Hazards and Earth System Science, 11, 1327-1339, 2011 Author(s): N. Mazarakis, V. Kotroni, K. Lagouvardos, A. A. Argiriou, and C. J. Anderson The sensitivity of quantitative precipitation forecasts to various modifications of the Kain-Fritsch (KF) convective parameterization scheme (CPS) is examined for twenty selected cases characterized by intense convective activity and widespread precipitation over Greece, during the warm period of 2005–2007. The study is conducted using the MM5 model with a two nested domains strategy, with horizontal grid increments of 24 and 8 km, respectively. Five modifications to the KF CPS, each designed to test the sensitivity of the model to the convective scheme formulation, are discussed. The modifications include: (i) the maximization of the convective scheme precipitation efficiency, (ii) the change of the convective time step, (iii) the forcing of the convective scheme to produce more/less cloud material, (iv) changes to the trigger function and (v) the alteration of the vertical profile of updraft mass flux detrainment. The simulated precipitation from the 8-km grid is verified against raingauge measurements. Although skill scores vary widely among the cases and the precipitation thresholds, model results using the modifications of the convective scheme show improvements in 6-h precipitation totals compared to simulations generated using the unmodified convective scheme. In general, forcing the model to produce less cloud material improves the precipitation forecast for the moderate and high precipitation amounts, while the same modification and the change of the convective time step to 1 min has the same result for the high precipitation thresholds. The increase of convective time step to 15 min, the maximization of precipitation efficiency and the changes to the trigger function give similar results for medium and high precipitation. On the other hand, the forecast for the light precipitation is improved by forcing the model to produce more cloud material as well as by the alteration of the vertical profile of updraft mass flux detrainment.

Description: Future changes in European winter storm losses and extreme wind speeds inferred from GCM and RCM multi-model simulations Natural Hazards and Earth System Science, 11, 1351-1370, 2011 Author(s): M. G. Donat, G. C. Leckebusch, S. Wild, and U. Ulbrich Extreme wind speeds and related storm loss potential in Europe have been investigated using multi-model simulations from global (GCM) and regional (RCM) climate models. Potential future changes due to anthropogenic climate change have been analysed from these simulations following the IPCC SRES A1B scenario. The large number of available simulations allows an estimation of the robustness of detected future changes. All the climate models reproduced the observed spatial patterns of wind speeds, although some models displayed systematic biases. A storm loss model was applied to the GCM and RCM simulated wind speeds, resulting in realistic mean loss amounts calculated from 20th century climate simulations, although the inter-annual variability of losses is generally underestimated. In future climate simulations, enhanced extreme wind speeds were found over northern parts of Central and Western Europe in most simulations and in the ensemble mean (up to 5%). As a consequence, the loss potential is also higher in these regions, particularly in Central Europe. Conversely, a decrease in extreme wind speeds was found in Southern Europe, as was an associated reduction in loss potential. There was considerable spread in the projected changes of individual ensemble members, with some indicating an opposite signature to the ensemble mean. Downscaling of the large-scale simulations with RCMs has been shown to be an important source of uncertainty. Even RCMs with identical boundary forcings can show a wide range of potential changes. The robustness of the projected changes was estimated using two different measures. First, the inter-model standard deviation was calculated; however, it is sensitive to outliers and thus displayed large uncertainty ranges. Second, a multi-model combinatorics approach considered all possible sub-ensembles from GCMs and RCMs, hence taking into account the arbitrariness of model selection for multi-model studies. Based on all available GCM and RCM simulations, for example, a 25% mean increase in risk of loss for Germany has been estimated for the end of the 21st century, with a 90% confidence range of +15 to +35%.

Description: Evidence of a previously unrecorded local tsunami, 13 April 2010, Cook Islands: implications for Pacific Island countries Natural Hazards and Earth System Science, 11, 1371-1379, 2011 Author(s): J. Goff Tsunami hazard assessments for Pacific Islands Countries (PICs) tend to focus on subduction zone sources. It is generally recognised that while volcanic-related tsunamigenic sources exist, they are probably only of minor relevance to the overall hazardscape of the Pacific. This paper outlines the evidence for a previously unrecorded local tsunami that struck the uninhabited south coast of Mangaia, Cook Islands, on 13 April 2010. The tsunami had a maximum inundation of 100 m inland and a runup of 12 m a.s.l. This event was most probably caused by a small submarine slope failure, the most recent of an unknown number of previous inundations. Since most PICs have a volcanic origin, it is suggested that current perceptions about the local and regional significance of such events is inaccurate. A review of volcanic-related tsunamigenic sources throughout the Pacific reveals a wealth of data concerning submarine slope failures in particular and a more general background of active volcanism. These sources are as relevant to PICs close to or far away from subduction zones. As populations grow and the coastlines of many PICs and those on the edge of the Pacific Ocean become increasing occupied, the likelihood for loss of life from these events increases.

Description: Active faults crossing trunk pipeline routes: some important steps to avoid disaster Natural Hazards and Earth System Science, 11, 1433-1436, 2011 Author(s): V. M. Besstrashnov and A. L. Strom Assessment of seismic strong motion hazard produced by earthquakes originating within causative fault zones allows rather low accuracy of localisation of these structures that can be provided by indirect evidence of fault activity. In contrast, the relevant accuracy of localisation and characterisation of active faults, capable of surface rupturing, can be achieved solely by the use of direct evidence of fault activity. This differentiation requires strict definition of what can be classified as "active fault" and the normalisation of methods used for identification and localisation of active faults crossing oil and natural gas trunk pipelines.

Description: Regional-scale analysis of lake outburst hazards in the southwestern Pamir, Tajikistan, based on remote sensing and GIS Natural Hazards and Earth System Science, 11, 1447-1462, 2011 Author(s): M. Mergili and J. F. Schneider This paper presents an analysis of the hazards emanating from the sudden drainage of alpine lakes in South-Western Tajik Pamir. In the last 40 yr, several new lakes have formed in the front of retreating glacier tongues, and existing lakes have grown. Other lakes are dammed by landslide deposits or older moraines. In 2002, sudden drainage of a glacial lake in the area triggered a catastrophic debris flow. Building on existing approaches, a rating scheme was devised allowing quick, regional-scale identification of potentially hazardous lakes and possible impact areas. This approach relies on GIS, remote sensing and empirical modelling, largely based on medium-resolution international datasets. Out of the 428 lakes mapped in the area, 6 were rated very hazardous and 34 hazardous. This classification was used for the selection of lakes requiring in-depth investigation. Selected cases are presented and discussed in order to understand the potentials and limitations of the approach used. Such an understanding is essential for the appropriate application of the methodology for risk mitigation purposes.

Description: Sensitivity of interglacial Greenland temperature and δ 18 O to orbital and CO 2 forcing: climate simulations and ice core data Climate of the Past Discussions, 7, 1585-1630, 2011 Author(s): V. Masson-Delmotte, P. Braconnot, G. Hoffmann, J. Jouzel, M. Kageyama, A. Landais, Q. Lejeune, C. Risi, L. Sime, J. Sjolte, D. Swingedouw, and B. Vinther The sensitivity of interglacial Greenland temperature to orbital and CO 2 forcing is investigated using the NorthGRIP ice core data and coupled ocean-atmosphere IPSL-CM4 model simulations. These simulations were conducted in response to different interglacial orbital configurations, and to increased CO 2 concentrations. These different forcings cause very distinct simulated seasonal and latitudinal temperature and water cycle changes, limiting the analogies between the last interglacial and future climate. However, the IPSL-CM4 model shows similar magnitudes of Arctic summer warming and climate feedbacks in response to 2 × CO 2 and orbital forcing of the last interglacial period (126 000 yr ago). The IPSL model produces a remarkably linear relationship between top of atmosphere incoming summer solar radiation and simulated changes in summer and annual mean central Greenland temperature. This contrasts with the stable isotope record from the Greenland ice cores, showing a multi-millennial lagged response to summer insolation. During the early part of interglacials, the observed lags may be explained by ice sheet-ocean feedbacks linked with changes in ice sheet elevation and the impact of meltwater on ocean circulation, as investigated with sensitivity studies. A quantitative comparison between ice core data and climate simulations requires to explore the stability of the stable isotope – temperature relationship. Atmospheric simulations including water stable isotopes have been conducted with the LMDZiso model under different boundary conditions. This set of simulations allows to calculate a temporal Greenland isotope-temperature slope (0.3–0.4 ‰ per °C) during warmer than present Arctic climates, in response to increased CO 2 , increased ocean temperature and orbital forcing. This temporal slope appears twice as small as the modern spatial gradient and is consistent with other ice core estimates. A preliminary comparison with other model results implies that other mechanisms could also play a role. This suggests that further simulations and detailed inter-model comparisons are also likely to be of benefit. Comparisons with Greenland ice core stable isotope data reveals that IPSL/LMDZiso simulations strongly underestimate the amplitude of the ice core signal during the last interglacial, which could reach +8–10 °C at fixed-elevation. While the model-data mismatch may result from missing positive feedbacks (e.g. vegetation), it could also be explained by a reduced elevation of the central Greenland ice sheet surface by 300–400 m.

Description: Perturbation of convection-permitting NWP forecasts for flash-flood ensemble forecasting Natural Hazards and Earth System Science, 11, 1529-1544, 2011 Author(s): B. Vincendon, V. Ducrocq, O. Nuissier, and B. Vié Mediterranean intense weather events often lead to devastating flash-floods. Extending the forecasting lead times further than the watershed response times, implies the use of numerical weather prediction (NWP) to drive hydrological models. However, the nature of the precipitating events and the temporal and spatial scales of the watershed response make them difficult to forecast, even using a high-resolution convection-permitting NWP deterministic forecasting. This study proposes a new method to sample the uncertainties of high-resolution NWP precipitation forecasts in order to quantify the predictability of the streamflow forecasts. We have developed a perturbation method based on convection-permitting NWP-model error statistics. It produces short-term precipitation ensemble forecasts from single-value meteorological forecasts. These rainfall ensemble forecasts are then fed into a hydrological model dedicated to flash-flood forecasting to produce ensemble streamflow forecasts. The verification on two flash-flood events shows that this forecasting ensemble performs better than the deterministic forecast. The performance of the precipitation perturbation method has also been found to be broadly as good as that obtained using a state-of-the-art research convection-permitting NWP ensemble, while requiring less computing time.

Description: Geomorphological and geochemical characterization of the 11 August 2008 mud volcano eruption at S. Barbara village (Sicily, Italy) and its possible relationship with seismic activity Natural Hazards and Earth System Science, 11, 1545-1557, 2011 Author(s): P. Madonia, F. Grassa, M. Cangemi, and C. Musumeci On 11 August 2008 a paroxysmal eruption occurred at Santa Barbara mud volcano (MV), located close to Caltanissetta, one of the most densely populated cities of Sicily (Italy). An associated minor event took place on August 2009. Both the events caused severe damage to civil infrastructures located within a range of about 2 km from the eruptive vent. Geomorphological, geochemical, and seismological investigations were carried out for framing the events in the appropriate geodynamic context. Geomorphological surveys recognized, in the immediate surrounding of the main emission point, two different families of processes and landforms: (i) ground deformations and (ii) changes in morphology and number of the fluid emitting vents. These processes were associated to a wider network of fractures, seemingly generated by the shock wave produced by the gas blast that occurred at the main paroxysm. Geochemical characterization allowed an estimation of the source of the fluids, or at least their last standing, at about 3 km depth. Finally, the close time relationships observed between anomalous increments of seismic activity and the two main paroxysmal events accounted for a possible common trigger for both the phenomena, even with different timing due to the very different initial conditions and characteristics of the two processes, i.e. seismogenesis and gas overloading.

Description: Moraine-dammed lake failures in Patagonia and assessment of outburst susceptibility in the Baker Basin Natural Hazards and Earth System Science, 14, 3243-3259, 2014 Author(s): P. Iribarren Anacona, K.P. Norton, and A. Mackintosh Glacier retreat since the Little Ice Age has resulted in the development or expansion of hundreds of glacial lakes in Patagonia. Some of these lakes have produced large (≥ 10 6 m 3 ) Glacial Lake Outburst Floods (GLOFs) damaging inhabited areas. GLOF hazard studies in Patagonia have been mainly based on the analysis of short-term series (≤ 50 years) of flood data and until now no attempt has been made to identify the relative susceptibility of lakes to failure. Power schemes and associated infrastructure are planned for Patagonian basins that have historically been affected by GLOFs, and we now require a thorough understanding of the characteristics of dangerous lakes in order to assist with hazard assessment and planning. In this paper, the conditioning factors of 16 outbursts from moraine-dammed lakes in Patagonia were analysed. These data were used to develop a classification scheme designed to assess outburst susceptibility, based on image classification techniques, flow routine algorithms and the Analytical Hierarchy Process. This scheme was applied to the Baker Basin, Chile, where at least seven moraine-dammed lakes have failed in historic time. We identified 386 moraine-dammed lakes in the Baker Basin of which 28 were classified with high or very high outburst susceptibility. Commonly, lakes with high outburst susceptibility are in contact with glaciers and have moderate (〉 8°) to steep (〉 15°) dam outlet slopes, akin to failed lakes in Patagonia. The proposed classification scheme is suitable for first-order GLOF hazard assessments in this region. However, rapidly changing glaciers in Patagonia make detailed analysis and monitoring of hazardous lakes and glaciated areas upstream from inhabited areas or critical infrastructure necessary, in order to better prepare for hazards emerging from an evolving cryosphere.

Description: Epic landslide erosion from mountain roads in Yunnan, China – challenges for sustainable development Natural Hazards and Earth System Science, 14, 3093-3104, 2014 Author(s): R. C. Sidle, M. Ghestem, and A. Stokes Expanding systems of mountain roads in developing countries have significantly increased the risk of landslides and sedimentation, and have created vulnerabilities for residents and aquatic resources. We measured landslide erosion along seven road segments in steep terrain in the upper Salween River basin, Yunnan, China and estimated sediment delivery to channels. Landslide erosion rates along the roads ranged from 2780 to 48 235 Mg ha −1 yr −1 , the upper end of this range being the highest rate ever reported along mountain roads. The two roads with the highest landslide erosion (FG1 = 12 966 Mg ha −1 yr −1 ; DXD = 48 235 Mg ha −1 yr −1 ) had some of the highest sediment delivery rates to channels (about 80 and 86%, respectively). Overall, 3 times more landslides occurred along cut slopes compared to fill slopes, but fill slope failures had a combined mass 〉 1.3 times that of cut slope failures. Many small landslides occurred along road cuts, but these were often trapped on the road surface. Given the magnitude of the landslide problem and the lack of attention to this issue, a more sustainable approach for mountain road development is outlined based on an analysis of landslide susceptibility and how thresholds for landslide trigger mechanisms would be modified by road location and different construction techniques.

Description: Evaluation of forest fire models on a large observation database Natural Hazards and Earth System Science, 14, 3077-3091, 2014 Author(s): J. B. Filippi, V. Mallet, and B. Nader This paper presents the evaluation of several fire propagation models using a large set of observed fires. The observation base is composed of 80 Mediterranean fire cases of different sizes, which come with the limited information available in an operational context (burned surface and approximative ignition point). Simulations for all cases are carried out with four different front velocity models. The results are compared with several error scoring methods applied to each of the 320 simulations. All tasks are performed in a fully automated manner, with simulations run as first guesses with no tuning for any of the models or cases. This approach leads to a wide range of simulation performance, including some of the bad simulation results to be expected in an operational context. Disregarding the quality of the input data, it is found that the models can be ranked based on their performance and that the most complex models outperform the more empirical ones. Data and source codes used for this paper are freely available to the community.

Description: The sharp decline of East Asian summer monsoon at mid-Holocene indicated by the lake-wetland transition in the Sanjiang Plain, northeastern China Climate of the Past Discussions, 10, 4595-4622, 2014 Author(s): Z. Q. Zhang, G. P. Wang, X. G. Lv, H. J. Jia, and Q. H. Xu The timing of the waxing and wining of the East Asian summer monsoon during the Holocene is still under debate. In present study, we present the high-resolution grain-size and LOI records from a well-dated mud/peat profile to reveal the lake-wetland transition in the Sanjiang Plain and discuss its significance to Holocene monsoon evolutions. The results show that the shallow-water lakes have developed in low-lying areas of the plain before 4600 yr BP, corresponding to the Holocene monsoon maximum. Thereafter, the wetlands began to initiate with the extinction of the paleolakes, marking a lake-shrinking stage with the relative dry climate. Considering the prevalent monsoon climate in the Sanjiang Plain, we suggest the lake-wetland transition at 4600 yr BP indicate a sharp decline of the summer monsoon rather than the basin infilling process. Such a remarkable monsoon weakening event has been widely documented in northern China, and we associated it with the ocean–atmosphere interacting processes in low-latitude regions.

Description: A two thousand year annual record of snow accumulation rates for Law Dome, East Antarctica Climate of the Past Discussions, 10, 4469-4497, 2014 Author(s): J. Roberts, C. Plummer, T. Vance, T. van Ommen, A. Moy, S. Poynter, A. Treverrow, M. Curran, and S. George Accurate high resolution records of snow accumulation rates in Antarctica are crucial for estimating ice sheet mass balance and subsequent sea level change. Snowfall rates at Law Dome, East Antarctica, have been linked with regional atmospheric circulation to mid-latitudes as well as regional Antarctic snowfall. Here, we extend the Law Dome accumulation record from 750 to 2035 years, using recent annual layer dating that extends to AD −22. Accumulation rates were calculated as the ratio of measured to modelled layer thicknesses, multiplied by the long term mean accumulation rate. The modelled layer thicknesses were based on a power law vertical strain rate profile fitted to observed annual layer thickness. The periods AD 380–442, AD 727–783 and AD 1970–2009 have above average snow accumulation rates, while AD 663–704, AD 933–975 and AD 1429–1468 were below average. The calculated snow accumulation rates show good correlation with atmospheric reanalysis estimates, and significant spatial correlation over a wide expanse of East Antarctica, demonstrating that the Law Dome record captures larger scale variability across a large region of East Antarctica well beyond the immediate vicinity of the Law Dome summit. Spectral analysis reveals periodicities in the snow accumulation record which may be related to ENSO and Interdecadal Pacific Oscillation frequencies.

Description: Technical Note: Probabilistically constraining proxy age–depth models within a Bayesian hierarchical reconstruction model Climate of the Past Discussions, 10, 4499-4533, 2014 Author(s): J. P. Werner and M. P. Tingley Reconstructions of late-Holocene climate rely heavily upon proxies that are assumed to be accurately dated by layer counting, such as measurement on tree rings, ice cores, and varved lake sediments. Considerable advances may be achievable if time uncertain proxies could be included within these multiproxy reconstructions, and if time uncertainties were recognized and correctly modeled for proxies commonly treated as free of age model errors. Current approaches to accounting for time uncertainty are generally limited to repeating the reconstruction using each of an ensemble of age models, thereby inflating the final estimated uncertainty – in effect, each possible age model is given equal weighting. Uncertainties can be reduced by exploiting the inferred space–time covariance structure of the climate to re-weight the possible age models. Here we demonstrate how Bayesian Hierarchical climate reconstruction models can be augmented to account for time uncertain proxies. Critically, while a priori all age models are given equal probability of being correct, the probabilities associated with the age models are formally updated within the Bayesian framework, thereby reducing uncertainties. Numerical experiments show that updating the age-model probabilities decreases uncertainty in the climate reconstruction, as compared with the current de-facto standard of sampling over all age models, provided there is sufficient information from other data sources in the region of the time-uncertain proxy. This approach can readily be generalized to non-layer counted proxies, such as those derived from marine sediments.

Description: High-resolution records of the beryllium-10 solar activity proxy in ice from Law Dome, East Antarctica: measurement, reproducibility and principal trends Climate of the Past Discussions, 7, 677-713, 2011 Author(s): J. B. Pedro, A. M. Smith, K. J. Simon, T. D. van Ommen, and M. A. J. Curran Three near-monthly resolution 10 Be records are presented from the Dome Summit South (DSS) ice core site, Law Dome, East Antarctica. The chemical preparation and Accelerator Mass Spectrometer (AMS) measurement of these records is described. The reproducibility of 10 Be records at DSS is assessed through intercomparison of the ice core data with data from two previously published and contemporaneous snow pits. We find generally good agreement between the five records, comparable to that observed between other trace chemical records from the site. This result allays concerns raised by a previous Antarctic study (Moraal et al., 2005) about poor reproducibility of ice core 10 Be records. A single composite series is constructed from the three ice cores providing a monthly-resolved record of 10 Be concentrations at DSS over the past decade (1999 to 2009). To our knowledge, this is the first published ice core data spanning the recent exceptional solar minimum of solar cycle 23. 10 Be concentrations are significantly correlated to the cosmic ray flux recorded by the McMurdo neutron monitor ( r xy = 0.64, with 95% CI of 0.53 to 0.71), suggesting that solar modulation of the atmospheric production rate may explain up to ~40% of the variance in 10 Be concentrations at DSS. Sharp concentration peaks occur in most years during the summer-to-autumn, possibly caused by stratospheric incursions. Our results underscore the presence of both production and meteorological signals in ice core 10 Be data.

Description: Tracking climate variability in the western Mediterranean during the Late Holocene: a multiproxy approach Climate of the Past Discussions, 7, 635-675, 2011 Author(s): V. Nieto-Moreno, F. Martínez-Ruiz, S. Giralt, F. Jiménez-Espejo, D. Gallego-Torres, M. Rodrigo-Gámiz, J. García-Orellana, M. Ortega-Huertas, and G. J. de Lange Climate variability in the western Mediterranean is reconstructed for the last 4000 yr using marine sediments recovered in the west Algerian-Balearic basin, near the Alboran basin. Fluctuations in chemical and mineralogical sediment composition as well as grain size distribution are linked to fluvial-eolian oscillations, changes in redox conditions and paleocurrent intensity. Multivariate analyses allowed us to characterize three main groups of geochemical and mineralogical proxies determining the sedimentary record of this region. These three statistical groups were applied to reconstruct paleoclimate conditions at high resolution during the Late Holocene. An increase in fluvial-derived elements (Rb/Al, Ba/Al, REE/Al, Si/Al, Ti/Al, Mg/Al and K/Al ratios), finer grain size, slower flows and oxygen-poor bottom waters – as suggested by sortable silt (10–63 μm), clays ( 〈 2 μm) and redox-sensitive elements (V/Al, Cr/Al, Ni/Al and Zn/Al ratios) – depict the Roman Humid Period (RHP) and the Little Ice Age (LIA), while drier environmental conditions are recognized during the Late Bronze Age-Iron Age (LBA-IA) and the Medieval Warm Period (MWP). Although no Ba excess was registered, other paleoproductivity indicators (total organic carbon content, Br/Al ratio, and organometallic ligands such as U and Cu) display the highest values during the RHP, this period exhibiting by far the most intense productivity of the last 4000 yr. These marine sediments evidence oscillations that support the link of the westernmost Mediterranean climate with the North Atlantic coupled ocean-atmosphere climatic system, pointing to solar irradiance and the North Atlantic Oscillation (NAO) variability as the main driving mechanisms behind natural climate variability over decadal to centennial time-scales for the last 4000 yr.

Description: Ionospheric variations before some large earthquakes over Sumatra Natural Hazards and Earth System Science, 11, 597-611, 2011 Author(s): A. M. Hasbi, M. A. Mohd Ali, and N. Misran The paper investigates the ionospheric variations before some large earthquakes that occurred during 2004–2007 in Sumatra using GPS and CHAMP data. The TEC shows the occurrence of positive and negative anomalies detected within a few hours to 6 days before the earthquakes. These anomalies mostly occur during the daytime hours between 4 and 17 LT. The TEC anomalies are mostly consistent with the CHAMP satellite electron density data. The electron density analysis over the 28 March 2005 earthquake epicenter shows that an equatorial anomaly modification took place a few days before the event. The modification took shape in the form of crest amplification during the daytime. The comparison between the TEC and electron density measurements during very quiet geomagnetic conditions is shown to be a useful indicator of a forthcoming earthquake.

Description: Critical analysis of the electrostatic turbulence enhancements observed by DEMETER over the Sichuan region during the earthquake preparation Natural Hazards and Earth System Science, 11, 561-570, 2011 Author(s): T. Onishi, J.-J. Berthelier, and M. Kamogawa In this paper, we report initial results from a detailed statistical study of the plasma waves observed by the DEMETER satellite over the Sichuan region during a period of 20 days encompassing the large earthquake of magnitude M =7.9 that occurred on 12 May 2008. The main objective of this paper is to present a statistical method to process and analyze plasma wave data and assist in detecting possible earthquake precursors among larger irregular disturbances arising from the natural variability of the ionized environment of the Earth. This method, presently used for dayside observations, involves two stages. First, VLF wave spectra are processed to recognize the various types of plasma waves usually observed at mid and low latitudes and derive a reduced number of parameters that fully characterize these emissions and may be conveniently used for a detailed statistical study. In a second stage, we perform a statistical analysis of the results by taking into account two "reference zones" displaced respectively 30 ° eastward and westward from the "epicentre zone". Plasma and wave disturbances possibly induced by earthquakes in preparation are likely to maximize close to the "epicentre zone", while natural disturbances associated, in particular, with the varying magnetic activity are rather uniform over a wider longitude sector, thus enabling the use of observations over the reference zones as a base line. The initial results of this study show a deviation of the power spectrum of electrostatic turbulence in the epicentre zone about 6 days prior to the earthquake but no significant anomalous variations can be observed on other characteristics of plasma waves. From the analysis of the data over the two reference~zones and using recently produced sector magnetic activity indices, we conclude that the enhancement of electrostatic turbulence is associated with magnetospheric processes rather than with pre-seismic activity.

Description: Impact of maximum borehole depths on inverted temperature histories in borehole paleoclimatology Climate of the Past Discussions, 7, 715-748, 2011 Author(s): H. Beltrami, J. E. Smerdon, G. S. Matharoo, and N. Nickerson A quantitative assessment is presented for the impact of the maximum depth of a temperature-depth profile on the estimate of the climatic transient and the resultant ground surface temperature (GST) reconstruction used in borehole paleoclimatology. The depth of the profile is important because the downwelling climatic signal must be separated from the quasi-steady state thermal regime established by the energy in the Earth's interior. This component of the signal is estimated as a linear increase in temperature with depth from the lower section of a borehole temperature profile, which is assumed to be unperturbed by recent changes in climate at the surface. The validity of this assumption is dependent on both the subsurface thermophysical properties and the character of the downwelling climatic signal. Such uncertainties can significantly impact the determination of the quasi-steady state thermal regime, and consequently the magnitude of the temperature anomaly interpreted as a climatic signal. The quantitative effects and uncertainties that arise from the analysis of temperature-depth profiles of different depths are presented. Results demonstrate that widely different GST histories can be derived from a single temperature profile truncated at different depths. Borehole temperature measurements approaching 500–600 m depths are shown to provide the most robust GST reconstructions spanning 500 to 1000 ybp. It is further shown that the bias introduced by a temperature profile of depths shallower than 500–600 m remains even if the time span of the reconstruction target is shortened.

Description: Evolution of extreme Total Water Levels along the northern coast of the Iberian Peninsula Natural Hazards and Earth System Science, 11, 613-625, 2011 Author(s): D. F. Rasilla Álvarez and J. C. García Codron This paper assesses the evolution of storminess along the northern coast of the Iberian Peninsula through the calculation of extreme (1%) Total Water Levels (eTWL) on both observed (tide gauge and buoy data) and hindcasted (SIMAR-44) data. Those events were first identified and then characterized in terms of oceanographic parameters and atmospheric circulation features. Additionally, an analysis of the long-term trends in both types of data was performed. Most of the events correspond to a rough wave climate and moderate storm surges, linked to extratropical disturbances following a northern track. While local atmospheric conditions seem to be evolving towards lesser storminess, their impact has been balanced by the favorable exposure of the northern coast of the Iberian Peninsula to the increasing frequency and strength of distant disturbances crossing the North Atlantic. This evolution is also correctly reproduced by the simulated long-term evolution of the forcing component (meteorological sea level residuals and wave run up) of the Total Water Level values calculated from the SIMAR 44 database, since sea level residuals have been experiencing a reduction while waves are arriving with longer periods. Finally, the addition of the rate of relative sea level trend to the temporal evolution of the atmospheric forcing component of the Total Water Level values is enough to simulate more frequent and persistent eTWL.

Description: Upper ocean climate of the Eastern Mediterranean Sea during the Holocene Insolation Maximum – a model study Climate of the Past Discussions, 7, 1457-1509, 2011 Author(s): F. Adloff, U. Mikolajewicz, M. Kucera, R. Grimm, E. Maier-Reimer, G. Schmiedl, and K. Emeis Nine thousand years ago, the Northern Hemisphere experienced enhanced seasonality caused by an orbital configuration with a minimum of the precession index. To assess the impact of the "Holocene Insolation Maximum" (HIM) on the Mediterranean Sea, we use a regional ocean general circulation model forced by atmospheric input derived from global simulations. A stronger seasonal cycle is simulated in the model, which shows a relatively homogeneous winter cooling and a summer warming with well-defined spatial patterns, in particular a subsurface warming in the Cretan and Western Levantine areas. The comparison between the SST simulated for the HIM and the reconstructions from planktonic foraminifera transfer functions shows a poor agreement, especially for summer, when the vertical temperature gradient is strong. However, a reinterpretation of the reconstructions is proposed, to consider the conditions throughout the upper water column. Such a depth-integrated approach accounts for the vertical range of preferred habitat depths of the foraminifera used for the reconstructions and strongly improves the agreement between modelled and reconstructed temperature signal. The subsurface warming is recorded by both model and proxies, with a light shift to the south in the model results. The mechanisms responsible for the peculiar subsurface pattern are found to be a combination of enhanced downwelling and wind mixing due to strengthened Etesian winds, and enhanced thermal forcing due to the stronger summer insolation in the Northern Hemisphere. Together, these processes induce a stronger heat transfer from the surface to the subsurface during late summer in the Western Levantine; this leads to an enhanced heat piracy in this region.