Land–Water Boundary Treatment for a Tsunami Model With Dimensional Splitting

The Method of Splitting Tsunamis (MOST) model adapted by National Oceanic and Atmospheric Administration (NOAA) for tsunami forecasting operations is praised for its computational efficiency, associated with the use of splitting technique. It will be shown, however, that splitting the computations between \(x\) and \(y\) directions results in specific sensitivity to the treatment of land–water boundary. Slight modification to the reflective boundary condition in MOST caused an appreciable difference in the results. This is demonstrated with simulations of the Tohoku-2011 tsunami from the source earthquake to Monterey Bay, California, and in southeast Alaska, followed by comparison with tide gage records. In the first case, the better representation of later waves (reflected from the coasts) by the modified model in a Pacific-wide simulation resulted in twice as long match between simulated and observed tsunami time histories at Monterey gage. In the second case, the modified model was able to propagate the tsunami wave and approach gage records at locations within narrow channels (Juneau, Ketchikan), to where MOST had difficulty propagating the wave. The modification was extended to include inundation computation. The resulting inundation algorithm (Cliffs) has been tested with the complete set of NOAA-recommended benchmark problems focused on inundation. The solutions are compared to the MOST solutions obtained with the version of the MOST model benchmarked for the National Tsunami Hazard Mitigation Program in 2011. In two tests, Cliffs and MOST results are very close, and in another two tests, the results are somewhat different. Very different regimes of generation/disposal of water by Cliffs and MOST inundation algorithms, which supposedly affected the benchmarking results, have been discussed.     

Expressions for the Global Gravimetric Moho Modeling in Spectral Domain

We apply a newly developed numerical method to improve the Moho geometry by the implementation of gravity data. This method utilizes expressions for the gravimetric forward and inverse modeling derived in a frequency domain. Methods for a spectral analysis and synthesis of the gravity field and crust density structures are applied in the gravimetric forward modeling of the consolidated crust-stripped gravity disturbances, which have a maximum correlation with the (a priori) Moho model. These gravity disturbances are obtained from the Earth’s gravity disturbances after applying the topographic and stripping gravity corrections of major known anomalous crust density structures; in the absence of a global mantle model, mantle density heterogeneities are disregarded. The isostatic scheme applied is based on a complete compensation of the crust relative to the upper mantle density. The functional relation is established between the (unknown) Moho depths and the complete crust-stripped isostatic gravity disturbances, which according to the adopted isostatic scheme have (theoretically) a minimum correlation with the Moho geometry. The system of observation equations, which describes the relation between spherical functions of the isostatic gravity field and the Moho geometry, is defined by means of a linearized Fredholm integral equation of the first kind. The Moho depths are determined based on solving the gravimetric inverse problem. The regularization is applied to stabilize the ill-posed solution. This numerical procedure is utilized to determine the Moho depths globally. The gravimetric result is presented and compared with the seismic Moho model. Our gravimetric result has a relatively good agreement with the CRUST2.0 Moho model by means of the RMS of differences (of 3.5 km). However, the gravimetric solution has a systematic bias. We explain this bias between the gravimetric and seismic Moho models by the unmodelled mantle heterogeneities and uncertainties in the CRUST2.0 global crustal model.     

Evidence for an Independent 26-s Microseismic Source near the Vanuatu Islands

The 26 s peak in the ambient seismic noise spectrum is persistently excited and observed at stations globally. Using noise cross-correlation functions (NCFs), the location suggests that the source could be situated in the Gulf of Guinea and Fiji Basin. However, the Fiji Basin was proposed to be the mirror site (near antipode) of the Gulf of Guinea source instead of an independent source, assuming that the surface waves more efficiently propagate along the major-arc paths of oceanic movements. To investigate the propagation of the Rayleigh waves along continental and oceanic paths, we analyzed the surface wave data recorded from an earthquake near the Gulf of Guinea and found that Rayleigh waves travel along continental minor-arc paths more efficiently than along oceanic major-arc paths. We then located the source in the western Pacific Ocean from group velocities measured with earthquake data by using the travel time misfit in NCFs after calibration and concluded that the source is in the Vanuatu Islands. Moreover, the temporal variation of the 26 s microseismic peak observed in the western Pacific seismic stations is very different from that in stations near the Gulf of Guinea, which suggests that they are excited by independent sources. Therefore, the Vanuatu source should be an independent microseismic source. As it is close to volcanoes in the Vanuatu islands, the Pacific 26 s microseismic source might be excited by magmatic processes, which are also responsible for very-long-period volcanic tremors.     

