Mineralogy and geochemistry of the Texeo Cu–Co mine site (NW Spain): screening tools for environmental assessment

The Cu–Co–Ni Texeo mine has been the most important source of Cu in NW Spain since Roman times and now, approximately 40,000 m³ of wastes from mine and metallurgical operations, containing average concentrations of 9,263 mg kg⁻¹ Cu, 1,100 mg kg⁻¹ As, 549 mg kg⁻¹ Co, and 840 mg kg⁻¹ Ni, remain on-site. Since the cessation of the activity, the abandoned works, facilities and waste piles have been posing a threat to the environment, derived from the release of toxic elements. In order to assess the potential environmental pollution caused by the mining operations, a sequential sampling strategy was undertaken in wastes, soil, surface and groundwater, and sediments. First, screening field tools were used to identify hotspots, before defining formal sampling strategies; so, in the areas where anomalies were detected in a first sampling stage, a second detailed sampling campaign was undertaken. Metal concentrations in the soils are highly above the local background, reaching up to 9,921 mg kg⁻¹ Cu, 1,373 mg kg⁻¹ As, 685 mg kg⁻¹ Co, and 1,040 mg kg⁻¹ Ni, among others. Copper concentrations downstream of the mine works reach values up to 1,869 μg l⁻¹ and 240 mg kg⁻¹ in surface water and stream sediments, respectively. Computer-based risk assessment for the site gives a carcinogenic risk associated with the presence of As in surface waters and soils, and a health risk for long exposures; so, trigger levels of these elements are high enough to warrant further investigation.     

Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Seismic Hazard Analysis — Quo vadis?

The paper is dedicated to the review of methods of seismic hazard analysis currently in use, analyzing the strengths and weaknesses of different approaches. The review is performed from the perspective of a user of the results of seismic hazard analysis for different applications such as the design of critical and general (non-critical) civil infrastructures, technical and financial risk analysis. A set of criteria is developed for and applied to an objective assessment of the capabilities of different analysis methods. It is demonstrated that traditional probabilistic seismic hazard analysis (PSHA) methods have significant deficiencies, thus limiting their practical applications. These deficiencies have their roots in the use of inadequate probabilistic models and insufficient understanding of modern concepts of risk analysis, as have been revealed in some recent large scale studies. These deficiencies result in the lack of ability of a correct treatment of dependencies between physical parameters and finally, in an incorrect treatment of uncertainties. As a consequence, results of PSHA studies have been found to be unrealistic in comparison with empirical information from the real world. The attempt to compensate these problems by a systematic use of expert elicitation has, so far, not resulted in any improvement of the situation. It is also shown that scenario-earthquakes developed by disaggregation from the results of a traditional PSHA may not be conservative with respect to energy conservation and should not be used for the design of critical infrastructures without validation. Because the assessment of technical as well as of financial risks associated with potential damages of earthquakes need a risk analysis, current method is based on a probabilistic approach with its unsolved deficiencies.

Traditional deterministic or scenario-based seismic hazard analysis methods provide a reliable and in general robust design basis for applications such as the design of critical infrastructures, especially with systematic sensitivity analyses based on validated phenomenological models. Deterministic seismic hazard analysis incorporates uncertainties in the safety factors. These factors are derived from experience as well as from expert judgment. Deterministic methods associated with high safety factors may lead to too conservative results, especially if applied for generally short-lived civil structures. Scenarios used in deterministic seismic hazard analysis have a clear physical basis. They are related to seismic sources discovered by geological, geomorphologic, geodetic and seismological investigations or derived from historical references. Scenario-based methods can be expanded for risk analysis applications with an extended data analysis providing the frequency of seismic events. Such an extension provides a better informed risk model that is suitable for risk-informed decision making.     

An operational multiscale system for hazards prediction, mapping, and response

By definition, a crisis is a situation that requires assistance to be managed. Hence, response to a crisis involves the merging of local and non-local emergency response personnel. In this situation, it is critical that each participant: (1) know the roles and responsibilities of each of the other participants; (2) know the capabilities of each of the participants; and (3) have a common basis for action. For many types of natural disasters, this entails having a common operational picture of the unfolding events, including detailed information on the weather, both current and forecasted, that may impact on either the emergency itself or on response activities. The Consequences Assessment Tool Set (CATS) is a comprehensive package of hazard prediction models and casualty and damage assessment tools that provides a linkage between a modeled or observed effect and the attendant consequences for populations, infrastructure, and resources, and, hence, provides the common operational picture for emergency response. The Operational Multiscale Environment model with Grid Adaptivity (OMEGA) is an atmospheric simulation system that links the latest methods in computational fluid dynamics and high-resolution gridding technologies with numerical weather prediction to provide specific weather analysis and forecast capability that can be merged into the geographic information system framework of CATS. This paper documents the problem of emergency response as an end-to-end system and presents the integrated CATS–OMEGA system as a prototype of such a system that has been used successfully in a number of different situations.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

