Earthquake energy dissipation in a fracture mechanics framework

Earthquakes are rupture-like processes that propagate along tectonic faults and cause seismic waves. The propagation speed and final area of the rupture, which determine an earthquake’s potential impact, are directly related to the nature and quantity of the energy dissipation involved in the rupture process. Here, we present the challenges associated with defining and measuring the energy dissipation in laboratory and natural earthquakes across many scales. We discuss the importance and implications of distinguishing between energy dissipation that occurs close to and far behind the rupture tip, and we identify open scientific questions related to a consistent modeling framework for earthquake physics that extends beyond classical Linear Elastic Fracture Mechanics.

Earthquakes are one of the most damaging natural hazards facing humankind.Improvements in understanding the fundamental physics of earthquakes could have dramatic consequences for our ability to plan and react to catastrophic earthquakes in densely populated areas.Seismological observations show that earthquakes comprise a rupture front propagating along a fault and leaving behind slip and stress drop, which is a form of fracture propagation.Thus the field of fracture mechanics has played a fundamental role in shaping what we know about earthquake physics.Classical models describe an earthquake as a shear crack and define, for example, the relationship between earthquake rupture area, propagation speed, and the spectral characteristics of radiated seismic waves [1][2][3] .While the overall energy budget that compares states before and after an earthquake has a wellestablished theoretical basis 4 , key aspects of the instantaneous energy balance governing the behavior of the earthquake rupture remain poorly understood.Fracture mechanics theory predicts that rupture growth is a balancing act involving three components: 1) energy dissipated to extend the crack, either by creating new surface area or generating frictional heat, 2) energy radiated as seismic waves, and 3) release of stored elastic energy from the surrounding rock.This view has been confirmed by a broad range of laboratory experiments and codified in the theory of Linear Elastic Fracture Mechanics (LEFM) 5 .However, the complexity of earthquake faults far exceeds that of typical laboratory setups, raising significant questions about the applicability and predictive power of LEFM in natural conditions 6,7 .One of our goals here is to paint a picture of the state of the art in understanding earthquake rupture and, in particular, the extent to which LEFM can further our understanding of the earthquake energy budget in cases where the fault zone has finite width and rupture propagation involves branching, off-fault fracturing and other processes that go beyond the simple assumptions underlying LEFM.
While there is transformative potential in extending knowledge and connections between the fields of fracture mechanics and earthquake physics, there are several key impediments.These include: 1) a lack of fundamental understanding of how various dissipative processes at different spatial and temporal scales contribute to the mechanics of earthquakes, 2) extreme discrepancies (of many orders of magnitude) between values of fracture energy measured in laboratory experiments 8 and inferred from natural earthquakes, and 3) vastly different terminology between the communities.In this Perspective, we aim to review the mechanics and energy dissipation in earthquake ruptures, define a clear terminology, and discuss the capabilities and limitations of current observations and measurement techniques and how they affect the observed discrepancies.Finally, we propose a path forward in the form of key outstanding questions, clear scientific objectives for future work, and suggestions to overcome the limitations of LEFM as applied to earthquake faulting.