PALEOMAGIA: A PHP/MYSQL database of the Precambrian paleomagnetic data

Most paleomagnetic applications require a precise, rationally organized and up-todate catalogue or database of paleomagnetic results worldwide. These include reconstructions of continents, calculations of the Apparent Polar Wander Paths (APWPs) or paleolatitude drift curves, testing the Geocentric Axial Dipole (GAD) model, studies of geomagnetic paleosecular variation or reversal asymmetries, comparison of coeval results obtained from different types of rocks, estimation of inclination shallowing in sedimentary rocks and understanding the delay in remanence acquisition caused by slow cooling in large intrusions. For this purpose, various databases, such as the Global Paleomagnetic Database (GPMDB), and the Magnetics Information Consortium Database (MagIC) have been generated. This paper presents a new relational database (PALEOMAGIA) where 3278 entries of Precambrian data have been split geographically, sorted according to age and rock types and ranked using a revised version of the Van der Voo grading scheme. The latest geochronologic information is included wherever available. Significant effort has been put to the retrieval and archiving of data published in the last decade, which are virtually nonexistent in GPMDB. Here we present the database and its browser-based user interface from a scientific and a technical point of view.     

Assessment of genotoxic and humoral immune system alterations in silica exposed workers from pottery industries in South India

Chemicals, especially silica, have been suspected to cause genetic alterations in pottery industry workers. The present study aims to analyze the frequency of chromosomal aberrations (CA), micronucleus (MN) and DNA damage (comet assay) in the peripheral blood lymphocytes and the immunological alterations workers chronically exposed to silica and in control subjects. In our study, 50 silica exposed workers and 35 control subjects were recruited and silica level was measured by respirable dust and respirable quartz concentrations of badge dosimeter. The serum immunoglobulins (IgM, IgG, and IgA) of all the subjects were measured by using ELISA method. The individuals exposed to silica have a significant increase in the frequency of CA, MN and the total DNA damage (p < 0.05). Immunoglobulin elevation on silica exposed workers was statistically significant (p < 0.05) on comparison with their respective controls. Investigation of the smoking and alcohol habitats coupled with silica exposure in exposed and control subjects represents alcohol consumption and smoking as additional risk factors and must be avoided. Multiple linear regression analysis obtained for CA, MN and comet assay confirm these tests as biomarkers for silica exposed pottery workers. Some confounding factors also showed significant influence on exposed subjects. These results indicate the mutagenic risk in the working environment has a high probability of association with the silica dust exposure in pottery industries. Nevertheless, the present study will create awareness and public concern not only among the silica exposed workers but also to the welfare of their progeny.     

FEM × DEM modelling of cohesive granular materials: Numerical homogenisation and multi-scale simulations

The article presents a multi-scale modelling approach of cohesive granular materials, its numerical implementation and its results. At microscopic level, Discrete Element Method (DEM) is used to model dense grains packing. At the macroscopic level, the numerical solution is obtained by a Finite Element Method (FEM). In order to bridge the micro- and macro-scales, the concept of Representative Elementary Volume (REV) is applied, in which the average REV stress and the consistent tangent operators are obtained in each macroscopic integration point as the results of DEM’s simulation. In this way, the numerical constitutive law is determined through the detailed modelling of the microstructure, taking into account the nature of granular materials. We first elaborate the principle of the computation homogenisation (FEM × DEM), then demonstrate the features of our multiscale computation in terms of a biaxial compression test. Macroscopic strain location is observed and discussed.     

The Beni-Ilmane (North-Central Algeria) Earthquake Sequence of May 2010

