scholarly journals Origin of calderas: discriminating between collapses and explosions

2017 ◽  
Vol 59 (6) ◽  
Author(s):  
Izumi Yokoyama

<p>Origins of calderas may differ according to their subsurface structure that may be characterized by high or low density deposits that may be observed as high or low gravity anomalies, respectively. In the Introduction, the pioneering work of Fouqué[1879] on Santorini caldera is referred to in relation to definition of calderas. First, our discussion is focused on four calderas that were seen forming during the period from 1815 (the Tambora eruption) to 1991 (the Pinatubo eruption). Coincidently, these four calderas are all low-gravity-anomaly type. Their formation processes and subsurface structure are summarized by the existing data analyzed by various authors. These results are confirmed by results of drillings at some other calderas. Then, caldera formation of both types is discussed: High-gravity-anomaly-type calderas are expected to originate from subsidence of high-density ejecta into the summit magma reservoir. On the calderas of this type, the genetic eruption<span style="text-decoration: line-through;">s</span> believed to be accompanied by subsidences were not actually observed, and consequently three examples are mentioned only briefly. The low-gravity-anomaly-type calderas are discussed from standpoint of both the models of collapses and explosions. It is also emphasized that dynamic pressure ofexplosions is an important factor in the caldera formation, not only volume of the ejecta. To confirm the possibility that volcanic ejecta and edifices collapse into magma reservoirs, we discuss stress propagation from a depleted reservoir upward towards the Earth surface. Formation mechanisms of large calderas of this type are speculated; large calderas measuring about 20 km across may develop by successive merging of component calderas over a long period of times. A Kamchatka caldera under enlargement during the Holocene period is interpreted by successive merging of five component calderas.</p>

2015 ◽  
Vol 58 (2) ◽  
Author(s):  
Izumi Yokoyama

<p>Firstly the volume of pyroclastic ejecta during the 1883 eruption of Krakatau is re-examined. To revise the volume of flow deposits, the author basically follows Verbeek’s observation while to estimate the fall deposits, as the last resort, the author assumes that volume ratios fall / flow are common to similar caldera eruptions, and the ratios determined by the caldera- forming eruptions of Novarupta and Pinatubo are applied to the Krakatau eruption. Verbeek’s estimation of the total volume of ejecta, 12 km<sup>3</sup> is revised to 19 km<sup>3</sup>. This is significantly different from the volume of disrupted volcano edifice, 8 km<sup>3</sup>. Such a result does not support the predecessors’ hypothesis that calderas are formed by collapses of volcano edifices into magma reservoirs in replacement of the total ejecta. Through the discussion on the volume estimation of volcanic ejecta on and around Krakatau, the author recognizes that such estimation should be originally very difficult to attain enough accuracy. Much importance of “caldera deposits” to post-eruption settlements of the ejecta is emphasized. In relation to caldera formation, mechanical stability of a cavity in the crust is discussed. Lastly, upon the basis of subsurface structure, especially caldera deposits, a structural image of Krakatau caldera is presented.</p>


2014 ◽  
Vol 23 (2) ◽  
pp. 213-229 ◽  
Author(s):  
Cangqi Zhou ◽  
Qianchuan Zhao

AbstractMining time series data is of great significance in various areas. To efficiently find representative patterns in these data, this article focuses on the definition of a valid dissimilarity measure and the acceleration of partitioning clustering, a common group of techniques used to discover typical shapes of time series. Dissimilarity measure is a crucial component in clustering. It is required, by some particular applications, to be invariant to specific transformations. The rationale for using the angle between two time series to define a dissimilarity is analyzed. Moreover, our proposed measure satisfies the triangle inequality with specific restrictions. This property can be employed to accelerate clustering. An integrated algorithm is proposed. The experiments show that angle-based dissimilarity captures the essence of time series patterns that are invariant to amplitude scaling. In addition, the accelerated algorithm outperforms the standard one as redundancies are pruned. Our approach has been applied to discover typical patterns of information diffusion in an online social network. Analyses revealed the formation mechanisms of different patterns.


