Wim Simons

Bio: Wim Simons is an academic researcher from Delft University of Technology. The author has contributed to research in topics: Sea level & Tide gauge. The author has an hindex of 17, co-authored 40 publications receiving 2071 citations. Previous affiliations of Wim Simons include Chulalongkorn University & Royal Thai Survey Department.

A decade of GPS in Southeast Asia: Resolving Sundaland motion and boundaries

Abstract: A unique GPS velocity field that spans the entire Southeast Asia region is presented. It is based on 10 years (1994–2004) of GPS data at more than 100 sites in Indonesia, Malaysia, Thailand, Myanmar, the Philippines, and Vietnam. The majority of the horizontal velocity vectors have a demonstrated global accuracy of ?1 mm/yr (at 95% confidence level). The results have been used to (better) characterize the Sundaland block boundaries and to derive a new geokinematic model for the region. The rotation pole of the undeformed core of the Sundaland block is located at 49.0°N–94.2°E, with a clockwise rotation rate of 0.34°/Myr. With respect to both geodetically and geophysically defined Eurasia plate models, Sundaland moves eastward at a velocity of 6 ± 1 to 10 ± 1 mm/yr from south to north, respectively. Contrary to previous studies, Sundaland is shown to move independently with respect to South China, the eastern part of Java, the island of Sulawesi, and the northern tip of Borneo. The Red River fault in South China and Vietnam is still active and accommodates a strike?slip motion of ?2 mm/yr. Although Sundaland internal deformation is general very small (less than 7 nanostrain/yr), important accumulation of elastic deformation occurs along its boundaries with fast?moving neighboring plates. In particular in northern Sumatra and Malaysia, inland?pointing trench?perpendicular residual velocities were detected prior to the megathrust earthquake of 26 December 2004. Earlier studies in Sumatra already showed this but underestimated the extent of the deformation zone, which reaches more than 600 km away from the trench. This study shows that only a regional Southeast Asia network spanning thousands of kilometers can provide a reference frame solid enough to analyze intraplate and interplate deformation in detail.

Insight into the 2004 Sumatra–Andaman earthquake from GPS measurements in southeast Asia

Abstract: Data collected at approximately 60 Global Positioning System (GPS) sites in southeast Asia show the crustal deformation caused by the 26 December 2004 Sumatra-Andaman earthquake at an unprecedented large scale. Small but significant co-seismic jumps are clearly detected more than 3,000 km from the earthquake epicentre. The nearest sites, still more than 400 km away, show displacements of 10 cm or more. Here we show that the rupture plane for this earthquake must have been at least 1,000 km long and that non-homogeneous slip is required to fit the large displacement gradients revealed by the GPS measurements. Our kinematic analysis of the GPS recordings indicates that the centroid of released deformation is located at least 200 km north of the seismological epicentre. It also provides evidence that the rupture propagated northward sufficiently fast for stations in northern Thailand to have reached their final positions less than 10 min after the earthquake, hence ruling out the hypothesis of a silent slow aseismic rupture.

Insight into the 2004 Sumatra–Andamanearthquake from GPS measurements insoutheast Asia

Abstract: Data collected at ,60 Global Positioning System (GPS) sites in southeast Asia show the crustal deformation caused by the 26 December 2004 Sumatra–Andaman earthquake at an unprecedented large scale. Small but significant co-seismic jumps are clearly detected more than 3,000 km from the earthquake epicentre. The nearest sites, still more than 400 km away, show displacements of 10 cm or more. Here we show that the rupture plane for this earthquake must have been at least 1,000 km long and that non-homogeneous slip is required to fit the large displacement gradients revealed by the GPS measurements. Our kinematic analysis of the GPS recordings indicates that the centroid of released deformation islocated at least 200 km north of the seismological epicentre. It also provides evidence that the rupture propagated northward sufficiently fast for stations in northern Thailand to have reached their final positions less than 10 min after the earthquake, hence ruling out the hypothesis of a silent slow aseismic rupture.

India and Sunda plates motion and deformation along their boundary in Myanmar determined by GPS

