PSHA of Van province for performance assessment using spectrally matched strong ground motion records

Within the framework of the performance based earthquake engineering, site specific earthquake spectra for Van province has been obtained. It is noteworthy that, in probabilistic seismic hazard assessment, as a first stage data from geological studies and records from the instrumental period were compiled to make a seismic source characterization for the study region. The probabilistic seismic hazard curves were developed based on selected appropriate attenuation relationships, at rock sites, with a probability of exceedance 2, 10 and 50% in 50 yrs period. The obtained results are compared with the spectral responses proposed for seismic evaluation and retrofit of building structure in Turkish Earthquake Code (2007), section 7. The acceleration response spectrums obtained from probabilistic seismic hazard analysis are matched to adjust earthquake accelerograms recorded during the 2011 Van earthquakes by using SeismoMatch v2.0 software. The aim of this procedure is to obtain a set of reasonable earthquake input motions for the seismic evaluation of existing buildings.     

Assessment of Liquefaction Potential of Guwahati City: A Case Study

The liquefaction potential of saturated cohesionless deposits in Guwahati city, Assam, was evaluated. The critical cyclic stress ratio required to cause liquefaction and the cyclic stress ratio induced by an earthquake were obtained using the simplified empirical method developed by Seed and Idriss (J soil Mech Found Eng ASCE 97(SM9):1249–1273, 1971, Ground motions and soil liquefaction during earthquakes. Earthquake Engineering Research Institute, Berkeley, CA, 1982) and Seed et al. (J Geotech Eng ASCE 109(3):458–483, 1983, J Geotech Eng ASCE 111(12):1425–1445, 1985) and the Idriss and Boulanger (2004) method. Critical cyclic stress ratio was based on the empirical relationship between standard penetration resistance and cyclic stress ratio. The liquefaction potential was evaluated by determining factor of safety against liquefaction with depth for areas in the city. A soil database from 200 boreholes covering an area of 262 km² was used for the purpose. A design peak ground acceleration of 0.36 g was used since Guwahati falls in zone V according to the seismic zoning map of India. The results show that 48 sites in Guwahati are vulnerable to liquefaction according to the Seed and Idriss method and 49 sites are vulnerable to liquefaction according to the Idriss and Boulanger method. Results are presented as maps showing zones of levels of risk of liquefaction.     

Seismic hazard analysis of Mersin Province,Turkey using probabilistic and statistical methods

Nearly 108-km lengths of Mersin shores are composed of natural beaches. The region is located between major tourist centers. In the future, this region is thought to be built with a great number of tourist facilities. Turkey’s largest seaport, Ata? refinery (Mersin International Port) is located in Mersin. Recently, Mersin is becoming of great importance to Turkey as the latter plans to construct its second nuclear power plant in the region. Therefore, as nuclear power plants are built to withstand environmental hazards, it is very important to analyze the seismic risk of the areas where the nuclear power plant will be constructed. The region is located between the East Anatolian Fault Zone and Center Anatolian Fault Zone. Based on the Turkey Earthquake Regions Map, Mersin is divided into second-, third-, and fourth-degree earthquake regions. In this study, we sampled earthquakes of magnitude of 4.0 or greater between 01 Jan 1900 and 31 Dec 2010 in the area; seismic hazard of Mersin province was estimated with probabilistic and statistical methods. The study area was selected as the coordinates between 36.03° and 37.42° North and 32.57° and 35.16° East. On the study area, different scaled magnitude values in the last 110 years converted to a common scale (Mw) and earthquake catalog was re-compiled and also seismic sources that may affect the area was determined. In this study, the seismic hazards of the region were obtained using the methods of probability and statistics. This study used three different attenuation relationships. Using the attenuation relationships suggested by Boore et al. (Seismol Res Lett 68(1):128–153, 1997) and Kalkan and Gülkan (Earthquake Spectra 20:1111–1138, 2004), the largest ground acceleration which corresponds to a recurrence period of 475 years was found as 0.08–0.09 g and Akkar and Ça?nan (Bull Seismol Soc Am 100 6:2978–2995, 2010), 0.04 g for bedrock at the central district. When computing for seismic hazard curves, Mut district appears to have a greater seismic hazard compared with other districts. Moreover, according to the attenuation relationships, seismic hazard curves corresponding to a recurrence period of 475 years were obtained for the Mersin Central, Mut, Erdemli, Çaml?yayla, and Tarsus districts.     

Assessment of liquefaction potential of Erzincan Province and its vicinity,Turkey

