W.J.F. Simons
Please Note
20 records found
1
The provinces of Bangkok, Samut Prakan, Samut Sakhon, and Nakhon Pathom in Thailand are experiencing subsidence caused by land subsidence, tectonic activity, and sea-level rise. INSAR result from 2015-2022 show that Bangkok and nearby provinces subsided up to 3 cm/yr in the past 20 years. GNSS results show absolute subsidence rates (below 20 m) up to 5 mm/yr in the past 25 years. According to satellite altimetry data, Bangkok is currently experiencing a sea-level rise of up to 5 mm per year in the Gulf of Thailand. Ground water pumping also play an important role on land subsidence.
A Tsunami Generated by a Strike-Slip Event
Constraints From GPS and SAR Data on the 2018 Palu Earthquake
A devastating tsunami struck Palu Bay in the wake of the 28 September 2018 Mw = 7.5 Palu earthquake (Sulawesi, Indonesia). With a predominantly strike-slip mechanism, the question remains whether this unexpected tsunami was generated by the earthquake itself, or rather by earthquake-induced landslides. In this study we examine the tsunami potential of the co-seismic deformation. To this end, we present a novel geodetic data set of Global Positioning System and multiple Synthetic Aperture Radar-derived displacement fields to estimate a 3D co-seismic surface deformation field. The data reveal a number of fault bends, conforming to our interpretation of the tectonic setting as a transtensional basin. Using a Bayesian framework, we provide robust finite fault solutions of the co-seismic slip distribution, incorporating several scenarios of tectonically feasible fault orientations below the bay. These finite fault scenarios involve large co-seismic uplift (>2 m) below the bay due to thrusting on a restraining fault bend that connects the offshore continuation of two parallel onshore fault segments. With the co-seismic displacement estimates as input we simulate a number of tsunami cases. For most locations for which video-derived tsunami waveforms are available our models provide a qualitative fit to leading wave arrival times and polarity. The modeled tsunamis explain most of the observed runup. We conclude that co-seismic deformation was the main driver behind the tsunami that followed the Palu earthquake. Our unique geodetic data set constrains vertical motions of the sea floor, and sheds new light on the tsunamigenesis of strike-slip faults in transtensional basins.
Hazardous tsunamis are known to be generated predominantly at subduction zones. However, the 2018 Mw 7.5 Palu (Indonesia) earthquake on a strike-slip fault generated a tsunami that devastated the city of Palu. The mechanism by which this tsunami originated from such an earthquake is being debated. Here we present near-field ground motion (GPS) data confirming that the earthquake attained supershear speed, i.e. a rupture speed greater than the shear wave speed of the host medium. We subsequently study the effect of this supershear rupture on tsunami generation by coupling the ground motion to a 1-D non-linear shallow-water wave model accounting for both time-dependent bathymetric displacement and velocity. With the local bathymetric profile of Palu bay around a tidal station, our simulations reproduce the tsunami arrival and motions observed by CCTV cameras. We conclude that Mach (shock) fronts, generated by the supershear speed, interacted with the bathymetry and contributed to the tsunami.
Copernicus Sentinel–1 is a C-Band Synthetic Aperture Radar (SAR) satellite mission within the European Copernicus Programme. The two satellites Sentinel-1A and -1B were launched in April 2014 and 2016, respectively. The Copernicus POD (Precise Orbit Determination) Service is responsible for the determination of orbital and auxiliary products required by the Payload Data Ground Segment (PDGS). Precise orbits are determined based on the dual-frequency GPS (Global Positioning System) data delivered by dedicated geodetic-grade GPS receivers on-board the satellites. Several updates in the operational orbit determination were done during the years including an update of the GPS antenna reference point coordinates. The switch to GPS carrier phase ambiguity-fixing was a major improvement. A reprocessing of the entire mission span of both satellites became necessary to provide a consistent orbit time series for the mission based on state-of-the-art models and processing settings. Due to the lack of independent observation techniques, the Sentinel-1 orbit quality has been assessed by analysing processing metrics, orbit overlaps and orbit comparisons. For this purpose, members of the Copernicus POD Quality Working Group (QWG) provided reprocessed Sentinel-1 orbit time series based on their software packages and their orbit determination settings. A weighted average of all five delivered solutions - a combined orbit - serves as reference for the comparisons. The quality and reliability of this reference orbit depends among others on the number of available orbit solutions and whether a manoeuvre has been performed during the processed day or not. The mean orbit consistency between all orbit solutions is below 1 cm in 3D RMS for the entire mission time interval for both satellites. Only few days show inferior quality due to data gaps or orbit manoeuvres. Following this sophisticated validation process, the reprocessed Sentinel-1 orbits from the Copernicus POD Service have been made available to the user community.
The Peninsular Malaysia Geodetic Vertical Datum 2000 (PMGVD2000) inherited several deficiencies due to offsets between local datums used, levelling error propagations, land subsidence, sea level rise, and sea level slopes along the southern half of the Malacca Strait on the west coast and the South China Sea in the east coast of the Peninsular relative to the Port Klang (PTK) datum point. To cater for a more reliable elevation-based assessment of both sea level rise and coastal flooding exposure, a new epoch-based height reference system PMGVD2022 has been developed. We have undertaken the processing of more than 30 years of sea level data from twelve tide gauge (TG) stations along the Peninsular Malaysia coast for the determination of the relative mean sea level (RMSL) at epoch 2022.0 with their respective trends and incorporates the quantification of the local vertical land motion (VLM) impact. PMGVD2022 is based on a new gravimetric geoid (PMGeoid2022) fitted to the RMSL at PTK. The orthometric height is realised through the GNSS levelling concept H = hGNSS–Nfit_PTK–NRMDT, where NRMDT is a constant offset due to the relative mean dynamic ocean topography (RMDT) between the fitted geoid at PTK and the local MSL datums along the Peninsular Malaysia coast. PMGVD2022 will become a single height reference system with absolute accuracies of better than ±3 cm and ±10 cm across most of the land/coastal area and the continental shelf of Peninsular Malaysia, respectively.
Altimetry for the future
Building on 25 years of progress
The 2018 (Formula presented.) Palu earthquake is a remarkable strike-slip event due to its nature as a shallow supershear fault rupture across several segments and a destructive tsunami that followed coseismic deformation. GPS offsets in the wake of the 2018 earthquake display a transient in the surface motions of northwest Sulawesi. A Bayesian approach identifies (predominantly aseismic) deep afterslip on and below the coseismic rupture plane as the dominant physical mechanism causing the cumulative, postseismic, surface displacements whereas viscous relaxation of the lower crust and poro-elastic rebound contribute negligibly. We confirm a correlation between shallow supershear rupture and postseismic surface transients with afterslip activity in the zone below an interseismically locked fault plane where the slip rate tapers from zero to creeping.
Malaysia is located at the stable part of the tec-tonic Sundaland platelet in SE Asia. The platelet is surrounded in almost every direction by tectonically active convergent boundaries, at which the Philippine Sea, the Australian and the Indian Plates are subducting respectively from the East, South and West.The current Malaysia geodetic reference frame called MGRF2000 is a static reference frame and hence did not incorporate the effects of plate motion and the ensuing deformation from (megath-rust) earthquakes. To prevent degradation of Continuously Operating Reference Station (CORS) coordinates, a new time-dependent national reference frame was developed. Taking advantage of the availability of the GNSS data of the CORS network in Malaysia, notably the Malaysia Active GPS System (MASS) and Malaysia Real-Time Kinematic GNSS Network (MyRTKnet), a more accurate and robust Malaysian geodetic reference frame was determined, fully aligned and compatible with ITRF2014. The cumulative solution obtained from stacking Malaysian CORS position time series formed the basis of the new MGRF2020 realization. It consists of 100+ station positions at epoch 2020.0, station velocities and Post-Seismic Deformation (PSD) parametric models for stations subjected to major earthquakes. The (1999-2018) position time series exhibit Weighted Mean Root Square (WRMS) values of 3.0, 3.2 and 7.6 mm in respectively the East, North and Vertical components. A new semi-kinematic geodetic datum (GDM2020) for Malaysia, useable for GIS, mapping and cadastre applications is proposed to replace the existing static datum (GDM2000). A transformation suite to convert the spatial databases from GDM2000 to GDM2020 was also developed.
Vertical motion of Phuket Island (1994–2018) due to the Sumatra-Andaman mega-thrust earthquake cycle
Impact on sea-level and consequences for coral reefs
Macro-tidal coral reefs are particularly sensitive to medium to long-term changes in sea-level. Vertical motions of the seabed contribute to both lower or higher relative sea-level changes, particularly in tectonic plate boundary deformation zones along active subduction trenches. Phuket Island in Southern Thailand is subject to both horizontal and vertical land deformations during the seismic cycle of mega thrust earthquakes along the Sumatra and Andaman trenches. The relative sea-level changes in this region were historically monitored using the space geodetic techniques GPS and satellite altimetry alongside the traditional tide-gauge measurements over a period of almost 25 years. The GPS results show that the south of Thailand is still undergoing post-seismic deformations from the 2004 Mw 9.2 Sumatra-Andaman earthquake, after a significant change in the vertical motion of Phuket: from stable quasi-linear uplift at 2.5 ± 0.2 mm/yr, to (temporary) non-linear subsidence rates of 5–10 mm/yr in the past 14 quake aftermath years. The satellite altimetry data estimates the absolute sea-level rise in the Andaman Sea around the island at 3.9 ± 0.5 mm/yr. Therefore relative sea-level changes in Phuket appear to have been slightly positive until the end of 2004, followed by a significant increase averaging ~9 mm/yr that accumulated in 12 ± 1 cm by the end of 2018. Tide-gauge data gives a very similar result. The increased sea-level rise from 2005 onwards, due to tectonic land subsidence, correlates with the enhanced shallow coral reef growth and swift recoveries from bleaching events reported by long-term coral studies in this area.
North Sulawesi Trench or Minahasa subduction area is a subduction zone between the oceanic crust of Sulawesi Sea and the North Sulawesi Arm located at the triple junction in Eastern Indonesia. This subduction activity causes the North Arm of Sulawesi as an earthquake-prone area. Tectonic activities in the region can be studied through geodetic monitoring using GNSS GPS observations and by physical modelling from the rate of geodetic geometric results. Yearly GNSS GPS campaign have been conducted in the region from 1997 to 2008 and continuously observed by BIG from 2008 to 2016 using permanent GNSS GPS stasion. The coordinates of monitoring stations realized in ITRF-2008 provide residual RMS values of 3.13 mm, 4.15 mm and 7.26 mm for the northern, eastern and vertical components, where this indicates a high degree of accuracy. A simple estimation profile using GNSS GPS data based on the Okada elastic equation for the subduction zone shows a subduction movement ranging from 4 to 5 cm/yr with a locking depth of about 50 km, a dip 300 and ending in the post-seismic phase due to the sequence of earthquakes occurring in Minahasa since January 1, 1996 Mw 7.9 to 16 June 2002 Mw 5.9.
Sulawesi microplate Island is located at famous triple junction area of the Eurasian, India-Australian, and Philippine Sea plates. Under the influence of the northward moving Australian plate and the westward motion of the Philippine plate, the island at Eastern part of Indonesia is collide and with the Eurasian plate and Sunda Block. Those recent microplate tectonic motions can be quantitatively determine by GNSS-GPS measurement. We use combine GNSS-GPS observation types (campaign type and continuous type) from 1997 to 2015 to derive newly velocity field of the area. Several strategies are applied and tested to get the optimum result, and finally we choose regional strategy to reduce error propagation contribution from global multi baseline processing using GAMIT/GLOBK 10.5. Velocity field are analyzed in global reference frame ITRF 2008 and local reference frame by fixing with respect alternatively to Eurasian plate - Sunda block, India-Australian plate and Philippine Sea plates. Newly results show dense distribution of velocity field. This information is useful for tectonic deformation studying in geospatial era.
One of methods that can be used to determine the tectonic deformation status is rate estimation from geometric rotation and strain using quantitative velocity data from GPS observations. Microplate Sulawesi region located in the zone of triple junction (Eurasia, Australia and Philippine Sea Plates) has very complex tectonic and seismic condition, which is why become very important to know its recent deformation status in order to study neotectonic and disaster mitigation. Deformation rate quantification is estimated in two parameters: rotation and geodetic strain rate of each GPS station Delaunay triangle in the study area. The analysis in this study is not done using the grids since there is no rheological information at location that can be used as the interpolation-extrapolation constraints. Our analysis reveals that Sulawesi is characterized by rapid rotation in several different domains and compression-strain pattern that varies depending on the type and boundary conditions of microplate. This information is useful for studying neo tectonic deformation status and earthquake disaster mitigation.
The existence of intra-plate deformation of the Sundaland platelet along its eastern edge in North Borneo, South-East Asia, makes it an interesting area that still is relatively understudied. In addition, the motion of the coastal area of North-West Borneo is directed toward a frontal fold-and-thrust belt and has been fueling a long debate on the possible geophysical sources behind it. At present this fold-and-thrust belt is not generating significant seismic activity and may also not be entirely active due to a decreasing shelfal extension from south to north. Two sets of Global Positioning System (GPS) data have been used in this study; the first covering a time period from 1999 until 2004 (ending just before the Giant Sumatra–Andaman earthquake) to determine the continuous Sundaland tectonic plate motion, and the second from 2009 until 2011 to investigate the current deformations of North Borneo. Both absolute and relative positioning methods were carried out to investigate horizontal and vertical displacements. Analysis of the GPS results indicates a clear trend of extension along coastal regions of Sarawak and Brunei in North Borneo. On the contrary strain rate tensors in Sabah reveal that only insignificant and inconsistent extension and compression occurs throughout North-West Borneo. Moreover, station velocities and rotation rate tensors on the northern part of North Borneo suggest a clockwise (micro-block) rotation. The first analysis of vertical displacements recorded by GPS in North-West Borneo points to low subsidence rates along the western coastal regions of Sabah and inconsistent trends between the Crocker and Trusmadi mountain ranges. These results have not been able to either confirm or reject the hypothesis that gravity sliding is the main driving force behind the local motions in North Borneo. The ongoing Sundaland–Philippine Sea plate convergence may also still play an active role in the present-day deformation (crustal shortening) in North Borneo and the possible clockwise rotation of the northern part of North Borneo as a micro-block. However, more observations need to be collected to determine if the northern part of North Borneo indeed is (slowly) moving independently.
Sentinel-1A is the first satellite of the European Copernicus programme. Equipped with a Synthetic Aperture Radar (SAR) instrument the satellite was launched on April 3, 2014. Operational since October 2014 the satellite delivers valuable data for more than two years. The orbit accuracy requirements are given as 5. cm in 3D. In order to fulfill this stringent requirement the precise orbit determination (POD) is based on the dual-frequency GPS observations delivered by an eight-channel GPS receiver.The Copernicus POD (CPOD) Service is in charge of providing the orbital and auxiliary products required by the PDGS (Payload Data Ground Segment). External orbit validation is regularly performed by comparing the CPOD Service orbits to orbit solutions provided by POD expert members of the Copernicus POD Quality Working Group (QWG). The orbit comparisons revealed systematic orbit offsets mainly in radial direction (approx. 3. cm). Although no independent observation technique (e.g. DORIS, SLR) is available to validate the GPS-derived orbit solutions, comparisons between the different antenna phase center variations and different reduced-dynamic orbit determination approaches used in the various software packages helped to detect the cause of the systematic offset. An error in the given geometry information about the satellite has been found. After correction of the geometry the orbit validation shows a significant reduction of the radial offset to below 5. mm. The 5. cm orbit accuracy requirement in 3D is fulfilled according to the results of the orbit comparisons between the different orbit solutions from the QWG.