Fundamentals of theoretical earthquake mechanics
The mechanics of earthquakes are complex and arguably equally challenging to measure as to theoretically describe.The goal of a theoretical earthquake model is to describe the essential processes with viable equations and tools that allow field observation to be interpreted.Hence, it builds on common observations, which show that earthquake ruptures begin in a localized region of a fault known as the hypocenter (see Fig. 1), which is the location where initial shear stresses (τ 0 ) are sufficient to overcome frictional strength and the fault motion begins to accelerate.The initial stress level is one of the most difficult parameters to constrain.It can be spatially variable, and it can have a large impact on rupture style and speed.Unlike static frameworks, such as slip-tendency analysis popular in structural geology, faults can be stressed well below strength almost everywhere and yet rupture spontaneously: only a small portion of the fault needs to reach its strength to nucleate an earthquake.For ordinary earthquakes with fast rupture speeds, the hypocenter is the location where seismic waves are first radiated.From there, earthquake rupture expands along the fault, causing the fault surfaces to begin to slip.The transition region between slipping and unslipped sections of the fault is called the rupture front (see Fig. 1).For fast ruptures, moving at speeds of several km/s, the slip rate of fault surfaces accelerates from below 1 μm/s to above 1 m/s over timescales of less than a second.Ahead of the rupture front, the shear stress on the fault increases in what is called a dynamic stress concentration.At the rupture tip, a rapid transition occurs.The slip rate increases as the shear stress drops rapidly to a dynamic level τ d that is below the initial level τ 0 .This drop in stress causes a release of stored strain energy, which drives the rupture.Simultaneously, part of this energy is dissipated through various processes, including fracture of the surrounding rock, comminution (the production of rock powders), heating, and possibly melting of rocks.The remaining energy is radiated as seismic waves that transport kinetic energy far from the source and cause ground shaking, ultimately to be dissipated as heat.Eventually, the rupture ceases to grow, and all sections of the earthquake rupture area arrest (or evolve to a very slow and long quasi-static front of postseismic slip).This can occur because continued slip necessitates a disproportionately large amount of energy dissipation or because the rupture front propagates into unfavorably stressed regions (i.e., τ 0 < τ d ).

Earthquakes as a rupture process described by fracture mechanics
The earthquake rupture process of propagation and arrest shares many features of a crack propagating through a solid material.Thus, the theoretical framework of LEFM, as summarized in Box 1, has been adapted to describe the mechanics of earthquake rupture [e.g., [9][10][11] ].Specifically, under suitable assumptions, LEFM provides an energy balance that governs rupture growth: where G is the energy release rate and Γ tot the fracture energy (precise definitions are provided in Box 1).LEFM allows rupture speed prediction based on a few assumptions given in Box 1 where C f is the rupture speed, C R is the Rayleigh wave speed, and G 0 is the static energy release rate.However, the application of LEFM to earthquake ruptures requires a few key adaptions from its classical form.For instance, LEFM was initially developed for cracks that were traction-free behind the rupture front (Box 1).While valid for opening cracks, this is clearly not valid for surfaces in frictional contact.However, for many of the equations commonly used to model frictional ruptures, including rate-and-state friction equations [e.g., 12,13 ], the shear stress well behind the rupture front varies only modestly with time and position.For these cases, one can approximate the dynamic friction as reaching a constant residual level τ r for slip larger than a characteristic slip distance d c .While this is inconsistent with expectations for complex natural faults, it provides a starting place for the application of LEFM and an opportunity to investigate the earthquake energy budget and rupture dynamics.
In the attempt to measure energy dissipation, the work done by the frictional stresses above the minimal stress through slip, commonly known as the "breakdown work" W b , is measured/computed (see Box 2).If residual friction is constant and d c is small enough, then the breakdown work is spatially localized in the process zone near the rupture tip.This would ensure the above-mentioned separation of scale that allows the application of LEFM because the energy dissipation associated with fracture propagation is determined entirely within the crack tip region and, therefore, independent of other processes on the rest of the rupture surface.In this particular case, the breakdown work is exactly equal to the associated fracture energy, i.e., Γ tot = W b (assuming that there are no other dissipative processes).Some important comments: 1. breakdown work ≠ fracture energy: breakdown work W b and associated fracture energy are not generally equal, as commonly assumed in some literature.The breakdown work is only equal to the fracture energy under specific conditions, i.e., when there is separation of scale (see Box 1).The precise limit of separation of scale remains unknown.Since only the localized part of the breakdown work is part of the fracture energy, Γ tot ≤ W b (in the absence of other dissipative processes).2. fracture energy and physical processes: the fracture energy Γ tot governing crack growth in Eqs. ( 1) and ( 2) is the cumulative quantity that may include the effects of many processes such as work done by frictional weakening, plastic dissipation, off-fault damage, pulverization within the fault zone, and others, i.e., Γ tot = ∑Γ processes .However, Γ tot only includes the part of the energy that is localized near the crack tip (on-or off-fault) to comply with separation of scale and to be consistent with LEFM.

