The emplacement of intermediate volume ignimbrites: A case study from Roccamonfina Volcano,Southern Italy

A model is presented for the emplacement of intermediate volume ignimbrites based on a study of two 6 km³ volume ignimbrites on Roccamonfina Volcano, Italy. The model considers that the flows were slow moving, and quickly deflated from turbulent to non-turbulent conditions. Yield strength and density increased whereas fluidisation decreased with time and runout of the pyroclastic flows. In proximal locations, on the caldera rim, heterogeneous exposures including discontinuous lithic breccias, stratified and cross-stratified units interbedded with massive ignimbrite suggest deposition from turbulent flows. In medial locations thick, massive ignimbrite occurs associated with three types of co-ignimbrite lithic breccia which we interpret as being emplaced by non-turbulent flows. Multiple grading of different breccia/lithic concentration types within single flow units indicates that internal shear occurred producing overriding or overlapping of the rear of the flow onto the slower-moving front part. This overriding of different parts of non-turbulent pyroclastic flows could be caused by at least two different mechanisms: (1) changes in flow regime, such as hydraulic jumps that may occur at breaks in slope; and (2) periods of increased discharge rate, possibly associated with caldera collapse, producing fresh pulses of lithic-rich material that sheared onto the slower-moving part of the flow in front.We propose that ground surge deposits enriched in pumice compared with their associated ignimbrite probably formed by a flow separation mechanism from the top and front of the pyroclastic flow. These turbulent clouds moved ahead of the non-turbulent lower part of the flow to form stratified pumice-rich deposits. In distal regions well-developed coarse, often clast-supported, pumice concentrations zones and coarse intra-flow-unit lithic concentrations occur within the massive ignimbrite. We suggest that the flows were non-turbulent, possessed a relatively high yield strength and may have moved by plug flow prior to emplacement.     

Potential cauchy-poisson waves generated by submarine eruptions of kick 'em Jenny volcano

Kick 'em Jenny is the only known currently active submarine volcano in the Lesser Antilles. The volcano has erupted at least 10 times since first being discovered in 1939 and the summit has shoaled from a depth of 232 m in 1962 to its present-day depth of 150 m. Kick 'em Jenny is located in a province of explosive volcanism, has a known history of explosive eruptions and erupts magma of an explosive type. Future eruptions are likely to become increasingly more violent as the effect of the overlying water pressure becomes less. A preliminary study (Smith and Shepherd, 1993) suggests that Kick 'em Jenny is a prime candidate for tsunamigenic eruptions on a potentially hazardous scale, possibly affecting the whole of the eastern Caribbean region.The classic approach to problems of water waves generated by sudden disturbances of the free surface makes use of the Cauchy-Poisson-Lamb theory. A large number of theoretical developments to this theory have been made for specific forms of surface disturbance. A development by Unoki and Nakano (1953a, b) considers both two- and three-dimensional Cauchy-Poisson waves generated by finite initial elevations and impulses applied to a free surface of infinitely deep water. Unoki and Nakano's results compared well to the wave systems recorded following submarine eruptions of the Myojinsho Reef volcano in 1952–53.Given the similarity of the two situations, Unoki and Nakano's theory is applied to Kick 'em Jenny to provide estimates of potential Cauchy-Poisson wave heights throughout the eastern Caribbean for a range of eruption magnitudes. The results show that, although the waves generated are unlikely to pose much of a threat to the eastern Caribbean as a whole, they should be considered a hazard to the islands immediately adjacent to the volcano including Grenada, the Grenadines, and St Vincent.     

Measurements of peroxyacetylnitrate in the marine boundary layer over the Atlantic

The atmospheric concentration of peroxyacetylnitrate (PAN) was measured during a cruise of the R.S. Polarstern from Bremerhaven (Germany) to Rio Grande do Sul (Brazil) in September/ October 1988. The measurements were made in-situ by a combination of electron capture gaschromatography with a cryogenic preconcentration step. The theoretical lower limit of detection (3) was 0.4 ppt. The mixing ratios of PAN varied by more than three orders of magnitude from 2000 ppt in the English Channel to less than 0.4 ppt south of the Azores (38° N). South of 35° N, PAN levels were below the detection limit, except at 30–31° S off the eastern coast of South America. Here, PAN mixing ratios of 10 to 100 ppt were detected in continentally influenced air masses. Detectable levels of PAN were mostly observed in air masses of continental or high northern origin. Changes in the wind directions were usually associated with substantial changes in the PAN mixing ratios.     

Spatial and temporal variability of Holocene temperature in the North Atlantic region

The early Holocene climate of the North Atlantic region was influenced by two boundary conditions that were fundamentally different from the present: the presence of the decaying Laurentide Ice Sheet (LIS) and higher than present summer solar insolation. In order to assess spatial and temporal patterns of Holocene climate evolution across this region, we collated quantitative paleotemperature records at sub-millennial resolution and synthesized their temporal variability using principal components analysis (PCA). The analysis reveals considerable spatial variability, most notably in the time-transgressive expression of the Holocene thermal maximum (HTM). Most of the region, but especially areas peripheral to the Labrador Sea and hence closest to the locus of LIS disintegration, experienced maximum Holocene temperatures that lagged peak summer insolation by 1000-3000 years. Many sites from the northeastern North Atlantic sector, including the Nordic Seas and Scandinavia, either warmed in phase with maximum summer insolation (11,000-9000 years ago) or were less strongly lagged than the Baffin Bay-Labrador Sea region. These spatially complex patterns of Holocene climate development, which are defined by the PCA, resulted from the interplay between final decay of the LIS and solar insolation forcing.     

Abrupt climate change: An alternative view

Hypotheses and inferences concerning the nature of abrupt climate change, exemplified by the Dansgaard-Oeschger (D-O) events, are reviewed. There is little concrete evidence that these events are more than a regional Greenland phenomenon. The partial coherence of ice core δ¹⁸O and CH₄ is a possible exception. Claims, however, of D-O presence in most remote locations cannot be distinguished from the hypothesis that many regions are just exhibiting temporal variability in climate proxies with approximately similar frequency content. Further suggestions that D-O events in Greenland are generated by shifts in the North Atlantic ocean circulation seem highly implausible, given the weak contribution of the high latitude ocean to the meridional flux of heat. A more likely scenario is that changes in the ocean circulation are a consequence of wind shifts. The disappearance of D-O events in the Holocene coincides with the disappearance also of the Laurentide and Fennoscandian ice sheets. It is thus suggested that D-O events are a consequence of interactions of the windfield with the continental ice sheets and that better understanding of the wind field in the glacial periods is the highest priority. Wind fields are capable of great volatility and very rapid global-scale teleconnections, and they are efficient generators of oceanic circulation changes and (more speculatively) of multiple states relative to great ice sheets. Connection of D-O events to the possibility of modern abrupt climate change rests on a very weak chain of assumptions.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.