2021 ◽  
Author(s):  
◽  
Alistair Stronach

<p><b>New Zealand’s capital city of Wellington lies in an area of high seismic risk, which is further increased by the sedimentary basin beneath the Central Business District (CBD). Ground motion data and damage patterns from the 2013 Cook Strait and 2016 Kaikōura earthquakes indicate that two- and three-dimensional amplification effects due to the Wellington sedimentary basin may be significant. These effects are not currently accounted for in the New Zealand Building Code. In order for this to be done, three-dimensional simulations of earthquake shaking need to be undertaken, which requires detailed knowledge of basin geometry. This is currently lacking, primarily because of a dearth of deep boreholes in the CBD area, particularly in Thorndon and Pipitea where sediment depths are estimated to be greatest.</b></p> <p>A new basin depth map for the Wellington CBD has been created by conducting a gravity survey using a modern Scintrex CG-6 gravity meter. Across the study area, 519 new high precision gravity measurements were made and a residual anomaly map created, showing a maximum amplitude anomaly of -6.2 mGal with uncertainties better than ±0.1 mGal. Thirteen two-dimensional geological profiles were modelled to fit the anomalies, then combined with existing borehole constraints to construct the basin depth map. </p> <p>Results indicate on average greater depths than in existing models, particularly in Pipitea where depths are interpreted to be as great as 450 m, a difference of 250 m. Within 1 km of shore depths are interpreted to increase further, to 600 m. The recently discovered basin bounding Aotea Fault is resolved in the gravity data, where the basement is offset by up to 13 m, gravity anomaly gradients up to 8 mGal/km are observed, and possible multiple fault strands identified. A secondary strand of the Wellington Fault is also identified in the north of Pipitea, where gravity anomaly gradients up to 18 mGal/km are observed.</p>


2021 ◽  
Vol 9 ◽  
Author(s):  
Shanwei Liu ◽  
Yinlong Li ◽  
Qinting Sun ◽  
Jianhua Wan ◽  
Yue Jiao ◽  
...  

The purpose of this paper is to analyze the influence of satellite altimetry data accuracy on the marine gravity anomaly accuracy. The data of 12 altimetry satellites in the research area (5°N–23°N, 105°E–118°E) were selected. These data were classified into three groups: A, B, and C, according to the track density, the accuracy of the altimetry satellites, and the differences of self-crossover. Group A contains CryoSat-2, group B includes Geosat, ERS-1, ERS-2, and Envisat, and group C comprises T/P, Jason-1/2/3, HY-2A, SARAL, and Sentinel-3A. In Experiment I, the 5′×5′ marine gravity anomalies were obtained based on the data of groups A, B, and C, respectively. Compared with the shipborne gravity data, the root mean square error (RMSE) of groups A, B, and C was 4.59 mGal, 4.61 mGal, and 4.51 mGal, respectively. The results show that high-precision satellite altimetry data can improve the calculation accuracy of gravity anomaly, and the single satellite CryoSat-2 enables achieving the same effect of multi-satellite joint processing. In Experiment II, the 2′×2′ marine gravity anomalies were acquired based on the data of groups A, A + B, and A + C, respectively. The root mean square error of the above three groups was, respectively, 4.29 mGal, 4.30 mGal, and 4.21 mGal, and the outcomes show that when the spatial resolution is satisfied, adding redundant low-precision altimetry data will add pressure to the calculation of marine gravity anomalies and will not improve the accuracy. An effective combination of multi-satellite data can improve the accuracy and spatial resolution of the marine gravity anomaly inversion.


2014 ◽  
Vol 57 (5) ◽  
Author(s):  
Izumi Yokoyama