Fundamentals of Linear Elastic Fracture Mechanics (LEFM)
Linear Elastic Fracture Mechanics (LEFM) is a theoretical framework to describe crack growth, originally developed with a focus on opening (Mode I) cracks [90][91][92] .The question of whether a crack grows or not is reduced to a comparison of two states (see Box Figure ): the current state with crack (half-)length ℓ and an incremental state with crack length ℓ + dℓ.The following assumptions are made: 1. the material surrounding the crack has predominantly linear elastic behavior, 2. dissipation is localized at the crack tip, 3. the crack surface is traction-free.
One can then express the comparison of states in terms of energies stating that the crack grows if G 0 > Γ tot , where G 0 is the static energy release rate and Γ tot the fracture energy 93,94 .The static energy release rate is the drop in the total potential energy Π per increment of crack length given as G 0 = − dΠ/dℓ, where dΠ ≈ (U 2 − U 1 ) − F with U i being the stored elastic strain energy for state i (see Box Figure ) and F the work done by external forces in between (e.g.gravitational).The fracture energy is the energy dissipated in the process of breaking the material per unit of crack length growth.This comparison of energies is only possible because the above-mentioned assumptions guarantee a "separation of scale" between the global quantity G 0 driving the crack and the local quantity Γ tot resisting crack growth.The key aspect of this approach is that the specific physical processes of how energy is dissipated (e.g., as decohesion, plastic work) are irrelevant as long as they are localized at the crack tip, satisfying assumption 2. They are all lumped into the parameter Γ tot .
To avoid infinite stresses at the crack tip, the dissipation is often smeared out in a still-localized-enough "process zone" of size s ≪ ℓan approach known as "small-scale yielding" and commonly implemented via cohesive zone models 9,10,95 .While the precise limit for small-scale yielding remains unknown, a process zone with s ≈ 0.4ℓ may, under some circumstances, still be enough localized [see pp.LEFM also predicts the speed at which the crack grows 5 through an energy balance, G = Γ tot , where the dynamic energy release rate G≤G 0 accounts for kinetic energy being radiated away from the rupture tip.Assuming time-invariant loading (i.e., no wave reflections) and 2D or circular 3D configurations, the dynamic energy release rate can be approximated by G ≈ (1 − C f /C R )G 0 , where C f = dℓ/dt is the rupture speed and C R the material Rayleigh wave speed.
Box Figure | Schematic of the fundamental principles for crack growth in Linear Elastic Fracture Mechanics.a The current state of a material system under remote applied load (gray arrows) containing a rupture of length ℓ with a process zone (marked green) at the rupture tip.b The incremental state is the same system as shown in a but the rupture length increased by an incremental length dℓ.
3. fracture energy variability: due to its multi-physical origin, the fracture energy may vary spatially along the fault and change with rupture speed.Hence, it may not be known a-priori.4. non-localized heat: the so-called "frictional heat," e.g., W H = R D 0 τ r dδ, which is the work of residual friction τ r on the fault, is not included in Γ tot as it is not localized.Consequently, the energy due to the residual stress is also excluded from G, and the Griffith energy balance (Eq.( 1)).We note that while the majority of the total dissipation goes into frictional heat, this does not detract from the importance of the fracture energy for rupture propagation.

Applicability of fracture mechanics to laboratory ruptures versus natural earthquakes
Frictional stick-slip events or fracture propagation on pre-existing surfaces represent the laboratory equivalent of earthquakes.While earthquakes generated in the laboratory (so-called "labquakes") share many features of natural earthquakes 14 on tectonic faults, the vast differences in scale raise important questions that include the application of LEFM to labquakes and earthquakes on natural faults.Yet, recent works provide a useful starting point because they include quantitative predictions of rupture speed 15,16 and arrest [17][18][19] for labquakes.These and similar experiments also measure the fracture energy Γ tot of labquakes with local dynamic shear stress measurements [20][21][22][23] and/or the stress-versus-slip relation 24,25 .Such experiments require a sample that is large compared to the process zone size and a critical length scale for rupture nucleation.For fault normal stresses of 1 − 10MPa this requires meter-scale rock samples 22,23 or 20 − 30cm sized samples composed of glassy polymer such as PMMA [e.g., 20,26 ].Dynamic rupture can also be studied on smaller samples, at higher normal stress levels (50 − 150MPa), if arrays of sensors are installed on the sample, inside a pressure vessel 27 and for cases where the fault zone contains sufficient wear material 28 .To infer rupture-related quantities, it is important to measure stress evolution on or near the fault as a rupture front propagates past the sensor location as opposed to sample-wide averages.
In summary, the simplifying assumptions of LEFM appear to be valid for large-scale laboratory experiments where LEFM quantitatively describes rupture speed and arrest.However, laboratory experiments differ from natural faults in several ways that must be accounted for to understand the limitations of LEFM and to develop appropriate extensions of that theory.First, the magnitude range of labquakes is relatively limited, which impedes a precise determination of earthquake scaling properties.Furthermore, laboratory experiments are often conducted at low stress levels ( ~5MPa) compared to the unconfined strength of the rock or polymer samples.This limits offfault damage or inelastic deformation that may strongly affect the rupture process and the overall energy dissipation.The experiments also typically employ simple fault geometries, while tectonic faults are much more complex.Under some conditions, the complexities of geological faults can be lumped into a single tip-localized parameter Γ tot ; however, for other cases, the framework of LEFM requires modification.A key question is if energy dissipation (aside from frictional heat) on natural faults with all its complexities (e.g., weakening, offfault inelasticity) truly is localized in the vicinity of the rupture tip, which would guarantee separation of scale and applicability of LEFM.Should this not be the case, is it enough if "most" of the energy is dissipated in a localized manner?The implications of these questions are important as the answers determine the extent to which LEFM can be applied to earthquakes in its current or modified forms.
While numerical simulations are a powerful tool to model earthquakes in complex systems, they also face significant challenges.As a result, they are not yet able to capture all processes at all space and time scales of the earthquake cycle and require significant further algorithmic and computational development to achieve fully realistic scenarios.Furthermore, theoretical tools such as LEFM are useful even if simulations can provide accurate results because they allow us to understand the "why" related to the obtained simulation results.This synergy can only be achieved if the theoretical model is correct and accounts for the relevant processes appropriately.

Tip or tail: spatiotemporal energy dissipation in earthquakes
Tectonic faulting complexity involves simultaneous dissipative processes during an earthquake, such as fracturing, comminution, heating, and possibly rock melting.These processes depend on the fault slip rate and the thickness of the shearing zone within the fault.Mechanisms such as flash heating [e.g., 42 ], melt lubrication [e.g., 43 ], thermal pressurization [e.g., 44 ], acoustic fluidization 45 , elastohydrodynamic lubrication 33,46 , off-fault deformation incurred during slow 47 or fast rupture [48][49][50][51][52] are among the various processes that have been proposed to explain energy dissipation during earthquakes.
These dissipative processes may influence the mechanics of earthquakes in different ways 8,53 .However, building on previous literature [53][54][55] , among others], we propose a conceptual picture (see Fig. 1) that distinguishes between the following key processes: 1. Tip processes dissipate energy near the rupture front and therefore contribute to the earthquake fracture energy that is equivalent to Γ tot utilized in LEFM.The rupture tip region is characterized by intense slip accelerations ( > 100m/s 2 ) and high slip velocities ( > 1m/s), but because it is highly transient, the associated slip is typically a relatively small fraction of the total coseismic slip.2. Tail processes occur behind the rupture tip where slip acceleration is much lower.However, slip velocities may remain relatively high ( ~1m/s) in the wake of the rupture tip, especially for crack-like ruptures that are characterized by wide-spread slip compared to pulse-like ruptures during which only small parts of the fault slide at a given time.Therefore, the slip accumulated in the rupture tail can be large if the rupture continues long enough and may result from complex secondary ruptures (as observed in labquakes [e.g., 22,56 ]).