Abstract: Using a regional GPS data set including ?190 stations in Asia, from Nepal to eastern Indonesia and spanning 11 years, we update the present?day relative motion between the Indian and Sundaland plates and discuss the deformation taking place between them in Myanmar. Revisiting measurements acquired on the Main Boundary Thrust in Nepal, it appears that points in southern Nepal exhibit negligible deformation with respect to mainland India. Including these points, using a longer time span than previous studies, and making an accurate geodetic mapping in the newest reference frame allows us to refine the present?day Indian motion. Our results confirm that the current motion of India is slower than predicted by the NUVEL?1A model, and in addition our India?Eurasia motion is significantly (?5 mm/yr) slower than previous geodetic determinations. This new Indian motion, combined with a refined determination of the Sundaland motion, gives way to a relative India?Sunda angular velocity of 20.2°N, 26.1°E, 0.370°/Myr in ITRF2000, predicting a relative motion of 35 mm/yr oriented N10° at the latitude of Myanmar. There, the Sagaing Fault accommodates only 18 mm/yr of right?lateral strike slip, only half of the shear component of motion. We present two models addressing how and where the remaining deformation may occur. A first model of distributed deformation implies convergence on the Arakan subduction (the northern continuation of the now famous Sumatra?Andaman Trench) and wrench faulting in the Arakan wedge. The second model uses localized deformation, where deformation observed west of the Sagaing Fault is entirely due to elastic loading on a faster and oblique Arakan subduction (23 mm/yr). This latter model predicts that a major earthquake of Mw 8.5 may occur every century on this segment of the subduction.

Microblock rotations and fault coupling in SE Asia triple junction (Sulawesi, Indonesia) from GPS and earthquake slip vector data

Abstract: The island of Sulawesi, eastern Indonesia, is located within the triple junction of the Australian, Philippine, and Sunda plates and accommodates the convergence of continental fragments with the Sunda margin. We quantify the kinematics of Sulawesi by modeling GPS velocities and earthquake slip vectors as a combination of rigid block rotations and elastic deformation around faults. We find that the deformation can be reasonably described by a small number of rapidly rotating crustal blocks. Relative to the Sunda Plate, the southwestern part of Sulawesi (Makassar Block) rotates anticlockwise at ?1.4°/Myr. The northeastern part of Sulawesi, the Bangai?Sula domain, comprises three blocks: the central North Sula Block moves toward the NNW and rotates clockwise at ?2.5°/Myr, the northeastern Manado Block rotates clockwise at ?3°/Myr about a nearby axis, and East Sulawesi is pinched between the North Sula and Makassar blocks. Along the boundary between the Makassar Block and the Sunda Plate, GPS measurements suggest that the trench accommodates ?15 mm/yr of slip within the Makassar Strait with current elastic strain accumulation. The tectonic boundary between North Sula and Manado blocks is the Gorontalo Fault, moving right laterally at about 11 mm/yr and accumulating elastic strain. The 42 mm/yr relative motion between North Sula and Makassar blocks is accommodated on the Palu?Koro left?lateral strike?slip fault zone. The data also indicate a pull?apart structure in Palu area, where the fault shows a transtensive motion and may have a complex geometry involving several active strands. Sulawesi provides a primary example of how collision can be accommodated by crustal block rotation instead of mountain building.

Geologically current plate motions

Abstract: SUMMARY We describe best-fitting angular velocities and MORVEL, a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface. Seafloor spreading rates and fault azimuths are used to determine the motions of 19 plates bordered by mid-ocean ridges, including all the major plates. Six smaller plates with little or no connection to the mid-ocean ridges are linked to MORVEL with GPS station velocities and azimuthal data. By design, almost no kinematic information is exchanged between the geologically determined and geodetically constrained subsets of the global circuit—MORVEL thus averages motion over geological intervals for all the major plates. Plate geometry changes relative to NUVEL-1A include the incorporation of Nubia, Lwandle and Somalia plates for the former Africa plate, Capricorn, Australia and Macquarie plates for the former Australia plate, and Sur and South America plates for the former South America plate. MORVEL also includes Amur, Philippine Sea, Sundaland and Yangtze plates, making it more useful than NUVEL-1A for studies of deformation in Asia and the western Pacific. Seafloor spreading rates are estimated over the past 0.78 Myr for intermediate and fast spreading centres and since 3.16 Ma for slow and ultraslow spreading centres. Rates are adjusted downward by 0.6–2.6 mm yr−1 to compensate for the several kilometre width of magnetic reversal zones. Nearly all the NUVEL-1A angular velocities differ significantly from the MORVEL angular velocities. The many new data, revised plate geometries, and correction for outward displacement thus significantly modify our knowledge of geologically current plate motions. MORVEL indicates significantly slower 0.78-Myr-average motion across the Nazca–Antarctic and Nazca–Pacific boundaries than does NUVEL-1A, consistent with a progressive slowdown in the eastward component of Nazca plate motion since 3.16 Ma. It also indicates that motions across the Caribbean–North America and Caribbean–South America plate boundaries are twice as fast as given by NUVEL-1A. Summed, least-squares differences between angular velocities estimated from GPS and those for MORVEL, NUVEL-1 and NUVEL-1A are, respectively, 260 per cent larger for NUVEL-1 and 50 per cent larger for NUVEL-1A than for MORVEL, suggesting that MORVEL more accurately describes historically current plate motions. Significant differences between geological and GPS estimates of Nazca plate motion and Arabia–Eurasia and India–Eurasia motion are reduced but not eliminated when using MORVEL instead of NUVEL-1A, possibly indicating that changes have occurred in those plate motions since 3.16 Ma. The MORVEL and GPS estimates of Pacific–North America plate motion in western North America differ by only 2.6 ± 1.7 mm yr−1, ≈25 per cent smaller than for NUVEL-1A. The remaining difference for this plate pair, assuming there are no unrecognized systematic errors and no measurable change in Pacific–North America motion over the past 1–3 Myr, indicates deformation of one or more plates in the global circuit. Tests for closure of six three-plate circuits indicate that two, Pacific–Cocos–Nazca and Sur–Nubia–Antarctic, fail closure, with respective linear velocities of non-closure of 14 ± 5 and 3 ± 1 mm yr−1 (95 per cent confidence limits) at their triple junctions. We conclude that the rigid plate approximation continues to be tremendously useful, but—absent any unrecognized systematic errors—the plates deform measurably, possibly by thermal contraction and wide plate boundaries with deformation rates near or beneath the level of noise in plate kinematic data.