<p> </p><p>One of the differences between volcanic craters and calderas is that the latter bottoms are flatways filled with caldera deposit with lower density in comparison to country rocks. The 1883 Krakatau eruption affords us important knowledge on caldera formation even if it was not observed with modern sophisticated instruments. First, volcanic activities of the Krakatau Islands before and after the 1883 eruption are reexamined: previous suppositions involving a caldera-forming eruption of the proto-Krakatau prior to 1883 proved to be unsupported by the bathymetric topographies and gravity anomalies on and around the Krakatau Islands. Then, Anak Krakatau is interpreted as a parasitic cone of the main Krakatau volcano. As supplementary knowledge to discussion of caldera deposits, the results of drillings at several calderas in Japan and Mexico are introduced. Mass deficiency of the caldera deposit at Krakatau caldera is estimated by the gravity anomaly observed there and converted to probable volume with suitable density. For quantitative examination of the subsurface structure beneath the Krakatau complex, spatial distributions of seismic S-wave attenuation and Vp/Vs ratios have been already studied by temporary seismological observations and their results have been published. The high ratios of Vp/Vs observed approximately at Krakatau caldera may be attributable to the caldera deposit that is low density and contain much water. As additional remarks, a zone having both the characters, S-wave attenuation and zones of relatively high Vp/Vs ratio, may be a probable magma reservoir centering at a depth of about 10 km.</p><p><span style="font-size: 12.000000pt; font-family: 'TimesNewRomanPSMT';"><br /></span></p>


2021 ◽  
Author(s):  
Mirko Scheinert ◽  
Philipp Zingerle ◽  
Theresa Schaller ◽  
Roland Pail ◽  
Martin Willberg

&lt;p&gt;In the frame of the IAG Subcommission 2.4f &amp;#8220;Gravity and Geoid in Antarctica&amp;#8221; (AntGG) a first Antarctic-wide grid of ground-based gravity anomalies was released in 2016 (Scheinert et al. 2016). That data set was provided with a grid space of 10 km and covered about 73% of the Antarctic continent. Since then a considerably amount of new data has been made available, mainly collected by means of airborne gravimetry. Regions which were formerly void of any terrestrial gravity observations and have now been surveyed include especially the polar data gap originating from GOCE satellite gravimetry. Thus, it is timely to come up with an updated and enhanced regional gravity field solution for Antarctica. For this, we aim to improve further aspects in comparison to the AntGG 2016 solution: The grid spacing will be enhanced to 5 km. Instead of providing gravity anomalies only for parts of Antarctica, now the entire continent should be covered. In addition to the gravity anomaly also a regional geoid solution should be provided along with further desirable functionals (e.g. gravity anomaly vs. disturbance, different height levels).&lt;/p&gt;&lt;p&gt;We will discuss the expanded AntGG data base which now includes terrestrial gravity data from Antarctic surveys conducted over the past 40 years. The methodology applied in the analysis is based on the remove-compute-restore technique. Here we utilize the newly developed combined spherical-harmonic gravity field model SATOP1 (Zingerle et al. 2019) which is based on the global satellite-only model GOCO05s and the high-resolution topographic model EARTH2014. We will demonstrate the feasibility to adequately reduce the original gravity data and, thus, to also cross-validate and evaluate the accuracy of the data especially where different data set overlap. For the compute step the recently developed partition-enhanced least-squares collocation (PE-LSC) has been used (Zingerle et al. 2021, in review; cf. the contribution of Zingerle et al. in the same session). This method allows to treat all data available in Antarctica in one single computation step in an efficient and fast way. Thus, it becomes feasible to iterate the computations within short time once any input data or parameters are changed, and to easily predict the desirable functionals also in regions void of terrestrial measurements as well as at any height level (e.g. gravity anomalies at the surface or gravity disturbances at constant height).&lt;/p&gt;&lt;p&gt;We will discuss the results and give an outlook on the data products which shall be finally provided to present the new regional gravity field solution for Antarctica. Furthermore, implications for further applications will be discussed e.g. with respect to geophysical modelling of the Earth&amp;#8217;s interior (cf. the contribution of Schaller et al. in session G4.3).&lt;/p&gt;


1980 ◽  
Vol 34 (3) ◽  
pp. 251-264 ◽  
Author(s):  
Gerard Lachapelle ◽  
K. P. Schwarz