Starting on 14 May 2010 and lasting several months, the village of Beni-Ilmane (Msila District, North-Central Algeria) and its surroundings were struck by an important seismic crisis marked by three successive moderate shocks (5.0 ≤ M _d ≤ 5.2). This sequence of events caused severe damage in the Beni-Ilmane village and in the epicentral area. The poor quality of masonry construction and the cumulative effects of the large number of aftershock events played a key role in the destruction. To follow this earthquake sequence, 11 temporary seismic stations, in addition to the permanent stations of the Algerian seismic network, were deployed in the region. A representative set of well located aftershocks in the period of maximum activity (lasting 18 days) were selected. The horizontal distribution of the aftershocks shows two main earthquake clusters located near Beni-Ilmane village, one cluster oriented E–W and the other oriented NNE–SSW, crossing the first cluster at its eastern tip. The aftershocks distribution suggests that the three main shocks ruptured two distinct and adjacent fault segments of about 8 km length. The focal mechanisms of the first and third events, located in the NNE–SSW cluster, show near-vertical left-lateral strike-slip fault planes. In the second cluster, oriented E–W, focal mechanisms show a high-angle reverse fault. A field survey, initiated immediately after the first main shock, identified surface fissures generated by the three largest events in the sequence. The fissures, concentrated in a narrow area at the western termination of the NE–SW Jebel Choukchot anticline (location of Beni-Ilmane village), showed several orientations which were mainly related to gravity instabilities. The 2010 Beni-Ilmane earthquake sequence, located in the Bibans–Hodna Mountains transition zone, demonstrates that the Tellian Atlas–High Plateaus border region is an active seismic zone marked by moderate and possibly strong earthquakes; thus, a reevaluation of the seismic hazard in the region is needed.     

Recent developments of the Middle East catalog

This article summarizes a recent study in the framework of the Global Earth model (GEM) and the Earthquake Model of the Middle East (EMME) project to establish the new catalog of seismicity for the Middle East, using all historical (pre-1900), early and modern instrumental events up to 2006. According to different seismicity, which depends on geophysical, geological, tectonic, and seismicity data, this region is subdivided to nine subregions, consisting of Alborz–Azerbaijan, Afghanistan–Pakistan, Saudi Arabia, Caucasus, Central Iran, Kopeh–Dagh, Makran, Zagros, and Turkey (Eastern Anatolia; after 30° E). After omitting the duplicate events, aftershocks, and foreshocks by using the Gruenthal method, and uniform all magnitude to Mw scale, 28,244 main events remain for the new catalog of Middle East from 1250 B.C. through 2006. The magnitude of completeness (Mc) was determined as 4.9 for five out of nine subregions, where the least values of Mc were found to be 4.2. The threshold of Mc is around 5.5, 5.0, 4.5, and 4.0, for the time after 1950, 1963, 1975, and 2000, respectively. The average of teleseismic depths in all regions is less than 15 km. Totally, majority of depth for Kopeh–Dagh and Central Iran, Zagros, and Alborz–Azerbaijan, approximately, is 15, 13, and 11 km and for Afghanistan–Pakistan, Caucasus, Makran, Turkey (after 30° E), and Saudi Arabia is about 9 km.     

Dynamic Site Characterization and Correlation of Shear Wave Velocity with Standard Penetration Test ‘N’ Values for the City of Agartala,Tripura State,India

