I will be filling this in over the next few days and wanted to start collating social media materials for this event.
There was a large volcanic eruption in the Tonga region. This eruption was observable from satellites and has generated a modest but observable tsunami from Australia to the United States.
This event is still unfolding and it will take months until we have a deeper understanding of the causes for the tsunami. We know it is related to the explosive volcanic eruption from Hunga Tonga-Hunga Ha’apai, about 55 kms (35 miles) northwest of the largest island of the Kingdom of Tonga, Tongatapu.
I will continue to fill in details. I am currently busy trying to manage our tsunami event response and am learning lots in the process. However, this delays my time available here.
This week, CGS sent teams to various harbors & beaches on the California coast to collect measurements, photos, & videos documenting the effects of the Jan. 15 #tsunami from Tonga. This info helps us understand how future tsunami might impact our coastal harbors & communities. pic.twitter.com/xGa8zmNNNs
After two weeks of work I can finally share my 3D reconstruction of the gigantic ash cloud from the January 15 Hunga #Tonga-Hunga-Ha'apai #eruption
Parts likely reached *close to 60km* according to my reconstruction, that's beyond the stratosphere and inside the mesosphere!🧵 pic.twitter.com/qMzPSjZj7P
It is incredible to see how the #Andes bounced back part of the #Tonga atmospheric Lamb wave on its first cross over South America!. Here depicted with signal processed IR data from the GOES 16 geostationary satellite pic.twitter.com/OLM0MD0neO
Below there are many tweets etc. and one may feel like they are scrolling forever. These tweets are loosely organized into several sections.
Background Material
Tsunami Notifications
Tsunami Education
Tsunami Observations
Tsunami Modeling
Volcano Eruption Observations
Fascinating Observations
Background Material
Excellent thread on how you can best help disaster relief efforts (and why) from someone who knows from professional experience. Donate $ (not 👠s) to reputable relief organizations. This is true for most disaster relief. https://t.co/zj6AIUuCZC
This is a truly excellent short article by @LoriDengler about yesterday’s #tsunami and eruption event. It makes it clear just how unique this was, and why the @NWS_NTWC folks had to improvise (brilliantly) to get the alert out. https://t.co/meypb1w2k9
This is a small volcanic island but below the ocean the volcano is huge at around 1.8 km high and 20km wide. Much of the 2014-15 cone is now gone. Read more about what we know by @scronin70: https://t.co/nMNQYGcLDwpic.twitter.com/Ejn3z6e1I4
Scientists said the volcano had been puffing away for about a month before rising magma, superheated to around 1,000 degrees Celsius (1,832 Fahrenheit), met with 20 degree (68 Fahrenheit) seawater on Saturday, causing an instantaneous and massive explosion pic.twitter.com/iNVC2tB3XM
An underwater volcano off Tonga erupted on Jan. 15, triggering tsunami warnings and evacuation orders in Japan and causing large waves in several South Pacific islands, where footage on social media showed waves crashing into coastal homes pic.twitter.com/L7uzK59jG7
Observations from Japan on why the far-field tsunami was likely triggered by air-sea coupling, not the standard shallow-water wave propagation from the source. https://t.co/70O5V095Xs
New preprint out on ESSOAr – "Tonga eruption triggered waves propagating globally from surface to edge of space" – in which we look at the deeply unusual atmospheric waves generated by the recent Hunga Tonga eurption. pic.twitter.com/IJfMgb1Bt9
*** California is in a NWS Tsunami "Advisory" – "If a tsunami does impact California, it is unlikely it will be a large tsunami but possibly in the Advisory range (0.3m to 1m) and it could arrive at about 0700AM Pacific time according to the NTWC." https://t.co/ioUmU0Yrd3https://t.co/0b8aE12CV5
#TsunamiAdvisory remains in effect for the #WestCoast, with peak tsunami wave heights of 1 to 3 feet possible. Main impacts expect to be strong rip currents and coastal flooding of some immediate coastal low-lying areas. #CAwx
#TsunamiAdvisory remains in effect for the #WestCoast, with peak tsunami wave heights of 1 to 2 feet possible. Main impacts expect to be strong rip currents, coastal flooding, and inundation of low lying areas is possible. Move to higher ground. #CAwx
A Tsunami Advisory has been issued for Humboldt Co. due to volcanic activity in the S. Pacific. A tsunami capable of producing strong currents that may be hazardous to swimmers, boats and coastal structures is expected beginning 7:30-8 a.m. Widespread inundation is not expected.
#SF remains in a Tsunami Advisory Strong, dangerous currents WILL be present. Stay out of water and away from coastal areas. pic.twitter.com/nMYm0wxvPA
— SAN FRANCISCO FIRE DEPARTMENT MEDIA (@SFFDPIO) January 15, 2022
We continue to monitor the #Tsunami closely. Please listen to local advisories — stay safe and off beaches that are under a warning. https://t.co/M1QdOhBleG
There has been another eruption from the volcano in Tonga. However, based on the latest information, there is NO TSUNAMI THREAT for American Samoa at this time. Because this volcano remains active, please stay tuned for further updates.
— National Weather Service WSO Pago Pago (@NWSPagoPago) January 16, 2022
Tonga Volcanic Eruption and Tsunami: World Bank Disaster Assessment Report Estimates Damages at US$90M https://t.co/bx6Iq5VTyb via @WorldBank
— Jason "Jay" R. Patton (@patton_cascadia) March 7, 2022
A #Tsunami Advisory means: a dangerous wave is on the way. Strong and unusual currents are expected along the coast, and in bays, marinas, and harbors. Move to high ground and away from the shore. More at https://t.co/npoUHxEZLS. pic.twitter.com/MCLDdN9qPp
yes, don't forget to check out the #Tsunami Hazard Areas (where people might want to evacuate from during a tsunami event) for #California are available at https://t.co/9LPXadKWEi
There are a number of #Tsunami HAZARDS that could directly impact harbors and boaters: > Strong and unpredictable currents, especially where there are narrow entrances, narrow openings, and other narrow parts of harbor.
If you're having a hard time interpreting news about this tsunami advisory — what the danger is and where — @DaveSnider breaks it down really well in this video on the warning center's facebook page.https://t.co/tHsstcS5xl
The Hunga Tonga-Hunga Ha'apai volcanic eruption was heard here in Alaska starting around 3:30 a.m. – 6,000 miles from the volcano! Infrasound measurements from the @alaska_avo confirm that it was indeed coincident with the volcanic pressure wave. Special thanks to Dr. David Fee. pic.twitter.com/Wp4tnwiaud
2. Tsunamis are not one wave. It's more like sloshing and that sloshing can continue for a day. Just because the first wave has passed, it is not time to go see the beach. 3. Much tsunami damage happens in ports because of the currents. Moving water has huge momentum. [2/2]
Wondering why you didn't get an alert about the #Tonga tsunami and the #TsunamiAdvisory in WA this morning? There are many ways to get tsunami alerts and it's best to be signed up for many kinds! Learn more about tsunami alerts (and other hazard alerts) at https://t.co/iU0UZFRnC2pic.twitter.com/aRkJ2CKB32
After checking with @LoriDengler, I think that today is the first-ever Pacific-wide #tsunami alert from a volcanic eruption. Interestingly, the warning center needed an earthquake magnitude to issue the bulletin. So they made one up: mag 0.1! (Now changed to mag 1) pic.twitter.com/znjfUa1o3S
This was not your run of the mill way a #tsunami is generated. Most often its from plates shifting abruptly on the sea floor. This tsunami was essentially caused by a massive underwater explosion of molten rock and lava that displaced the water above it. pic.twitter.com/JyLwZYHcRp
Not sure I have ever seen this before. My @noaaocean colleague just flagged this. The pressure wave caused by the #TongaVolcano is also actually causing a tsunami – in this case a meteotsunami of about 10cm in Puerto Rico. Wild. pic.twitter.com/EapEuNhjB5
Fascinating part of the tsunami warning process today: This wasn't an earthquake, so the tsunami warnings were sent out with a default magnitude of 1. Shows that our whole system is set up for earthquake tsunamis. 1/2 pic.twitter.com/pwsBuicyyX
Great questions here about hydrophones as a tool to monitor submarine volcanism. This is my favorite topic, so @syabilazriAS is going to get a longer answer than maybe was desired. 😆🧵 https://t.co/VcbgmUEnZn
High tide is around 9:30 AM along the Northwest California coast. Even though the initial wave is expected 7:30 AM to 8 AM, additional waves and strong currents will continue. As much as 3 feet of tsunami wave is expected. Stay back from the beaches and lowest-lying areas. #CAwxpic.twitter.com/BmoCVtJc2f
Monterey Tide Station #Tsunami The first surge may not be the largest. In other locations today, the largest surge came much later than the first arrival time. In Monterey, CA it may have 2.5 feet above high-tide conditions, similar to a King Tide event. pic.twitter.com/3h7cpXiPZr
Tsunami energy arriving at Del Monte Beach. Wave run up onto the beach is impressive. I'm safely distanced, but I saw other folks have to scramble when the waves unexpectedly reached them. Follow @NWSBayArea for #tsunami safety info. pic.twitter.com/JRl2XsTftI
Seeing some surges on the Port San Luis tsunami gauge. Reporting up to a 24 cm residual so far. That's 9.4 inches or about 19 inches from the bottom and top of the residual. More at https://t.co/SGd8WQoeji. #tsunamipic.twitter.com/bKYRRXuW4W
Made it to the high ground edge (6 feet+ above high water level) at Meeker Slough meets the SF Bay. Another observer tells me she has seen two surges already. Check the bridge piling for the high water mark from an earlier surge. pic.twitter.com/Vv8O00SLBF
So far, business as usual at the Ferry Building, which sits above the expected tsunami surge. High tide in SF is at 9:25. It is probably just my imagination that the bay looks more unsettled this morning. pic.twitter.com/GkkSUSbBgY
10 AM | Here are the latest observed tsunami wave heights from along the West Coast of the US. Generally tsunami wave heights have been around 1 foot or less along the Oregon and Washington coast. (1/2) pic.twitter.com/r94If9ODXM
High tide occurred through the morning, and total water levels are decreasing, lowering risk of coastal flooding, but rapid fluctuating surges of water onto and off of the coast/strong currents will continue thru the remainder of the day. pic.twitter.com/5ylsGGvr0m
Here are the latest maximum observed wave heights. Highlights for our area include 3.7 feet at both Crescent City and Arena Cove. The advisory continues for the West Coast, stay tuned to the latest updates on the advisory from @NWS_NTWCpic.twitter.com/UYwFWWLGWd
Tsunami surge up the mouth of San Luis Creek at Avila Beach (Port San Luis) #tsunami (L) Sat. morning Jan 15 (R) Thurs. afternoon Jan 20 pic.twitter.com/XAAnzYK8GH
Tide gage in American Samoa registering the tsunami created by the Tonga submarine volcanic eruption. Peak amplitude so far is 0.74 meter (~2.5 feet). pic.twitter.com/ATsmnP5clg
~ 2 Meter #Tsunami in Nuku'alofa, Tonga nach der erneuten heftigen Explosion des Hunga Tonga-Hunga Ha'apai Vulkans. Die ersten Wellen haben auch Fidschi erreicht. Ich hoffe, die seit gestern bestehende Tsunami-Warnung wurde ernst genommen https://t.co/akfdQqtrP7pic.twitter.com/b90hbS4oTW
While a 1m wave might not sound big, tsunami have much longer periods (the time between each wave) than wind waves so even a 1m wave can cause significant damage and flooding! pic.twitter.com/27DoahlJEb
Muy buena foto !! Interesante comparar los efectos con lo que observa el mareógrafo. Este indica 40 a 60cm de desviación respecto de valores medios. El tsunami sigue en desarrollo. Gracias !! pic.twitter.com/f4WzPs2sMn
Timelapse video of the #tsunami taken at Mogareeka inlet at 7-7:20am this morning. Tide is rising (flow left to right) but here is what happens as the waves come through. Mogareeka is usually very flat so the effects are amplified @anuearthscience @ourANU pic.twitter.com/y00Tj1iFb3
Urgente Marejadas destruyeron otro muelle en Las Coloradas sector Isla del Rey comuna de Corral, los ríos. El oleaje arrasó con al menos dos muelles uno que aparentemente sería de acceso público y otro privado, se mantiene alerta de tsunami en la zona.pic.twitter.com/ChYMvIM2wr
A #tsunami is occurring. Tsunami Advisories have been ended for portions of Southcentral and Southeast Alaska- they continue elsewhere. See https://t.co/npoUHxWBas for the latest. Port San Luis, CA: 4.3 ft King Cove, AK: 3.3 ft Crescent City, CA: 3.7 ft Point Reyes, CA: 2.9 ft pic.twitter.com/HeZJldZlxZ
https://t.co/Xw7BUtepv4 article with early news about what is happening inside Tonga. Communications are difficult. Tonga volcano: Photos and video give first glimpse of tsunami's impact https://t.co/j5LDZ3B6lP
1/n Lot of damages and very complicated aftermath in Tonga islands. Hopefully aerial views by New Zealand Defence Force suggest #tsunami height of only few meters and limited inundation distance. Here Nomuka island 70km NE of Hunga Tonga volcano. Google earth 2016 image follows pic.twitter.com/bUgwtsJ7f8
Offical announcement from the Government of Tonga. Tragically, 3 people confirmed to have died.
Also, while there are many satellite images coming out showing the awful eruption impact, please be compassionate and considerate to those impacted or awaiting news of loved ones. https://t.co/Vw2SfHEHjh
⬇️The effects of the explosion on Nomuka Island are visible when comparing the #Sentinel2 🇪🇺🛰️images of ↙️18 Dec. 2021 ↘️17 Jan. pic.twitter.com/rtrWOXgpC4
🇹🇴The miracle survival of a 57-year-old disabled man who survived in the ocean for 27 hours after being swept away by a tsunami wave is one of the first astonishing accounts to emerge from Tonga
This map, based on #alos2#sar satellite data, shows the devastating impact of the huge #Tonga volcanic eruption. Damaged areas are marked by yellow to red pixels, with red indicating the most damage. More info, GeoTIFF and KMZ files at https://t.co/uxKeZHJ0ixpic.twitter.com/2xSSZeYVbH
Quite the day at the beach sampling the January 2022 Hunga tsunami deposits- if you squint there are ?3 slightly different grainsize layers here – plus thin black layer of ash near the top, western Tongatapu pic.twitter.com/oZcEPGxeO6
Extracting oriented tubes of tsunami deposits will help colleagues figure out deposition histories using 3D magnetic and textural properties- Hunga tsunami deposits – western Tongatapu pic.twitter.com/ojcoi9kiSP
Tonga tsunami, quite impressive that it propagated all the way to stations in the coast of Mexico (Manzanillo over 1.4m) and California… Here is a order one attempt to model this tsunami @geosmx#geoclawpic.twitter.com/drzz5GMHJK
Truly incredible imagery from the Himawari meso sector of multiple shockwaves from a volcanic eruption on Tonga propagating through what was recently subtropical storm Cody. My back of the napkin math estimates put the shockwave speed somewhere in the vicinity of 500-600 mph. pic.twitter.com/lDyNXFpLbE
This volcano eruption is producing some sights we won't see on satellite for a while after. The explosive updraft powers right into the stratosphere, and the warming with height there makes the blob appear warm/shallow. There's also the shockwave and condensation pushing NE pic.twitter.com/6ud79uws9w
1.14.2021: Large volcanic eruption near Tonga (Hunga Tonga-Hunga Ha'apai volcano) today as seen from outer space. Shown on visible imagery using the Himawari satellite. #hiwx#tsunami#earthquakepic.twitter.com/zOTj6Qu1Wv
Fantastic #infrasound signals on the Australian IMS arrays from the #Tonga#volcano. Data courtesy of @GeoscienceAus and made available via @IRIS_EPO – these atmospheric waves travel much slower than seismic waves and over an hour separates the signals on mainland Australia. pic.twitter.com/7A8KUZdA22
The pressure wave from the Hunga Tonga-Hunga Ha'apai eruption arrived here in Anchorage at 3:30 a.m. AKST. This is exactly 7 hours after the eruption. The volcano is 5,820 miles away (9,360 km). That means it travelled at 830 mph (1,340 kmh). pic.twitter.com/R3rgzAbo6r
Seismogram from Monasavu, Fiji ~800km NW of Hunga Tonga-Hunga Ha'apai eruption. Assuming main eruption at ~04h15mUTC, shows P waves, seismic surface waves (Lq, Lr) and oceanic SOFAR acoustic waves (T). Signal from the eruption continues for 2+ hours.https://t.co/etU65z1wyqpic.twitter.com/2PublssOYm
A number of folks have posted similar obs, but here are two pressure traces showing the #HungaTongaHungaHaapai shock wave observed at UNR in Reno and UoU in Salt Lake City this morning. The SLC data seem to show subsequent oscillations (sloshing) in the valley cold pool. #UTwxpic.twitter.com/ekDRXUUrq8
Here's the latest eruption. Again ice detectable. This time ash signal appearing. Plume appears to be significantly stronger than the last one. I've set the colour bar and scale to try to highlight the tropospheric and stratospheric portions of the plume (top left plot). pic.twitter.com/zRaFvG2jvJ
You can see the atmospheric shock wave ripple out Pacific-wide after the initial #eruption in #Tonga (quite hard to see but it crosses #Hawaii and #Australia).
Una reciente imagen tomada por los satélites de #SkySat propiedad de Planet Labs. muestra que el corredor de tierra que existía desde 2014 entre las islas desapareció después de las erupciones cataclísmicas de ayer. pic.twitter.com/2fM8rwNPZS
That map, created with April 2016 data from the R/V Falkor, showed details of the overall volcanic edifice. Measured on this quick Google Earth overlay, the diameter of the caldera rim is ~6 km. Recent eruptions have been on the N and NE sides. Red * is 2009 vent. pic.twitter.com/oPdKMR44G3
Jumping on board pressure perturbation Twitter to share this animation of @okmesonet pressure data. Several pressure waves created by the Hunga Tonga–Hunga Ha'apai volcano eruption in Tonga passed from southwest to northeast across Oklahoma between 7 and 9 AM CST this morning. pic.twitter.com/aBpRXNbNeX
Tonga tsunami is arriving at BC! Last twelve hours of seafloor pressure data show: Series of waves, starting to arrive at @Ocean_Networks stations at 8:22 PST, and at the west coast about 9 AM PST. Height about 5 cm offshore. Coastal currents probably a bigger issue. #BCTsunamipic.twitter.com/VpEzXymHB3
Ionospheric total electron content (TEC) perturbations derived from a GNSS site on Samoa from the Hunga Tonga eruption were not small to say the least @IGSorgpic.twitter.com/bMa8MKCZ3o
15 minute pressure altimeter change via ASOS NWS/MADIS 5 minute interval data. Shows the shockwave from the #Tongaeruption , feel free to use as you wish. pic.twitter.com/P31Aq1SYku
Longwave infrared channel via #GOESWest of the #Tonga eruption.. one of the most incredible satellite animations I've ever seen. The relative warmth of the ash cloud atop the very cold tropospheric convective anvil. Waves upon waves. Simply incredible. pic.twitter.com/MoBcIxkblW
New data alert 🚨 We just overflew Hunga Tonga-Hunga Ha'apai volcano with #Sentinel1 🛰️🌋 The datatake didn't include #Tonga main island, where I desperately hope everyone is safe 🇹🇴❤️🩹Here's our last 3 passes over the volcano… pic.twitter.com/KD39030U5S
The evolution of the volcanic island of Hunga Tonga over time, with the last image having been taken just two hours before the massive eruption last night. Curious to see the scene after that… Images from Google Earth and @planet. pic.twitter.com/lOmca4Du7I
Looks to me like we see the seismic signal from the Tonga eruption at Weston, MA and Westport, CT. The spectrum shows signal at <0.06 Hz (>17 sec) & at 0.1 to 0.2 Hz (5 to 10 sec), which fits for ocean & continental Rayleigh waves.@jpulli@stevenjgibbons@KaseyAderholdpic.twitter.com/hJZJvO59PM
Before and after photos show that the island of Hunga Tonga and Hunga Ha'apai is essentially gone following the explosive eruption of the volcano last night. https://t.co/6GvI5nNGV2pic.twitter.com/wEjBhfmWFZ
A lot of talk about just how big the eruption at Hunga Tonga-Hunga Ha'apai was. It might be awhile before we know & we don't know if there is more to come. Questions abound about what caused the tsunami, why the eruption was so explosive, etc https://t.co/oXas1XevjC@DiscoverMag
With latest satellite imagery, we get a step closer to understand what happen with Hunga #Tonga leading to this ocean-wide #tsunami.
Latest @sentinel_hub imagery shows the loss of a majority of the volcano's emerged landmass. However, most of its structure lies underwater. pic.twitter.com/odz5VcNphl
Atmospheric wave response to Tonga eruption, from 4 UTC to 10:50 UTC. Slightly smoothed 10-minute change in GOES-17 band 13 (IR). Looks like some modest filtering would pull out a really clear signal. pic.twitter.com/CHZY7iv4HH
There are a lot of questions about VEI (Volcano Explosivity Index). I recognize that the want to compare this eruption is there, but there is so much information that we simply do not have. More important is how the people of Tonga are and is this over? We do not know. https://t.co/YiJ4vahkPu
Pressure wave #2 passed here last night at about 11 p.m. This is the wave travelling from the opposite direction. Still had an impressive magnitude. pic.twitter.com/1fG4aKhWqs
Trying to understand why the weather stations at Stornaway on the Outer Hebrides measured the blast before us. Basically because the blast came from the north. Unbelievably the shortest distance between here and the South Pacific is over the North Pole. 🤯 #TongaVolcanoEruptionpic.twitter.com/6MdZVJ4VWI
It appears that there is some minor activity ongoing at Hunga Tonga-Hunga Ha'apai, as expected. This is based purely on satellite data. No volcanic lightning detected.https://t.co/13uhR0353u
Why does it take awhile to get satellite images of what's going on in #Tonga?
We have to wait for satellites to fly over, or redirect them. Also, they "see" in different wavelengths, so night, clouds, and ash can obscure the view.
In Germany two main air pressure waves from the #Tonga eruption could be detection: The first wave traveled from north to south, while the second wave moved from south to north. The reason might be explained by the animation below, where I visualized an outgoing circular wave… pic.twitter.com/B57uRyy3ik
A news story about Tonga. Most of the news in the U.S. is very American-centric, so you have to poke around the int’l news scene to get any news about Tonga.https://t.co/GXZlAIlsrQ
Displacements measured at GPS/GNSS station TONG in Tonga about 70 km from #HungaTongaHungaHaapai shows large motion over about 10 minutes that returns close to previous position. Some kind of shock wave or seismic wave, probably. https://t.co/7Bpq2U8GsN
Copernicus Sentinel-1A radar imaged #HungaTongaHungaHaapai on 15 January 2022 after major eruption. Most of two islands and entire new cone was blown away, along with reef south of underwater caldera. @googleearth Engine HV radar polarization animation Aug-Jan by @TheHandwergerpic.twitter.com/gwnA52Q12e
So the Hunga Tonga and Hunga Ha'apai were two separate islands before an eruption in 2014-15, and they split up before the violent eruption last Saturday. And now, very little of the two islands are left. pic.twitter.com/5G7Zy3n0td
SNPP/OMPS limb-profiler (OMPS-LP) aerosol vertical profiles from Jan 16 shown below captured the stratospheric #volcanic aerosol cloud reaching altitudes up to ~30 km (in same location as the highest SO₂ columns). h/t @NASAGoddard Ozone & Air Quality teamhttps://t.co/PBUPJCgEtJpic.twitter.com/3eNIwouCfs
Don't know if anyone has plottet something similar already. Here is the pressure wave from the #HungaTongaHungaHaapai eruption in the infrasound section. The wave taking the path the other way around Earth is also clearly visible a few hours later. Dashes: speed of sound#Tongapic.twitter.com/ukkm1AXZCx
A THIRD pressure anomaly associated with the #TongaVolcano passed through #Miami on Sunday evening… the timing means that it was the first wave making a full trip around the globe! Absolutely mind-blowing power. pic.twitter.com/lpW9FY97Mw
Sure looks like we got a fourth passage of the Hunga Tonga – Hunga Ha'apai eruption shockwave in Utah. Timing is spot-on and signal is similar to previous passages. Interesting for sure but worried for the people of Tonga. pic.twitter.com/ZzZzTrIJcH
Revisiting the Tonga volcanic shockwave: Here's the latest Eureka barograph showing the 1st shockwave, another distinctive shockwave just after midnight Sat night (the other side of the initial shockwave), & another possible shockwave just after noon today. #CAwx#Tongaeruptionpic.twitter.com/c4bXjPcZrY
The @EOS_SG blog post on the Hunga Tonga-Hunga Ha'apai eruption is up! Featuring the infrasound signal from the eruption as recorded in Singapore https://t.co/YMD56WZpqY
I took a quick dive into the science of the volcanic eruption in Tonga. It was an extraordinary event that will keep researchers busy for a while… #TongaVolcanohttps://t.co/K0K3aQUv3s
Pressure waves from #HungaTongaHungaHaapai have travelled 3 times around the globe as of this morning. This pressure graph is from Iceland. It shows 6 peaks rather than just 3 because.. 1/n pic.twitter.com/W4WMOYswSV
(1/4) #CTBTO continues to analyze data from the Hunga Tonga-Hunga Ha'apai volcanic eruption. In terms of infrasound technology only, this is the largest event ever recorded by the #IMS infrasound network; much larger than the Chelyabinsk meteor in 2013. pic.twitter.com/T7y9Nk9Rhm
The low-frequency signal from the Hunga Tonga-Hunga Ha'apai vulcanic eruption that generated the tsunami (15/01) was well captured by >400 #BMKG broadband seismic stations over the Indonesia region. The vertical record from five selected stations (bandpass filtered 0.01-0.05 Hz): pic.twitter.com/k2D3M6zClL
Island nation of #Tonga is completely offline following a #tsunami triggered by a massive volcanic eruption in the Pacific Ocean.
According to @kentikinc data, traffic volumes began to drop around 4:30 UTC (5:30pm local) before finally going to zero at 5:40 UTC (6:40pm local). pic.twitter.com/g4QZilBrd5
Shockwave from Hunga Tonga-Hunga Haʻapai eruption plume, seen by pressure change at UK sites. Wave moves southward down the country 18-20Z 15Jan. The same wave, but travelling the other way around the globe, moves northward up the country 01-03Z 16Jan. @RoostWeather@Silkstinihopic.twitter.com/2jXaWwyzih
"The waves are red because of airglow, an aurora-like phenomenon caused by chemical reactions in the upper atmosphere. Airglow is usually too faint to see, but gravity waves from the volcano boosted the reaction rates." https://t.co/nntTFiBOrEhttps://t.co/qWinu9InSV
Tonga volcano : This is shockwave as measured with the IASI satellite mission (temperature perturbation between the day of the eruption and the next day). First time we see this ! pic.twitter.com/7uTCwP3tNP
Oopps! When a volcano erupts with such intensity in one part of the world and causes an #oilspill 10.000 km away. The Mare Doricum vessel was offloading at La Pampilla (Peru) refinery when the tsunami hit. Now, one of the worst oil disasters in the region. https://t.co/7nBNtWgX01pic.twitter.com/MVPFGyIFWe
In the middle of the night (my time) I got a notification from the EMSC earthquake notification service. I encourage everyone to download and use this app.
There was an intermediate depth magnitude M 7.5 earthquake in Peru. The tectonics in this region of the world are dominated by the convergent plate boundary, a subduction zone formed by the convergence of the oceanic Nazca and continental South America plates.
As the Nazca plate subducts, it dips below the South America plate at different dip angles. In this region of Peru, the dip angle is shallow and we term this flat-slab subduction.
This M 7.5 earthquake occurred in the downgoing Nazca plate, so was not a subduction zone megathrust event, but a “slab” event (for being in the Nazca slab).
I prepared a much more extensive report for a M 8.0 earthquake in a nearby location that happened on 26 May 2019. Read more about the tectonics of this region in that report here.
Was this M 7.5 an aftershock of the M 8.0? Probably not, based on the USGS M 8.0 slip model. However this M 7.5 could have been triggered by changes in static coulomb stress following the M 8.0.
I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.
Below is my interpretive poster for this earthquake
I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
In the upper left corner is a large scale plate tectonic map showing the major plate boundary faults.
In the lower left center is a map showing how the Nazca slab is configured in different locations (Ramos and Folguera, 2009).
In the left center is a cross section showing seismicity in this region (Kirby et al., 1995). The source area for this plot is designated by a dashed yellow box on the map.
In the upper right corner is a pair of maps that show the landslide probability (left) and the liquefaction susceptibility (right) for this M 7.5 earthquake. I spend more time describing these types of data here. Read more about these maps here.
In the lower right corner I plot the USGS modeled intensity (Modified Mercalli Intensity scale, MMI) and the USGS “Did You Feel It?” observations (labeled in yellow). Above the map is a plot showing these same data plotted relative to distance from the earthquake. Read more about what these data sets are and what they represent in the report here.
Here is the map with 3 month’s seismicity plotted.
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
Specific References
Antonijevic, S.K., et a;l., 2015. The role of ridges in the formation and longevity of flat slabs in Nature, v. 524, p. 212-215, doi:10.1038/nature14648
Bishop, B.T., Beck, S.L., Zandt, G., Wagner, L., Long, M., Knezevic Antonijevic, S., Kumar, A., and Tavera, H., 2017, Causes and consequences of flat-slab subduction in southern Peru: Geosphere, v. 13, no. 5, p. 1392–1407, doi:10.1130/GES01440.1.
Rhea, S., Hayes, G., Villaseñor, A., Furlong, K.P., Tarr, A.C., and Benz, H.M., 2010. Seismicity of the earth 1900–2007, Nazca Plate and South America: U.S. Geological Survey Open-File Report 2010–1083-E, 1 sheet, scale 1:12,000,000.
Villegas-Lanza, J. C., M. Chlieh, O. Cavalié, H. Tavera, P. Baby, J. Chire-Chira, and J.-M. Nocquet (2016), Active tectonics of Peru: Heterogeneous interseismic coupling along the Nazca megathrust, rigid motion of the Peruvian Sliver, and Subandean shortening accommodation, J. Geophys. Res. Solid Earth, 121, 7371–7394, https://doi.org/10.1002/2016JB013080.
Wagner, L.S., and Okal, E.A., 2019. The Pucallpa Nest and its constraints on the geometry of the Peruvian Flat Slab in Tectonophysics, v. 762, p. 97-108, https://doi.org/10.1016/j.tecto.2019.04.021
Yepes,H., L. Audin, A. Alvarado, C. Beauval, J. Aguilar, Y. Font, and F. Cotton (2016), A new view for the geodynamics of Ecuador: Implication in seismogenic source definition and seismic hazard assessment, Tectonics, 35, 1249–1279, https://doi.org/10.1002/2015TC003941.
I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.
Last afternoon (my time) there was an M 7.0 earthquake near Acapulco, Mexico. This event generated a tsunami, landslides, building damage, casualties (one fatality as I write this), and many emotions. https://earthquake.usgs.gov/earthquakes/eventpage/us7000f93v/executive
I present my interpretive poster and a few figures. Read more about the tectonics of this region here, in a report for an M 7.4 earthquake in 2020.
Below is my interpretive poster for this earthquake
I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 7.0 in one version.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
In the upper left corner is a small scale map showing the major plate boundaries.
Below the plate tectonic map is a plot showing the tide gage data from Acapulco, Mexico. Note the clear tsunami signal.
To the right of the plate tectonic map is a large scale map showing aftershocks in the region of the M 7.1 mainshock. Note that these aftershocks are from the Servicio Sismológico Nacional (SSN) Catálogo de sismos and that there are two mainshock locations (USGS M 7.0 and SSN M 7.1).
In the lower right corner is a map that shows a comparison of earthquake intensity between the USGS models and the Did You Feel It observations.
Above the intensity comparison map is a plot showing these same data, intensity is on the vertical axis an distance from the earthquake [Hypocenter] is on the horizontal axis.
In the upper right corner is a map that shows the results of an earthquake induced liquefaction model. Read more about this model here.
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.
On 14 August ’21 there was a magnitude M 7.2 oblique strike-slip earthquake in Haiti. This earthquake was along the Enriquillo-Plantain Garden fault zone, which also ruptured in 2010. Here is my report for the 2010 Haiti earthquake (see more about the tectonics of this region of the world). https://earthquake.usgs.gov/earthquakes/eventpage/us6000f65h/executive
Below is my interpretive poster for this earthquake
I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 7.0 in one version.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
in the lower right corner is a small scale plate tectonic map.
Above that map is a plot showing the USGS finite fault slip model. This shows the location of the fault and color represents how much the fault slipped during the earthquake.
In the upper right corner is a map that compares the USGS earthquake intensity models (the contoured lines) with the USGS Did You Feel It? observations from real people.
In the lower center is a map that shows the aftershocks from the M 7.2 earthquake and from the 2010 M 7.0 earthquake.
In the upper left are two maps that show models of earthquake triggered landslides and earthquake induced liquefaction for this M 7.2 event. Read more about these models here.
Here is the map with 3 month’s seismicity plotted.
Earthquake Aftershocks
Below a map showing the aftershocks from the 2021 M 7.2 and 2010 M 7.0 Haiti earthquakes.
Potential for Ground Failure
Below are a series of maps that show the potential for landslides and liquefaction. These are all USGS data products.
There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:
FOS = Resisting Force / Driving Force
When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.
Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.
Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.
Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides.
Below are maps showing a comparison between the USGS modeled earthquake triggered landslides and liquefaction potential with the Centre Nationale De Information Géospatiale (CNIGS) probabilistic models of ground failure.
— ANU Earth Sciences 🌏 (@anuearthscience) July 14, 2022
Given the larger magnitude and farther west location of this 2021 Haiti quake relative to the 2010 earthquake, it is also worth noting the short time elapsed between the last historical sequence of large quakes on this fault. Fig from our 2012 paperhttps://t.co/zLQ8wSbtWrpic.twitter.com/ux17G6bz2F
Why was there a gap between the 2010 and 2021 Haiti earthquakes? Because a sequence of moderate quakes in 1860 released strain in the gap! https://t.co/rpuXfEw78G
I was excited to see & have a chance to comment on a study published in Science yesterday, discussed here. The response to the 2021 Nippes, Haiti, earthquake was very different from the response in 2010, 1/https://t.co/y6l3jjsAYQ
Saturday's M7.2 earthquake in Haiti was close to the 2010 M7.0 earthquake. Both events are devastating on their own but compounded by ongoing problems the region faces. Compare @IRIS_EPO's Teachable Moments:
— Southern California Earthquake Center (@SCEC) August 16, 2021
Here is a comparison of Peak Ground Acceleration (perceived shaking) for the January 2010 M7.0 and August 2021 M7.2* events. Note the notable difference at Port-au-Prince. #Haiti#Earthquake
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events. https://earthquake.usgs.gov/earthquakes/eventpage/us6000f53e/executive
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.
Below is my interpretive poster for this earthquake
I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here. https://earthquake.usgs.gov/earthquakes/eventpage/us6000f48v/executive
Below is my interpretive poster for this earthquake
I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
A few days ago, I was passed out on my couch (sleep apnea) and for some reason I awoke and noticed that I had gotten a CSEM notification of a large earthquake offshore of Alaska. Well, after looking into that, I sent my boss, Rick, a text message: “8.2.” https://earthquake.usgs.gov/earthquakes/eventpage/us6000f02w/executive
Rick Wilson runs the tsunami program at the California Geological Survey (CGS) and works with the California Governor’s Office of Emergency Services (Cal OES) to use official forecasts of tsunami size from the National Tsunami Warning Center (NTWC) to alert coastal emergency managers about the level of potential evacuation that they may want to act upon.
More about this process can be found here. Take a look at the CGS Special Report 236 to learn about the Tsunami Playbooks and the “FASTER” approach for tsunami evacuation guidance. Evacuation is something that is done at the local level, so CGS and Cal OES can only provide recommendations.
Needless to say, we were both at the ready to respond. Rick has hourly phone calls with the NTWC and follows up with phone calls and emails to specific interested parties (e.g. the emergency managers). We each went into tsunami response mode. I manage the Tsunami Event Response Team, which may be activated to collect observations of tsunami inundation or ocean currents.
I started looking at tide gage and DART Buoy data to see how large the tsunami was in the epicentral region. The M 8.2 was in the region of the 1938 M 8.2 earthquake which generated a transoceanic tsunami. I also looked into the literature about the 1938 tsunami, to see what size that tsunami was. The 1938 tsunami had a decimeter scale wave height (peak to trough) for gages in Alaska and in California (Johnson and Satake, 1994). Jeff Freymueller et al. (2021) had also recently worked on the 1938 earthquake source area and tsunami modeling as well.
The nearest tide gage for this 2021 event is at Sand Point, but the nearest gage in 1938 was in Unalaska. So, in order to get a modest comparison between 1938 and 2021, I felt a need to wait for the Unalaska data to trickle in. This may give us some idea whether the 1938 tsunami recorded in Crescent City and San Francisco might be a decent analogue. Of course, we need to get the official forecast from the NTWC prior to sending out any information. But, that process can take hours (over 3 hours in this case). So, we need to get our minds wrangled around the possibilities in the absence of more information. Earthquake and Tectonic Background:
The plate boundary in the north Pacific is a convergent (pushing together) plate boundary where the Pacific plate on the south ‘subducts’ northwards beneath the North America plate on the north. The Alaska-Aleutian subduction zone forms a deep sea trench which can be seen in maps of the region. The subduction zone fault dips into the Earth, getting deeper to the north.
Between earthquakes (the interseismic period), the megathrust fault is seismogenically coupled (i.e. ‘locked’) just like velcro has the ability to hold together one’s wallet. The plates are always moving towards each other. Because the fault is locked, the crust surrounding the fault bends elastically to accommodate this convergent motion.
As the crust bends and flexes, it stores energy (i.e. tectonic strain). The part of the fault closest to the seafloor (the southernmost part of this subduction zone fault) gets pulled downwards, while the part of the crust further to the north flexes upwards.
The materials along the earthquake fault have properties that resist motion (like the velcro). But, as the plates converge and increase the amount of energy stored, the forces on the fault may exceed the strength of the fault. At this time, the fault slips, causing an earthquake.
The part of the fault that was being pulled downwards gets pushed upwards during the earthquake (the coseismic period), while the crust that was being flexed upwards between earthquakes thus subsides downwards during the earthquake.
The Alaska-Aleutian subduction zone has a history of subduction zone earthquakes and tsunami, plus there exists a prehistory of earthquakes and tsunami in some parts of this plate boundary. Geologists are often asked to determine the potential hazard of future earthquakes and tsunami and their answers are based on what we know from the past (using both historic and prehistoric data).
The 2021 M 8.2 earthquake happened in the same location as a 1938 M 8.2 earthquake, just to the east of a sequence of earthquakes from last year (22 July and 19 October 2020). Tsunami:
When the earthquake fault slips, and the upper plate deforms, the vertical motion of the plate can elevate (or lower) the overlying ocean water. After the water changes position, it seeks to return to sea-level (an equipotential surface). If elevated, the water drops downwards and then oscillates up and down. This is the process that generates waves that radiate from the area with seafloor deformed by the earthquake.
Things that make a tsunami larger are [generally]:
More vertical land motion (possibly from larger slip on the fault, e.g. from a larger magnitude earthquake)
Deeper water (deeper water = more volume of water moving = more energy to create larger tsunami waves)
So let’s take a look at the things that may have affected the size of the tsunami from this 2021 M 8.2 earthquake.
First of all, based on the earthquake slip models (estimates of how the earthquake slipped, in meters, and how that slip varied along the fault) suggest that a majority of the largest slip happened beneath the continental shelf. The water depth on the shelf is similar to many shelfs worldwide, shallower than about 200 meters. How does this affect the size of the tsunami?
Well, I guess that is the main point, the ground deformation that generated the tsunami was beneath shallow water.
These slip models are based on a variety of data and most of the data are seismic data. Some tsunami are generated by slow slip (not generating seismic waves) on the shallow part of the fault. These are called tsunami earthquakes.
Because tsunami earthquakes may be generated by slip in this way, slip models using seismic data cannot resolve the location of the slip on the fault that created these tsunami. However, the tsunami from this 2021 M 8.2 earthquake were small. Therefore the updip part of the fault probably did not contribute significantly to the tsunamigenic ground deformation.
Below is my interpretive poster for this earthquake
I plot the seismicity from the past 3 months, with diameter representing magnitude (see legend). I also include earthquake epicenters from 1921-2021 with magnitudes M ≥ 7.5.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
I include outlines of the historic subduction zone earthquakes as prepared by Peter Haeussler from the USGS in Anchorage. He appears in the video about the 1964 earthquake below.
Some of the tide gage and DART buoy locations are labeled.
Note how there are still aftershocks from the 2018 M 7.9 earthquake sequence.
I include some inset figures. Some of the same figures are located in different places on the larger scale map below. I present 3 posters, each with slightly different information.
This is the first poster I prepared.
In the upper center is a low-angle oblique view of the plate boundary. Note the oceanic Pacific plate is subducting beneath the continental North America plate. As the plate goes down, the water embedded in the rocks and sediment are released into the overlying mantle wedge. This water causes the mantle to melt, which rises, erupts as lava and forms the volcanic chain we call the Aleutian Islands. I place a green star in the “epicentral” location of the 2021 M 8.2 earthquake.
In the upper left corner is part of a figure from Witter et al. (2019) that shows sections of the megathrust fault relative to how much the fault is thought to be locked. This is called the coupling ratio. For a fault that is fully coupled (or locked), the ratio is 1.0. For a fault that is slipping about 50% and accumulating about 50% of the plate motion rate, the coupling ratio is 0.5. Many subduction zones have low coupling ratios of 0.2-0.6. The region of the fault west of the 1938 and 2021 M 8.2 earthquakes is called the Shumagin Gap, thought to be possibly aseismic (with a coupling ratio closer to 0). But the 2020 sequence of M 7.8 and 7.9 earthquakes filled much of this gap.
In the upper right corner is a plot showing the earthquake shaking intensity using the Modified Mercalli Intensity Scale (MMI). This is a USGS model based on observations of intensity from thousands of earthquakes. Read more about MMI here.
In the center right is a plot showing the aftershocks within a couple hours of the mainshock
In the lower right corner is the initial record of the tsunami at the Sand Point tide gage (see map for gage location).
I labeled the USGS slab 2.0 slab contours (Hayes et al., 2018). These depth contours represent the depth of the megathrust fault at these locations. The M 8.2 hypocentral depth is 32.2 km and the slab2 depth is about 35 km. Nice!
Here is the map with 3 month’s seismicity plotted. There are 3 posters. The first one is something I put together around 2 hours after I awoke on the couch (abt 2am my time). I prepared the 2nd poster an hour later, which includes some information about tsunami prehistory. I prepared the 3rd poster late Sunday evening, about 3 days after the earthquake.
This is the second map I prepared and some figures are the same as in the first poster.
Below the low-angle oblique map is a slip model from the USGS. The color represents the amount of slip on the fault. Note that the maximum slip is close to the epicenter. This is not always the case, as for the 1938 event, it appears that the maximum slip was not where the mainshock epicenter was.
In the upper left corner is a map from Nelson et al. (2015). Those authors studied the prehistoric tsunami records at Chrikof Island, an island about 200 km to the east of the 2021 M 8.2 epicenter. The lower map shows GPS derived plate motion rates.
In the lower right corner is also from Nelson et al. (2015). On this plot, the vertical axis represents time with “today” at the top and over 5000 years ago at the bottom. The horizontal axis is space, west to east from left to right. Each colored symbol represents the time of a prehistoric tsunami. The vertical size of these symbols represents the uncertainty (or “error”) associated with those chronologic data. We can take the number of earthquakes or tsunami over a period of time to estimate how frequently those process happen over time.
To the left is a more updated version of the Sand Point tide gage, showing a wave height (peak to trough) of about 45 cm. We cannot compare this to the 1938 tsunami as there was not a tide gage at Sand Point in 1938
I prepared a 3rd poster, but updated it to this 4th poster.
In the Intensity Data area, I added USGS “Did You Feel It?” data, which come from reports from real people. Learn more about dyfi here. The model data are the colored lines labeled in white and the dyfi data are colored polygons labeled in yellow.
In the aftershocks plot, I added epicenters from the several days after the mainshock. I also added a transparent overlay of the USGS finite fault model (the slip model). Compare the overlap, or non-overlap, of the slip region and the aftershocks. Why do you think that they are not completely overlapping?
In the lower right section are tide gage records from gages in the area included in the poster. I plot the tidal forecast (dark blue), the tide gage observed water surface elevation (medium blue), and the difference between these data (in light blue) which is a record of the tsunami (and other waves, like wind waves). I made a rough approximation estimate of the maximum wave height and labeled this in yellow. The San Point tide gage has a mx wave height of about 0.8 m!
I also plot the data from the DART buoy 46403, which is the closest DART buoy to the mainshock epicenter. The DART buoy network is used to help calibrate tsunami forecast models during tsunami events. These are basically pressure transducers on the seafloor that measure changes in pressure caused by waves and atmospheric processes. The data plotted here are not tsunami data, but seismic wave data. One reason we know that this is not a tsunami is that the waveform initiated about 3 minutes after the earthquake. A tsunami would take longer to get to the buoy.
In the upper left corner is a pair of maps that show USGS earthquake induced ground failure models. The map on the right shows what areas have likelihood of having landslides triggered by the 2021 M 8.2 earthquake. The panel on the right shows the possibility that areas might experience liquefaction induced by the earthquake.
I added aftershocks associated with the 2020 M 7.8/7.5 sequence that filled the Shumagin Gap (green circles) and outlined the aftershock region for both 2020 and 2021 sequences. The 2021 sequence is not yet over. The largest aftershock so far has only been M 6.1. The 1938 M 8.2 event had a M~7 event 5 days after the mainshock. Stay tuned?
Tectonic Overview
Below is an educational video from the USGS that presents material about subduction zones and the 1964 earthquake and tsunami in particular. Youtube Source IRIS mp4 file for downloading.
Credits:
Animation & graphics by Jenda Johnson, geologist
Directed by Robert F. Butler, University of Portland
U.S. Geological Survey consultants: Robert C. Witter, Alaska Science Center Peter J. Haeussler, Alaska Science Center
Narrated by Roger Groom, Mount Tabor Middle School
This is a map from Haeussler et al. (2014). The region in red shows the area that subsided and the area in blue shows the region that uplifted during the earthquake. These regions were originally measured in the field by George Plafker and published in several documents, including this USGS Professional Paper (Plafker, 1969).
Here is a cross section showing the differences of vertical deformation between the coseismic (during the earthquake) and interseismic (between earthquakes).
This figure, from Atwater et al. (2005) shows the earthquake deformation cycle and includes the aspect that the uplift deformation of the seafloor can cause a tsunami.
Here is a figure recently published in the 5th International Conference of IGCP 588 by the Division of Geological and Geophysical Surveys, Dept. of Natural Resources, State of Alaska (State of Alaska, 2015). This is derived from a figure published originally by Plafker (1969). There is a cross section included that shows how the slip was distributed along upper plate faults (e.g. the Patton Bay and Middleton Island faults).
Here is a graphic showing the sediment-stratigraphic evidence of earthquakes in Cascadia, but the analogy works for Alaska also. Atwater et al., 2005. There are 3 panels on the left, showing times of (1) prior to earthquake, (2) several years following the earthquake, and (3) centuries after the earthquake. Before the earthquake, the ground is sufficiently above sea level that trees can grow without fear of being inundated with salt water. During the earthquake, the ground subsides (lowers) so that the area is now inundated during high tides. The salt water kills the trees and other plants. Tidal sediment (like mud) starts to be deposited above the pre-earthquake ground surface. This sediment has organisms within it that reflect the tidal environment. Eventually, the sediment builds up and the crust deforms interseismically until the ground surface is again above sea level. Now plants that can survive in this environment start growing again. There are stumps and tree snags that were rooted in the pre-earthquake soil that can be used to estimate the age of the earthquake using radiocarbon age determinations. The tree snags form “ghost forests.
This is a photo that I took along the Seward HWY 1, that runs east of Anchorage along the Turnagain Arm. I attended the 2014 Seismological Society of America Meeting that was located in Anchorage to commemorate the anniversary of the Good Friday Earthquake. This is a ghost forest of trees that perished as a result of coseismic subsidence during the earthquake. Copyright Jason R. Patton (2014). This region subsided coseismically during the 1964 earthquake. Here are some photos from the paleoseismology field trip. (Please contact me for a higher resolution version of this image: quakejay at gmail.com)
This is another video about the 1964 Good Friday Earthquake and how we learned about what happened.
Here is a map that shows historic earthquake slip regions as pink polygons (Peter Haeussler, USGS). Dr. Haeussler also plotted the magnetic anomalies (grey regions), the arc volcanoes (black diamonds), and the plate motion vectors (mm/yr, NAP vs PP).
Here is the figure from Sykes et al. (1980) that shows the space time relations for historic earthquakes in relation to the map.
Above: Rupture zones of earthquakes of magnitude M > 7.4 from 1925-1971 as delineated by their aftershocks along plate boundary in Aleutians, southern Alaska and offshore British Columbia [after Sykes, 1971]. Contours in fathoms. Various symbols denote individual aftershock sequences as follows: crosses, 1949, 1957 and 1964; squares, 1938, 1958 and 1965; open triangles, 1946; solid triangles, 1948; solid circles, 1929, 1972. Larger symbols denote more precise locations. C = Chirikof Island. Below: Space-time diagram showing lengths of rupture zones, magnitudes [Richter, 1958; Kanamori, 1977 b; Kondorskay and Shebalin, 1977; Kanamori and Abe, 1979; Perez and Jacob, 1980] and locations of mainshocks for known events of M > 7.4 from 1784 to 1980. Dashes denote uncertainties in size of rupture zones. Magnitudes pertain to surface wave scale, M unless otherwise indicated. M is ultra-long period magnitude of Kanamori 1977 b; Mt is tsunami magnitude of Abe[ 1979]. Large shocks 1929 and 1965 that involve normal faulting in trench and were not located along plate interface are omitted. Absence of shocks before 1898 along several portions of plate boundary reflects lack of an historic record of earthquakes for those areas.
Here is a great illustration that shows how forearc sliver faults form due to oblique convergence at a subduction zone (Lange et al., 2008). Strain is partitioned into fault normal faults (the subduction zone) and fault parallel faults (the forearc sliver faults, which are strike-slip). This figure is for southern Chile, but is applicable globally.
Proposed tectonic model for southern Chile. Partitioning of the oblique convergence vector between the Nazca plate and South American plate results in a dextral strike-slip fault zone in the magmatic arc and a northward moving forearc sliver. Modified after Lavenu and Cembrano (1999).
In 2016, there was an earthquake along the Alaska Peninsula, a M 7.1 on 2016.01.24. Here is my earthquake report for this earthquake. Here is a map for the earthquakes of magnitude greater than or equal to M 7.0 between 1900 and today. This is the USGS query that I used to make this map. One may locate the USGS web pages for all the earthquakes on this map by following that link.
The tidal forecasts are shown as a dark blue line.
The actual observed water surface elevation is plotted in medium blue.
By removing (subtracting) the tide forecast from the observed data, we get the signal from wind waves, tsunami, and atmospheric phenomena. This residual is plotted in light blue.
The scale for the tsunami wave height is on the right side of the chart.
Note the all tsunami wave height plots are the same vertical scale, except for Sand Point.
I measured the largest wave heights for each site, displayed in yellow. Alaska
Here are the data from the DART buoy nearest the M 8.2. People often mistake these data for tsunami data, but this is generated by seismic waves.
One way to test one’s hypothesis about whether these buoy data are seismic waves or tsunami waves, one simply need to take a look at the time that the wave begins to be recorded by the DART buoy.
Seismic waves travel through water at about 1.5 kms per second. While tsunami wave velocity (based on the shallow water wave equation) for depths ranging from 200-4000 meters is between ~0.02 to 0.2 kms per second, much slower than seismic waves.
Surface Deformation
Below are surface deformation data generated by the USGS based on their finite fault model. The three panels show surface deformation in the north, east, and vertical directions.
North, East, and Up are positive (blue) while South, West, and Down are negative (red).
Note the upper panel and how the Pacific plate is moving to the north and the North America is moving south. Does this make sense?
The middle panel is interesting too, but skip to the lower panel, vertical. The accretionary prism (forming the continental slope), directly above the aftershocks and mainshock, rises up during the earthquake. The upper North America plate landward of the slip patch subsides. Does this make sense?
Earlier in this report we took a look at the geologic evidence for megathrust subduction zone earthquakes, evidence that records this “coseismic” subsidence.
Shaking Intensity and Potential for Ground Failure
Below are a series of maps that show the shaking intensity and potential for landslides and liquefaction. These are all USGS data products.
There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:
FOS = Resisting Force / Driving Force
When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.
Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.
Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.
Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides.
Johnson and Satake (1994) studied tsunami waveforms from the 10 November 1938 Alaska M 8.2 earthquake. Their analysis was designed to estimate the source for the tsunami. Below are some figures from their paper, with figure captions beneath each figure.
This first plot shows the tsunami records from tide gages. This is the plot I used to consider the potential impact to the coast from the 2021 M 8.2 tsunami.
Digitized marigrams from 1938 Alaskan earthquake recorded in Crescent City, San Diego, and San Francisco. The tidal componenht asn ot beenr emoved.S tartt ime listedf or each record is the time in minutes from the origin time of the earthquaketo the startt ime of the digitizedr ecord.
Here is a map that shows the fault model that they used, as well as the amount of slip that they used for each fault element.
Location of subfaults used in inversion of tsunami waveforms. Graph shows slip distribution in meters.
This is a figure comparing their model results (synthetic = dashed) compared to the tide gage records (solid lines).
Observed and synthetic waveforms from inversion for four subfaults. Start time of each record is different. The arrows indicate the parts of the waveforms used for the inversion.
Freymueller et al. (2021) also studied the 1938 M 8.2 event, seeking to resolve the slip on the fault using tsunami modeling.
Below are figures with their captions in blockquote.
Here are some maps showing 2 of the slip distrubutions that they used for their modeling.
Example slip distributions for two of the slip models, shallow eastern and shallow far eastern. For each model the slip is the product of a function f(x) representing the along-strike variation and g(y) representing the downdip variation, and then scaled to a constant magnitude MW 8.25. The functions f(x) and g(y) are based on relations in Freund and Barnett [1976]. For the central and western models, the rupture area is the same as for the eastern model, but the area of higher slip is shifted to the west. For the mid-depth and deep models, the main area of high slip is shifted downdip.
Here are some maps showing vertical seafloor displacements for some of their tsunami scenarios.
Vertical seafloor displacements caused by representative slip scenarios. On the left side, the slip is concentrated in the east and the deep, mid-depth and shallow slip distribution scenarios are shown. On the right, the Western, Central and Far Eastern slip distribution scenarios are shown assuming the shallow rupture. Displacements are in meters. Red contours show depth to the plate interface from 0 to 80 km with a 10 km increment.
Here are plots that show some results of their modeling. The tide gage data are plotted in black and their simulated waves are plotted with red and blue lines.
Tide gauge data and model predictions for the eastern and far eastern source models.
Here is an animation from one of the Ferymueller et al. (2021) models for the 1938 M 8.2 tsunami.
Nelson et al. (2015) presented their evidence for prehistoric tsunami on Chirikof Island, an island in the forearc in the eastern part of the 1938 earthquake slip patch.
They found evidence for many tsunami over a timespan from before 5000 years ago.
Below are some figures from their paper, with figure captions in blockquote.
This figure shows the tectonic setting and the area of their field study.
A) Location of Chirikof Island within the plate tectonic setting of the Alaska-Aleutian subduction zone. Rupture areas for great twentieth century earthquakes on the megathrust are in pink. (B) Velocity field of the Alaska Peninsula and the eastern Aleutian Islands observed by global positioning system (GPS) (Fournier and Freymueller, 2007). Colors show inferred rupture areas for earthquakes in 1788 (green) and 1938 (orange). Both A and B are modified from Witter et al. (2014). The section of the megathrust between Kodiak Island and the Shumagin Islands has been referred to as the Semidi segment (e.g., Shennan et al., 2014b). (C) Physiography of Chirikof Island (Google Earth image, 2012) showing the location of our study area at Southwest Anchorage, a prominent moraine, a fault scarp (facing southeast) that probably records the 1880 earthquake, the New Ranch valley reconnaissance core site, and UNAVCO GPS station AC13 (http:// pbo .unavco .org /station /overview /AC13). In the eighteenth and nineteenth centuries, Chirikof Island was known to native Alutiiq and Russians as Ukamuk Island.
Here is a plot that shows the timing for the prehistoric tsunami inferred by these authors. The vertical axis is the time scale, with “today” at the top. Each colored pattern represents the age range for a tsunami deposit.
These data are plotted left to right, west to east, so we can compare tsunami records at different locations along the margin. These comparisons are important so that we can test different hypotheses about how subduciton faults may slip over time. In the 2021 case, the slip area was close to the 1938 earthquake. But, did has this always occured here?
Age probability distributions for probable (red) and possible (orange) tsunami deposits at Southwest Anchorage (labels as in Fig. 11) compared with age distributions for possible tsunami deposits at Sitkinak Island (Briggs et al., 2014a) and with age estimates for great earthquakes and tsunamis on Kodiak Island (from studies referenced on this figure;
Fig. 1). Dotted horizontal lines show our correlation of evidence for some younger earthquakes and tsunamis. Times of great earthquakes inferred from episodes of village abandonment determined from archaeological stratigraphy in the eastern Alaska-Aleutian megathrust region are also shown (Hutchinson and Crowell, 2007).
— Jason "Jay" R. Patton (@patton_cascadia) July 29, 2021
From @BNONews BREAKING: Tsunami sirens sound in Kodiak, Alaska after a major magnitude 8.2 earthquake struck off the coast; risk being evaluated for the Pacific pic.twitter.com/amxpLGX70s
Here's the #USGS MT for the recent M 8.2 on Fig. 1 of Freymueller et al. 2021 (https://t.co/FN8owbDqEY). Orange outline is aftershocks of the 1938 M 8.2. Red lines are 1 m contours of 1938 slip models. Grey is slip deficit inferred from geodesy. Obvious similarities 1938 -> 2021! pic.twitter.com/DIUh4YVhXc
UPDATE: The timing and form of this signal looks like it is the DART response to the seismic waves directly from the earthquake, NOT to a tsunami wave. pic.twitter.com/bxeF5TPqjv
Small tsunami waves continue arriving at Sand Point & other coastal areas of Alaska. Tomorrow these waves will create swirly currents in boat harbors up & down the west coast, so tie up your boats real good. pic.twitter.com/nofvKqJoU5
And since I have a drone workshop to attend tomorrow, I will bow out now and get some sleep. My initial guess that today's event may have been similar to the 1938 M8.2 earthquake still looks like it has some merit. Follow https://t.co/A1MNRg1WKF for updates on tsunami warnings. pic.twitter.com/g4qME2w0SI
@NOAA Tsunami Warning System has issued a tsunami watch for the West Coast. The warning for Hawaii has been cancelled, because the waves are focused east of Hawaii and the event isn't that large. @NWS_NTWCpic.twitter.com/h2KRBmOKNL
Records of tsunami deposits show significant tsunamis in 1788, 1880 and 1938 (https://t.co/NsFfTuqigs), indicating recurrence intervals of large earthquakes in the Semidi segment every 58-92 years. We are now 83 years since 1938, so that seems roughly consistent. pic.twitter.com/CGIM40Fv0g
— Jason "Jay" R. Patton (@patton_cascadia) July 29, 2021
Preliminary finite fault for this morning's M8.2 earthquake is available. Rupture primarily to the NE of the hypocenter, away from the Shumagin Gap.https://t.co/dVkYuR2kPCpic.twitter.com/idGBqxRhbX
Tonight's M8.2 event occurred close to the rupture area of the 2020 M7.8 earthquake and was the largest U.S. earthquake in 50 years. We'll continue to update as this sequence unfolds, but here is a short piece on our website with what we know so far. https://t.co/PzHaaQ8Zblpic.twitter.com/vcM8fq9IV7
— Alaska Earthquake Center (@AKearthquake) July 29, 2021
Some Perryville M 8.2 thoughts: One of the arresting things about Chirikof coastal geology is that the island is clearly sinking like a stone today, evident in geodesy and coastal geology. Figure from Nelson et al. 2015 https://t.co/vGKDp0WYuN
The two closest sites to the M8.2 Alaska earthquake today show some decent surface wave signals. There are several other closer sites that should give us better insight. @UNAVCOpic.twitter.com/lN22i7arEP
8.2 Earthquake is the largest in Alaska since 1965. I was sitting in the upper wheelhouse of my 125' steel schooner ALEUTIAN EXPRESS at Chignik Harbor and the whole boat bounced and vibrated for about a minute. 14' range of gradual Tsunami one foot every 4 minutes both directions pic.twitter.com/IlYox48ejg
Clear NE propagation from the M8.2 in Alaska, but look at 102 sec- action to E way updip by the trench. Early aftershock or where rupture finally expired? It's small amplitude, but coherent and seen by 4 very different arrays. I await better analyses.https://t.co/6hFfZ64Elwpic.twitter.com/hacucSnOil
Good morning all! The tsunami waves are still bouncing around the Aleutian Islands in Alaska (max height measured was ~2 feet). The tsunami turned out not to be very big & all @NWS_NTWC alerts for the US west coast are CANCELLED. 🚨NO alerts for CA, OR, WA. #earthquakepic.twitter.com/uEHSdzzvv9
Good morning PNW- ICYMI, last night there was a M8.2 earthquake off the Alaska Peninsula. Here, you can see waves from it (bottom) compared to a nearby Alaskan M6.8 (top, similar to our 2001 Nisqually M6.8) at station LEBA near the SW Washington coast. pic.twitter.com/GCCAjbYpII
Here I show the cross-section through the Alaska seismicity with projected mechanisms. The largest two events are yesterday's M8.2 quake and last year's M7.8, both subduction interface events. For reference, a cartoon of the shallow subduction zone from https://t.co/Gces1m71C8pic.twitter.com/bgchTpPT5n
You may not have felt it, but a groundwater well in Washington County, Maryland did! An 8.2 magnitude earthquake rocked southern Alaska overnight and the water level in our well sloshed almost a foot. https://t.co/kafsMsaaph. For more real-time well data https://t.co/w56ACDNk4hpic.twitter.com/87fVSz0MLz
15 second sample rate data for AB13 is now available for the M8.2 Alaska earthquake, we see a pretty appreciable SE offset with 10 cm of subsidence. The event started NE of AC12 and ruptured to the NE, so this site is in the middle of it all. @UNAVCOpic.twitter.com/8nTTZstIDX
The largest earthquake to hit the U.S. in the last few decades took place in Alaska yesterday. The Mw 8.2 quake broke the Aleutian megathrust in the Shumagin seismic gap. The rupture did not propagate to the trench, causing only a minor tsunami. Figure by @QQtecGeodesypic.twitter.com/02ylFet8P6
14+ hours after the #alaska earthquake and there is still a tsunami bouncing around at the closest tide gauge (small tsunami) pic.twitter.com/3xtjS4hhge
— Dr Janine Krippner (@janinekrippner) July 29, 2021
As of 12 hours following the M8.2 we've located ~140 aftershocks. The locations and magnitudes are subject to change upon further review, but look to be occurring to the east of 2020 sequence. The map here shows 2020 in gray and the recent aftershocks in red. pic.twitter.com/hQ93k7HVUZ
— Alaska Earthquake Center (@AKearthquake) July 29, 2021
Last night's magnitude 8.2 earthquake serves as a powerful reminder of the restlessness of our planet's surface—and it presents an exciting opportunity to peer deeper at our planet’s inner workings.
(1/3) The "Lame Monster": Today's largest US earthquake in >50 years did not make a large tsunami. Why? B/c most of elastic energy was released deeper in the Earth.
— California Department of Conservation (@CalConservation) July 30, 2021
GPS receivers can be used as seismometers. In blue are the 5 Hz velocities recovered on Kodiak Island with the variometric approach for the M8.2 earthquake yesterday. In red, the collocated accelerometer, S19K, downsampled to 5 Hz. pic.twitter.com/jxJR1v7Fj1
A notable characteristic of the M8.2 Alaska earthquake is that it was relatively deep and doesn’t appear to have ruptured the shallow plate boundary. Could overpressured sediments on the shallow plate boundary inhibited shallow slip? Check out this seismic image updip of event. pic.twitter.com/HRQEPxrAZk
While the M8.2 was the largest earthquake in the U.S. in 50 years, Alaska has experienced some significantly sized events during that time. The plot here shows the largest Alaska earthquake magnitude each year since 1964. Since 2000, we're experienced at least a M6.4 annually. pic.twitter.com/Iq9pFanPqi
— Alaska Earthquake Center (@AKearthquake) July 30, 2021
Whopper M8.2 earthquake in Alaska moved GPS stations, revealing the broad pattern and extent of deformation. Stations near Denali NP ~900 km moved a few mm… See https://t.co/4zpOW4m1pJ for more info and data. pic.twitter.com/j1GSIazVfJ
— Jason "Jay" R. Patton (@patton_cascadia) August 2, 2021
#Sentinel1 co-seismic interferograms (ascending track) over western Alaska, show ground deformation towards the southern coast, above the main M8.2 #earthquake fault rupture. Aftershock epicenters (yellow) from USGS. pic.twitter.com/RtavuJZGSZ
The M 8.2 Chignik earthquake that occurred off the Alaskan Peninsula on July 28 was the largest US earthquake in 50 years. This 2013 simulation from the same region shows how a hypothetical M 9.1 (almost 30x stronger!) earthquake can create a far-reaching tsunami. @USGS_Quakespic.twitter.com/tLtxWxoal7
— USGS Coastal Change (@USGSCoastChange) July 30, 2021
— Jason "Jay" R. Patton (@patton_cascadia) August 4, 2021
References:
Basic & General References
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
Specific References
Abe, K., 1972. Lithospheric Normal Faulting Beneath the Aleutian Trench in Phys. Earth Planet. Interiors, v. 5, p. 1990-198.
Benz, H.M., Tarr, A.C., Hayes, G.P., Villaseñor, Antonio, Hayes, G.P., Furlong, K.P., Dart, R.L., and Rhea, Susan, 2011. Seismicity of the Earth 1900–2010 Aleutian arc and vicinity: U.S. Geological Survey Open-File Report 2010–1083-B, scale 1:5,000,000.
Brown, J.R., Prejan, S.G., Beroza, G.C., Gomberg, J.S., and Hauessler,m P.J., 2013. Deep low-frequency earthquakes in tectonic tremor along the Alaska-Aleutian subduction zone in JGR Solid Earth, v. 118, p. 1079-1090, doi:10.1029/2012JB009459
Freymueller, J.T., Suleimani, E.N., and Nocolski, D.J., 2021. Constraints on the Slip Distribution of the 1938 MW 8.3 Alaska Peninsula Earthquake From Tsunami Modeling in GRL, v. 48, no. 9, https://doi.org/10.1029/2021GL092812
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Gaedicke, C., Baranov, B., Seliverstov, N., Alexeiev, D., Tsukanov, N., Freitag, R., 2000. Structure of an active arc-continent collision area: the Aleutian-Kamchatka junction. Tectonophysics 325, 63–85
Johnson, J.M. and Satake, K., 1994. Rupture extent of the 1938 Alaskan earthquake as inferred from tsunami waveforms in GRL, v. 21, no. 8, p. 733-736
Koulakov, I.Y., Dobretsov, N.L., Bushenkova, N.A., and Yakovlev, A.V., 2011. Slab shape in subduction zones beneath the Kurile–Kamchatka and Aleutian arcs based on regional tomography results in Russian Geology and Geophysics, v. 52, p. 650-667.
Konstantnovskaia, 2001. Arc-continent collision and subduction reversal in the Cenozoic evolution of the Northwest Pacific: an example from Kamchatka (NE Russia) in Tectonophysics, v. 333, p. 75-94.
Lay, T., Ye, L., Bai, Y., Cheung, K. F., Kanamori, H., Freymueller, J., … Kogan, M. G. (2017). Rupture along 400 km of the Bering fracture zone in the Komandorsky Islands earthquake (MW 7.8) of 17 July 2017. Geophysical Research Letters, 44, 12,161–12,169. https://doi.org/10.1002/2017GL076148
Nelson, A.R., Briggs, R.W., Dura, T., Engelhart, S.E., Gelfenbaum, G., Bradley, L.-A., Forman, S.L., Vane, C.H., and Kelley, K.A., 2015, Tsunami recurrence in the eastern Alaska-Aleutian arc: A Holocene stratigraphic record from Chirikof Island, Alaska: Geosphere, v. 11, no. 4, p. 1172–1203, doi:10.1130/GES01108.1.
Plafker, G., 1969. Tectonics of the March 27, 1964 Alaska earthquake: U.S. Geological Survey Professional Paper 543–I, 74 p., 2 sheets, scales 1:2,000,000 and 1:500,000, http://pubs.usgs.gov/pp/0543i/.
Shevchenko, V.I., Lukk, A.A., and Prilepin, M.T., 2006. The Sumatra Earthquake of December 26, 2004, as an Event Unrelated to the Plate-Tectonic Process in the Lithosphere in Physics of the Solid Earth, v. 42, no. 12, p. 1018–1037.
Stauder, W., 1968. Mechanism of the Rat Island Earthquake Sequence of February 4, 1965, with Relation to Island Arcs and Sea-Floor Spreading in JGR, v. 73, no. 12, p. 3847-3858
Sykes, L.R., Kissinger, J.B>, House, L., Davies, J.N>, and Jacob, K.H., 1980. Rupture Zones and Repeat Times of Great Earthquakes Along the Alaska-Aleutian Arc, 1784-1980, in Maurice Ewing Series, Earthquake Prediction, An International Review, AGU
Torsvik, T. H. et al., 2017. Pacific plate motion change caused the Hawaiian-Emperor Bend in Nat. Commun., v. 8, doi: 10.1038/ncomms15660
Wilson, J. Tuzo, 1963. “A possible origin of the Hawaiian Islands” in Canadian Journal of Physics. v. 41, p. 863–870 doi:10.1139/p63-094.
I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here. https://earthquake.usgs.gov/earthquakes/eventpage/us6000exs5/executive
Below is my interpretive poster for this earthquake
I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
Central China here, part of the Tibetan Plateau, is dominated by the plate tectonics and climate.
To the south is the India plate that is moving northwards, pummeling into the southern part of the Eurasia plate, at a rate of about 25 to 50 mm per year, depending upon the reference frame (Pusok and Stegman, 2020). The India plate began moving northward from Antarctica since before 80 million years ago.
The details of the story has changed as more geological information is interpreted. But the general story is that the India plate moved away from Antarctica as an oceanic spreading center formed between these plates. The India plate moved towards Asia.
Prior to about 45-50 Ma, there was oceanic crust between India and Eurasia. But at this time, the continental crust of the India and Eurasia plates collided. This collision would eventually cause uplift of the Himalaya mountain range and the Tibetan Plateau to the north of the Himalaya.
There are marine fossils on the top of the mountains in the Himalaya! (this is how we know there was ocean between these plates in the past)
Here is a time series showing the convergence of these two plates modified from Pusok & Stegman (2020) and the USGS.
A part of the tectonic story is told by one of the rock stars of plate tectonics, Dr. Paul Tapponier. Tapponier conducted experiments that showed how north-south convergence, like that of India and Asia, coupled with a backstop (something that is more difficult to move) from the west, would lead to some crust to squish out to the east.
This is called extrusion tectonics as the crust of eastern Asia is being extruded to the east, like a watermelon seed is extruded from between one’s fingers when they squeeze on the wet seed.
Below is a color version of the results from Tapponier’s experiment. Compare this with maps showing the GPS motion of the crust in this region.
Note that the plastic has numerous faults develop as part of the extrusion. We can see how the blue and yellow lines show lateral offset along these faults.
Many of these faults are left-lateral strike-slip faults. Strike-slip means that the curst moves side-by-side when looking down on the crust from outer space (or an airplane, or Google Earth).
Left-lateral means that, when standing on one side of the fault, looking across the fault at something, that thin one is looking at is moving to the left during an earthquake. More about tectonic fundamentals here.
In the map on the right ^^^ we can see that these left-lateral strike-slip faults that are mapped in the region are just like the faults in the blue-yellow plastic.
One of the major left-lateral strike-slip faults along the Tibetan Plateau is the Kunlun fault, which has been well studied in places. A tectonic history of the region, and how the Kunlun fault fits into this history, is presented by Staisch et al., 2020.
Earthquake!
The earthquake yesterday was a magnitude M 7.3 earthquake that happened along the Kunlun fault system. The earthquake mechanism is that of an oblique strike-slip fault. Earthquake mechanisms (focal mechanisms or moment tensors) are derived from seismic wave data surrounding the earthquake.
These mechanisms are visual depictions of the orientation of the fault that slipped and there are always two possible fault planes shown on these mechanisms (also called beach balls).
Without additional information, we don’t know which potential fault plane is the correct one. Because the mapped faults in the region are left-lateral strike-slip faults, and the mechanism is oriented (relative to north) to these fault lines, one can use these faults as a basis to interpret which fault is the correct fault on the mechanism and the type of earthquake that happened. Which type of earthquake do you think slipped yesterday?
I interpret (as have many others) that this M 7.3 was a left-lateral strike-slip earthquake. There is a splay fault off the Kunlun fault system called the Maduo Grande fault (e.g. Chen Xia et al., 2011).
The aftershocks align with this mapped fault, supporting my interpretation. The fault is mapped for about 150 km and the aftershocks span about 160 km. Using Wells and Coppersmith empirical fault relations (i.e. using the relations between earthquake fault surface rupture length and earthquake magnitude), we can take the length of the fault that has cut through to the surface to estimate the earthquake magnitude.
Using these relations, taking 150-160 km would produce a M 7.5-7.6 earthquake. This is much larger than an M 7.3. However, the actual surface rupture length is probably less than the span of the aftershocks. regardless, this magnitude and mapped fault are compatible with each other (the length and the magnitude of the earthquake are a moderate match).
There are probably secondary effects from this M 7.3, like maybe landslides or liquefaction. The epicenter is in a part of the river valley where the river has a braided form. Basically, there is more sediment that the water available to move it. In places like this, there are lots of wetlands and swamps. The waterlogged subsurface can tend to promote liquefaction, so there is a high possibility that evidence of liquefaction will be found. Because this is an arid region, InSAR satellite imagery analysis may reveal evidence of liquefaction (in addition to revealing tectonic deformation).
There may have even been a seiche or tsunami in the nearby lake.
Of note is a M 6.1 earthquake that occurred hundreds of kilometers to the south before the M 7.3. While these two earthquakes are unlikely to be related (because of the distance), it is interesting because they are both evidence supporting the extrusion tectonics interpretation.
Below is my interpretive poster for this earthquake
I plot the seismicity from the past 3 months, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 7.0.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
I include some inset figures. Some of the same figures are located in different places on the larger scale map below.
In the upper left corner is a plate tectonic map that shows the major faults, plates, and plate boundaries. The Gray arrows show the ways the plates are moving (compare with the Extrusion Tectonics figure).
Below that tectonic map is the photo of the Tapponier et al. (1982) experiment on the left and a map showing the faults that really exist that match faults from the model.
On the right is a diagram that shows the evolution of faulting in the region that are the result of changes in driving forces (e.g. broad scale tectonic plate motions and how their relative motions cause different types of faults at different times; also that these older faults remain and may be reactivated as different types of faults later (Staisch et al., 2020).
In the lower left center is a map that has red and blue arrows (Tayler and Lin, 2009). These arrows represent places that are moving as measured by GPS observations. The orientation of the arrow shows the direction the plate is moving relative to “stable” Eurasia (data are from Zhang et al., 2004). I place a yellow star in the location of the M 7.5 epicenter.
In the upper left-center there is a map that shows some of the mapping of faults in the area (Chen Xia et al., 2011). The Kunlun fault is the main fault line with a thicker width. I place a yellow star in the location of the M 7.5 epicenter.
In the upper right corner are maps that display data from the USGS event page. On the left is the modeled landslide probability for this earthquake. On the right is a map that shows the susceptibility (chance of) liquefaction from this M 7.3 earthquake. The headwaters of the Yellow River feed the lakes to the west of the earthquake and this river valley is where the Maduo Grande fault is mapped.
In the lower right corner is a map that shows the earthquake shaking intensity modeled for this earthquake. The colors use the Modified Mercalli Intensity scale (MMI), which is listed in the upper right margin (compare the numbers on the map with the data in the table). I also plot the aftershocks here to show how they relate to the mapped fault.
Here is the map with 3 month’s seismicity plotted.
Other Report Pages
Some Relevant Discussion and Figures
Global and Regional Tectonic Faults
Here is an excellent overview of the faults in the region from Taylor and Yin (2009).
I especially like this figure as it helps us understand what the fault patterns mean by looking at the fault types. For example, I thought that the faults in the northwest Tarim Bain would have been strike-slip, but they appear to be predominantly thrust or reverse faults (e.g. the Southern Tian Shan thrust).
A color-shaded relief map with active to recently active faults related to the Indo-Asian collision zone and surrounding regions. (The paper lists the sources of their fault data, but are “augmented by our own kinematic interpretations.”)
Thrust faults have barbs on the upper plate, normal faults have bar and ball on the hanging wall, arrows indicate direction of horizontal motion for strike-slip faults. Dashed white lines are Mesozoic suture zones: IYS—Indus Yalu suture zone; BNS—Bangong Nujiang suture zone; JS—Jinsha suture zone; SSZ—Shyok suture zone; TS—Tanymas suture zone; AMS—Anyimaqen-Kunlun-Muztagh suture zone.
Here is the map from ChenXia et al, (2011) that shows the faults in teh area of the M 7.3 earthquake (located near the MD-GD F fault label. See the interpretive poster.
Geometric distribution map of the mid-eastern part of the Kunlun Fault [5, 9, 12], Diebu-Wudu Fault and Awancang Fault [12]. Other faults are drawn according to ref. [21], the areal geological map with scale 1:200000 [20] and interpretation of ETM satellite imagery. Topography is generated from Shuttle Radar Topography Mission (SRTM) data (http://srtm.csi.cgiar.org/SELECTION/ listImages.asp). HDSW, Huaideshuiwai; ELS F, Elashan Fault; RYS F, Riyueshan Fault; LT-TC F, Lintan-Tanchang Fault; DB-WD F, Diebu-Wudu Fault; TZ F, Tazang Fault; MD-GDF, Maduo-Gande Fault; MDSF, the southern Maduo Fault; Dari F, Dari Fault; LRB F, Longriba Fault; MJ F, Minjiang Fault; XGZ, Xigongzhou; TZ, Tazang.
Here is a map that shows some earthquake mechanisms (centroid moment tensors) for earthquakes from 1977 to 2009 (Taylor and Yin, 2009).
Compare this map with the one above and the mechanisms match pretty well with the types of faults mapped in the above map.
Color-shaded relief map overlain with Harvard centroid moment tensor (CMT) earthquake focal mechanisms from 1 January 1977 to 1 January 2009 and background seismicity from Engdahl and Villasenor (2002) with events >M5.5 for both data sets. Green, purple, and light-blue earthquake focal mechanisms are locations of 2008 western Kunlun, Nima, and Wenchuan events, respectively.
Here is a map from Zhu et al. (2019) that shows the slip rates for some of the faults. Paleoseismologists (geologists who study the prehistoric record of earthquakes) dig holes into the Earth to excavate and expose earthquake faults. They are in search of evidence of past earthquakes, generally in the form of discrete offsets of geologic materials (like rocks, soils, and other stuff).
The timing of earthquakes can be interpreted and numerical ages (like radiocarbon ages) can be used as a basis for this timing. If there are ways to measure how much the fault is offset during these earthquakes, the distance and time information can be used to calculate a “slip rate” for the fault
The slip rate tells us how fast the fault slips over time. The map below shows the slip rates calculated from these studies. Also, geologists can take older geologic units, measure the offset of these units, and calculate a long term geologic slip rate. The map below presents long term geologic slip rates.
For example, The Kunlun fault (the thick red line) has a slip rate of 10+-1.05 mm per year near the 2001 M 8.1 earthquake (labeled M 7.8 on the interpretive poster).
Structural background of the East Kunlun Fault, seismic activity, and long-term geologic slip rate.
Dr. Lydia Staisch and her colleagues (Staisch et al., 2020) use geologic and fault mapping to interpret the tectonic geologic history of this region of central Asia. The next few figures walk us through this story, but read their paper for more information. I will try to update these figures with better quality versions when i get a good copy of the paper (this one is from sci-hub).
This first map shows their interpretation of when each fault system had initiated through time.
Map of major active strike-slip and normal faults in the Tibetan Plateau, adapted from Styron et al. (2010) and Fu et al. (2011). The estimated initiation age of each fault system is denoted, along with abbreviated names. Abbreviations for strike-slip faults are as follows; cKF: central Kunlun Fault, eKF: east Kunlun Fault , wKF: west Kunlun Fault, wHF: west Haiyuan Fault, EF: Elashan Fault, RF: Riyueshan Fault, ATF: Altyn Tagh Fault, JF: Jiali Fault, RRF: Red River Fault; KF: Karakoram Fault, XF: Xianshuihe Fault, RCF: Riganpei Co Fault, GCF: Gyaring Co Fault, BCFS: Bue Co Fault System. Abbreviations for normal faults are as follows; TG: Thakkola Graben, ADR: Ama Drime Rift, KCR: Kung Co Rift, GMR: Gurla Mandhata Rift, RG: Ringbung Graben, YR: Yadong Rift, LR: Lunggar Rift, NR: Nyainqentanghlah Rift, LKR: Lopukangri Rift, TYCR: Tangra Yum Co Rift, PXR: Pumqu-Xainza Rift, GR: Gulu Rift, SHG: Shuang Hu Graben.
Here is their low-angle oblique view of the cut-away view (cross-section) of the fault bound tectonic blocks in the region of the Dongdatan Valley (a valley formed by the Kunlun fault system).
In the Late Triassic, the dominant tectonic forcing was North-South shortening.
Over time, the tectonic forcing shifted to include lateral motion (i.e. strike-slip faulting) in the Miocene.
Schematic block diagrams showing the evolution of deformation within the Dongdatan Valley. (A) Field observations of Permian carbonates thrust over Triassic metapelites and regional dating of plutons and metamorphic cooling episodes (Mock et al., 1999; Liu et al., 2005; Wu et al., 2019) suggest that the East Kunlun Shan experienced late Triassic north-south oriented compression from the northward accretion of the Qiangtang block. (B) North-south compression may have been reactivated during Jurassic – Cretaceous accretion of the Lhasa block based on the timing of Yangqu Group deposition and a separate regional cooling event documented in the West and East Kunlun Shan (Liu et al., 2005; Li et al., 2019). Permian marbles and Triassic metapelites were exposed at the time of Yangqu deposition. (C) Late Cretaceous to Eocene shortening from collision between India and Eurasia resulted in thrust faulting along the Wenquan Hu thrust fault and burial of terrestrial strata in the footwall. Shortening and exhumation may have continued elsewhere in the East Kunlun Shan into late Eocene time, but ceased along the Wenquan Hu thrust fault by 43 Ma. (D) East-west oriented strike-slip faulting locally causes exhumation and erosion. Thermochronologic modeling suggest that strike-slip faulting initiated by ~20 Ma. Miocene – present strike-slip faulting results in basin subsidence in the East Wenquan basin, deposition of terrestrial strata, and juxtaposition of Jurassic – Cretaceous and Cenozoic strata.
Here they present a more detailed cross section of the region through time.
Schematic orogen-scale cross section of the geodynamic evolution of the Tibetan Plateau. (a) The onset of crustal shortening and thickening in northern Tibet soon after the Indo-Asian collision continued into late Oligocene time with moderate elevation gain. (b) Shortening in northern Tibet continues as a southward sweep of magmatism suggests the onset of slab rollback in southern Tibet. (c) Shortening ceases within the Hoh Xil Basin by 27 Ma and likely by 24 Ma in the East Kunlun Shan. Surface uplift may have continued due to crustal thickening via lower crustal flow in northern Tibet. Continued slab rollback may have destabilized the northern Tibetan mantle root by removing is southern buttress. (d) The onset of strike-slip faulting in the East Kunlun Shan between 23 and 20 Ma is coincident with proposed slab breakoff and elevation gain in southern Tibet, and proposed mantle root loss, surface uplift and increased magmatism within the northern Tibet. (e) After 20 Ma, strike-slip and normal faulting expanded throughout the Tibetan Plateau, coincident with the proposed onset of eastward-directed lower crustal flow.
This is an excellent figure that shows how Haibing et al. (2005) used fluvial geomorphologic features (shapes in the landscape formed by fluvial or other physical processes) to derive slip rates.
Basically, there are features in the landscape that are formed over time and their formation is abandoned at some time. Also, these “geomorphic” features (like river banks) can be offset by the earthquake fault during earthquakes.
If we know something about the age of a feature (when it was abandoned or when it was formed) and we know the amount that feature was offset during an earthquake, we can calculate the slip rate of that fault. The more slip rate calculations that can be made provides us with knowledge about how the slip rate can vary through time.
This figure shows a Satellite image of an area with the Kunlun fault curring through river banks. Can you see the fault? Hint, the river banks are on the left and are generally running north-south. The fault cuts across these river banks (and other features)
The lower part of the figure is their interpretation of these features and the ages of them. Harkins et al. (2010) has an excellent compilation of similar features for the eastern Kunlun fault.
(A) Enlarged Ikonos image of eastern part of Hongshui Gou confluence with (B) corresponding geomorphic field map, showing different terrace levels offset by the Kunlun fault. Black dots indicate locations of sampling pits for thermo-luminescence dating on T2 and T3, with corresponding ages
Here Haibing et al. (2005) present the results from their analysis (offsets at different times) and these data all support an interpretation of 10 mm per year of slip on the western Kunlun fault.
(A) Enlarged Ikonos image of eastern part of Hongshui Gou confluence with (B) corresponding geomorphic field map, showing different terrace levels offset by the Kunlun fault. Black dots indicate locations of sampling pits for thermo-luminescence dating on T2 and T3, with corresponding ages
Geodetic Analyses
Geodesy is the study of the motion of the Earth. The data used to measure how the Earth moves can be from GPS data, tide gage data, benchmark surveys, and satellite remote sensing data (e.g. InSAR, LiDAR, etc.). The motion can be partitioned into different directions (e.g. horizontal, vertical, and rotational).
The maps below are from Taylor and Yin (2009, lower map) and Zhang et al. (2004, upper map) shows the velocity (speed) of locations where GPS locations (positions) have been collected over a period of time (probably decades). The arrows called vectors represent the direction of motion and the rate of motion (speed or velocity). The GPS sites are where the dots are and the uncertainty (sometimes called error) of the velocity calculation is represented by the ellipse at the tip of the arrow.
These arrows represent motion relative to stable Eurasia. So, arrows that are pointing to the north tell us that that GPS site is moving north relative to Eurasia. Unfortunately there is no scale, but based on the Zhang et al. (2004) paper (also shown below), the most southwest GPS site (in India) has a velocity of about 25 mm per year (mm/yr)
We can make some simple observations and interpretations from these data. Look at the lower map with the red and blue arrows (vectors).
GPS sites in northern India (in the lower left (southwest) part of the map) show that this region is moving north-northeast relative to Eurasia. This matches the long term motion of the India plate we discussed in the introduction to this report above.
GPS sites in the Tian Basin (the low, green colored area in the central upper left (northwest) part of the map) are also moving north relative to Asia. However, they are moving to the north more slowly than the sites in India
Because the GPS sites (and the crust in that location) are moving slower in the north, north of the Himalaya and Tibetan Plateau. This tells us that the crust is slowing down between India and the Tarim Basin. Why is this?
The crust is slowing down because the crust is deforming, either elastically where the deformation of the crust buldges up or flexes sideways, or anelastically where the deformation is accommodated by fault slip on tectonic compressional faults (e.g. reverse or thrust faults).
Note how these GPS plate motion vectors (the red and blue arrows) change whether they are slightly to the east or slightly to the west of North. In northeast India, the motion is slightly to the northeast and in the Tarim Basin some of them are moving slightly to the west. If this difference is larger than the error ellipses, it would tell us that the crust may also be experiencing changes in lateral motion through this region. This type of lateral motion may be accommodated by strike-slip faults (see the arrow shaped figure from Taylor and Yin (2009) below).
Caption from Zhang et al., 2004)
Global positioning system (GPS) velocities (mm/yr) in and around Tibetan Plateau with respect to stable Eurasia, plotted on shaded relief map using oblique Mercator projection. Ellipses denote 1s errors. Blue polygons show locations of GPS velocity profiles in Figures 3 and DR1 (see footnote 1). Dashed yellow polygons show regions that we used to calculate dilatational strain rates. Yellow numbers 1–7 represent regions of Himalaya, Altyn Tagh, Qilian Shan, Qaidam Basin, Longmen Shan, Tibet, and Sichuan and Yunnan, respectively.
Caption from Taylor and Yin (2009), Figure 4 is the arrow shaped figure below.
Color-shaded relief map of the Indo-Asian collision zone with global positioning system (GPS) velocities (arrows) from the Zhang et al. (2004) compilation. Blue arrows indicate data used in Figure 4A and line indicates data used in Figure 4B.
Here is the cross section of the India-Eurasia plate convergence through time from Pusok and Stegman (2020) shown above. They use numerical modeling of the mantle convection to try to interpret the data presented in this cross section.
Proposed convergence history between India and Eurasia. (A) Schematic diagram for the evolution of Neo-Tethys Ocean as a sequence of stages with approximate times for important events. (B) Hypothesized convergence rate versus age, with numbered segments corresponding to stages in (A), superimposed on the relative plate motions between India-Africa and India-Antarctica (dotted blue and green lines, respectively) (4).
Here Pusock and Stegman (2020) present the long term geodetic data (geologic rates) and compare the results of their modeling with these source data.
India-Eurasia convergence data and results of numerical models. (A) Observed spreading rates versus age for the CIR (black lines) and SEIR (gray lines). Solid lines use the geomagnetic polarity time scale GTS04 and a three-plate algorithm to calculate Euler rotations (2, 4), while dotted lines use the GTS12 time scale and a two-plate rotation algorithm (17). Blue line shows computed convergence rate between the left plate and the overriding plate in the reference model (ConvIndia35). Colored markers correspond to model snapshots in Fig. 3. (B) Speedup over time: Observational data are scaled to a velocity of 70 mm/year, representative of the present-day global subduction. Numerical data (blue line) are scaled to velocities in single subduction experiments (ConvIndia31). (C) Convergence velocity versus the DSF (see Materials and Methods and the Supplementary Materials) for all simulations. Time-averaged values given for single subduction (black), periods of plume push (white), and periods of free double subduction (blue) including reference model (red). Maximum/minimum values during model evolution shown with gray bars. (D) Speedup versus the DSF for all numerical simulations [colors same as (C)]. Least-squares fit (dashed lines) for periods of forced convergence (white markers) and periods of free double subduction (blue markers).
Here Pusok and Stegman (2020) present the details from their numerical modeling. What do you think about this?
Time evolution of the reference model ConvIndia35. The system is composed of three plates representing India (left plate), an intraoceanic plate (middle plate), and Eurasia (overriding plate). This model was performed with plume-like influx boundary conditions imposed for 9 Ma (Influx BC06) and a mantle density and viscosity profile (25, RLLB15-Mean_norm; see details in Materials and Methods and the Supplementary Materials). Frames correspond to colored markers in Fig. 2A, and full-time evolution is available in movie S1. The dynamic pressure (background field) indicates when each regime dominates. (A) During the plume push regime, the entire system is dominated by the influx velocity, seen by positive dynamic pressures throughout the domain. (B) During the double subduction regime, the mantle flow and plate convergence are dominated by interaction between the two slabs, as indicated by pressure increasing within the region confined between them. Changes in driving forces are also seen by the mantle velocity field (uniformly spaced arrows). The influx boundary conditions, which mimic the stages of a plume push, gradual arrival (1), peak activity (2), and decline (3), are sufficient to initiate a secondary subduction at a weak zone. As the plume push force wanes, the secondary subduction becomes self-sustaining, transitioning the system into a double subduction–driven system (4). The pull from both slabs drives fast convergence and then slows when the middle plate is consumed during arc-continent collision (5) and continental collision (6). Relative motion between the red markers is used to calculate convergence rates.
Here is a map that shows the Zhang et al. (2004) interpretation of how the GPS motion is accommodated in different areas.
Simplified tectonic map showing active faults and movements of Tibetan Plateau and its margins constrained by global positioning system measurements. Numbers are rates of movement (mm/yr). Map covers same region as Figure 1. Bold black lines show active faults. Bold purple arrows indicate N208E shortening across interior of plateau. Blue arrows indicate shortening perpendicular to margins of plateau. Green arrows indicate extension in western and central Tibet. Red arrows indicate strikeslip faulting. Open black arrows denote relative motion with respect to stable Eurasia. Colored numbers show respective amounts
Here Zhang et al (2004) show how they interpret these GPS data to be accommodated along strike-slip fault systems.
Components velocities of global positioning system (GPS) control points along profiles perpendicular to northeastern and eastern margins of Tibet with respect to stable Eurasia. A: Velocity components parallel to N208E vs. distance (in km) from south to north (left to right) along each of four profiles in Figure 1. Blue—profile A-A9, red— B-B9, green—C-C9, and orange—D-D9. Squares show GPS stations on Indian plate; triangles— within Himalaya; diamonds— in interior of Tibetan Plateau; stars— in Qaidam and Qilian Shan; dots—in Tarim or Gobi Alashan regions north of Tibetan Plateau. B: Velocity components parallel to N1108E at stations between 31.58N and 348N and between 798E and 1008E (profile E- 9). Yellow bars denote ranges of velocity measurements. C: Velocity components normal to N208E, inferred India-Eurasia relative motion. Symbols as in A. Curves of different colors sketch approximately eastward flow of each respective profile. D: N108E trending profile orthogonal to Kunlun fault showing components of east-west motion (profile FF9,
Fig. 1). Yellow bars mark range of measurement (~1 standard deviation). Thick red line marks Kunlun fault. Red arrows denote locations of other active faults within this ~400 km zone.
This is the Taylor and Yin (2009) plot that shows their interpretation of the GPS data.
Note how they further note that the strike-slip relative motion can be seen as right-lateral in the south and left-lateral in the north.
Zhu et al. (2019) attempted to use GPS data to model the faults to see where they are accumulating tectonic strain. This map shows how they separated the Earth into different rigid tectonic blocks.
Observed crustal deformation and optimal inversion values. a The fit conditions of the 1999–2013 GPS data. The three inserts are the frequency distributions of north–south and east–west residuals, respectively. Black arrows are observed values, and blue arrows are fitted GPS values. b The location of leveling profile. EP, WP, and MP indicate the best fit of 1970–2013 long-term leveling data. Blue dots show the original observed leveling data, and red dots are fit data from west to east along the three leveling lines. c, d The best fit velocity of 2013–2015 and 2015–2017 GPS data, respectively
Here is a result from Zhu et al. (2019) that shows their estimate of the slip rate on the East Kunlun fault. This is based on an “inversion” of the GPS data along with their block model.
Their paper includes additional analyses to incorporate knowledge about earthquakes that happened during the time GPS data had been recorded.
Slip velocities along the East Kunlun Fault calculated using 1999–2013 GPS velocity data. The 3D drawings represent the results obtained using the inversion strategy proposed in this paper. P1–P4 are the long-term slip velocities and locking depths calculated using the arc tangent method (Savage and Burford 1973). In P1–P4, the blue curves are the fitted curves. The red dots are the observed slip rates parallel to the fault (including errors)
Kirby et al. (2007) use GPS velocities to also interpret strike-slip motion through the region.
The upper map shows the source data for cross sections of those data in the lower plots. Note that there is abundant evidence for strike-slip shear in the western cross sections, but this appears to diminish to the east. Do you agree with this interpretation?
Geodetic velocities across the Kunlun fault. (a) Map of velocity vectors relative to a stable Eurasian reference frame. Data compiled from Shen et al. [2005] and Zhang et al. [2004]. (b–d) Fault-parallel velocities (110 azimuth) in transects located across the central portion of the fault (Figure 6b), the eastern fault tip (Figure 6c), and the Tazang
fault (Figure 6d). Solid lines represent velocities predicted from a dislocation model in an elastic half-space with a locking depth of 15 km and serve to illustrate the range of reasonable slip rates permitted by the data. Position of Kunlun fault shown as vertical line.
Here Kirby et al. (2007) compare their results with the slip-rate estimates from other studies.
Compilation of slip rate estimates and recent seismicity along the Kunlun fault. Open and shaded boxes represent slip rate estimates from previous studies, dark circles new slip rates presented in this work, and heavy vertical lines represent the allowable range of slip rates from geodetic velocities. The lateral extent of previous ruptures along the Kunlun fault are depicted for earthquakes in 1937 [Li et al., 2005; Tapponnier and Molnar, 1977] and 1963 [Tapponnier and Molnar, 1977] and for the 2001 Kokoxili event [Lin et al., 2002; Van der Woerd et al., 2002a]. Note that slip rates (gray shading) decrease markedly toward the eastern tip of the fault, coincident with an absence of historic seismicity and imply strain accumulation in the surrounding Tibetan Plateau.
InSAR Analyses (stay tuned for more)
Shaking Intensity and Potential for Ground Failure
Below are a series of maps that show the shaking intensity and potential for landslides and liquefaction. These are all USGS data products.
There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:
FOS = Resisting Force / Driving Force
When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.
Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.
Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.
Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides. We can see that there is a low probability for landslides. However, we have already seen photographic evidence for landslides and the lower limit for earthquake triggered landslides is magnitude M 5.5 (from Keefer 1984)
— Jason "Jay" R. Patton (@patton_cascadia) May 21, 2021
#EarthquakeReport updated #Earthquake interpretive poster for M7.3 in Central China left-lateral strike slip, along Maduo-Grande fault part of well studied Kunlun fault result of extrusion tectonics historical analogies (2001 M7.8 1997 M 7.5) aftershocks suggest 160km length pic.twitter.com/jK1umjZW1h
— Jason "Jay" R. Patton (@patton_cascadia) May 22, 2021
Quick static stress change modeling caused by the 21/05/2021 Mw7.3 Qinghai Earthquake on parallèle fault-ruptures with Strike/Dip/Rake=91.6°/67.5°/-40° at depth =08km. The red color indicate increase hazard on choosed fault planes. pic.twitter.com/MOobuzUr91
— Jugurtha Kariche (@JkaricheKariche) May 22, 2021
The epicenter corresponds to an unnamed strike-slip fault in the HimalTibet-1.0 compilation by Taylor and Yin (2009). https://t.co/WFMo752CVX
Both this earthquake, and the one that struck 1000 km to the south four hours earlier, are manifestations of the extrusion tectonics so beautifully modelled in 1982 by Tapponnier et al. – crustal blocks escaping to the east in response to the collision of India into Asia. pic.twitter.com/RKJCFa9sRU
— Dr. Judith Hubbard (@JudithGeology) May 22, 2021
Shallow 7.3 Mw EQ.(Qinghai, #CHINA 🇨🇳), probably ruptured along the left-lateral Maduo-Gande Fault, a active fault within the Bayan Har Block.
The data link for LSA is not currently working, here is a similarly impressive seismogram from Kakani, Nepal, in the Himalaya, just at the S edge of the Tibetan Plateau. https://t.co/xmSPjXRWqYpic.twitter.com/ARCvf8JQ16
A second powerful earthquake shook parts of central China on Friday night causing significant damage to buildings and infrastructure. https://t.co/VcKoZqkxd4
Watch the waves from the M 7.3 #earthquake in Southern Qinghai, #China roll across seismic stations in North America. P waves (pale blue) are very subtle, S waves (mostly purple, some orange) are clearer, surface waves are a bit messy on this one. pic.twitter.com/4kNrLsdJJ5
Sand Liquefaction near the epicenter area of the M7.3 Qinghai Earthquake Central China (Credits: Prof. Xiaochen Zhou from China Earthquake Administration) pic.twitter.com/dNc64u1qJY
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
Specific References
Harkins, N., Kirby, E., Shi, X., Burbank, D., and Chun, F., 2010. Millennial slip rates along the eastern Kunlun fault: Implications for the dynamics of intracontinental deformation in Asia in Lithosphere, v. 2, no. 4, p. 247-266, doi: 10.1130/L85.1
Kirby, E., Harkings, N., Wang, E., Shi, X., Fan, C., and Burbank, D., 2007. Slip rate gradients along the eastern Kunlun fault in Tectonics, v. 26, TC2010, doi:10.1029/2006TC002033
Pusok, A.E. and Stegman, D.R., 2020. The convergence history of India-Eurasia records multiple subduction dynamics processes in Science Advances, v. 6 : eeaz8681, 8 p., DOI: 10.1126/sciadv.aaz8681
Staisch, L.M., Niemi, N.A., Clark, M.K., and Hong, C., 2020. The Cenozoic evolution of crustal shortening and left-lateral shear in the central East
Kunlun Shan: implications for the uplift history of the Tibetan Plateau in Tectonics, v. 39, no. 9, https://doi.org/10.1029/2020TC006065
Tapponier, P., Peltzer, G., Le Dain, A.Y., Armijo, R., and Cobbold, P., 1982. Propagating extrusion tectonics in Asia: New insights from simple experiments with plasticine in Geology, v. 10, p. 611-616, doi: 10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2
Taylor, M. and Yin, A., 2009. Active structures of the Himalayan-Tibetan orogen and their relationships to earthquake distribution, contemporary strain field, and Cenozoic volcanism in Geosphere, v. 5, no. 3, doi: 10.1130/GES00217.1
Zhang, P-Z., Shen, Z., Wang, M., Gan, W., Bürgmann, R., Molnar, P., Wang, Q., Niu, Z., Sun, J., Wu, J., Hanrong, S., Xinzhao, Y, 2004. Continuous deformation of the Tibetan Plateau from global positioning system data in Geology, v. 23, no. 9, p. 809-812, https://doi.org/10.1130/G20554.1
This year we look back and remember what happened ten years ago in Japan and across the entire Pacific Basin.
There are numerous web experiences focused on this type of reflection. Here is a short list, some of which I have been involved in.
The California Geological Survey (CGS) prepared several websites devoted to the Tōhoku-oki Earthquake and Tsunami, focusing on how the tsunami impacted California and what has changed in the state since 2011. The entire tsunami program contributed to these expoeriences. Learn more about this program at the CGS Tsunami Program website.
Here is a fantastic web page that Nick Graehl put together using an ArcGIS Story Map system. Don’t miss the two multimedia tours, one for observations in Japan and one for observations in California.
Here are all the pages for this earthquake and tsunami:
I focus mostly on new material I prepared for the following report.
Updated Interpretive Poster
I plot the seismicity from the year after the M 9.1, as well as large events from the past century, with diameter representing magnitude (see legend).
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
In the upper left corner is a small scale plate tectonic map showing the plate boundary faults with the magnetic anomalies overlain in transparency. There is an inset low angle oblique illustrative map showing how these plates interact in the subsurface (Lin et al., 2016).
In the lower right corner is a map that shows a comparison between the USGS modeled earthquake intensity and the USGS Did You Feel It? observations. These data are also included in a web map lower down in this update.
To the left of the intensity map are two tide gage plots that show a tsunami record. The upper plot is from Crescent City, California. The lower plot is from Naha, a location southwest of the earthquake, labeled on tectonic map. These and other tide gage records are viewable in the tide gage web map below.
In the upper right corner are two maps displaying the results from ground failure models from the USGS. The map on the left shows the potential for landslides triggered by the M 9.1 earthquake. The map on the right shows the chance that an area may have experienced liquefaction. These are included in a web map below.
Here is the map with a year’s and century’s seismicity plotted.
Seismicity
Web Map
Use this map to see the magnitudes of different earthquakes experienced in Japan. The map shows earthquake epicenters for large-magnitude historic events of the past century. It also includes epicenters for all aftershocks and triggered earthquakes for a year after the M 9.1 earthquake, and an outline of the aftershocks, which illustrates the area of the fault that slipped during the M 9.1 earthquake.
If you want to see this map in a larger window, click here.
Earthquake Intensity
Earthquake intensity is a measure of how strongly earthquake shaking is felt by people and objects. The further away from the epicenter, the lower the earthquake intensity. Seismologists use computer models to estimate what the intensity will be from an earthquake. The U.S. Geological Survey uses its “Did You Feel It?” (DYFI) system to collect observations about how strongly people in different places felt an earthquake.
Here is a figure that shows a more comparison between the modeled intensity and the reported intensity. Both data use the same color scale, the Modified Mercalli Intensity Scale (MMI). More about this can be found here. The colors and contours on the map are results from the USGS modeled intensity. The DYFI data are plotted as colored dots (color = MMI, diameter = number of reports).
The 3 panels, from left to right, show the USGS Shakemap (the model estimate), the DYFI reports, and an overlay comparing both of these data.
Web Map
Use this map to see the level of intensity people felt in different parts of Japan. The map displays the USGS intensity model for the M 9.1 earthquake as transparent colors. The map also shows, as colored circles, the “Did You Feel It?” report results from people who experienced shaking from this earthquake.
If you want to see this map in a larger window, click here.
Tsunami
Tsunami can be caused by a variety of processes, including earthquakes, volcanic eruptions, landslides, and meteorological phenomena. Earthquakes, eruptions, and landslides cause tsunami when these processes displace water in some way. We may typically associate tsunami with subduction zone earthquakes because these earthquakes are the type that generate vertical land motion along the sea floor.
Here is a great illustration of how a subduction zone earthquake can generate a tsunami (Atwater et al., 2005).
We think that the earthquake slipped at least 50 meters (165 feet) during several minutes. This is the largest coseismic measurement of any subduction zone earthquake (so far).
When the fault slipped, it caused the seafloor to deform and move. This motion also displaced the overlying water column.
As the water column is elevated, it gains potential energy. As this uplifted water expends this energy by oscillating up and down, it radiates energy in the form of tsunami waves.
Tsunami were observed across the entire Pacific Basin, causing extensive damage and casualties in Japan, but also in other places too. There was about $100 million damage to coastal infrastructure in California alone.
This is an animated model of the Great East Japan tsunami of ten years ago. The warmer the colors, the larger the wave. The first surges reached the closest Japan coasts in about 25 minutes. The first surges reached Crescent City in 9.5 hours. (modified text from Dr. Lori Dengler)
This is the same map used as an overlay in the web map below.
Here is the tide gage record from Crescent City, California, USA.
Time is represented by the horizontal axis and elevation is represented on the vertical axis. The darker blue line in this image represents NOAA’s tidal forecast. The data recorded by the tide gage are represented by the light blue colored lines. Wave height is the distance measured between the wave crest and trough. Wave amplitude is the level of water above sea level.
Some of these data came from the IOC sea level monitoring website.
Web Map
Use this map to see tsunami wave data as recorded by tide gages across the entire Pacific Basin. Click on a white triangle and there is a link to open the tide gage data as a graphic.
There is an overlay of color that represents the size of the tsunami as it travelled across the ocean. Learn more about these data here.
If you want to see this map in a larger window, click here.
Ground Failure
Below are a series of maps that show the shaking intensity and potential for landslides and liquefaction. These are all USGS data products.
There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:
FOS = Resisting Force / Driving Force
When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.
Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.
Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.
Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides. We can see that there is a low probability for landslides. However, we have already seen photographic evidence for landslides and the lower limit for earthquake triggered landslides is magnitude M 5.5 (from Keefer 1984)
Use this map to see the magnitudes of different earthquakes experienced in Japan. The map shows earthquake epicenters for large-magnitude historic events of the past century. It also includes epicenters for all aftershocks and triggered earthquakes for a year after the M 9.1 earthquake, and an outline of the aftershocks, which illustrates the area of the fault that slipped during the M9.1 earthquake.
Web Map
If you want to see this map in a larger window, click here.
— Jason "Jay" R. Patton (@patton_cascadia) March 11, 2018
References:
Basic & General References
Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
Hirose et al., 2011. Outline of the 2011 off the Pacific coast of Tohoku Earthquake (Mw 9.0) Seismicity: foreshocks, mainshock, aftershocks, and induced activity in Earth Planets Space, v. 63, p. 655-658
Ikuta et al., 2012. A small persistent locked area associated with the 2011 Mw9.0 Tohoku-Oki earthquake, deduced from GPS data in Journal of Geophysical Research, v. 117, DOI: 10.1029/2012JB009335
Ito et al., 2011. Slip distribution of the 2011 off the Pacific coast of Tohoku Earthquake inferred from geodetic data in Earth Planets Space, v. 63, p. 627-630
Koper et al., 2011. Frequency-dependent rupture process of the 2011 Mw 9.0 Tohoku Earthquake: Comparison of short-period P wave back projection images and broadband seismic rupture models in Earth Planets Space, v. 63, p. 599-602.
Kosuga et al, 2011. Seismic activity around the northern neighbor of the 2011 off the Pacific coast of Tohoku Earthquake with emphasis on a potentially large aftershock in the area in Earth Planets Space, v. 63, p. 719-723.
Lay et al., 2011 a. The 2011 Mw 9.0 off the Pacific coast of Tohoku Earthquake: Comparison of deep-water tsunami signals with finite-fault rupture model predictions in Earth Planets Space, v. 63, p. 797-801.
Lay et al., 2011 b. Possible large near-trench slip during the 2011 Mw 9.0 off the Pacific coast of Tohoku Earthquake in Earth Planets Space, v. 63, p. 687-692.
Lay et al., 2011 c. Outer trench-slope faulting and the 2011 Mw 9.0 off the Pacific coast of Tohoku Earthquake in Earth Planets Space, v. 63, p. 713-718.
Lee et al., 2011. Evidence of large scale repeating slip during the 2011 Tohoku‐Oki earthquake in Geophysical Research Letters, v. 38, DOI: 10.1029/2011GL049580.
Newman et al., 2011. Hidden depths in Nature, v. 474, p. 441-443.
Nishimura et al., 2011. The 2011 off the Pacific coast of Tohoku Earthquake and its aftershocks observed by GEONET in Earth Planets Space, v. 63, p. 631-636.
Orzawa et al., 2011. Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake in Nature, v. 000, p. 1-4.
Satake et al., 2013. Time and Space Distribution of Coseismic Slip of the 2011 Tohoku Earthquake as Inferred from Tsunami Waveform Data in Bulletin of the Seismological Society of America, v. 1032, p. 1473-1492.
Simons et al., 2011. The 2011 Magnitude 9.0 Tohoku-Oki Earthquake: Mosaicking the Megathrust from Seconds to Centuries in Science, v. 332, p. 1421-1425.
Terakawa et al., 2013. Changes in seismic activity following the 2011 Tohoku-oki earthquake: Effects of pore fluid pressure in Earth and Planetary Science Letters, v. 365, p. 17-24.
Toda et al., 2011. Using the 2011 Mw 9.0 off the Pacific coast of Tohoku Earthquake to test the Coulomb stress triggering hypothesis and to calculate faults brought closer to failure in Earth Planets Space, v. 63, p. 725-730.
Uchida and Matsuzawa, 2011. Coupling coefficient, hierarchical structure, and earthquake cycle for the source area of the 2011 off the Pacific coast of Tohoku earthquake inferred from small repeating earthquake data in Earth Planets Space, v. 63, p. 675-679.
Wang et al., 2013. The 2011 Mw 9.0 Tohoku Earthquake: Comparison of GPS and Strong-Motion Data in Bulletin of the Seismological Society of America, v. 103, p. 1336-1347.
Yagi and Fukahata, 2011. Rupture process of the 2011 Tohoku‐oki earthquake and absolute elastic strain release in Geophysical Research Letters, v. 38, DOI: 10.1029/2011GL048701
Yamazaki et al., 2011. Modeling near‐field tsunami observations to improve finite‐fault slip models for the 11 March 2011 Tohoku earthquake in Geophysical Research Letter,s v. 38, DOI: 10.1029/2011GL049130
Yomogida et al., 2011. Along-dip segmentation of the 2011 off the Pacific coast of Tohoku Earthquake and comparison with other megathrust earthquakes in Earth Planets Space, v. 63, p. 697-701.
Yue and Lay, 2013. Source Rupture Models for the Mw 9.0 2011 Tohoku Earthquake from Joint Inversions of High-Rate Geodetic and Seismic Data in Bulletin of the Seismological Society of America, v. 103, p. 1242-1255.