Surface Anomalies Prior to Earthquakes Habibeh Valizadeh, Shattri B. Mansor Husaini Omar and Farid Azad Department of Civil Engineering Universiti Putra Malaysia Serdang, Selangor Malaysia shattri@eng.upm.edu.my Introduction Recently, new theories on underground geophysical and geochemical interactions occur during preparation stages of earthquakes and the resultant measurable variations have been put into test and some warning factors were suggested as earthquake precursors. In case of oceanic and coastal earthquakes, with thinner crust, these pre-earthquake activities may be detected through secondary oceanic and atmospheric phenomenon. Earthquake Precursor Earthquake Prediction 1
Earthquakes? Are they really predictable? Vibrations in the earth are caused by sudden release of energy. This energy is produced somewhere within the crust. Its formation and existence produce phenomena under, on and above the ground. Satellite-based measurements and ground observation networks can be specialized to monitor the earthquakes-related changes. Earthquake Precursors Temperature anomalies SLHF (higher atmosphere-surface energy exchange) Chl-a concentration Radon gas emission Crust Deformations Strange cloud formation Seismic pattern 2
Frequency and distribution( 1973-212) Earthquake Precursors Short Term Hours/Days/Weeks Long Term Months/Years/Decade Seismology Remote Sensing In-situ measurement Earthquake event Our concern is: Remote Sensing Earthquake precursors prior to the event (Hours/Days/Weeks) 3
Monitoring the Precursors Possibilities Free available remote sensing data covering large scales allow monitoring the earth s surface. Data providers produce high-quality and trustable data using in-situ measurement networks and validation models. Problems Attenuation of received signals. Significant seasonal and natural changes on the surface masking the earthquake-related anomalies. Anomalies due to human activities. Limited knowledge on the local fault regime; earthquake formation site and places of vibration. Low resolution remote sensing data and insufficient number of ground stations. Workflow Identifying the available maps, remote sensing data, reanalysis information, fault distribution maps and geological setting Seismographs Descriptive Analysis Analytical Analysis Spacebased data Historical shakes Fault maps and tectonic information Statistical analysis, visual inspection, abnormality detection, mapping the spatial distribution of variations Detection of the concurrent precursors, determination of active faulting Assessing the extents of the earthquake area and the possibility of monitoring seismic activity from satellite data for the case study area Evaluation of the available earthquake preparation theories by concurrencies of RS-based precursors and seismic records Area Characterization Expected results Earthquake characteristics Minor shake mapping Seismic gaps Statistical analysis Long-term prediction Determining the suitable precursors Recognizing of the Suspected Area for future quakes 4
Data Surface Latent Heat Flux from NCEP Surface Temperature from ASTER, AVHRR or AMSR-E Data Chlorophyll-a from MODIS Upwelling Indices from PFEL NCEP: National Center for Environmental Prediction ASTER: Advanced Spaceborn Thermal Emission and Reflection Radiometer AVHRR: Advanced Very High Resolusion Radiometer AMSR: Advanced Microwave Scaning Radiometer MODIS: Moderate Resolution Imaging Spectrodiometer PFEL: Pacific Fisheries Environmental Laboratory Oceanic Case Studies 1 5
Earthquake of California 25 Time series of thermal anomalies at the epicenter of the California earthquake showing high values a month before the main event. dashed line is the 1- year average of SST for the region. SST deg.c 2 25 1994-23 Sigma 2 Sigma 18 16 14 12 1 1-May 16-May 31-May 15-Jun 3-Jun 15-Jul 3-Jul 14-Aug 29-Aug Surface Heat Flux (w.m-2) 14 12 1 8 6 4 2 25 2-24 1-May 16-May 31-May 15-Jun 3-Jun 15-Jul 3-Jul 14-Aug 29-Aug Temporal variation in SLHF of the California earthquake covering the epicenter pixel showing increase in some occasions prior to the main event; dashed line is the 5- last-year average of SLHF for the region. 11 Changes in SLHF Spatio-temporal variation in SLHF prior and after the main event of the Northern California earthquake. 12 6
SST Anomalies 1 Jun 2 Jun 3 Jun 4 Jun 5 Jun 6 Jun 7 Jun 8 Jun 9 Jun 1 Jun 11 Jun 12 Jun 13 Jun 14 Jun 15 Jun 16 Jun 17 Jun 18 Jun Chl-a concentration in the ocean are intimately linked with the SST. Sudden changes in Chl-a distribution arises from sudden changes of sea thermal structure. 13 Variation of Chl-a 1 Jun 2 Jun 3 Jun 4 Jun 5 Jun 6 Jun 7 Jun 9 Jun 1 Jun 11 Jun 12 Jun 13 Jun 14 Jun 15 Jun 16 Jun 18 Jun 14 7
Upwelling Index Two major factors which cause rising in Chl-a concentration are ocean upwelling and sea surface temperature both of which are pre seismic indicators. 6 Daily averaged upwelling index for Northern California earthquake showing maximum rise some days prior to the main event; dashed line is the 1-year average of upwelling index for the region. 5 25 26-211 Chl 8d Composite 4 3 2 1 Upwelling Index (m3/s/1m) 4 2-2 -4 1-May 16-May 31-May 15-Jun 3-Jun 15-Jul 3-Jul 14-Aug 29-Aug 25 1995-24 -6 1-May 16-May 31-May 15-Jun 3-Jun 15-Jul 3-Jul 14-Aug 29-Aug 8-day averaged Chl-a for Northern California earthquake showing some high Chl-a matched the upwelling in terms of location and time; dashed line is the 6- year average of Chl-a for the region. 15 Earthquake of California 24 Surface Temperature deg.c 35 3 25 2 15 24 1994-23 Sigma 2 Sigma 1 1-Aug 16-Aug 31-Aug 15-Sep 3-Sep 15-Oct 3-Oct The anomalous SLHF values before and during the earthquake of September 28, 24; Red bar indicates the day of the main event. Surface Heat Flux (w.m-2) 1 8 6 4 2 24 1-Aug 16-Aug 31-Aug 15-Sep 3-Sep 15-Oct 3-Oct Time series of surface temperature; shows several anomalies during the preparation stage and sudden fall after the main event. Upwelling Index (m3/s/1m) 3 25 2 15 1 5-5 24 1994-23 Sigma -1 1-Aug 16-Aug 31-Aug 15-Sep 3-Sep 15-Oct 3-Oct Daily averaged upwelling index, showing rises before the main event. 16 8
Variation of Chl-a 2 Sep 5 Sep 8 Sep 14 Sep 2 Sep 21 Sep 22 Sep 28 Sep 29 Sep 1 Oct 4 Oct 7 Oct Temporal distribution of Chl-a concentration ; the increasing trend to the day of the event and general decrease in the area afterwards is obvious. 17 Earthquake of California 23 Surface Temperature deg.c Upwelling Index (m3/s/1m) 18 17 16 15 14 23 1994-23 Sigma 2 Sigma 13 1-Nov 16-Nov 1-Dec 16-Dec 31-Dec 15-Jan 3-Jan 1-1 -2-3 23-4 1993- -5 22-6 1-Nov 16-Nov 1-Dec 16-Dec 31-Dec 15-Jan 3-Jan Surface Heat Flux (w.m-2) 2 16 12 8 4 23 1998-22 1-Nov 16-Nov 1-Dec 16-Dec 31-Dec 15-Jan 3-Jan The SST and SLHF time series of the epicentral oceanic water of the earthquake of December 22, 23 generally higher before the earthquake. The effect of aftershocks during the second half of December is also shown. 3 1.8 2 1.6 23 24-26 Sigma 2 Sigma Chl-a 8d Composite 1.4 1.2 1.8.6.4.2 1-Nov 16-Nov 1-Dec 16-Dec 31-Dec 15-Jan 3-Jan Daily averaged upwelling index and 8-day composite Chlorophyll-a time series, showing anomalies in some occasions from one month before the main event followed by a sudden downwelling and Chl-a decrease immediately after that. 1 8 9
Thermal and Heat Surface Heat Flux (w.m-2) 1 24 8 1999-23 6 4 2 1-Aug 16-Aug 31-Aug 15-Sep 3-Sep 15-Oct 3-Oct F LH = LvCeUa (qs qa) FLH: Surface evaporation Lv: Latent heat of condensation Ce: Surface exchange coefficient for moisture Ua: Surface wind speed, Qs: Saturated specific humidity at ocean surface Qa: Air specific humidity at 2 m above the surface. Surface Temperature deg.c 35 3 25 2 15 24 1994-23 Sigma 2 Sigma 1 1-Aug 16-Aug 31-Aug 15-Sep 3-Sep 15-Oct 3-Oct Ts: LST ε: Surface emissivity γ and δ: two parameters dependent on the Planck s function ψ1, ψ2, and ψ3 : Referred to as atmospheric functions (AFs) Increased Chl-a on ocean surface 1
Chl-a & Upwelling Chl-a 8d Composite 1.8 1.6 1-Yr Average Sigma 2 Sigma 1.4 1.2 1.8.6.4.2 1-Nov 16-Nov 1-Dec 16-Dec 31-Dec 15-Jan 3-Jan Upwelling Index (m3/s/1m) 2-2 -4-6 1-Nov 16-Nov 1-Dec 16-Dec 31-Dec 15-Jan 3-Jan Increased SLHF & active faults zones 25 2 28 2-27 Sigma 2 Sigma SLHF (W/m2) 15 1 5 1-Jan 16-Jan 31-Jan 15-Feb 1-Mar 16-Mar 31-Mar 11
Earthquakes of Indonesia No. Place Date Longitude Latitude Magnitude Focal Depth (km) 1 Simeulue, Indonesia Feb 2, 28 95.978 E 2.778 N 7.4 35 2 Kepulauan, Indonesia Feb 25, 28 1.18 E 2.351 S 7.2 35 SSH and SLHF Anomalies Images of SSH retrieved from AMSR-E in the Indian Ocean during the Simeulue and Kepulauan earthquakes of February, 28 showing significant rises near epicenters one week before and during the earthquake events. Sharp rises in SLHF values of the pixels covering the epicenter of 25 th Feb, 28 earthquake is observable from the end of January to few days before the main event. Red bar is the day of the main event. SLHF (W/m2) 25 28 2-27 Sigma 2 Sigma 2 15 1 5 1-Jan 16-Jan 31-Jan 15-Feb 1-Mar 16-Mar 31-Mar 12
Seismic study micro-shake detection using seismograph interpretation: Evaluating the shaking rate before the main events Understanding the possible hidden fault pattern and local faulting activity by statistical analyses of the various information, related to foreshocks and aftershocks. Discovering the time and intensity frames of the possible correlation between seismic and remote sensing precursors. Findings The systematic patterns of SLHF along earthquake origins. Relative humidity, surface and air temperature values are warning signals of an impending earthquake (2-3 weeks prior to the main event). 2-3 weeks before the earthquakes the productivity rate of the open ocean water exceeded the average values. 13
Benefits Remote sensing techniques allow monitoring the earthquake precursory factors anomalies over large areas to detect tectonic activity and understand the mechanism of earthquake preparation processes to provide possibilities of a reliable prediction of these potential precursors in different parts of the world. Q & A?? For further information please contact: Shattri Mansor Remote Sensing & GIS Research Centre Faculty of Engineering Universiti Putra Malaysia 434 UPM Serdang Phone: +619-2244333 E-mail: shattri@eng.upm.edu.my 14