Seismic site characterization is the basic requirement for seismic microzonation and site response studies of an area. Site characterization helps to gauge the average dynamic properties of soil deposits and thus helps to evaluate the surface level response. This paper presents a seismic site characterization of Agartala city, the capital of Tripura state, in the northeast of India. Seismically, Agartala city is situated in the Bengal Basin zone which is classified as a highly active seismic zone, assigned by Indian seismic code BIS-1893, Indian Standard Criteria for Earthquake Resistant Design of Structures, Part-1 General Provisions and Buildings. According to the Bureau of Indian Standards, New Delhi (2002), it is the highest seismic level (zone-V) in the country. The city is very close to the Sylhet fault (Bangladesh) where two major earthquakes (M _w > 7) have occurred in the past and affected severely this city and the whole of northeast India. In order to perform site response evaluation, a series of geophysical tests at 27 locations were conducted using the multichannel analysis of surface waves (MASW) technique, which is an advanced method for obtaining shear wave velocity (V _s) profiles from in situ measurements. Similarly, standard penetration test (SPT-N) bore log data sets have been obtained from the Urban Development Department, Govt. of Tripura. In the collected data sets, out of 50 bore logs, 27 were selected which are close to the MASW test locations and used for further study. Both the data sets (V _s profiles with depth and SPT-N bore log profiles) have been used to calculate the average shear wave velocity (V _s30) and average SPT-N values for the upper 30 m depth of the subsurface soil profiles. These were used for site classification of the study area recommended by the National Earthquake Hazard Reduction Program (NEHRP) manual. The average V _s30 and SPT-N classified the study area as seismic site class D and E categories, indicating that the city is susceptible to site effects and liquefaction. Further, the different data set combinations between V _s and SPT-N (corrected and uncorrected) values have been used to develop site-specific correlation equations by statistical regression, as ‘V _s’ is a function of SPT-N value (corrected and uncorrected), considered with or without depth. However, after considering the data set pairs, a probabilistic approach has also been presented to develop a correlation using a quantile–quantile (Q–Q) plot. A comparison has also been made with the well known published correlations (for all soils) available in the literature. The present correlations closely agree with the other equations, but, comparatively, the correlation of shear wave velocity with the variation of depth and uncorrected SPT-N values provides a more suitable predicting model. Also the Q–Q plot agrees with all the other equations. In the absence of in situ measurements, the present correlations could be used to measure V _s profiles of the study area for site response studies.     

Determination of Dips and Depths of Near Surface Layers by Radon Transform

The refracted arrivals on seismic shot records have long been recognized as an efficient tool for the computation of detailed near-surface information. In this paper, a new concept of refraction static, which is based on the Radon transform and avoids the tedious process of picking first arrival times, is proposed. This method is particularly suitable when a rough near-surface problem necessitates the utilization of numerous shallow refraction data for the one reflector case. Quasi-linearity of refractors and a constant velocity medium are assumed within the shooting range. Synthetic and real cases have been tested to evaluate the performance of the method. The result is revealed to be satisfactory. Comparison of the synthetic model with the results obtained through the Radon transform reveals a very good accuracy for the proposed method.     

Explosive felsic volcanism on El Hierro (Canary Islands)

The Canary Islands consist of seven basaltic shield volcanoes whose submerged portion is much more voluminous than the subaerial part of each island. Like so many other volcanic oceanic islands, the indicative deposits of explosive felsic volcanism are not a common feature on the Canary archipelago. Hitherto, they have only been documented from the central islands of Gran Canaria and Tenerife, which are the largest volcanic complexes of the islands. On the other Canary Islands, the presence of felsic rocks is mostly restricted to intrusions and a few lava flows, generally within the succession in the oldest parts of individual islands. In this paper, we present a detailed stratigraphic, lithological and sedimentological study of a significant felsic pumice deposit on the island of El Hierro, referred here as the Malpaso Member, which represents the only explosive episode of felsic volcanism found on the Canary Islands (outside of Gran Canaria and Tenerife). The products of the eruption indicate a single eruptive event and cover an area of about 15 km². This work provides a detailed stratigraphic and chronological framework for El Hierro, and four subunits are identified within the member on the basis of lithological and granulometric characteristics. The results of this study demonstrate the importance of an explosive eruption in a setting where the activity is typified by effusive basaltic events. Given the style and the spatial distribution of the Malpaso eruption and its products, a future event with similar characteristics could have a serious impact on the population, infrastructure and economy of the island of El Hierro.     

The gravity field and crustal thickness of Venus