Methodology for preliminary assessment of Natech risk in urban areas

Concern for natural hazard-triggered technological disasters (Natech disasters) in densely populated and industrialized areas is growing. Residents living in urban areas subject to high natural hazard risk are often unaware of the potential for secondary disasters such as hazardous materials releases from neighboring industrial facilities, chemical storage warehouses or other establishments housing hazardous materials. Lessons from previous disasters, such as the Natech disaster during the Kocaeli earthquake in Turkey in 1999 call for the need to manage low frequency/high consequence events, particularly in today’s densely populated areas. However, there is little guidance available on how local governments and communities can assess Natech risk. To add to the problem, local governments often do not have the human or economic resources or expertise to carry out detailed risk assessments. In this article, we propose a methodology for preliminary assessment of Natech risk in urban areas. The proposed methodology is intended for use by local government officials in consultation with the public. The methodology considers possible interactions between the various systems in the urban environment: the physical infrastructure (e.g., industrial plants, lifeline systems, critical facilities), the community (e.g., population exposed), the natural environment (e.g., delicate ecosystems, river basins), and the risk and emergency management systems (e.g., structural and nonstructural measures). Factors related to vulnerability and hazard are analyzed and qualitative measures are recommended. Data from hazardous materials releases during the Kocaeli, Turkey earthquake of August 17, 1999 are used as a case study to demonstrate the applicability of the methodology. Limitations of the proposed methodology are discussed as well as future research needs.

Reducing hazard vulnerability through local government engagement and action

The concept of a natural hazard is a human construct. It is the interaction with human communities and settlements that defines a natural phenomenon as a natural hazard. Thus the end point of hazard mitigation and hazard vulnerability assessment must involve an attempt to reduce, or mitigate, the impact of the natural hazard on human communities. The responsibility to mitigate hazard impact falls primarily upon governments and closely connected non-government and private institutional agencies. In particular, it is most often local government that takes the responsibility for safeguarding its own communities, infrastructure and people. Hazard vulnerability of specific local communities is best assessed by the local government or council, which then faces the responsibility to translate that assessment into community education and infrastructural safeguards for hazard mitigation. This paper illustrates the process of local government engagement in hazard mitigation in Australia, through the Natural Disaster Risk Management Studies, as a first step towards natural disaster reduction.     

Nonlinear Analysis of Local Site Effects on Seismic Ground Response in the Bam Earthquake

The influence of local geologic and soil conditions on the intensity of ground shaking is addressed in this study. The amplification of the ground motion due to local site effects resulted in severe damage to dwellings in the Bam area during the 2003 Bam Earthquake. A unique set of strong motion acceleration recordings was obtained at the Bam accelerograph station. Although the highest peak ground acceleration recorded was the vertical component (nearly 1 g), the longitudinal component (fault-parallel motion) clearly had the largest maximum velocity as well as maximum ground displacement. Subsurface geotechnical and geophysical (down-hole) data in two different sites have been obtained and used to estimate the local site condition on earthquake ground motion in the area. The ground response analyses have been conducted considering the nonlinear behavior of the soil deposits using both equivalent linear and nonlinear approaches. The fully nonlinear method embodied in FLAC was used to evaluate the nonlinear soil properties on earthquake wave propagation through the soil layer, and compare with the response from the equivalent linear approach. It is shown that thick alluvium deposits amplified the ground motion and resulted in significant damage in residential buildings in the earthquake stricken region. The comparison of results indicated similar response spectra of the motions for both equivalent and nonlinear analyses, showing peaks in the period range of 0.3–1.5 s. However, the amplification levels of nonlinear analysis were less than the equivalent linear method especially in long periods. The observed response spectra are shown to be above the NEHRP building code design requirements, especially at high frequencies.     

乌鲁木齐河流域地下水水质监测网设计

文章运用地下水易污性编图及污染源分布图法进行了乌鲁木齐河流域地下水水质监测网设计。共设计了130监测孔,现有46个监测孔,另需要84个新的监测孔。按监测类型分为面源监测点22个,点源监测点87个,重点水源地与泉水监测点21个。按监测运行分长期监测点55个,流域普查监测点75个。普查监测点监测频率为1次/5年,长期监测点监测频率为1次/年。首期有针对性地在污染严重的柴窝堡新化厂排污区、乌鲁木齐河谷老排污区、米泉污灌区、米泉工业污染区、老龙河污染区取了25个污染水样测试分析,结果显示地下水已经受到严重污染。     

浙江省小流域山洪灾害临界雨量确定方法分析

小流域山洪灾害具有突发性,预测预警难度较大.本文结合浙江省小流域山洪灾害防御的实践经验,研究提出了以水位反推法计算临界雨量,简单实用,具有较好的实践价值和推广意义.