Global trends in wind speed and wave height over the past 25 years

Abstract: Wind speeds over the world’s oceans have increased over the past two decades, as have wave heights. Studies of climate change typically consider measurements or predictions of temperature over extended periods of time. Climate, however, is much more than temperature. Over the oceans, changes in wind speed and the surface gravity waves generated by such winds play an important role. We used a 23-year database of calibrated and validated satellite altimeter measurements to investigate global changes in oceanic wind speed and wave height over this period. We find a general global trend of increasing values of wind speed and, to a lesser degree, wave height, over this period. The rate of increase is greater for extreme events as compared to the mean condition.

Global variations of stress drop for moderate to large earthquakes

Abstract: [1] We investigate the global variation of earthquake stress drops using spectra of about 2000 events of mb ≥ 5.5 between 1990 and 2007. We use an iterative least squares method to isolate source displacement spectra from travel path and receiver contributions, based on a convolutional model. The observed P wave source spectra are corrected with a globally averaged empirical correction spectrum and estimates of near-source attenuation. Assuming a Brune-type source model, we estimate corner frequencies and compute stress drops. Stress drop estimates for individual earthquakes range from about 0.3 to 50 MPa, but the median stress drop of about 4 MPa does not vary with moment, implying earthquake self-similarity over the Mw = 5.2 to 8.3 range of our data. A comparison of our results with previous studies confirms this observation over most of the instrumentally observable magnitude range. While the absolute values of our estimated stress drops depend upon the assumed source model, we identify relative regional variations of stress drop that are robust with respect to the processing parameters and modeling assumptions, which includes an inherent assumption of constant rupture velocity. We find a dependence of median stress drop on focal mechanism, with a factor of 3–5 times higher stress drops for strike-slip earthquakes and also find a factor of 2 times higher stress drops for intraplate earthquakes compared to interplate earthquakes.

A geodetic plate motion and Global Strain Rate Model

Abstract: We present a new global model of plate motions and strain rates in plate boundary zones constrained by horizontal geodetic velocities. This Global Strain Rate Model (GSRM v.2.1) is a vast improvement over its predecessor both in terms of amount of data input as in an increase in spatial model resolution by factor of ∼2.5 in areas with dense data coverage. We determined 6739 velocities from time series of (mostly) continuous GPS measurements; i.e., by far the largest global velocity solution to date. We transformed 15,772 velocities from 233 (mostly) published studies onto our core solution to obtain 22,511 velocities in the same reference frame. Care is taken to not use velocities from stations (or time periods) that are affected by transient phenomena; i.e., this data set consists of velocities best representing the interseismic plate velocity. About 14% of the Earth is allowed to deform in 145,086 deforming grid cells (0.25° longitude by 0.2° latitude in dimension). The remainder of the Earth's surface is modeled as rigid spherical caps representing 50 tectonic plates. For 36 plates we present new GPS-derived angular velocities. For all the plates that can be compared with the most recent geologic plate motion model, we find that the difference in angular velocity is significant. The rigid-body rotations are used as boundary conditions in the strain rate calculations. The strain rate field is modeled using the Haines and Holt method, which uses splines to obtain an self-consistent interpolated velocity gradient tensor field, from which strain rates, vorticity rates, and expected velocities are derived. We also present expected faulting orientations in areas with significant vorticity, and update the no-net rotation reference frame associated with our global velocity gradient field. Finally, we present a global map of recurrence times for Mw=7.5 characteristic earthquakes.