An evaluation of the empirical gravity anomaly covariance function using over 95 000 surface gravity anomalies in the North American Western Cordillera was carried out. A regression analysis of the data exhibits a strong and quasi-linear correlation of free air gravity anomalies with heights. This height correlation is removed from the free air anomalies prior to the numerical evaluation of the gravity anomaly covariance function. This covariance function agrees well with that evaluated previously by the authors for the remainder of Canada. A possible use for such a covariance function of ‘height independent’ gravity anomalies in mountainous areas is described. First, the height independent gravity anomaly at a point of known height is evaluated by least squares prediction using neighboring measured height independent gravity anomalies. Secondly, the part caused by the height correlation is calculated using linear regression parameters estimated previously and added to the predicted height independent gravity anomaly to obtain a predicted standard free air anomaly. This technique can be used to densify the coverage of free air anomalies for subsequent use in integral formulas of physical geodesy, e.g., those of Stokes and Vening Meinesz. This method requires that point topographic heights be given on a grid.


2014 ◽  
Vol 4 (1) ◽  
Author(s):  
L.E. Sjöberg

AbstractTraditional gravimetric geoid determination relies on Stokes’ formula with removal and restoration of the topographic effects. It is shown that this solution is in error of the order of the quasigeoid-to-geoid difference, which is mainly due to incomplete downward continuation (dwc) of gravity from the Earth’s surface to the geoid. A slightly improved estimator, based on the surface Bouguer gravity anomaly, is also biased due to the imperfect harmonic dwc the Bouguer anomaly. Only the third estimator,which uses the (harmonic) surface no-topography gravity anomaly, is consistent with the boundary condition and Stokes’ formula, providing a theoretically correct geoid height. The difference between the Bouguer and no-topography gravity anomalies (on the geoid or in space) is the “secondary indirect topographic effect”, which is a necessary correction in removing all topographic signals.


Author(s):  
Alessandro Bianchini ◽  
Giulia Andreini ◽  
Giovanni Ferrara ◽  
Lorenzo Ferrari ◽  
Dante Tommaso Rubino

Recent studies showed that a prompt detection of the stall inception, connected with a specific model to predict its associated aerodynamic force, could provide room for an extension of the left margin of the operating curve of high-pressure centrifugal compressors. In industrial machines working in the field, however, robust procedures to detect and identify the phenomenon are still missing, i.e. the operating curve is almost ever cut preliminary by the manufacturer by a proper safety margin; moreover, no agreement is found in the literature about a well-defined threshold to define the onset of the stall. In particular, in some cases the intensity of the arising subsynchronous frequency is compared to the revolution frequency, while in many other ones it is compared to the blade passage frequency. A large experience in experimental stall analyses collected by the authors revealed that in some cases unexpected spikes could make this direct comparison not reliable for a robust automatic detection. To this end, a new criterion was developed based on an integral analysis of the area subtended to the entire subsynchronous spectrum of the dynamic pressure signal of probes positioned just outside the impeller exit. A dimensionless parameter was then defined to account for the spectrum area increase in proximity to stall inception. This new parameter enabled the definition of a reference threshold to highlight the arising of stall conditions, whose validity and increased robustness was here verified based on a set of experimental analyses of different types of full-stage test cases of industrial centrifugal compressors at the test rig.


Geophysics ◽  
1977 ◽  
Vol 42 (7) ◽  
pp. 1458-1461 ◽  
Author(s):  
Bijon Sharma ◽  
T. K. Bose

The application of the method of the Fourier transform in interpreting gravity anomalies of faults has so far been based upon the Fourier transform of the gravity anomaly due to a single semi‐infinite block cut by a vertical fault. A general expression for the Fourier transform of the fault anomaly is here derived which is valid for an arbitrary angle of inclination of the fault plane. For deriving the general expression, the gravity anomaly of the fault is first separated into a constant and a variable term. The transforms of the two terms are calculated separately and then added to give the general expression for the Fourier transform of the fault anomaly.


Sign in / Sign up

Export Citation Format

Share Document