The gravity and topography of Venus obtained from observations of the Magellan mission, as well as the gravity and topography from our numerical mantle convection model, are discussed in this paper. We used the hypothesis that the geoid of degrees 2–40 is produced by sublithospheric mantle density anomalies that are associated with dynamical process within the mantle. We obtained the model dynamical admittance(the geoid topography ratio based on a convection model) by a numerical simulation of the Venusian mantle convection, and used it to correct the dynamical effect in the calculation of crustal thickness. After deducting the dynamical effect, the thickness of the Venusian crust is presented. The results show that the gravity and topography are strongly correlated with the Venusian mantle convection and the Venusian crust has a significant influence on the topography. The Venusian crustal thickness varies from 28 to 70 km. Ishtar Terra, and Ovda Regio and Thetis Regio in western Aphrodite Terra have the highest crustal thickness(larger than 50 km). The high topography of these areas is thought to be supported by crustal compensation and our results are consistent with the hypothesis that these areas are remnants of ancient continents. The crustal thickness in the Beta, Themis, Dione, Eistla, Bell, and Lada regiones is thin and shows less correlation with the topography, especially in the Atla and Imdr regiones in the eastern part of Aphrodite Terra. This is consistent with the hypothesis that these highlands are mainly supported by mantle plumes. Compared with the crustal thickness calculated with the dynamical effect, our results are more consistent with the crust evolution and internal dynamical process of Venus.     

Indian Summer Monsoon Rainfall Characteristics During Contrasting Monsoon Years

The present paper presents a diagnostic study of two recent monsoon years, of which one is dry monsoon year (2009) and the other is wet monsoon year (2010). The study utilized the IMD gridded rainfall data set in addition to the Reynolds SST, NCEP-NCAR reanalysis wind and temperature products, and NOAA OLR. The study revealed that the months July and August are the most crucial months to decide whether the ISMR is wet or dry. However, during July 2009, most of the Indian subcontinent received more than 60 % in the central and western coastal regions. In a wet monsoon year, about 35–45 % of rainfall is contributed during June and July in most parts of India. During these years, the influence of features in the Pacific Ocean played vital role on the Indian summer monsoon rainfall. During 2009, Pacific SST was above normal in nino regions, characteristic of the El Nino structure; however, during 2010, the nino regions were clearly below normal temperature, indicating the La Nina pattern. The associated atmospheric general circulation through equatorial Walker and regional Hadley circulation modulates the tropospheric temperature, and hence the organized convective cloud bands. These cloud bands show different characteristics in northward propagation during dry and wet years of ISMR. During a dry year, the propagation speed and magnitudes are considerably higher than during a wet monsoon year.     

Microseismic monitoring of columnar jointed basalt fracture activity: a trial at the Baihetan Hydropower Station,China

Severe stress release has occurred to the surrounding rocks of the typically columnar jointed basalt after excavation at the Baihetan Hydropower Station, Jinsha River, China, where cracking, collapse, and other types of failure may take place occasionally due to relaxation fracture. In order to understand the relaxation fracture characteristics of the columnar jointed basalt in the entire excavation process at the diversion tunnel of the Baihetan Hydropower Station, real-time microseismic monitoring tests were performed. First, the applicability of a geophone and accelerometer was analyzed in the columnar jointed basalt tunnel, and the results show that the accelerometer was more applicable to the cracking monitoring of the columnar jointed basalt. Next, the waveform characteristics of the microseismic signals were analyzed, and the microseismic signals were identified as follows: rock fracture signal, drilling signal, electrical signal, heavy vehicle passing signal, and blast signal. Then, the attenuation characteristics of the microseismic signals in the columnar jointed basalt tunnel were studied, as well as the types and characteristics of the columnar jointed basalt fracture. Finally, location analysis was conducted on the strong rock fracture events, in which four or more sensors were triggered, to obtain the temporal and spatial evolution characteristics and laws of the columnar jointed basalt relaxation fracture after excavation. The test results are not only of important reference value to the excavation and support of diversion tunnel at the Baihetan Hydropower Station, but also of great referential significance and value to the conduction of similar tests.     

近地层能量闭合度对陆面过程模式影响

