MODELS AT MERAPI VOLCANO DERIVED FROM GPS AND GRAVITY DATA A SUMMARY
|
|
- Mervin Hamilton
- 6 years ago
- Views:
Transcription
1 Tiede et al. 117 MODELS AT MERAPI VOLCANO DERIVED FROM GPS AND GRAVITY DATA A SUMMARY C. Tiede 1, K.F. Tiampo 2, J. Fernández 3 1 Max-Planck Institut für Astronomie, Heidelberg, Germany 2 Department of Earth Sciences, University of Western Ontario, London, ON, Canada. 3 Instituto de Astronomía y Geodesia (CSIC-UCM), Facultad de Ciencias Matemáticas,Ciudad Universitaria, Madrid, Spain Abstract Merapi, located at Central Java, Indonesia is one of the most active volcanoes in the entire world. Within the joint Indonesian-German research project several models have been generated to understand the volcano, its structure and behavior. We describe here a model based on the joint interpretation of the GPS and gravimetry measurements as well as a sensitivity analysis approach for network configuration purposes. 1. Introduction Merapi volcano (φ = ; λ = ; h = 2969 m (ellipsoidal, based on WGS84, measured in 2000)) located on Central Java, Indonesia, is one of the approximately 129 active volcanoes in the Indonesian islands. It is one of the most active and hazardous volcanoes worldwide as a result of the approximately 2 million people living in the nearby surroundings of this high-risk volcano who are permanently threatened by its activity. The high rate of volcanic and seismic activity along the entire Sunda arc between Sumatra and Flores is caused by the subduction zone which forms the 6000 to 7000 m deep Java trench approximately 250 km south of the coastline of Indonesia. An interdisciplinary research cooperation called MERAPI (Mechanism Evaluation, Risk Assessment, Prediction Improvement) between the Geoforschungszentrum Potsdam and the Volcanological Survey of Indonesia as coordinating parties studied the region from 1997 through An overview of all 13 projects of MERAPI is given by Zschau et al. (1998) whereby the project parts can be distinguished via their different targets into History: Geological work on magma evolution and eruption history Monitoring: Continuous measurements of seismology, gas analysis and deformation for monitoring purposes Structure analysis: Geophysical measurements of DC geoelectrics, gravity, Long Offset Transient Electromagnetics (LOTEM), magnetotellurics and seismics for exploration of the structure of the volcano. The group under the lead of Prof. C. Gerstenecker (Darmstadt University of Technology) was responsible for the project Untersuchungen säkularer Schwereänderungen am Merapi, Java: Ursachen und Wirkungen, Ge 381/ in the time period between 1997 and Results of the measurements, carried out within this subproject are shown in this article. For a detailed
2 118 Tiede et al. overview about the measurements, the circumstances and the details behind the observations, see Läufer et al. (2007), this issue. 2. Gravity and deformation network In 1997 a gravity and deformation network at Merapi volcano consisting of 23 stations was established, see Figure 1, to measure deformation and gravity changes related to the volcanic activity. Gravity changes dgi and three-dimensional deformations in east, north and height uxi, uyi, uzi with i=1, 2,..., 23, have been observed between 1997 and 2002 during 6 field campaigns, using 4 relative spring gravimeters (LaCoste & Romberg) and up to 11 Trimble GPS receivers. Figure 1: Gravity and three-dimensional deformation network organized in a loop configuration around Merapi volcano, Java, Indonesia. The numbers of the measured points correspond to the points given in Figure 2 Prior to the modelling of the gravity and deformation measurements a data preprocessing and analysis was carried out which required a detailed digital elevation model. This had been generated by the fusion of SAR-data and photogrammetric images, Läufer (2003) and validated by comparison with further satellite based digital elevation models Gerstenecker et al. (2005). Figure 2 shows the actual changes between observations in 2000 and Observed dg, ux, uy and uz are given as a straight, bold line including 1-σ error bars for each of the 20 observation points for which data was available in both epochs (2000 and 2002).
3 Tiede et al. 119 Figure 2: Gravity changes (dg) and three-dimensional deformation (ux, uy, uz) measured between 2000 and 2002 at Merapi volcano. The measured data are shown as thick solid line including 1-σ error bars. Points are sorted with respect to the distance to the top of the volcano (distance given in the right legend). The dashed line shows the effects of the combined fault, the thin solid line shows the superposed effects of the combined fault and the elasticgravitational model, Tiede et al. (2007) Only very small gravity changes and deformations at stations of distances larger than 2 km from the crater rim were measured (Figure 2, stations point number 8-20), and the near field stations show significant gravity changes and three-dimensional deformations. These local effects at the crater are modelled and explained by the superposition of the effects produced by a) a local, extensional, reverse fault in the summit region, b) a spherical magma intrusion, modeled as a point source in an elastic-gravitational medium. To model deformation and gravity changes produced by both of the sources, we consider the theoretical model developed by Manshina et al. (1971) for the source (a), and by Rundle (1982), Fernández et al. (1997; 2006) for source (b). The unknown model parameters have been determined by using global optimization techniques. For details see Tiede (2005a) and Tiede et al. (2007). The results of the model computation are given in Figure 2. The dotted line shows the effect of the reverse fault zone, the straight line shows the superposition of effects produced by both sources. Here it can be seen that the model explains those observations which are near the summit (Figure 2, point number 1-5) as well as those located at larger distances. A visualization of the superimposed models is shown in Figure 3. Both plots show the structural map after Beauducel et al. (2000) in which fracture zones are displayed as black straight lines. Again, the derived fault model fits very well to the existing and displayed faults from Beauducel et al. (2000).
4 120 Tiede et al. Figure 3: Plan view and three-dimensional view of Merapi. The plan view shows the superimposed model given by an elastic-gravitational source (black sphere) in the east as well as a normal fault (black area). The structural map of Beauducel et al. (2000) including fracture zones as solid lines is merged on the digital elevation model derived by Läufer (2003). In the three-dimensional view the dip angle of the fault is shown, Tiede et al. (2007).
5 Tiede et al. 121 The fault is explained by a normal fault with a horizontal opening of the fracture to the north. The elastic-gravitational source is located near the fumarole fields Woro and Gendol. Due to the very small mass of the source it is unlikely that this model represents a magmatic chamber, but is instead related to the fumarole fields. It is only necessary to consider both types of source, a fault and magma intrusion, in the modeling for the area closest to the crater, therefore the question can be raised which of the observation points are of interest for deriving this specific physical model. Or equivalent to that: If one can assume an a-priori physical model (e.g. the elastic-gravitational source model) for the region of interest, where should the sensors be placed in order to best determine the unknown model parameters and which sensor placements which can be neglected for the correct determination of the physical model of interest. In order to follow up this idea a global sensitivity analysis (Tiede et al., 2005a, 2005b) was applied to a synthetic generated sensor grid in the area of interest. The analysis has been computed for the different measurements (dg, ux, uy, uz) and for the unknown parameters given by the elastic-gravitational model on each of the grid points. 3. Sensitivity analysis Before the determination of the elastic-gravitational source parameters within the optimization approach a global sensitivity analysis has been applied, Tiede et al. (2005b) with the purpose of computing how and to what extent each unknown model parameter contributes to the variance in the observed values (dg, ux, uy, uz). Sensitivity analysis are used for answering questions like how much does each input parameter contribute to the variance in the output values are there any interactions between the input parameters is there a region of input parameter range for which the model variation is the largest are there insignificant model input parameters Sensitivity analysis can be divided into local and global approaches with respect to how they treat the input factors, Saltelli et al. (2004). Local sensitivity analyses are based on computing partial derivatives of the output function. It is practicable if the variation around the midpoint of the input factors is small. Generally the relation between input and output variables is assumed to be linear. Global sensitivity analyses explain the uncertainty in the output variables by the uncertainty of each input factor. The input factors are given with a probability density distribution. They are varied simultaneously and the sensitivity is measured over the whole range of each input factor. One example of a global sensitivity analysis method is given by the variance-based techniques which are based on Monte Carlo sampling, Saltelli et al. (2000). The basic idea of such types of sensitivity analyses is given by the expression of the sensitivity of the variance sigma. It is determined how much the variance of the input parameters contributes to the variance of each output value. When the model is non-linear and various input variables are affected by uncertainties of different orders of magnitude, a global sensitivity analysis should be used. Therefore, in this computation the Sobol variance-based global sensitivity analysis is used, Sobol (1993).
6 122 Tiede et al. Sensitivity analysis computation For the sensitivity analysis computation a synthetic grid consisting of 100 points was generated covering the second and third loop of the gravity and deformation network, see Figure 4. (East m, North m). Figure 4 (following page): Location of the synthetic point grid which forms the basis for carrying out the sensitivity analysis. The probability density functions are all assumed to be unique. The unknown input parameters for the elastic-gravitational model are given by the threedimensional position of the source (ξ, ψ, ζ ), its mass (m) and its energy (e) as the product of the cubic radial component and a pressure term. A Monte Carlo sampling for these unknown input parameters with samples is investigated within the boundaries given in Table 1. For each of these samples the forward model is applied (given by the elastic-gravitational source model) and the four different output values dg, ux, uy and uz are computed. Afterwards the variation of these output values is expressed by the variation of the input parameters of the elastic-gravitational model. The resulting sensitivity indices for each of the four types of measurements are summed up and normalized. Figure 5 shows the sensitivity map whereby the sensitivities are given for all unknown parameters in common and per each point (here one point covers the measurements of dg, ux, uy, uz).
7 Tiede et al. 123 Table 1. Upper and lower limits on the unknown parameters of the elastic-gravitational model which are used for the Monte Carlo sampling parameter lower limit upper limit ξ [m] ψ [m] ζ [m] m [1012 kg] 0 2 e [1014N m] Figure 5: Sensitivity map for the synthetic grid shown in Figure 4. The color coding is done regarding the sensitivity at each point for dg, ux, uy, uz regarding the unknown parameters ξ, ψ, ζ, m, e of the elastic-gravitational source model Color coding is done with respect to the sensitivity. A bright color indicates highly sensitive regions around Merapi for the measurement of dg, ux, uy and uz and the estimation of the parameters of the elastic-gravitational source model. By comparing these results with the network which was installed at Merapi (the grey line in Figure 5 indicates the second loop around Merapi; the small square shows the region of the third loop around the crater) it can be seen that the theoretical sensitivity suggests that only the points very near to the summit are of interest for carrying out the physical model. These are the points 1-10 which are within a radius of 3000 m around the crater where a larger sensitivity occurs.
8 124 Tiede et al. 4. Summary Gravity and three-dimensional deformation data were used in order to model the summit region of Merapi volcano. Our measurements cannot be explained by a shallow or deep magma chamber as anticipated in previous modelling approaches by Beauducel and Cornet (1999), Ratdomopurbo and Poupinet (2000), Camus et al. (2000) or Setiawan (2002). Instead of the magma source model, we could explain our measurements by a normal fault with an opening of the fracture to the north in the crater region of Merapi which coincides well with the results of Beauducel et al. (2000). The attempts to model an elastic-gravitational source result in a very shallow chamber with very small mass located near the fumarole fields Gendol and Woro. A magma source cannot be the source but it may reflect changes in the fumaroles, Tiede et al. (2007). Furthermore, we showed an approach using sensitivity analysis for sensor network configuration for those cases in which the physical model is known prior to the actual measurements. The results of the approach were compared with the actual model results generated by applying the elastic-gravitational source model to the gravity residuals and three-dimensional deformations. The comparison shows that before establishing the actual measurement locations, those points which are sensitive to variation in the physical model parameters can be selected. With the use of synthetic sensitivity maps it is possible to put the sensors in these regions where they are most sensitive to the changes in the unknown parameters of the expected physical model. Results show that just the inner part of the installed network at Merapi contributes to the determination of the physical model parameters. As the computed grid is very broad, the sensitivity map does not look smooth. In future attempts the sensitivities will be carried out for a denser grid. Acknowledgements The authors want to thank C. Gerstenecker for the very fruitful collaboration in recent years. Carola Tiede would like to thank C. Gerstenecker for his patient guidance, encouragement and advice that he provided throughout the time she was a PhD student at Darmstadt University of Technology - it has been a fantastic time. Kristy Tiampo is grateful for the opportunity to have collaborated with C. Gerstenecker in a very interesting and exciting part of the world, and for his excellent judgment in evaluating both scientific opportunities and young scientists. José Fernández has been happy to cooperate with C. Gerstenecker during the last years developing new and innovative scientific ideas and research and on the training of young scientists in geodesy and volcanism, thus having the opportunity to know great colleagues, good researchers and better friends. In addition we thank B. Lühr for reviewing the article. This research has been carried out in the frame of the projects MERAPI, Untersuchungen säkularer Schwereänderungen am Merapi, Java: Ursachen und Wirkungen, Ge 381/12 1-4; Spanish MEC CGL C02/BTE and an NSERC Discovery Grant. References Beauducel, F. and F. Cornet(1999): Collection and three-dimensional modeling of GPS and tilt data at Merapi volcano, Java, J. Geophys. Res. 104,
9 Tiede et al. 125 Beauducel, F., F. Cornet, E. Suhanto, Duquesnoy and T., M. Kasser (2000): Constraints on magma flux from displacements data at Merapi volcano, Java, Indonesia, J. Geophys. Res. 105, Fernández J., J. Rundle, R. Granell and T.-T. Yu (1997): Programs to compute deformation due to a magma intrusion in elastic-gravitational layered earth models, Computer and Geosciences, 23: Camus, G., A. Gourgaud, P.-C. Mossand-Berthommier and P.-M. Vincent (2000): Merapi (Central Java, Indonesia): An outline of the structural and magmatological evolution, with a special emphasis to the major pyroclastic events, J. Volcanol. and Geotherm. Res. 100, Fernández, J., M. Charco, J.B. Rundle and K.F. Tiampo (2006): A revision of the FORTRAN codes GRAVW to compute deformation produced by a point magma intrusion in elasticgravitational layered Earth models, Computers & Geosciences, 32/2, doi: /j.cageo Gerstenecker, C., G. Läufer, D. Steineck, C. Tiede and B. Wrobel (2005), Validation of digital elevation models around Merapi Volcano, Java, Indonesia, NHESS, 5, Läufer G. (2003): Erzeugung hybrider digitaler Hoehenmodelle aktiver Vulkane am Beispiel des Merapi, Indonesien, PhD thesis, Darmstadt University of Technology, ISBN Manshina, L., D. Smylie and D. Orphal (1971): The displacement fields of inclined faults, Bulletin of the Seismological Society of America, 61: Ratdomopurbo, A. and G. Poupinet (2000): An overview of the seismicity of Merapi volcano (Java Indonesia) , J. Volcanol. Geother. Res. 100, Rundle, J. B. (1982), Deformation, gravity, and potential changes due to volcanic loading of the crust, J. Geophys. Res., 87(B12), 10,729 10,744. Saltelli, A., K. Chan, E. Scott (Eds.) (2000): Sensitivity analysis, John Wiley and Sons, LTD, Chichester, New York, Weinheim, Brisbane, Singapore, Toronto Setiawan, A. (2002): Modeling of gravity changes on Merapi volcano observed between , PhD thesis, Darmstadt University of Technology, Darmstadt, Germany Sobol, I. (1993): Sensitivity analysis for non linear mathematical models, Math. Model. Comput. Exp., 1: Tiede, C., J. Fernández, C. Gerstenecker and K. Tiampo (2007): A hybrid model for the summit region of Merapi volcano, Java, Indonesia, derived from gravity changes and deformation measured between 2000 and 2002, Pure and Applied Geophysics, 164, Tiede (2005a): Integration of optimization algorithms with sensitivity analysis, with application to volcanic regions, Ph.D. Thesis, Online publication at Universitäts- und Landesbibliothek Darmstadt:
10 126 Tiede et al. Tiede C., K. Tiampo, J. Fernández and C. Gerstenecker (2005b): Deeper understanding of nonlinear geodetic data inversion using a quantitative sensitivity analysis, Nonlinear Processes in Geophysics, 12: Zschau, J., R. Sukhyar, M. Purbawinata, B. Lühr and M. Westerhaus (1998): Projekt MERAPI Interdisciplinary research at a high-risk volcano in Decade Volcanoes under Investigation, Dt. Geophys. Gesellschaft, 3: 3-8, edited by J. Zschau and M. Westerhaus Dt. Geophys. Gesellschaft
MEASUREMENT OF SURFACE DEFORMATION ASSOCIATED WITH HYDROFRACTURE
Proceedings, 11 th FIG Symposium on Deformation Measurements, Santorini, Greece, 2003. MEASUREMENT OF SURFACE DEFORMATION ASSOCIATED WITH HYDROFRACTURE Bassam Saleh 1, Balqies Sadoun 1 and Pierre Antoine
More informationUsing topographic signatures to classify internally and externally driven tilt anomalies at Merapi Volcano, Java, Indonesia
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L05310, doi:10.1029/2007gl032262, 2008 Using topographic signatures to classify internally and externally driven tilt anomalies at Merapi Volcano, Java, Indonesia
More informationStructural investigation of Mt. Merapi by an active seismic experiment
Structural investigation of Mt. Merapi by an active seismic experiment N. Maercklin 1,2, C. Riedel 2, W. Rabbel 2, U. Wegler 1, B.-G. Lühr 1, and J. Zschau 1 1 GeoForschungsZentrum Potsdam, Germany; 2
More informationWhat is the Relationship between Pressure & Volume Change in a Magma Chamber and Surface Deformation at Active Volcanoes?
SSAC-pv2007.QE522.PL1.1 What is the Relationship between Pressure & Volume Change in a Magma Chamber and Surface Deformation at Active Volcanoes? What factors control the magnitude of surface deformation?
More informationVolcano-tectonic interaction at Soufrière Hills volcano, Montserrat (W.I.), constrained by dynamic gravity data
IOP Conference Series: Earth and Environmental Science Volcano-tectonic interaction at Soufrière Hills volcano, Montserrat (W.I.), constrained by dynamic gravity data To cite this article: Stefanie Hautmann
More informationMalte Westerhaus 1 & Wolfgang Welle 2
1 of 10 2/22/2011 12:52 PM Malte Westerhaus 1 & Wolfgang Welle 2 1 Geodetic Institute, Karlsruhe University, Englerstr. 7, D-76128 Karlsruhe 2 GeoForschungsZentrum Potsdam, Telegrafenberg, D- 14473 Potsdam
More informationThe Santorini Inflation Episode, Monitored by InSAR and GPS
The Santorini Inflation Episode, Monitored by InSAR and GPS Ioannis Papoutsis 1,2, Xanthos Papanikolaou 2, Michael Floyd 3, Kang Hyeun Ji 4, Charalampos Kontoes 1, Demitris Paradissis 2, Demitris Anastasiou
More informationThe Afar Rift Consortium
The Afar Rift Consortium University of Bristol University of Oxford University of Cambridge British Geological Survey University of Rochester University of Auckland, NZ Opportunistic programme following
More informationInternal structure of Mount Merapi, Indonesia, derived from long-offset transient electromagnetic data
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B9, 2187, doi:10.1029/2001jb000148, 2002 Internal structure of Mount Merapi, Indonesia, derived from long-offset transient electromagnetic data M. Müller,
More informationRESISTIVITY AND REFRACTION SEISMIC MAPPING IN PRAMBANAN TEMPLE AREA
Sismanto, et.al., Resistivity and Refraction Seismic Mapping RESISTIVITY AND REFRACTION SEISMIC MAPPING IN PRAMBANAN TEMPLE AREA Sismanto 1, Eddy Hartantyo 1, and Agus Waluyo 2 1 Geophysics Laboratory,
More informationGEO-VIII November Geohazard Supersites and Natural Laboratories Progress Report. Document 9
GEO-VIII 16-17 November 2011 Geohazard Supersites and Natural Laboratories Progress Report Document 9 This document is submitted to GEO-VIII for information. Geohazard Supersites and Natural Laboratories
More informationMAGMATIC, ERUPTIVE AND TECTONIC PROCESSES IN THE ALEUTIAN ARC, ALASKA
MAGMATIC, ERUPTIVE AND TECTONIC PROCESSES IN THE ALEUTIAN ARC, ALASKA Introduction The Aleutian Arc contains roughly ten percent of the world s active volcanoes. Hardly a year goes by without a major eruption
More informationDETERMINATION OF SLIP DISTRIBUTION OF THE 28 MARCH 2005 NIAS EARTHQUAKE USING JOINT INVERSION OF TSUNAMI WAVEFORM AND GPS DATA
Synopses of Master Papers Bulletin of IISEE, 47, 115-10, 013 DETERMINATION OF SLIP DISTRIBUTION OF THE 8 MARCH 005 NIAS EARTHQUAKE USING JOINT INVERSION OF TSUNAMI WAVEFORM AND GPS DATA Tatok Yatimantoro
More informationReport for 15th PCGIAP Meeting at 18th UNRCC-AP Working Group 1 Regional Geodesy
Report for 15th PCGIAP Meeting at 18th UNRCC-AP Working Group 1 Regional Geodesy Chairman Shigeru Matsuzaka, Japan - shigeru@gsi.go.jp Vice Chairman Prof. Pengfei Cheng, China - chengpf@casm.ac.cn John
More informationENV-5004B/ENVK5005B. Figure 6. Student Registration No. ENV-5004B/ENVK5005B Version 2
ENV-5004B/ENVK5005B Figure 6 Student Registration No UNIVERSITY OF EAST ANGLIA School of Environmental Sciences Main Series UG Examination 014-15 SOLID EARTH GEOPHYSICS SOLID EARTH GEOPHYSICS WITH FIELDCOURSE
More informationANOTHER MEXICAN EARTHQUAKE! Magnitude 7.1, Tuesday Sept. 19, 2017
ANOTHER MEXICAN EARTHQUAKE! Magnitude 7.1, Tuesday Sept. 19, 2017 Why is there no oceanic crust older than 200 million years? SUBDUCTION If new oceanic crust is being continuously created along the earth
More informationPreliminary Microgravity Measurement Of Paka Geothermal Prospect
Preliminary Microgravity Measurement Of Paka Geothermal Prospect Presented by: Levi Shako and Calistus Ndongoli Date: 4/November/2016 Venue: Addis Ababa, Ethiopia Introduction Microgravity is a geophysical
More informationPlate Tectonics. entirely rock both and rock
Plate Tectonics I. Tectonics A. Tectonic Forces are forces generated from within Earth causing rock to become. B. 1. The study of the origin and arrangement of Earth surface including mountain belts, continents,
More informationA) B) C) D) 4. Which diagram below best represents the pattern of magnetic orientation in the seafloor on the west (left) side of the ocean ridge?
1. Crustal formation, which may cause the widening of an ocean, is most likely occurring at the boundary between the A) African Plate and the Eurasian Plate B) Pacific Plate and the Philippine Plate C)
More informationGRAVIMETRIC MAP OF CHILE. Services Mining and Engineering Geophysics Instituto Geográfico Militar
GRAVIMETRIC MAP OF CHILE Araneda M., Rivas L., Avendaño M, Sottolichio G., Rubio W. Services Mining and Engineering Geophysics segmi@netexpress.cl) Instituto Geográfico Militar (lrivas@igm.cl) Abstract
More informationREFERENCE: The Blue Planet An Introduction to Earth System Science. Brian J. Skinner and Barbara W. Murck (2011) Third Edition. John Wiley and Sons
REFERENCE: The Blue Planet An Introduction to Earth System Science. Brian J. Skinner and Barbara W. Murck (2011) Third Edition. John Wiley and Sons Inc. PLATE BOUNDARIES OCEAN FLOOR SEISMIC ACTIVITY WORLD'S
More informationSemi-Dynamic Datum of Indonesia
Technical Seminar Reference Frame in Practice, Semi-Dynamic Datum of Indonesia H.Z. Abidin, A.B. Wijanarto, S.T. Wibowo, Susilo, J. Efendi Geospatial Information Agency, Indonesia (BIG) Sponsors: Page
More informationUniversità degli Studi di Napoli Federico II, Dipartimento di Scienze della Terra Napoli, Italy
Three years of monitoring using leveling and hybrid gravimetry applied to geothermal sites of Soultz-sous-Forêts and Rittershoffen, Rhine Graben, France Gilbert FERHAT 1, Nolwenn PORTIER 1, Jacques HINDERER
More informationSylvain Charbonnier. PASI Workshop About 60% of Indonesians live around 16 active volcanoes on the island of Java
Numerical modelling of pyroclastic flows: a case study from the recent activity of Merapi Volcano, Central Java, Indonesia Sylvain Charbonnier PASI Workshop 2011 Case study: Merapi Volcano! About 60% of
More informationCOMPOSITION and PHYSICAL PROPERTIES GENERAL SUBJECTS. GEODESY and GRAVITY
COMPOSITION and PHYSICAL PROPERTIES Composition and structure of the continental crust Composition and structure of the core Composition and structure of the mantle Composition and structure of the oceanic
More informationDynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 2016 Mw 7.0 Kumamoto Earthquake
Dynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 016 Mw 7.0 Kumamoto Earthquake Heng-Yi Su 1 *, Aitaro Kato 1 Department of Earth Sciences, National Central University, Taoyuan
More informationEAS 116 Earthquakes and Volcanoes
EAS 116 Earthquakes and Volcanoes J. Haase Forecasting Volcanic Eruptions Assessment of Volcanic Hazard Is that volcano active? Mount Lassen: 12000 BP and 1915 Santorini, IT: 180,000 BP, 70,000 BP, 21000
More informationSupporting the response to the 2018 lower East Rift Zone and summit collapse at Kīlauea Volcano, Hawaiʻi
Hawaiʻi Supersite success story Supporting the response to the 2018 lower East Rift Zone and summit collapse at Kīlauea Volcano, Hawaiʻi Since 1983, Kīlauea Volcano, on the Island of Hawaiʻi, has actively
More informationRecent activity. Current episode 12 years Transition between styles 4 periods of dome growth Since 2003, Vulcanian explosions (4 25 per day)
Recent activity Current episode 12 years Transition between styles 4 periods of dome growth Since 2003, Vulcanian explosions (4 25 per day) 02 December 09 Explosions originating from dome e.g. 10 Jan 2010
More informationGPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise
GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise Strain Analysis Introduction Name: The earthquake cycle can be viewed as a process of slow strain accumulation
More informationMonitoring Tidal and Non-tidal Tilt Variations in Lanzarote Island (Spain)
, Journal of the Geodetic Society of Japan Vol. 47, No. 1, (2001), pp. 456-462 Monitoring Tidal and Non-tidal Tilt Variations in Lanzarote Island (Spain) Jose Arnosol)*, Ricardo Vieiral), Emiho Velezl),
More informationFig. 1. Joint volcanological experiment on volcanic structure and magma supply system in Japan.
2. Joint Volcanological Experiment on Volcanic Structure and Magma Supply System Since 1994, joint experiments have been conducted in several volcanoes in Japan to reveal the structure and the magma supply
More informationJoint inversion of active and passive seismic data in Central Java
Geophys. J. Int. (2007) doi: 10.1111/j.1365-246X.2007.03435.x Joint inversion of active and passive seismic data in Central Java Diana Wagner, 1 I. Koulakov, 2,3 W. Rabbel, 1 B.-G. Luehr, 2 A. Wittwer,
More informationBeneath our Feet: The 4 Layers of the Earty by Kelly Hashway
Beneath our Feet: The 4 Layers of the Earty by Kelly Hashway The Earth is more than a giant ball made up of dirt, rocks, and minerals. The Earth may look like a giant ball from when looking at it from
More informationAmbient Noise Tomography of Merapi Complex, Central Java, Indonesia: A Preliminary Result
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Ambient Noise Tomography of Merapi Complex, Central Java, Indonesia: A Preliminary Result To cite this article: Massita Putriastuti
More informationFrom Punchbowl to Panum: Long Valley Volcanism and the Mono-Inyo Crater Chain
From Punchbowl to Panum: Leslie Schaffer E105 2002 Final Paper Long Valley Volcanism and the Mono-Inyo Crater Chain Figure 1. After a sequence of earthquakes during the late 1970 s to the early 1980 s
More informationDirected Reading. Section: The Theory of Plate Tectonics. to the development of plate tectonics, developed? HOW CONTINENTS MOVE
Skills Worksheet Directed Reading Section: The Theory of Plate Tectonics 1. The theory that explains why and how continents move is called. 2. By what time period was evidence supporting continental drift,
More informationSurface deformation studies of Tenerife Island, Spain from joint GPS-DInSAR observations
Surface deformation studies of Tenerife Island, Spain from joint GPS-DInSAR observations Sergey Samsonov, Kristy Tiampo, Pablo J. González, Juan Prieto, Antonio G. Camacho and José Fernández GNS Science,
More information4.Mashu. Summary. Latitude: 43 34'20" N, Longitude: '39" E, Elevation: 857 m (Kamuinupuri) (Elevation Point) (4. Mashu)
4.Mashu Latitude: 43 34'20" N, Longitude: 144 33'39" E, Elevation: 857 m (Kamuinupuri) (Elevation Point) Overview of Mashu taken from 3 rd Observation Platform on west side on October 16, 2012 by the Japan
More informationAnswers: Internal Processes and Structures (Isostasy)
Answers: Internal Processes and Structures (Isostasy) 1. Analyse the adjustment of the crust to changes in loads associated with volcanism, mountain building, erosion, and glaciation by using the concept
More informationvolcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1
volcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1 Università della Calabria, Cosenza, Italy 2 Istituto Nazionale di Geofisica e Vulcanologia Osservatorio Vesuviano,
More informationStudents will be able, using GIS, to locate the largest and most destructive earthquakes;
Exploring Earthquake Hazards with GIS Laboratory for Introduction to Geophysics, Spring 2007 Prof. Constantin Cranganu Brooklyn College of the City University of New York Learning Objectives: Students
More informationSupplementary Material
1 Supplementary Material 2 3 4 Interseismic, megathrust earthquakes and seismic swarms along the Chilean subduction zone (38-18 S) 5 6 7 8 9 11 12 13 14 1 GPS data set We combined in a single data set
More informationMount Stromlo Gravity Station Gravimetry at the ANU
ANU COLLEGE OF SCIENCE RESEARCH SCHOOL OF EARTH SCIENCES Canberra ACT 0200 Australia Mount Stromlo Gravity Station Gravimetry at the ANU Background Terrestrial observing has been conducted at Mount Stromlo
More informationInternational Research Collaboration in Indonesia: LIPI as a special reference
Workshop on Science & Technology Cooperation With Developing Countries on Global Issues. Tokyo, 8 October 2008 International Research Collaboration in Indonesia: LIPI as a special reference Hery Harjono
More informationIn order to obtain a long term monitoring result for the Kilauea Volcano, ALOS PALSAR images taken on Track 287, Frame 38, ascending orbit with 21.5 d
ALOS PALSAR OBSERVATION OF KILAUEA VOLCANO ACTIVITIES FROM 2006 TO 2009 Zhe Hu, Linlin Ge, Xiaojing Li, Kui Zhang, Alex Hay-Man NG and Chris Rizos Cooperative Research Centre for Spatial Information &
More informationIntegrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya
More informationHow summit calderas collapse on basaltic volcanoes: new insights from the April 2007 caldera collapse of Piton de la Fournaise volcano
IOP Conference Series: Earth and Environmental Science How summit calderas collapse on basaltic volcanoes: new insights from the April 2007 caldera collapse of Piton de la Fournaise volcano To cite this
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi: 10.1038/ngeo739 Supplementary Information to variability and distributed deformation in the Marmara Sea fault system Tobias Hergert 1 and Oliver Heidbach 1,* 1 Geophysical
More informationValidation of Digital Elevation Models around Merapi Volcano, Java, Indonesia
Natural Hazards and Earth System Sciences, 5, 863 876, 2005 SRef-ID: 1684-9981/nhess/2005-5-863 European Geosciences Union 2005 Author(s). This work is licensed under a Creative Commons License. Natural
More informationVolcanic Benefits & Forecasting
Volcanic Benefits & Forecasting Review: https://www.youtube.com/watch?v=ydy28qtdyjy 1. Based on what we know about volcanoes, predict where you might you expect to see volcanoes in and around New Zealand?
More informationLatitude: 42 49'36" N, Longitude: '41" E, Elevation: 1,898 m (Ezo-Fuji) (Elevation Point)
16.Yoteizan Latitude: 42 49'36" N, Longitude: 140 48'41" E, Elevation: 1,898 m (Ezo-Fuji) (Elevation Point) Overview of Yoteizan taken from northwest side on May 18, 2003 by the Japan Meteorological Agency
More informationSupplementary Figure 1. Distribution of seismic event locations determined using the final 3-D velocity model. We separate the crust-related
Supplementary Figure 1. Distribution of seismic event locations determined using the final 3-D velocity model. We separate the crust-related seismicity at depths of less than 40 km (panel A) from the deeper
More informationLECTURE #11: Volcanoes: Monitoring & Mitigation
GEOL 0820 Ramsey Natural Disasters Spring, 2018 LECTURE #11: Volcanoes: Monitoring & Mitigation Date: 15 February 2018 I. What is volcanic monitoring? the continuous collection of one or more data sources
More information19. Esan Continuously Monitored by JMA
19. Esan Continuously Monitored by JMA Latitude: 41 48'17" N, Longitude: 141 09'58" E, Elevation: 618 m (Esan) (Triangulation Point) Overview of Esan, taken from east side on March 13, 2009 by the Japan
More informationStudents will be able, using GIS, to locate the largest and most destructive earthquakes;
Exploring Earthquake Hazards with GIS Laboratory for Introduction to Geophysics, Spring 2007 Prof. Constantin Cranganu Brooklyn College of the City University of New York Learning Objectives: Students
More informationAxial Symmetric Crustal Deformation Model for Long Valley Caldera, California
Axial Symmetric Crustal Deformation Model for Long Valley Caldera, California V. Magni 1*, M. Battaglia 1, P. Tizzani 2, A. Manconi 3 and T. Walter 3 1 Dept of Earth Sciences, University of Rome La Sapienza,
More informationOn Integrated Geodetic Monitoring for Sinkhole-Induced Surface Deformation and Mass Dislocation
On Integrated Geodetic Monitoring for Sinkhole-Induced Surface Deformation and Mass Dislocation - JISDM2016 - Joint International Symposium on Deformation Monitoring - Session - Multi-Sensor-Systems for
More informationUSING LIDAR MEASUREMENTS FOR IMPROVING OR UP-DATING A DEM
JAES_1(14)_3_2011 VAIS M. et. all, pp.123-128 SECTION Geodesic Engineering USING LIDAR MEASUREMENTS FOR IMPROVING OR UP-DATING A DEM VAIS Manuel*, IOSIF Gheorghe, Bucharest University, *e-mail: manuel.vais@sipg.ro
More informationLab 1: Plate Tectonics April 2, 2009
Name: Lab 1: Plate Tectonics April 2, 2009 Objective: Students will be introduced to the theory of plate tectonics and different styles of plate margins and interactions. Introduction The planet can be
More informationDo Now: Vocabulary: Objectives. Vocabulary: 1/5/2016. Wegener? (Can they move?) the idea that continents have moved over time?
Do Now: 1. Who was Alfred Wegener? 2. What was Pangaea? 3. Are continents fixed? (Can they move?) 4. What evidence supports the idea that continents have moved over time? Objectives What evidence suggests
More informationTopic 12 Review Book Earth s Dynamic Crust and Interior
Topic 12 Review Book Earth s Dynamic Crust and Interior Define the Vocabulary 1. asthenosphere 2. continental crust 3. Convection current 4. Convergent plate boundary 5. Divergent plate boundary 6. earthquake
More informationEXPLORE PLATE TECTONICS & MORE THROUGH GPS DATA. Shelley Olds, UNAVCO April 12, 2018 NGSS Webinar
EXPLORE PLATE TECTONICS & MORE THROUGH GPS DATA Shelley Olds, UNAVCO April 12, 2018 NGSS Webinar Today s Outline Central question: How do we know the tectonic plates are still moving? -- What evidence
More informationGEODETIC EVIDENCE FOR SLOW INFLATION OF THE SANTORINI CALDERA
Proceedings, 11 th FIG Symposium on Deformation Measurements, Santorini, Greece, 2003. GEODETIC EVIDENCE FOR SLOW INFLATION OF THE SANTORINI CALDERA Stathis Stiros 1, Aris Chasapis 2 and Villy Kontogianni
More informationMineral Element Compound Atom Crystal. Silicate Nonsilicate Inorganic Luster Streak. Cleavage Fracture Hardness Density Reclamation
6 th Grade Science Midterm Study Guide Due February 5 th * Midterm Date: February 9th (Tuesday) * Midterm will consist of 50 multiple choice questions (about 10 per chapter) YOU MUST DEFINE THE KEY TERMS
More informationApplication of differential SAR interferometry for studying eruptive event of 22 July 1998 at Mt. Etna. Abstract
Application of differential SAR interferometry for studying eruptive event of 22 July 1998 at Mt. Etna Coltelli M. 1, Puglisi G. 1, Guglielmino F. 1, Palano M. 2 1 Istituto Nazionale di Geofisica e Vulcanologia,
More informationMagnetic and Gravity Methods for Geothermal Exploration
Magnetic and Gravity Methods for Geothermal Exploration Dr. Hendra Grandis Geophysics - ITB method and survey procedure Aero- or ground magnetic (covers a large area) Schlumberger resistivity mapping and
More informationVolcanic Plumes. JOHN WILEY & SONS Chichester New York Weinheim Brisbane Singapore Toronto
Volcanic Plumes R. S. J. SPARKS University of Bristol, UK M. I. BURSIK State University of New York, USA S. N. CAREY University of Rhode Island, USA J. S. GILBERT Lancaster University, UK L. S. GLAZE NASA/Goddard
More informationStatus after Six Years
German-Indonesian Tsunami Early Warning System () Status after Six Years Joern Lauterjung and the -Team Challenge Indonesia Short Early Warning Time Tsunamis are generated at the trench (orange line) Traveltimes
More informationVolcanic Disaster Mitigation for Living in Harmony with Active Volcanoes
Volcanic Disaster Mitigation for Living in Harmony with Active Volcanoes SURONO Expert Staf of Ministry of Energy and Mineral Resources, Republic of Indonesia Presented on : The 9 th International Graduate
More informationOverview of Indonesian Geothermal System
Overview of Indonesian Geothermal System The Resources Beyond High Enthalpy Magmatic/Volcanic Systems By Dr.Eng. Suryantini ITB 4/3/2017 1 Outline Current Condition of Indonesian Geothermal Resources and
More informationOhio s State Tests PRACTICE TEST GRADE 8 SCIENCE. Student Name
Ohio s State Tests PRACTICE TEST GRADE 8 SCIENCE Student Name The Ohio Department of Education does not discriminate on the basis of race, color, national origin, sex, religion, age, or disability in employment
More informationSlow Deformation of Mt. Baekdu Stratovolcano Observed by Satellite Radar Interferometry
Slow Deformation of Mt. Baekdu Stratovolcano Observed by Satellite Radar Interferometry Sang-Wan Kim and Joong-Sun Won Department of Earth System Sciences, Yonsei University 134 Shinchon-dong, Seodaemun-gu,
More informationTAKE HOME EXAM 8R - Geology
Name Period Date TAKE HOME EXAM 8R - Geology PART 1 - Multiple Choice 1. A volcanic cone made up of alternating layers of lava and rock particles is a cone. a. cinder b. lava c. shield d. composite 2.
More information4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution
4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution PI: Lijun Liu Department of Geology, University of Illinois at Urbana-Champaign Corresponding author: Lijun Liu, ljliu@illinois.edu
More informationM. Charco (1), J. Fernández (1), M. J. Sevilla (1) and J. B. Rundle (2)
ISSN: 0214-4557 Modeling magmatic intrusion s effects on the geoid and vertical deflection. Application to Lanzarote, Canary Islands, and Long Valley Caldera, California. M. Charco (1), J. Fernández (1),
More informationGeophysical surveys Anomaly maps 2D modeling Discussion Conclusion
Introduction Geophysical surveys Anomaly maps 2D modeling Discussion Conclusion General context Growth of volcanic domes is a multiphase process: extrusion phases 2005 lava dome growth inside the Mount
More informationLecture 2 The Earth. I. The Interior. Earth vital statistics. Planet size & mass determination. Orbital semi-major axis x 10 8 km (1.
Lecture 2 The Earth. I. The Interior Vital statistics Size & mass of a planet Gravity & magnetic mapping Seismology Structure of Earth s interior Plate tectonics Reminder: lectures at http://www.star.le.ac.uk/~pto/planets.html
More informationDiverse deformation patterns of Aleutian volcanoes from InSAR
Diverse deformation patterns of Aleutian volcanoes from InSAR Zhong Lu 1, Dan Dzurisin 1, Chuck Wicks 2, and John Power 3 U.S. Geological Survey 1 Cascades Volcano Observatory, Vancouver, Washington 2
More informationModule 1, Investigation 3: Predicting Eruptions
Module 1, Investigation 3: Predicting Eruptions Introduction Welcome! Volcanoes are either "active" or "extinct". Active means that the volcano has erupted during the past 10,000 years. It can also mean
More informationBathymetry Measures the vertical distance from the ocean surface to mountains, valleys, plains, and other sea floor features
1 2 3 4 5 6 7 8 9 10 11 CHAPTER 3 Marine Provinces Chapter Overview The study of bathymetry determines ocean depths and ocean floor topography. Echo sounding and satellites are efficient bathymetric tools.
More informationGeophysical Surveys of The Geothermal System of The Lakes District Rift, Ethiopia
Geophysical Surveys of The Geothermal System of The Lakes District Rift, Ethiopia By: Befekadu Oluma By: Geophysics Department Geological Survey of Ethiopia The formation of the rift was preceded by a
More informationReconstruction of Subsurface Structures from Ambient Seismic Noise in Jakarta Area
Reconstruction of Subsurface Structures from Ambient Seismic Noise in Jakarta Area Sri Widiyantoro 1 Andri Dian Nugraha 1 Zulfakriza 2 Bayu Pranata 3 1 Global Geophysics Research Group, Faculty of Mining
More informationGravimetric Tide observation at Lake Nasser Region, Aswan, Egypt
Gravimetric Tide observation at Lake Nasser Region, Aswan, Egypt R.M. Hassan, E.M. Abdelrahman, A. Tealeb, K.H. Zahran and G. Jentzsch ABSTRACT The LaCoste and Romberg gravimeter D-218 of the National
More informationOhio s State Tests PRACTICE TEST LARGE PRINT GRADE 8 SCIENCE. Student Name
Ohio s State Tests PRACTICE TEST LARGE PRINT GRADE 8 SCIENCE Student Name The Ohio Department of Education does not discriminate on the basis of race, color, national origin, sex, religion, age, or disability
More informationContinental Drift and Plate Tectonics
Continental Drift and Plate Tectonics Continental Drift Wegener s continental drift hypothesis stated that the continents had once been joined to form a single supercontinent. Wegener proposed that the
More informationC5 Magnetic exploration methods data analysis techniques
C5 Magnetic exploration methods data analysis techniques C5.1 Data processing and corrections After magnetic field data have been collected a number of corrections are applied to simplify the interpretation.
More informationTHE INTERNAL STRUCTURE OF THE EARTH
UNIT 1 THE INTERNAL STRUCTURE OF THE EARTH 1.1.Earth s interior layers The interior of the Earth can be divided into layers according to: -Composition layers ( organized in order of increasing density
More informationCrustal Shear Wave Velocities Inferred From Love Wave Group Velocity Dispersion Between Mt. Merapi and Mt. Lawu, Central Java
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:12 No:06 38 Crustal Shear Wave Velocities Inferred From Love Wave Group Velocity Dispersion Between Mt. Merapi and Mt. Lawu, Central Java
More informationEarth. Properties of Earth. Earth's Interior. Earth is the planet that we know best
Earth Earth is the planet that we know best Our observations of the other planets are all compared to the conditions and properties of Earth Properties of Earth Semimajor Axis Orbital Eccentricity Orbital
More informationVolcano Seismicity and Tremor. Geodetic + Seismic
Volcano Seismicity and Tremor Seismic Imaging Geodetic + Seismic Model based joint inversion Geodetic Monitoring How is magma stored in the crust? geometry, volume and physical state of crustal melts.
More informationASSESSMENT OF GEOTHERMAL POTENTIAL AT UNGARAN VOLCANO, INDONESIA DEDUCED FROM NUMERICAL ANALYSIS
PROCEEDINGS, Thirty-Fourth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 9-11, 2009 SGP-TR-187 ASSESSMENT OF GEOTHERMAL POTENTIAL AT UNGARAN VOLCANO,
More informationWas the February 2008 Bukavu seismic sequence associated with magma intrusion?
Was the February 2008 Bukavu seismic sequence associated with magma intrusion? N. d Oreye, C. Lucas National Museum of Natural History, Luxembourg nicolas.doreye@ecgs.lu P. Gonzalez, J. Fernández Univ.
More informationInquiry: Sumatran earthquakes with GPS Earth Science Education
Inquiry: Sumatran earthquakes with GPS Earth Science Education www.earthobservatory.sg Preparation: Before doing this investigation, complete two introductory investigations using GPS data from UNAVCO
More informationHow to Build a Mountain and other Geologic Structures. But first a short review
How to Build a Mountain and other Geologic Structures But first a short review Where do we see deep earthquakes? What is happening there? What can happen at a plate boundary? 1. Plates can move apart
More informationThe Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets
The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets P. Banerjee, 1 F. F. Pollitz, 2 R. Bürgmann 3 * 1 Wadia Institute of Himalayan Geology, Dehra Dun, 248001, India. 2
More informationChapter Introduction Lesson 1 Earthquakes Lesson 2 Volcanoes Chapter Wrap-Up
Chapter Introduction Lesson 1 Earthquakes Lesson 2 Volcanoes Chapter Wrap-Up What causes earthquakes and volcanic eruptions? What do you think? Before you begin, decide if you agree or disagree with each
More informationLatitude: 34 31'13" N, Longitude: '45" E, Elevation: 508 m (Miyatsukayama) (Triangulation Point - Toshima)
59. Toshima Latitude: 34 31'13" N, Longitude: 139 16'45" E, Elevation: 508 m (Miyatsukayama) (Triangulation Point - Toshima) Overview of Toshima taken from northwest side on October 30, 2002 by the Japan
More informationOn the Development and Implementations of the New Semi- Dynamic Datum for Indonesia
On the Development and Implementations of the New Semi- Dynamic Datum for Indonesia SusiloSUSILO 1, HasanuddinZ. ABIDIN 2, IrwanMEILANO 2, Benyamin SAPEII 2 1) Center for Geodetic Control Network and Geodynamic,
More informationHow to Build a Mountain and other Geologic Structures. But first, questions
How to Build a Mountain and other Geologic Structures But first, questions Questions your students might ask How were Montana s mountains formed? How old are the mountains? What are the different ways
More information