{"id":16256,"date":"2026-07-17T09:06:23","date_gmt":"2026-07-17T07:06:23","guid":{"rendered":"https:\/\/www.proefschriftmaken.nl\/portfolio\/celine-marsman\/"},"modified":"2026-07-17T09:06:31","modified_gmt":"2026-07-17T07:06:31","slug":"celine-marsman","status":"publish","type":"us_portfolio","link":"https:\/\/www.proefschriftmaken.nl\/en\/portfolio\/celine-marsman\/","title":{"rendered":"Celine Marsman"},"content":{"rendered":"","protected":true},"excerpt":{"rendered":"","protected":true},"author":7,"featured_media":16257,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_acf_changed":false,"footnotes":""},"us_portfolio_category":[45],"class_list":["post-16256","us_portfolio","type-us_portfolio","status-publish","post-password-required","hentry","us_portfolio_category-new-template"],"acf":{"naam_van_het_proefschift":"The Tohoku Megathrust Earthquake Cycle","samenvatting":"Aardbevingen langs zogenaamde subductiezones behoren tot de krachtigste en meest verwoestende natuurrampen op aarde. Gedurende de aardbevingscyclus, waaronder de interseismische, coseismische en postseismische fases, veranderen het landschap (met name langs de kust) en de relatieve zeespiegel voortdurend. In de afgelopen eeuw hebben vooruitgang in seismologische en geodetische waarnemingen, samen met gedetailleerde datasets van recente grote aardbevingen, zoals de Sumatra-Andaman aardbeving van 2004 in Indonesi\u00eb, de Maule aardbeving van 2010 in Chili en de Tohoku-Oki aardbeving van 2011 in Japan, bijgedragen aan een beter fysisch begrip van de aardbevingscyclus in subductiezones. Tegelijkertijd zijn numerieke modellen steeds geavanceerder geworden, dankzij de toegenomen rekenkracht en de integratie van aanvullende fysische processen. Tegenwoordig weten we dat interseismische koppeling, co- en postseismische verplaatsingen langs de plaatgrensbreuk, en visco-elastische relaxatie in de mantel, gezamenlijk verantwoordelijk zijn voor de waargenomen oppervlakverplaatsingen. Verschillende combinaties van deze processen kunnen echter vergelijkbare verplaatsingspatronen reproduceren op de beschikbare observatielocaties. Dit leidt tot een probleem dat niet uniek is. Het is nog steeds onduidelijk of dit het gevolg is van verschillende modelleerstrategie\u00ebn of te weinig waarnemingen die gebruikt worden om modelparameters te bepalen die de onderliggende geodynamische processen sturen.\n\nHet doel van dit proefschrift is om bij te dragen aan een beter fysisch begrip van de aardbevingscyclus in subductiezones door de fysische processen te kwantificeren die de waargenomen oppervlakverplaatsingen veroorzaken. Data-assimilatie technieken worden toegepast, waarbij geodetische waarnemingen aan het aardoppervlak worden ge\u00efntegreerd met fysisch onderbouwde modellen van subductiezones. Deze aanpak maakt het mogelijk modelparameters systematisch te schatten en onderlinge afhankelijkheden te analyseren, terwijl onzekerheden in zowel de modellen als de observaties expliciet worden meegenomen. In dit proefschrift staat de Tohoku-Oki aardbeving van 2011 in Japan centraal, vanwege de grote hoeveelheid beschikbare waarnemingen v\u00f3\u00f3r, tijdens en na de aardbeving, waaronder landmetingen en zeebodemmetingen nabij de trog. Een belangrijk aspect van dit proefschrift is dat de interseismische, coseismische en postseismische fases binnen \u00e9\u00e9n rheologisch raamwerk worden beschreven in het model. Hierdoor worden spanningsopbouw en -ontlading over de volledige seismische cyclus op consistente wijze beschreven.\n\nData-assimilatie technieken zijn tot op heden slechts beperkt toegepast binnen dit vakgebied en zijn daardoor relatief onbekend. Om die reden begint het onderzoek in dit proefschrift met het verkennen van het potentieel van data-assimilatie technieken voor het schatten van modelparameters, het kwantificeren van probabilistische onzekerheden en het identificeren van onderlinge afhankelijkheden tussen parameters. Hiervoor gebruiken wij een particle-methode om verticale oppervlakverplaatsingen te assimileren in een elastisch flexuurmodel van de overrijdende plaat tijdens de interseismische periode. Met synthetische experimenten kunnen we analyseren hoe de plaatsing van stations de resultaten be\u00efnvloedt en daarmee de datagevoeligheid. Synthetische experimenten laten zien dat zeebodemmetingen nabij de trog essentieel zijn indien de plaat met lange golflengtes vervormd, terwijl metingen op enkel het land voldoende zijn om kortgolvige deformatie goed te beschrijven. Toepassing van de methode op waarnemingen voorafgaand aan de Tohoku-Oki aardbeving van 2011 resulteert echter in onrealistisch lage schattingen van de effectieve plaatdikte, wat erop wijst dat het vereenvoudigde model onvoldoende realistisch is. Daarnaast blijken 104 tot 105 ensembleleden noodzakelijk om robuuste parameterschattingen te verkrijgen. Ondanks deze beperkingen blijkt de particle-methode waardevol voor het analyseren van datagevoeligheid, het identificeren van parameterafhankelijkheden en het schatten van onzekerheden. Om data-assimilatie technieken echter succesvol toe te passen op waarnemingen uit de praktijk, is een realistischer modelopzet noodzakelijk.\n\nOm de interseismische, coseismische en postseismische waarnemingen van de Tohoku-Oki aardbeving van 2011 te reproduceren, wordt vervolgens een realistischer tweedimensionaal (2-D) aardbevingscyclusmodel ontwikkeld. Hierbij wordt gebruikgemaakt van een raamwerk waarin spanningsopbouw en -ontlading consistent worden beschreven. De geometrie van de subducerende plaat en het verplaatsingsveld langs de plaatgrensbreuk worden als a priori informatie in het model verwerkt. De viscositeitsstructuur is afhankelijk van de temperatuur en spanning. Met behulp van de ensemble smoother met multiple data assimilation (ESMDA) worden modelparameters geschat die de viscositeitsstructuur van de mantel en aseismische verplaatsingen langs de plaatgrensbreuk bepalen. ESMDA maakt het mogelijk deze parameters te schatten met slechts 100 ensembleleden.\n\nData-assimilatie experimenten met het 2-D model, waarin verschillende keuzes voor de mechanische eigenschappen van de mantel worden onderzocht, benadrukken het belang van niet-lineaire dislocatiekruip in de asthenosfeer, in overeenstemming met laboratoriumexperimenten, evenals van ondiepe visco-elastische relaxatie (30\u201360 km diepte) in de continentale mantel om voldoende flexuur te reproduceren. Na de aardbeving vinden postseismische verplaatsingen vooral plaats dieper op het subductiecontact. Deze verplaatsingen zijn te beschrijven met een niet-lineaire kruip, waarbij diepteafhankelijke viscositeitsvariaties wijzen op thermische afhankelijkheid. De resultaten tonen bovendien aan dat rheologische keuzes sterk bepalen in hoeverre postseismische verplaatsingen langs de plaatgrensbreuk en visco-elastische relaxatie in de mantel van elkaar kunnen worden onderscheiden. Het opnemen van niet-lineaire rheologie vermindert deze ambigu\u00efteit en maakt het mogelijk de parameters van beide mechanismen beter te bepalen, ook al kunnen we de individuele signalen aan het oppervlak niet expliciet bepalen vanwege de niet-lineariteit van het model. Met name ruimtelijke patronen van verticale verplaatsing en de vroege postseismische ontwikkeling van horizontale bewegingen zijn belangrijk voor het schatten van de afzonderlijke parameters die postseismische verplaatsingen langs de plaatgrensbreuk en visco-elastische relaxatie sturen.\n\nDe postseismische 2-D simulaties laten verder zien dat versterkte landwaartse beweging nabij de trog kan worden verklaard door een combinatie van twee mechanismen: (1) relaxatie onder de subducerende plaat en (2) herstelbeweging, die de respons van de overrijdende plaat beschrijft ten opzichte van een referentiesimulatie zonder koppeling tussen de platen. Mechanische koppeling tussen de overrijdende plaat en de subducerende plaat is noodzakelijk om de effecten van relaxatie in de oceanische mantel over te dragen op de overrijdende plaat. Deze herstelbeweging weerspiegelt het geleidelijke herstel van de interseismische verkorting van de oceanische plaat tijdens de vroege postseismische fase en wordt pas zichtbaar wanneer de volledige aardbevingscyclus wordt gesimuleerd. Dislocatiekruip onder de subducerende plaat, met viscositeiten van 7.0 \u2212 8.6 \u00d7 1018 Pa\u00b7s gedurende het eerste jaar na de aardbeving, draagt bij aan de waargenomen landwaartse beweging nabij de trog. Wij laten zien dat een expliciete laag met lage viscositeit onder de subducerende plaat niet noodzakelijk is.\n\nWij breiden het onderzoek uit naar een driedimensionaal (3-D) model om rekening te houden met de laterale variabiliteit die wordt waargenomen in de data binnen Honshu, evenals met variaties in de geometrie van de subducerende plaat en de trog. Het 3-D aardbevingscyclusmodel maakt het mogelijk zowel diepteafhankelijke als laterale variaties in koppeling, co- en postseismische verplaatsingen langs de plaatgrensbreuk, spanningen en rheologie te beschrijven. Het doel is de heterogene viscositeitsstructuur van de mantel en van aseismische verplaatsingen langs de plaatgrensbreuk in kaart te brengen, en daarmee geodetische waarnemingen in Honshu en Hokkaido te kunnen verklaren. Mantelstroming en aseismische verplaatsingen langs de plaatgrensbreuk worden gemodelleerd met een samengestelde rheologie, bestaande uit een spanningsafhankelijke viscositeit en een minimale achtergrondviscositeit, waardoor zowel vroege als langdurige deformatieprocessen kunnen worden weergegeven. Met behulp van ESMDA worden de rheologische parameters geschat die deze samengestelde stromingswet beschrijven. Interseismische en postseismische geodetische waarnemingen uit Honshu worden geassimileerd om deze parameters te schatten.\n\nDe data-assimilatie resultaten van de 3-D experimenten laten eveneens zien dat er geen significante afhankelijkheden bestaan tussen parameters die aseismische verplaatsingen langs de plaatgrensbreuk en visco-elastische relaxatie sturen. Dit toont aan dat beide processen elk een duidelijk te onderscheiden effect hebben op de gecombineerde interseismische (3.4 jaar) en postseismische (10 jaar) verplaatsingstijdreeksen. Aseismische verplaatsingen langs de plaatgrensbreuk treden voornamelijk op onder het coseismisch ontkoppelde gedeelte, op dieptes tussen ongeveer 40 en 80 km, met maximale verplaatsingen rond 50 km diepte. Het grootste deel van deze verplaatsingen vindt plaats in de eerste vijf jaar na de aardbeving en bereikt ongeveer 5 m na 10 jaar. Laterale regio\u2019s vertonen daarentegen kleinere aseismische verplaatsingen, met maximaal 1.6 m na 10 jaar, maar nemen langzamer af in de tijd in vergelijking met de regio\u2019s direct onder het coseismisch ontkoppelde gebied.\n\nDe data-assimilatie experimenten laten zien dat het ondiepe gedeelte van de continentale mantel lage viscositeiten nodig heeft, met effectieve viscositeiten in de orde van 1019 Pa\u00b7s op dieptes van ongeveer 40\u201360 km, waar deformatie wordt gedomineerd door dislocatiekruip, in overeenstemming met de 2-D experimenten. Dieper in de continentale mantel is een lateraal uniforme achtergrondviscositeit van orde 1020 Pa\u00b7s vereist om mantelstroming in koudere gebieden mogelijk te maken. Dit benadrukt de noodzaak van een samengestelde rheologie die niet kan worden beschreven met een uitsluitend temperatuurafhankelijke stromingswet. Dislocatiekruip domineert de vroege postseismische deformatie, maar de bijdrage in de reksnelheid neemt af van circa 90% in de eerste 2.5 jaar tot ongeveer 30% in de periode 7.5\u201310 jaar na de aardbeving. Onder de subducerende plaat is dislocatiekruip bepalend voor het postseismische deformatiepatroon, met effectieve viscositeiten van orde 5 \u00d7 1018 Pa\u00b7s op dieptes van 100\u2013200 km gedurende de eerste tien jaar na de aardbeving.\n\nDe mechanismen die postseismische landwaartse beweging nabij de trog verklaren, zijn vergelijkbaar met die in de 2-D modellen en omvatten zowel het herstel van de plaat nabij de trog als viskeuze stroming onder de plaat. Hoewel enige aseismische slip is toegestaan in het ondiepe deel van de plaatgrensbreuk nabij de trog, in tegenstelling tot volledige koppeling, vereist dit versterkte stroming onder de plaat om de waargenomen oppervlakte-deformatie te reproduceren. Dislocatiekruip in de oceanische mantel blijft het dominante deformatiemechanisme gedurende de postseismische periode van 10 jaar. De bijbehorende viscositeitsvermindering is lateraal beperkt tot ongeveer 200 km ten noorden en zuiden van het centrum van de breuk.\n\nOnze 3-D modellen geven aan dat verbeterde weergaven van zowel koppeling als aseismische slip vereist zijn, met name langs de laterale randen van de breukzone. Ten noorden van de breukzone wordt aseismische slip waarschijnlijk onderschat, wat wordt ondersteund door vergelijkingen met aseismische slip afgeleid uit herhalende aardbevingen. Ondanks deze beperkingen reproduceert het 3-D model met succes het eerste-orde patroon van postseismische landwaartse beweging in Hokkaido, waarbij de omvang, richting en persistentie gedurende de postseismische fase worden vastgelegd. Resterende discrepanties zijn hoogstwaarschijnlijk toe te schrijven aan vereenvoudigde weergaven van aseismische slip en aan onzekerheden in de elastische parameters van de overrijdende plaat. Het onderliggende mechanisme dat verantwoordelijk is voor de waargenomen landwaartse beweging wordt verondersteld de in-het-vlak-buiging van de overrijdende plaat te zijn als reactie op postseismische relaxatie.\n\nAl met al draagt dit proefschrift bij aan een beter fysisch begrip van de aardbevingscyclus langs subductiezones. Wij laten zien hoe data-assimilatie technieken kunnen worden ingezet om geodetische waarnemingen te integreren met de geofysische processen die oppervlakverplaatsingen bepalen. Onze toepassing op de Tohoku-Oki aardbeving van 2011 levert nieuwe inzichten op. Ten eerste blijkt dat de beschikbare geodetische data voldoende informatie bevatten om de parameters die aseismische verplaatsingen langs de plaatgrensbreuk en de visco-elastische structuur van de mantel bepalen. Wij tonen dit aan met modellen waarbij spanningsopbouw en -ontlading consistent worden beschreven gedurende de volledige aardbeving- scyclus. Ten tweede tonen de resultaten aan dat aseismische verplaatsingen langs de plaatgrensbreuk en visco-elastische relaxatie elk een duidelijk te onderscheiden effect hebben op de waargenomen oppervlakverplaatsingen en daarom kunnen worden onderscheiden wanneer realistische modellen worden gebruikt. Ten derde tonen we aan dat een zone met lage viscositeit onder de subducerende plaat op natuurlijke wijze ontstaat en geen bijzondere aannames vereist, zoals recent is gesuggereerd. Ten vierde kan de postseismische landwaartse beweging in Hokkaido succesvol worden gereproduceerd met realistische modellen. Dit onderzoek vormt een basis voor toekomstige studies waarin steeds grotere hoeveelheden geodetische data, waaronder zeebodemwaarnemingen, worden ge\u00efntegreerd. Daarmee kan het fysisch begrip van grote subductieaardbevingen verder worden verdiept.","summary":"Megathrust earthquakes are among the most powerful and destructive natural events on Earth. Throughout the megathrust earthquake cycle, including interseismic, coseismic and postseismic phases, the coastal landscape and relative sea level continuously evolve. During the past century, progress in seismological and geodetic observations, together with unprecedented datasets from recent great earthquakes such as the 2004 Sumatra-Andaman earthquake in Indonesia, 2010 Maule earthquake in Chile, and 2011 Tohoku-Oki earthquake in Japan, have contributed to our understanding of processes driving surface deformation. In parallel, numerical models have become more sophisticated through advances in computing resources and the inclusion of additional physical mechanisms. It is now recognized that interseismic locking, coseismic slip, afterslip, and viscoelastic relaxation all shape the observed deformation, yet different combinations of these processes can reproduce similar displacement patterns at the available data locations, underscoring the non-uniqueness of the inverse problem. It is still unclear whether the non-uniqueness comes from different modelling strategies or how well observational data constrain model parameters that control the underlying geodynamic processes.\n\nThe aim of this thesis is to advance our understanding of the earthquake cycle at subduction zones by quantifying the physical processes that drive observed surface deformation. To achieve this, data assimilation techniques are used to integrate geodetic surface observations with physics-based subduction zone models. This approach enables systematic estimation of model parameters and their trade-offs, while accounting for uncertainties in both models and data. The 2011 Tohoku-Oki earthquake serves as a central study case due to extraordinary observational coverage before, during and after the earthquake, including both onshore and offshore observations near the trench. An important element of this work is the use of models that treat interseismic, coseismic, and postseismic phases within a single rheological framework, ensuring a consistent description of strain accumulation and release throughout the entire seismic cycle.\n\nAs data assimilation techniques are relatively new to earthquake cycle modelling, the research begins by investigating their potential to estimate model parameters, quantify their probabilistic uncertainties and identify parameter trade-offs. A particle method is applied to assimilate vertical surface displacements into an elastic flexural model of the overriding plate for interseismic deformation. Synthetic experiments show that seafloor data near the trench are essential to constrain long-wavelength deformation, while onshore data suffice for shorter wavelengths. When applied to real observations prior to the 2011 Tohoku-Oki event, the method yields unrealistically low estimates of the effective elastic plate thickness, highlighting the limitations of the simplified forward model. Moreover ensemble sizes of the order 104 to 105 were required to obtain robust model estimates. Despite these shortcomings, the particle method proves useful for exploring data sensitivity, identifying parameter trade-offs, and providing uncertainty estimates, thereby demonstrating its potential for constraining geodynamic parameters. Its application to real observations highlights the need for a more realistic model setup.\n\nTo reproduce interseismic, coseismic and postseismic observations of the 2011 Tohoku-Oki earthquake, a more realistic two-dimensional (2-D) earthquake cycle model is developed within a physically consistent framework for strain accumulation and release. This setup incorporates slab geometry, an a priori coseismic slip distribution, and a temperature- and stress-dependent viscosity structure. We estimate model parameters that determine the viscosity structure of the mantle and the shear zone using the ensemble smoother with multiple data assimilation (ESMDA). ESMDA made it possible to estimate model parameters with as little as 100 ensemble members.\n\nData assimilation experiments with the 2-D model, exploring different rheological formulations of the mantle and shear zone, highlight the role of nonlinear dislocation creep in the asthenosphere, consistent with laboratory constraints, and of shallow viscoelastic relaxation (30\u201360 km depth) in the continental mantle to reproduce plate bending. Afterslip is localized downdip of the coseismic rupture and requires nonlinear creep in the shear zone, with depth-dependent viscosity variations reflecting thermal control. Additionally, the data assimilation results show that rheological choices strongly influence the degree to which afterslip and viscoelastic relaxation can be distinguished. Incorporating nonlinear rheology reduces ambiguity and allows the controlling parameters of these processes to be more tightly constrained, even if their surface signals could not be fully separated due to the nonlinear model. In particular, spatial vertical displacement patterns and the early postseismic evolution of horizontal motion provide important constraints for estimating the individual parameters that control afterslip and viscoelastic relaxation.\n\nThe 2-D postseismic simulations further reveal that enhanced near-trench landward motion are reproduced by a combination two processes: 1) sub-slab relaxation and 2) recovery motion, representing the upper plate response relative to a free-slip reference simulation. Full mechanical coupling between the overriding plate and the slab is required to transfer the effects of sub-slab relaxation to the overriding plate. This recovery motion reflects the gradual recovery of interseismic shortening of the oceanic plate during the early postseismic phase, which becomes evident only when simulating the full earthquake cycle. Dislocation creep below the slab, with viscosities of 7.0\u22128.6\u00d71018 Pa\u00b7s during the first year after the earthquake, contributes to the observed near-trench landward motion without requiring a low-viscosity lithosphere\u2013asthenosphere boundary.\n\nWe extend our investigation to a fully three-dimensional (3-D) framework to account for the pronounced lateral variability observed across Honshu, as well as along-strike variations in slab geometry and trench configuration. The 3-D earthquake cycle model enables us to resolve both depth dependent and lateral variations in locking, coseismic slip, stress, rheology, and fault slip. Our objective is to constrain the heterogeneous viscosity structure of the mantle and the spatiotemporal variability of aseismic slip along the megathrust, thereby explaining geodetic observations across both Honshu and Hokkaido. Mantle flow and afterslip are represented using a composite rheology consisting of a stress-dependent power-law component and a background constant viscosity component, allowing us to capture both transient and steady-state deformation processes. We apply ESMDA to estimate rheological parameters governing the composite flow law in both the mantle and the shear zone. Interseismic and postseismic geodetic observations from Honshu are assimilated to inform these estimates.\n\nData assimilation results from the 3-D experiments indicate that there are no significant trade-offs between parameters controlling afterslip and those governing viscoelastic relaxation, demonstrating that each process leaves a distinct and identifiable signature in the combined 3.4- year interseismic and 10-year postseismic displacement time series. Afterslip is inferred to occur predominantly downdip of the coseismic rupture, at depths between approximately 40 and 80 km, with peak slip concentrated near 50 km depth. The majority of this downdip afterslip is accommodated within the first five years following the earthquake, reaching cumulative slip amplitudes of up to \u223c5 m after 10 years. In contrast, lateral regions experience smaller amounts of afterslip, with peak cumulative values of approximately 1.6 m after 10 years, but exhibit a slower temporal decay compared to the downdip region.\n\nMantle wedge viscosities inferred from the data assimilation indicate a low-viscosity mantle wedge, with effective viscosities on the order of 1019 Pa\u00b7s at depths of approximately 40\u201360 km, where deformation is dominated by dislocation creep, consistent with the 2-D model experiments. At greater depths within the mantle wedge, a laterally uniform background viscosity of order 1020 Pa\u00b7s is required to permit mantle flow in colder regions, highlighting the necessity of a composite rheology that cannot be represented by a purely temperature-dependent power-law formulation. Beyond depths of approximately 100 km, dislocation creep dominates the early postseismic deformation but its relative contribution to the total strain rate decreases with time, from about 90% in the first 0\u20132.5 years to roughly 30% after 7.5\u201310 years. Beneath the slab, dislocation creep plays a key role in shaping the observed postseismic deformation pattern, with effective viscosities on the order of 5 \u00d7 1018 Pa\u00b7s at depths of 100\u2013200 km over the first 10 years following the earthquake.\n\nThe mechanisms driving postseismic landward motion near the trench are similar to those identified in the 2-D models and include unshortening of the slab near the trench as well as viscous flow beneath the slab. While some aseismic slip is permitted in the shallow part of the megathrust near the trench, in contrast to complete locking, this requires enhanced sub-slab flow to reproduce the observed surface deformation. Dislocation creep in the oceanic mantle remains the dominant deformation mechanism throughout the 10-year postseismic period. The associated viscosity reduction is laterally confined to within approximately 200 km north and south of the rupture centre.\n\nOur 3-D models indicate that improved representations of both locking and afterslip are required, particularly along the lateral edges of the rupture zone. North from the rupture zone, afterslip is likely underestimated, as supported by comparisons with aseismic slip inferred from repeating earthquakes. Despite these limitations, the 3-D model successfully reproduces the first- order pattern of postseismic landward motion in Hokkaido, capturing its magnitude, direction, and persistence throughout the postseismic phase. Remaining discrepancies are most likely attributable to simplified representations of afterslip and to uncertainties in the elastic parameters of the overriding plate. The underlying mechanism responsible for the observed landward motion is inferred to be in-plane bending of the overriding plate in response to postseismic relaxation.\n\nOverall, this thesis advances our understanding of the megathrust earthquake cycle by demonstrating how data assimilation can be used to link geodetic observations to the geophysical processes governing surface deformation. Application to the 2011 Tohoku-Oki earthquake brings new insights, 1) the spatiotemporal evolution of surface deformation contains sufficient information to constrain the parameters controlling aseismic slip along the megathrust and the viscoelastic structure of the mantle, within a modelling framework that provides a consistent formulation of strain accumulation and release throughout the earthquake cycle. 2) The results demonstrate that afterslip and viscoelastic relaxation have a unique footprint on the observed surface displacements and can therefore be unravelled when realistic models are used. 3) A low-viscosity zone emerges naturally beneath the slab and does not require special circumstances as recently advocated. 4) Postseismic landward motion in Hokkaido can be successfully reproduced when realistic models are used. This work lays the groundwork for future studies that integrate increasingly dense geodetic data, including seafloor observations, ultimately improving our ability to interpret deformation signals and to understand the mechanics of megathrust earthquakes.","auteur":"Celine Marsman","auteur_slug":"celine-marsman","publicatiedatum":"22 september 2026","taal":"EN","url_flipbook":"https:\/\/ebook.proefschriftmaken.nl\/ebook\/celinemarsman?iframe=true","url_download_pdf":"https:\/\/ebook.proefschriftmaken.nl\/download\/a473b873-e45e-49c3-a57b-6221091ebc2d\/optimized","url_epub":"","ordernummer":"19292","isbn":"978-90-6266-750-5","doi_nummer":"","naam_universiteit":"Universiteit Utrecht","afbeeldingen":16258,"naam_student:":"","binnenwerk":"","universiteit":"Universiteit Utrecht","cover":"","afwerking":"","cover_afwerking":"","design":""},"_links":{"self":[{"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/us_portfolio\/16256","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/us_portfolio"}],"about":[{"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/types\/us_portfolio"}],"author":[{"embeddable":true,"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/users\/7"}],"replies":[{"embeddable":true,"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/comments?post=16256"}],"version-history":[{"count":1,"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/us_portfolio\/16256\/revisions"}],"predecessor-version":[{"id":16259,"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/us_portfolio\/16256\/revisions\/16259"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/media\/16257"}],"wp:attachment":[{"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/media?parent=16256"}],"wp:term":[{"taxonomy":"us_portfolio_category","embeddable":true,"href":"https:\/\/www.proefschriftmaken.nl\/en\/wp-json\/wp\/v2\/us_portfolio_category?post=16256"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}