One of the most important causes of damages after the earthquakes is the soil liquefaction. Liquefaction can be defined as temporary loss in strength of saturated sandy and silty deposits under transient and cyclic loadings due to excess pore water pressure. This study includes determination of liquefaction potential in Erzincan city center and its vicinity. Due to the proximity of the North Anatolian Fault Zone, in a probable earthquake, Erzincan Province is thought to be affected. In this context, the earthquake scenarios were produced using the empirical expressions. Liquefaction potential for different earthquake magnitudes was determined. These earthquake magnitudes were selected as 6.0, 6.5, 7.0, 7.5, respectively. Liquefaction potential was investigated using standard penetration test (SPT) data. The first stage of the study, 63 boreholes in different locations was drilled and SPT was performed. Disturbed and undisturbed soil samples were taken from these boreholes. Laboratory testing was performed to determine physical properties of soil samples, and liquefaction potential analyses were examined using three methods, namely Seed and Idriss (J Soil Mech Found Div ASCE 97(9):1249–1273, 1971), Tokimatsu and Yoshimi (Soil Found 23(4):56–74, 1983), Iwasaki et al. (Soil liquefaction potential evaluation with use of the simplified procedure. International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics, St. Louis, pp 209–214, 1981). In order to complete liquefaction analysis within a short time, MATLAB program was prepared. Liquefaction potential analyses were carried out with the MATLAB program. At the final stage of this study, liquefaction potential maps were prepared for different earthquake magnitudes. The expected results will be shared with the local authorities and important engineering remedial measurements will be proposed to prevent further life losses and to mitigate property losses.     

Seismic hazard analysis of Sinop province,Turkey using probabilistic and statistical methods

Using 4.0 and greater magnitude earthquakes which occurred between 1 January 1900 and 31 Dec 2008 in the Sinop province of Turkey this study presents a seismic hazard analysis based on the probabilistic and statistical methods. According to the earthquake zonation map, Sinop is divided into first, second, third and fourth-degree earthquake regions. Our study area covered the coordinates between 40.66°– 42.82°N and 32.20°– 36.55°E. The different magnitudes of the earthquakes during the last 108 years recorded on varied scales were converted to a common scale (Mw). The earthquake catalog was then recompiled to evaluate the potential seismic sources in the aforesaid province. Using the attenuation relationships given by Boore et al. (1997) and Kalkan and Gülkan (2004), the largest ground accelerations corresponding to a recurrence period of 475 years are found to be 0.14 g for bedrock at the central district. Comparing the seismic hazard curves, we show the spatial variations of seismic hazard potential in this province, enumerating the recurrence period in the order of 475 years.     

Convergence of iterative coupling for coupled flow and geomechanics

In this paper, we study solving iteratively the coupling of flow and mechanics. We demonstrate the stability and convergence of two widely used schemes: the undrained split method and the fixed stress split method. To our knowledge, this is the first time that such results have been rigorously obtained and published in the scientific literature. In addition, we propose a new stress split method, with faster convergence rate than known schemes. These results are specially important today due to the interest in hydraulic fracturing (Dean and Schmidt SPE J. 14:707–714, 2009; Ji et al. SPE J. 14:423–430, 2009; Samier and De Gennaro 2007; Settari and Maurits SPE J. 3:219–226, 1998), in oil and gas shale reservoirs.     

Chromite alteration processes within Vourinos ophiolite

The renewed interest in chromite ore deposits is directly related to the increase in Cr price ruled by international market trends. Chromite, an accessory mineral in peridotites, is considered to be a petrogenetic indicator because its composition reflects the degree of partial melting that the mantle experienced while producing the chromium spinel-bearing rock (Burkhard in Geochim Cosmochim Acta 57:1297–1306, 1993). However, the understanding of chromite alteration and metamorphic modification is still controversial (e.g. Evans and Frost in Geochim Cosmochim Acta 39:959–972, 1975; Burkhard in Geochim Cosmochim Acta 57:1297–1306, 1993; Oze et al. in Am J Sci 304:67–101, 2004). Metamorphic alteration leads to major changes in chromite chemistry and to the growth of secondary phases such as ferritchromite and chlorite. In this study, we investigate the Vourinos complex chromitites (from the mines of Rizo, Aetoraches, Xerolivado and Potamia) with respect to textural and chemical analyses in order to highlight the most important trend of alteration related to chromite transformation. The present study has been partially funded by the Aliakmon project in collaboration between the Public Power Corporation of Greece and Institute of Geology and Mineral Exploration of Kozani.     