大量近地层观测试验表明,利用涡动相关法观测的湍流通量小于近地层可利用能量,即近地层能量是不闭合的,这种不闭合度一般为20%甚至更高.而陆面过程模式是基于地气间能量平衡建立,并且模式中的湍流边界层参数化方案通常根据实际观测的湍流通量来确定,因此能量不闭合必将对陆面过程模式造成一定的影响.本文利用2007年春季SACOL站的近地层观测资料,依据能量守恒将能量不闭合中的残余能量通过波文比分配到观测的湍流通量中,即修正涡动相关法观测的湍流通量使得近地层能量达到平衡;之后分别利用观测和修正的湍流通量,建立了能量不闭合和闭合情形下的湍流参数化方案,借助陆面过程模式SHAW,通过数值模拟和对比分析方法考察近地层能量闭合度对陆面过程模式的影响.研究结果表明近地层能量闭合对陆面过程模式有显著的影响:在陆面过程数值模拟中,当应用近地层能量不闭合的湍流通量形成的湍流参数化方案时,陆面过程模式会明显高估地表长波辐射及土壤温度;但当应用修正湍流通量使得近地层能量达到闭合形成的湍流参数化方案后,在不改变任何地表土壤物理生化属性的情况下,陆面过程模式能较好地模拟地表长波辐射和土壤温度.     

On the ability of Moho reflections to affect the ground motion in northeastern Italy: a case study of the 2012 Emilia seismic sequence

It has been observed that post-critically reflected S-waves and multiples from the Moho discontinuity could play a relevant role on the ground motion due to medium to strong size earthquakes away from the source. Although some studies investigated the correlation between the Moho reflections amplitudes and the damage in the far field, little attention was given to the frequency content of these specific phases and their scaling with magnitude. The 2012 Emilia seismic sequence in northern Italy, recorded by velocimetric and accelerometric networks, is here exploited to investigate Moho reflections and multiples (SmSM). A single station method for group velocity-period estimation, based on the multiple filter technique, is applied to strong motion data to detect SmSM. Amplitude and frequency scaling with magnitude is defined for earthquakes from \(\hbox {Mw}=3.9\) to \(\hbox {Mw}=5.9\) . Finally, the ability of SmSM to affect the ground motion for a maximum credible earthquake within the Po plain is investigated by extrapolating observed engineering parameters. Data analysis shows that high amplitude SmSM can be recognized within the Po plain, and at the boundaries between the Po plain and the Alpine chain, at epicentral distances larger than 80 km, in the period range from 0.25 to 3 s and in the group velocity window from about 2.6 to 3.2 km/s. 5 % damped pseudo-spectral accelerations at different periods (0.3, 1.0 and 2.0 s), and Housner intensities, are obtained from data characterized by large amplitude SmSM. A scaling relationship for both pseudo-spectral accelerations and Housner intensities is found for the earthquakes of the 2012 Emilia seismic sequence. \(\hbox {I}_{\mathrm{MCS}}\) from VII to VIII is estimated, as a result of SmSM amplitude enhancement, at about 100 km for a maximum credible earthquake ( \(\hbox {Mw}=6.7\) ) in the Po plain, showing that moderate to high damage cloud be caused by these specific phases.     

Estimation of key surface parameters in semi-arid region and their impacts on improvement of surface fluxes simulation

Uncertainties in some key parameters in land surface models severely restrict the improvement of model capacity for successful simulation of surface-atmosphere interaction. These key parameters are related to soil moisture and heat transfer and physical processes in the vegetation canopy as well as other important aerodynamic processes. In the present study, measurements of surface-atmosphere interaction at two observation stations that are located in the typical semi-arid region of China, Tongyu Station in Jilin Province and Yuzhong Station in Gansu Province, are combined with the planetary boundary layer theory to estimate the value of two key aerodynamic parameters, i.e., surface roughness length z_0m and excess resistance κB^-1. Multiple parameterization schemes have been used in the study to obtain values for surface roughness length and excess resistance κB^-1 at the two stations. Results indicate that z_0m has distinct seasonal and inter-annual variability. For the type of surface with low-height vegetation, there is a large difference between the default value of z_0m in the land surface model and that obtained from this study. κB^-1 demonstrates a significant diurnal variation and seasonal variability. Using the modified scheme for the estimation of z_0m and κB^-1 in the land surface model, it is found that simulations of sensible heat flux over the semi-arid region have been greatly improved. These results suggest that it is necessary to further evaluate the default values of various parameters used in land surface models based on field measurements. The approach to combine field measurements with atmospheric boundary layer theory to retrieve realistic values for key parameters in land surface models presents a great potential in the improvement of modeling studies of surface-atmosphere interaction.     