BOX 2
The breakdown work The breakdown work 67,96 is the measurable portion of the frictional work density, which, when integrated on the fault surface, gives an estimate of the irreversible part of the total strain energy change that does not go into radiated energy 6,8 .Ref. 67 defined the breakdown work W b as the excess of work over the minimum shear stress achieved during slip τ min : where _ δðtÞ is the slip rate, τ(t) the shear stress, and t b is the time at which τ min and the critical slip distance d c are reached 97 .We note that a sharp boundary between the tip and tail processes likely does not exist and that some dissipative processes are affected by both tip and tail [e.g., 57 ].The tip and tail terminology is not limited to interface processes (i.e., fault processes) but may also include dissipation in the bulk material (i.e., host rocks), consistent with the initial formulation of LEFM.The "tip" and "tail" terminology becomes specifically useful when discussing how different dissipative processes may affect different aspects of earthquake rupture propagation and arrest.For example, flash heating may be a weakening mechanism that is active as a tip process.In contrast, thermal pressurization will likely only occur as a tail process after sufficient slip has occurred [e.g., 58 ].
Laboratory experiments offer valuable insights into this problem, albeit with notable limitations.Most small-scale experiments, for instance, cannot achieve slip acceleration that is fast enough to fully emulate the loading conditions of a dynamic rupture front (i.e., tip processes), while large-scale rupture experiments do not exhibit enough slip for tail processes to become dominant.It is far from trivial to set up laboratory experiments capable of reproducing the slip values, velocities, and accelerations under realistic loading conditions representative of a propagating earthquake rupture.These challenges can be addressed with numerical simulations, which allow us to assess tip processes under non-trivial friction conditions or assess contributions by other dissipative mechanisms under limited conditions [e.g., 48,59 ].Aside from limitations of laboratory and theoretical work, we also note that the theoretical definition of where the tip ends and the tail starts is not well defined and is likely model-dependent 57 and hence requires further investigation.
Finally, how do the tip and tail processes influence the mechanics of earthquakes?This is one of the key open questions in earthquake physics and is at the center of this Perspective.Without tail processes, LEFM shows, as outlined in Sec I, that the dissipative energy in the rupture tip, together with the energy release rate, controls rupture speed and arrest.Whether this is equally true for systems with significant tail processes and a "fuzzy" transition from tip to tail processes remains to be shown.Recent results 54,60,61 suggest that the tip processes dictate the rupture growth even in the presence of nonnegligible tail dissipation.This suggests that the fracture energy is generally well defined by the tip processes and that in the considered cases 54,60,61 , the tail processes do not significantly affect the energy release rate.However, the tail may become important when the earthquake propagates slowlypossibly during a slow arrestpropagates as multiple fronts, as a self-healing slip pulse 62 , or in multiple sliding episodes across rough fault surfaces [e.g., 63 ].It may also become important when considering how earthquakes prepare the fault for subsequent events.Here again, numerical simulations provide a tool to systematically study the link between tip and tail processes and the mechanics of earthquake ruptures.
In summary, defining tip and tail processes and how various dissipative processes contribute to them is crucial to a better understanding of how earthquakes propagate, arrest, and prepare the fault for subsequent events.The size of the yielding zone near the tip of a propagating rupture front is also an open question, which affects the values of inferred fracture energy, as we will discuss in the next section.Laboratory experiments and numerical simulations, best in synergistic combination, may provide crucial insight into these processes but require further development.Finally, the proposed framework needs to be applied to natural earthquakes but this requires a precise understanding of how field observations are linked to these tip and tail processes, which is another important open question.