Investigation of the variability of strong ground motions from Vrancea earthquakes

A systematic investigation of the applicability of several ground motion prediction models for Vrancea intermediate-depth seismic source is conducted in this research. Two ground motion prediction models recommended by previous evaluations (Vacareanu et al. in Bull Earthq Eng 11(6):1867–1884, 2013a; Pavel et al. in Earthq Struct 6(1):1–18, 2014), as well as two new state-of-the-art ground motion prediction equations (Vacareanu et al. in J Earthq Eng, 2013b; Earthq Struct 6(2):141–161, 2014) are tested using an increased strong ground motion database consisting of 150 recordings from Vrancea subcrustal earthquakes. The evaluation is performed by using several goodness-of-fit parameters from the literature. Moreover, the applicability of the single-station sigma method is also investigated by using the same strong ground motion database recorded in 30 seismic stations from southern and eastern Romania. The influence of the soil conditions on the numerical results obtained in this study is investigated and discussed using the results provided by the analysis of variance method. The impact of the single-station standard deviation on the levels of seismic hazard is also assessed in this study, and the results show, in the analyzed cases, significant reductions of the hazard levels.     

Probabilistic seismic hazard analysis of Tripura and Mizoram states

A probabilistic seismic hazard analysis for the states of Tripura and Mizoram in North East India is presented in this paper to evaluate the ground motion at bedrock level. Analyses were performed considering the available earthquake catalogs collected from different sources since 1731–2010 within a distance of 500 km from the political boundaries of the states. Earthquake data were declustered to remove the foreshocks and aftershocks in time and space window and then statistical analysis was carried out for data completeness. Based on seismicity, tectonic features and fault rupture mechanism, this region was divided into six major seismogenic zones and subsequently seismicity parameters (a and b) were calculated using Gutenberg–Richter (G–R) relationship. Faults data were extracted from SEISAT (Seismotectonic atlas of India, Geological Survey of India, New Delhi, 2000) published by Geological Survey of India and also from satellite images. The study area was divided into small grids of size 0.05° × 0.05° (approximately 5 km × 5 km), and the hazard parameters (rock level peak horizontal acceleration and spectral accelerations) were calculated at the center of each of these grid cells considering all the seismic sources within a radius of 500 km. Probabilistic seismic hazard analyses were carried out for Tripura and Mizoram states using the predictive ground motion equations given by Atkinson and Boore (Bull Seismol Soc Am 93:1703–1729, 2003) and Gupta (Soil Dyn Earthq Eng 30:368–377, 2010) for subduction belt. Attenuation relations were validated with the observed PGA values. Results are presented in the form of hazard curve, peak ground acceleration (PGA) and uniform hazard spectra for Agartala and Aizawl city (respective capital cities of Tripura and Mizoram states). Spatial variation of PGA at bedrock level with 2 and 10 % probability of exceedance in 50 years has been presented in the paper.     

Impedance Functions of Slab Foundations with Rigid Piles

The earthquake resistant design of structures depends upon the soil-structure interaction during seismic excitation. The dynamic behavior of surface foundations and deep foundations has been investigated by several authors (Gazetas in Soil Dyn Earthq Eng 2(1):2–42, 1983; Pecker in Dynamique des sols. Presses de l’Ecole Nationale des Ponts et Chaussées, Paris, 1984; Sieffert and Cevaer in Manuel des fonctions d’impédance. Ouest-Editions, Nantes, 1992, etc.) but the dynamic behavior of pile-reinforced soils has not been sufficiently studied yet. The design of pile-reinforced soils comprises the design of the rigid piles and the design of the earth-platform. The foundation system studied in this article consists of an earth-platform over a soft ground reinforced by deep piles. In order to understand the dynamic behavior of a rigid pile, a series of dynamic tests is conducted on an experimental site. The vertical and horizontal impedances of the slab foundation are obtained with and without rigid piles. The numerical models are developed to interpret these dynamic tests. The numerical and experimental impedance functions are compared in both the vertical and the horizontal directions. A sensitivity analysis on the influence of the physical and geometrical properties of rigid piles on the impedance functions is presented.     

Effect of frequency-dependent radiation pattern in the strong motion simulation of the 2011 Tohoku earthquake,Japan, using modified semi-empirical method