Numerical Simulation of the Sensitivity of Summer Monsoon Circulation and Rainfall over India to Land Surface Processes

—The influence of soil moisture and vegetation variation on simulation of monsoon circulation and rainfall is investigated. For this purpose a simple land surface parameterization scheme is incorporated in a three-dimensional regional high resolution nested grid atmospheric model. Based on the land surface parameterization scheme, latent heat and sensible heat fluxes are explicitly estimated over the entire domain of the model. Two sensitivity studies are conducted; one with bare dry soil conditions (no latent heat flux from land surface) and the other with realistic representation of the land surface parameters such as soil moisture, vegetation cover and landuse patterns in the numerical simulation. The sensitivity of main monsoon features such as Somali jet, monsoon trough and tropical easterly jet to land surface processes are discussed.¶Results suggest the necessity of including a detailed land surface parameterization in the realistic short-range weather numerical predictions. An enhanced short-range prediction of hydrological cycle including precipitation was produced by the model, with land surface processes parameterized. This parameterization appears to simulate all the main circulation features associated with the summer monsoon in a realistic manner.     

Erosion and aggradation on persistently active volcanoes—a case study from Semeru Volcano,Indonesia

Erosion processes on active volcanoes in humid climates result in some of the highest sediment yields on Earth. Episodic sediment yields after large eruptions have been evaluated, but not the long-term and continuous patterns on persistently active volcanoes. We have used high-spatial resolution satellite imagery and DEMs/DSMs along with field-based geologic mapping to assess accurately sediment budgets for the active Semeru Volcano in Java, Indonesia. Patterns of aggradation and degradation on Semeru differ from that of other active volcanoes because (1) both episodic pyroclastic density currents (PDC) and continuous supplies of tephra generate pulses of sediment, (2) sediment is transferred via cycles of aggradation and degradation that continue for >15 years in river channels after each PDC-producing eruption, and (3) rain-triggered lahars remove much greater material than fluvial transport during long, intense rainfall events. The geomorphic response of two of Semeru’s rivers to volcanic sediment migration indicates that (1) each river experiences alternating aggradation and degradation cycles following PDC-producing eruptions and (2) spatial patterns of sediment transfer are governed by geomorphic characteristics of the river reaches. Usually high degradation in the steep source reach is followed by a long bypassing middle reach. Aggradation predominates in the depositional reaches further down valley on the ring plain. Average sediment yields (10³–10⁵ t/km²/year) at persistently active volcanoes are two to three orders of magnitude lower than sediment yields after large and infrequent eruptions, but the continuous and steady sediment transfer in rivers removes more sediment on a mid-term (10 years) to long-term (30 years) basis. In contrast to the trend observed on composite cones after large and infrequent eruptions, decay of sediment yields is not exponential and river channels do not fully recover at steadily active volcanoes as episodic inputs from BAF eruptions, superimposed on the background remobilization of daily tephra, have a greater cumulative effect.     

Volcanic ash hazard in the Central Mediterranean assessed from geological data

Volcanic ash produced during explosive eruptions can have very severe impacts on modern technological societies. Here, we use reconstructed patterns of fine ash dispersal recorded in terrestrial and marine geological archives to assess volcanic ash hazards. The ash-dispersal maps from nine Holocene explosive eruptions of Italian volcanoes have been used to construct frequency maps of distal ash deposition over a wide area, which encompasses central and southern Italy, the Adriatic and Tyrrhenian seas and the Balkans. The maps are presented as two cumulative-thickness isopach maps, one for nine eruptions from different volcanoes and one for six eruptions from Somma-Vesuvius. These maps represent the first use of distal ash layers to construct volcanic hazard maps, and the proposed methodology is easily applicable to other volcanic areas worldwide.