Observations of energy dissipation in natural, laboratory, and simulated earthquakes
Given these theoretical considerations, we explore and compare observations of energy dissipation in both labquakes and tectonic earthquakes, with a focus on what this reveals about applying LEFM in earthquake physics.Estimates of Γ tot and the total energy dissipation vary widely for labquakes, tectonic earthquakes, and numerical models of earthquake rupture.For instance, estimates by Abercrombie and Rice 64 , which are calculated from a combination of seismically derived parameters (see Box 3), have suggested that the average energy dissipation in natural earthquakes ranges across multiple orders of magnitude from 10 2 to 10 7 J/m 2 .A compilation of seismologically inferred energy dissipation extends this range to 10 −2 − 10 8 J/m 2 8,44 , but these estimates are all highly model-dependent and subject to large, and potentially systematic, uncertainties 65 .Pseudo-dynamic earthquake modeling that infers shear stress evolution based on slip history [e.g., 66 ] yields estimates of MJ/m 2 for magnitudes M larger than 5 67 but commonly shows large variability along the fault plane [67][68][69][70] .Near-fault observations are used to infer constitutive parameters, such as the critical slip distance 71 , which would imply large values of breakdown work or average energy dissipation.Other approaches based on dynamic models [72][73][74][75][76] that simulate spontaneous dynamic rupture and employ a frictional constitutive law yield estimates of total energy dissipation that range from 1 to 10MJ/m 2 .
Similar estimates of energy dissipation from acoustic emission spectra of labquakes yield values of 10 −6 − 1J/m 2 77 .In more traditional shear fracture experiments of intact rock, however, the fracture energy has been measured in the range 10 3 − 10 4 J/m 2 at the 100MPa pressures expected in much of the seismogenic crust [78][79][80][81][82] .This is the energy required to form a fault in intact rock and is supposed to be orders of magnitude higher than the fracture energy required to rupture an existing tectonic fault.Energy dissipation estimated from friction experiments on rough surfaces that are flat at long wavelength yield estimates of 10 −1 − 10 1 J/m 2 22-25 .Other experiments where surfaces experience large slip and concentrated shear heating show continued weakening of the interface up to 1m of slip 83,84 , which has been interpreted as energy dissipation up to 1MJ/m 2 .Still, other data come from mining-induced earthquakes where the faults intersect working faces.The fraction of energy dissipated during an M2.1 mining event was estimated to be about 1 − 9% of the total energy released 85,86 .More recently ref. 87 arrived at a similar fraction of < 1% of the total energy based on fault gouge analysis, which corresponds to an average dissipated energy per event of roughly 0.5MJ/m 2 .
The interpretation of the broad range of values inferred for energy dissipation (per unit area) requires careful analysis, as some may correspond to (tip-localized) fracture energy, and others correspond to energy dissipation within a broader region that is not localized and should, therefore, not be included in fracture energy.For example, quasi-static laboratory experiments, including rotary shear and other experiments with modest slip acceleration (5 m/s 2 ), produce conditions appropriate for tail processes, not tip processes.Thus, under these conditions, the inferred dissipation does not correspond to the fracture energy and W b ≠ Γ tot .Spectral seismological estimates [e.g., compilations of refs.44,64] use information from the entire earthquake rupture area, rather than from just the propagating rupture front, and can therefore also include 'tail' dissipation mechanisms.
It has also been shown 88 that the resolution of shear stress evolution inferred from pseudodynamic modeling [e.g., 67 ] is limited by the bandwidth of the input data, either from the band limits imposed on ground motions or from smoothing operators used to regularize the kinematic finite-fault inversions.However, integral quantities such as fracture energy or breakdown work are less affected by bandwidth limitations, and thus, they are considered more reliable measures.Scale dependence of the physical processes governing dynamic weakening might also explain both the broad range of values and the scaling of energy dissipation with slip [see ref.Several recent works have discussed the increase of fracture energy and breakdown work with total slip or earthquake size.A coherent interpretation of this scaling is still lacking, which implies that some caution is warranted.However, we note that this scaling is observed in seismological estimates from both natural earthquakes 64 and labquakes 77 as well as in numerical modeling studies, as shown in ref. 8. On the other hand, Ke et al. 89 proposed a numerical model to suggest that scaling of seismologically estimated energy dissipation with slip can result from stress overshoot (i.e., lower final stress than dynamic friction), rather than a true increase of fracture energy with rupture size and fault slip.While overshoot cannot be used to explain the scaling reported from pseudo-dynamic modeling 67 , that work raises important questions about how to reconcile the huge range of fracture energy measurements and the earthquake energy budget.This highlights the importance of further research to adjust current practices used to perform dynamic modeling of natural earthquakes and ground motion predictions.

Conclusion & outlook
This perspective discusses the earthquake energy budget and the potential of Linear Elastic Fracture Mechanics (LEFM) theory to describe earthquake rupture in the laboratory and nature.The key condition for LEFM applicability is that energy dissipation during rupture propagation (except for frictional heat) must be localized in a small-scale zone at the rupture tipa condition that is commonly satisfied in large-scale laboratory experiments but may not be fully met for tectonic faults with all of their complexity.This raises important questions about how to consistently and correctly describe the energy dissipation of natural earthquakes.We suggest distinguishing between tip processes that account for localized dissipation and tail processes that occur further away from the rupture tip.In this framework, tip processes govern earthquake rupture extension and propagation, while tail processes are more important for other measures of earthquake energy dissipation and the global energy budget.We also highlight the large range of measured or inferred energy dissipation from labquakes and earthquakes across many orders of magnitude and the possibility that this could result from comparing localized with non-localized dissipation.
While the proposed tip-versus-tail perspective provides a useful approach to discussing energy dissipation in earthquakes, it also opens important scientific questions that are to be addressed in future research.For instance, the boundary between the tip and tail, i.e., the localization of energy dissipation, is neither well defined nor known.Experiments and field observations with improved sensing are needed to measure the contributions to tip and tail energy dissipation.Here it is important to note that any physical process may contribute to dissipation in the tip and the tail concurrently, and hence, separating the contributions to each is required.It is also important to use precise and consistent terminology to avoid misinterpretation of data.Specifically, only rupture-tip energy dissipation should be termed "fracture energy"; and when there is no proof that energy dissipation is local to the rupture tip, more general terms, such as "breakdown work," should be used.
Another important open question concerns the effect of significant tail processes on earthquake propagation and arrest mechanisms.Do these processes affect the energy balance (Eq.( 1)) and, hence the rupture speed?Here, numerical simulations are particularly useful to systematically explore and isolate these effects and to update the fracture theory for the description of rupture growth in the presence of tail processes.Such simulations could also provide a tool to determine the link between fault properties (tip and tail processes) and averaged global observations as inferred from seismological data and hence support the correct interpretation of earthquake energy dissipation across scales.

BOX 3
Seismologically derived earthquake parameters Abercrombie and Rice 56 proposed a parameter G', here denoted W ' b to avoid confusion with energy release rate: where D is the average slip over the fault plane, Δσ the average stress drop, μ the rock shear modulus, E R the radiated energy, and M 0 the seismic moment.W ' b is theoretically equal to the breakdown work as long as the final stress τ E f is equal to the residual sliding strength of the fault τ r .This approach is entirely based on seismologically derived parameters and can be based on rupture averages from simple source models or derived from finite-fault modeling with spatially and temporally varying slip.The most reliable seismologically derived parameter is the seismic moment where A is the rupture area.In finite-fault models, the stress drop can be determined from the spatially varying slip, but for smaller earthquakes, it is typically determined by assuming a simple circular source model 1 : For circular ruptures, r can be estimated from the corner frequency of the spectrum of teleseismic waves (f 0 ) or the reciprocal of the pulse duration in time domain modeling, where k is a geometrical constant that varies widely for commonly used source models 1,3,89 , and C S is the shear wave speed.These seismologically determined source parameters are subject to large systematic and random uncertainties and should be interpreted and modeled with extreme caution 57 .
In conclusion, as a community, we need to synergistically combine field observations, laboratory experiments, and numerical simulations to determine the degree of rupture-tip localization of various energy dissipative processes and the effect of non-localized dissipation on rupture mechanics to build a consistent model for earthquake physics.

Fig. 1 |
Fig. 1 | Schematic of an earthquake with its rupture front and seismic waves.a Earthquake fault complexity includes fault geometry, depth variation in geological units, wear and gouge formation within the fault zone, and fault branching.The red zone near the rupture front indicates regions with high fault slip rates.Note seismic waves (blue) radiate from the fault zone.b The evolution of the fault zone shear stress (blue) and slip rate (red) are shown for the region around the rupture tip and toward the hypocenter where the fault has already slipped.Seismic waves are indicated by gray lines and the associated ground motions are illustrated by a bright blue line.Processes within the rupture tip and in the tail behind the tip, socalled tip-and-tail processes, are a major focus of this Perspective.