We perform a strong ground motion simulation using a modified semi-empirical technique (Midorikawa in Tectonophysics 218:287–295, 1993), with frequency-dependent radiation pattern model. Joshi et al. (Nat Hazards 71:587–609, 2014) have modified the semi-empirical technique to incorporate the modeling of strong motion generation areas (SMGAs). A frequency-dependent radiation pattern model is applied to simulate high-frequency ground motion more precisely. Identified SMGAs (Kurahashi and Irikura in Earth Planets Space 63:571–576, 2011) of the 2011 off the Pacific coast of Tohoku earthquake (M _w  = 9.0) were modeled using this modified technique. We analyzed the effect of changing seismic moment values of SMGAs on the simulated acceleration time series. Final selection of the moment values of SMGAs is based on the root-mean-square error (RMSE) of waveform comparison. Records are simulated for both frequency-dependent and constant radiation pattern function. Simulated records for both cases are compared with observed records in terms of peak ground acceleration, peak ground velocity and pseudo-acceleration response spectra at different stations. Comparison of simulated and observed records in terms of RMSE suggests that the method is capable of simulating record, which matches in a wide frequency range for this earthquake and bears realistic appearance in terms of shape and strong motion parameters. The results confirm the efficacy and suitability of rupture model defined by five SMGAs for the developed modified technique.     

Applicability of different ground-motion prediction models for northern Iran

A total of 163 free-field acceleration time histories recorded at epicentral distances of up to 200 km from 32 earthquakes with moment magnitudes ranging from M _w 4.9 to 7.4 have been used to investigate the predictive capabilities of the local, regional, and next generation attenuation (NGA) ground-motion prediction equations and determine their applicability for northern Iran. Two different statistical approaches, namely the likelihood method (LH) of Scherbaum et al. (Bull Seismol Soc Am 94:341–348, 2004) and the average log-likelihood method (LLH) of Scherbaum et al. (Bull Seismol Soc Am 99:3234–3247, 2009), have been applied for evaluation of these models. The best-fitting models (considering both the LH and LLH results) over the entire frequency range of interest are those of Ghasemi et al. (Seismol 13:499–515, 2009a) and Soghrat et al. (Geophys J Int 188:645–679, 2012) among the local models, Abrahamson and Silva (Earthq Spectra 24:67–97, 2008) and Chiou and Youngs (Earthq Spectra 24:173–215, 2008) among the NGA models, and finally Akkar and Bommer (Seism Res Lett 81:195–206, 2010) among the regional models.     

Solution of the MoMaS reactive transport benchmark with MIN3P—model formulation and simulation results

This paper summarizes the governing equations as implemented in the MIN3P multicomponent flow and reactive transport code (Mayer et al., Water Resour Res 38:1174, 2002) and introduces the equations in discretized form. Linearization and solution methods are presented including adaptive time stepping and update modification schemes. Code-specific details for the implementation of the GdR MoMaS benchmark simulations (Carrayrou et al., Comput Geosci, 2009) are presented. The standard version of the MIN3P code was used to solve the Easy, Medium, and Hard Test Cases, in one and two spatial dimensions, for both advection- and diffusion-dominated conditions. An analysis of the sensitivity of the solution in relation to spatial and temporal discretization parameters is provided for the Easy Test Case, selected results are presented for the Medium and Hard Test Cases, and the performance of the code as a function of discretization parameters is evaluated for all test cases.     

Geology and geochemistry of Tertiary basalt in south Wadi Hodein area, South Eastern Desert, Egypt

Tertiary basalt is widespread in the area south of Wadi Hodein, south Eastern Desert, Egypt. It is the youngest unit in the basement rocks of the Central Eastern Desert classification of El Shazly (Proc 22nd Intl Geol Congr, New Delhi 10:88–101, 1964) and El Ramly (Ann Geol Surv Egypt II:1–17, 1972), traversed all the previous succession of the basement rocks as well as the Nubia Sandstone of Cretaceous age, forming sheets, small hills, ridges, and dikes. This Tertiary basalt is strongly associated with the opening of the Red Sea. Geologic, petrographic, and petrochemical studies as well as microprobe and X-ray analyses were performed on samples from Wadi Hodein Tertiary basalt. Field and petrographic studies classified the Tertiary basalt in south Wadi Hodein into porphyritic olivine basalt, plagiophyric basalt, and doleritic basalt. Opaque minerals (magnetite and ilmenite) constitute 6–7.5% of this basalt. Petrochemical studies and microprobe analyses reveal that they are low-TiO₂ basalt with low uranium and thorium contents, classified as being basaltic andesite to andesite, originated from calc-alkaline magma, and developed in within-plate tectonic environment. Scanning electron microscopy shows that magnetite and ilmenite are the prevalent opaque minerals in this Tertiary basalt. Field radiometric measurements of the Tertiary basalt in south Wadi Hodein reveals low uranium and thorium contents. Uranium contents range from 0.5 to 0.9 ppm, while thorium contents range from 1.2 to 3.2 ppm. Fractional crystallization and mass balance modeling indicate that the most-silica low-TiO₂ Tertiary basalt in south Wadi Hodein can be derived from the relatively less-silica low-TiO₂ Tertiary basalt of south Quseir and Gabal Qatrani through fractional crystallization of plagioclase, olivine, augite, and titanomagnetite oxides. Tertiary basalts in south Wadi Hodein and south Quseir have nearly the same age, 25 Ma (Sherif, The Fifth International Conference on the Geology of Africa, 2007), 24 Ma (Meneisy and Abdel Aal, Ain Shams Sci Bull 25(24B): 163–176, 1984), and 27 Ma (El Shazly et al., Egypt J Geol 1975), respectively. Finally, the fractionation modeling and geochemical characteristics of these basalts suggested their origination from one basaltic magma emplaced in late Oligocene.     

Comment on “A Comparison of Modified Fuzzy Weights of Evidence,Fuzzy Weights of Evidence,and Logistic Regression for Mapping Mineral Prospectivity” by Daojun Zhang,Frits Agterberg,Qiuming Cheng,and Renguang Zuo Math Geosci DOI 10.1007/s11004-013-9496-8

Despite a missing definition of equivalence of mathematical models or methods by Zhang et al. (Math Geosci, 2013), an “equivalence” (Zhang et al., Math Geosci, 2013, p. 6,7,8,14) of modified weights-of-evidence (Agterberg, Nat Resour Res 20:95–101, 2011) and logistic regression does not generally exist. Its alleged proof is based on a previously conjectured linear relationship between weights of evidence and logistic regression parameters (Deng, Nat Resour Res 18:249–258, 2009), which does not generally exist either (Schaeben and van den Boogaart, Nat Resour Res 20:401–406, 2011). In fact, an extremely simple linear relationship exists only if the predictor variables are conditionally independent given the target variable, in which case the contrasts, i.e., the differences of the weights, are equal to the logistic regression parameters. Thus, weights-of-evidence is the special case of logistic regression if the predictor variables are binary and conditionally independent given the target variable.     

Comprehensive seismic hazard assessment of Tripura and Mizoram states

Northeast India is one of the most highly seismically active regions in the world with more than seven earthquakes on an average per year of magnitude 5.0 and above. Reliable seismic hazard assessment could provide the necessary design inputs for earthquake resistant design of structures in this region. In this study, deterministic as well as probabilistic methods have been attempted for seismic hazard assessment of Tripura and Mizoram states at bedrock level condition. An updated earthquake catalogue was collected from various national and international seismological agencies for the period from 1731 to 2011. The homogenization, declustering and data completeness analysis of events have been carried out before hazard evaluation. Seismicity parameters have been estimated using G–R relationship for each source zone. Based on the seismicity, tectonic features and fault rupture mechanism, this region was divided into six major subzones. Region specific correlations were used for magnitude conversion for homogenization of earthquake size. Ground motion equations (Atkinson and Boore 2003; Gupta 2010) were validated with the observed PGA (peak ground acceleration) values before use in the hazard evaluation. In this study, the hazard is estimated using linear sources, identified in and around the study area. Results are presented in the form of PGA using both DSHA (deterministic seismic hazard analysis) and PSHA (probabilistic seismic hazard analysis) with 2 and 10% probability of exceedance in 50 years, and spectral acceleration (T = 0. 2 s, 1.0 s) for both the states (2% probability of exceedance in 50 years). The results are important to provide inputs for planning risk reduction strategies, for developing risk acceptance criteria and financial analysis for possible damages in the study area with a comprehensive analysis and higher resolution hazard mapping.     

Net earthquake hazard and elements at risk (NEaR) map creation for city of Istanbul via spatial multi-criteria decision analysis

The creation of earthquake hazard maps requires various datasets with selected attenuation relations. Based on the selected attenuation relation, the calculation time varies from half an hour to a couple of days. The length of time needed to create an earthquake hazard map also depends on the resolution of the resulting map. The time gets longer as the resolution of the resulting earthquake hazard map gets higher. The basic form of an attenuation relation requires complex calculation algorithms including geospatial information related to the region of interest. Nowadays, next-generation attenuation (NGA) models are introduced to generate more realistic earthquake hazard maps. However, the more complex the attenuation relation is, the longer time will be required to create a hazard map. This paper offers a new method to create high-resolution earthquake hazard maps, faster than using traditional attenuation relation methods, by using an analytic hierarchy process of spatial multi-criteria decision analysis and geographic information systems. This method has been generated and tested for the city of Istanbul. The resulting maps are compared with the earthquake hazard maps created for the city of Istanbul by using the NGA model of Boore and Atkinson (in Boore–Atkinson NGA ground motion relations for the geometric mean horizontal component of peak and spectral ground motion parameters (trans: Engineering Co, University of California B). Pacific Earthquake Engineering Research Center 2007). A second output of this paper is a map of the elements at risk (EaR) for the population and buildings of Istanbul, and the introduction of a new approach of net elements at risk (NEaR).     

Oxygen fugacity control in piston-cylinder experiments: a re-evaluation

Jakobsson (Contrib Miner Petrol 164(3):397–407, 2012) investigated a double capsule assembly for use in piston-cylinder experiments that would allow hydrous, high-temperature, and high-pressure experiments to be conducted under controlled oxygen fugacity conditions. Using a platinum outer capsule containing a metal oxide oxygen buffer (Ni–NiO or Co–CoO) and H₂O, with an inner gold–palladium capsule containing hydrous melt, this study was able to compare the oxygen fugacity imposed by the outer capsule oxygen buffer with an oxygen fugacity estimated by the AuPdFe ternary system calibrated by Barr and Grove (Contrib Miner Petrol 160(5):631–643, 2010). H₂O loss or gain, as well as iron loss to the capsule walls and carbon contamination, is often observed in piston-cylinder experiments and often go unexplained. Only a few have attempted to actually quantify various aspects of these changes (Brooker et al. in Am Miner 83(9–10):985–994, 1998; Truckenbrodt and Johannes in Am Miner 84:1333–1335, 1999). It was one of the goals of Jakobsson (Contrib Miner Petrol 164(3):397–407, 2012) to address these issues by using and testing the AuPdFe solution model of Barr and Grove (Contrib Miner Petrol 160(5):631–643, 2010), as well as to constrain the oxygen fugacity of the inner capsule. The oxygen fugacities of the analyzed melts were assumed to be equal to those of the solid Ni–NiO and Co–CoO buffers, which is incorrect since the melts are all undersaturated in H₂O and the oxygen fugacities should therefore be lower than that of the buffer by 2 log $a_{{{\text{H}}_{ 2} {\text{O}}}}$ .     

The 2012 Northern Italy Earthquakes: modelling human behaviour

The aim of this study was to investigate human behaviour during the 2012 Northern Italy Earthquakes. Furthermore, the current study used Kuligowski and Mileti’s (Fire Saf J 44:487–496, 2009) extended model and the Social Attachment Model as a framework to explain the behavioural responses. The study included 1839 participants who were affected by the earthquake occurred in the Emilia-Romagna region (Italy) on 20 May 2012. The most frequent behavioural responses during the earthquake were moving to another room of the house, escaping from home, and waiting in bed. According to Kuligowski and Mileti’s (Fire Saf J 44:487–496, 2009) extended model, perceived risk was associated with evacuation, and emergency preparedness was related to more efficient and effective responses during the earthquake. In line with the Social Attachment Model, affiliation behaviours were more frequent than flight behaviours, while, contrary to predictions, the social context did not influence emotional responses, evacuation behaviour, and search for protection.     

A phase-by-phase upstream scheme that converges to the vanishing capillarity solution for countercurrent two-phase flow in two-rock media

We discuss the convergence of the upstream phase-by-phase scheme (or upstream mobility scheme) towards the vanishing capillarity solution for immiscible incompressible two-phase flows in porous media made of several rock types. Troubles in the convergence were recently pointed out by Mishra and Jaffré (Comput. Geosci. 14, 105–124, 2010) and Tveit and Aavatsmark (Comput. Geosci. 16, 809–825, 2012). In this paper, we clarify the notion of vanishing capillarity solution, stressing the fact that the physically relevant notion of solution differs from the one inferred from the results of Kaasschieter (Comput. Geosci. 3, 23–48, 1999). In particular, we point out that the vanishing capillarity solution depends on the formally neglected capillary pressure curves, as it was recently proven in by Andreianov and Cancès (Comput. Geosci. 17, 551–572, 2013). Then, we propose a numerical procedure based on the hybridization of the interfaces that converges towards the vanishing capillarity solution. Numerical illustrations are provided.