{"id":"14a5ab78-4465-4528-938a-df4af45739ca","arxiv_id":"1908.08607","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper proposes and analyzes an iterative scheme for elastic-gravitational rupture dynamics with rate-and-state friction, claiming existence and uniqueness for the viscous regularized problem, but the proof relies on a false estimate.","lead":"This paper develops an iterative coupling scheme for dynamic earthquake ruptures in a self-gravitating elastic planet and claims to prove its convergence by adding artificial viscosity. The proof hinges on a contraction estimate that contains a false inequality, so the central well-posedness claim is not supported as written.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The contraction proof of Theorem 4.1 relies on the false inequality C_eps |s^k| <= |s^{eps,k}| after (4.14): for one slip rate zero and the other approaching zero the regularized difference is of order t^2/(2eps) while the unregularized difference is t, so no uniform positive C_eps exists.","rationale":"The reader's identified weakest assumption is precisely where the proof breaks. The paper's setup is clear, the function spaces are well presented, and the authors honestly state that the analysis is for the viscous regularized problem. However, the false monotonicity estimate is internal to the proof, not a matter of disagreeing with scientific consensus. Physically, the regularized friction law becomes flat as slip rate tends to zero, so the dissipation it provides vanishes; the initial condition (4.3) sets slip rate to zero, making this regime unavoidable. The contraction inequality (4.5) therefore does not follow, and the existence theorem built on it is unproven. The paper may be repairable, for example by using a regularization with linear growth near zero slip rate or by adding assumptions that prevent slip rate from approaching zero, but as written the central claim fails. I agree with the reader's REJECT verdict and recommend no change.","tokens_in":27487,"tokens_out":5795,"duration_ms":52186,"concrete_test":"Compute R(t) = |sqrt(t^2+eps^2) - eps| / t for j=0, i = t e with t/eps = 10^{-1}, 10^{-2}, ..., 10^{-6}; the asserted uniform lower bound fails because R(t) -> 0 as t -> 0. This single counterexample settles the estimate after eq. (4.14). Alternatively, keep the exact Taylor remainder in (4.15) and check whether the remainder can be dominated by the negative term; the paper supplies no such bound.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After eq. (4.14) the authors assert a uniform constant C_eps > 0 with C_eps ||i|-|j|| <= |sqrt(|i|^2+eps^2) - sqrt(|j|^2+eps^2)| for all slip-rate vectors i,j. This is false: with j=0 and i = t e, t -> 0, the right side equals sqrt(t^2+eps^2) - eps ~ t^2/(2eps), while the left side is C_eps t, so the ratio tends to 0. The only negative term in the energy estimates, -C_{F,s} C_eps ||eps_s^k||^2 in (4.15), depends on this constant; without it the coupling terms from normal stress and state error have no dissipation to counteract them. The proof also replaces nonlinear differences by first-order expansions marked with approximately equal signs in (4.15) and (4.20) without estimating remainders; F and G are only Lipschitz with bounded first derivatives (Assumptions 2.3 and 2.4), which does not justify the expansions at the level needed for a contraction. Since Theorem 4.1 feeds directly into Corollary 5.1, Theorem 5.2, and the discrete-time contraction Theorem 6.1, the central well-posedness and convergence claims are not supported as written.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops an iterative coupling scheme for the elastic-gravitational wave equation coupled to a rate-and-state friction law on an interior fault surface. Two regularizations are introduced: a Kelvin–Voigt viscosity coefficient gamma and a smoothed slip-rate functional Psi_eps(v) = sqrt(|v|^2 + eps^2) - eps. The authors formulate a weak problem (Problem 3.1), propose a Gauss–Seidel-type iteration (Problem 4.1), and claim in Theorem 4.1 that the iteration map is a contraction in V1 x V2 under conditions (4.4) on T and gamma. From the contraction they derive strong convergence (Corollary 5.1) and existence of a weak solution (Theorem 5.2), and they extend the contraction argument to backward-Euler time discretization in Theorem 6.1. The abstract states that viscosity guarantees convergence of the scheme for solutions of the regularized problems in both continuous and discrete time.","tokens_in":27691,"tokens_out":4653,"duration_ms":45502,"significance":"If the central results were correct, the paper would make a substantial contribution: a rigorous well-posedness framework for dynamic rupture with rate-and-state friction in a self-gravitating, prestressed elastic body, with explicit conditions on the artificial viscosity, the time step, and the contraction rate. The paper builds in a transparent way on established energy estimates of de Hoop, Holman and Pham and on the Martins–Oden approach, and the authors state the role of each assumption clearly. The discrete-time analysis is a useful addition. However, the main contraction proof rests on a quantitative monotonicity claim that is false, and this invalidates the central existence and convergence theorems as written.","major_comments":[{"comment":"The proof of Theorem 4.1 requires a uniform positive constant C_eps with C_eps |eps_s^k| <= |eps_s^{eps,k}| for all slip-rate vectors, where eps_s^{eps,k} is the difference of the regularized slip rates. This assertion is false. Taking j = 0 and i = t e with t > 0 small, the unregularized difference is t while the regularized difference is sqrt(t^2 + eps^2) - eps ~ t^2/(2 eps) as t -> 0, so the ratio tends to 0 and no positive constant C_eps independent of t can exist. The constant C_eps enters the only negative term -C_{F,s} C_eps ||eps_s^k||^2 in (4.15), which must absorb the C*_{F,sigma} and C*_{F,psi} coupling terms and the C*_{G,s} term in (4.19). Without this term the dissipation balance (4.21) and hence the contraction (4.5) do not follow. The same constant is reused in the discrete-time proof through (6.12), so Theorem 6.1 is affected equally. Since Theorem 4.1 feeds directly into Corollary 5.1, Theorem 5.2, and Theorem 6.1, the central well-posedness and convergence claims are not supported as written.","section":"§4.2, after Eq. (4.14)"},{"comment":"The proof replaces nonlinear differences by first-order Taylor expansions marked with approximate equality signs without estimating the remainder terms. Assumptions 2.3 and 2.4 only give Lipschitz continuity and bounded first derivatives of F and G; they do not justify pointwise equalities up to controlled remainders in the required L^2(Σ_f) norms. For a contraction argument with a strict rate lambda < 1, every remainder must be explicitly bounded and absorbed in the estimates. As written, these expansions constitute a second, independent gap in the proof of Theorem 4.1.","section":"§4.2, Eqs. (4.15) and (4.20)"}],"minor_comments":[{"comment":"There is an indexing inconsistency in the normal-stress error term: the first equality writes bar_sigma(eps^k_u), while the bound on the right is expressed in terms of eps^{k-1}_u. Since the iteration uses u^{k-1} in the normal stress, the error should be eps^{k-1}_u.","section":"§4.2, Eq. (4.16)"},{"comment":"There are several typographical errors: 'vatious' in the introduction, 'simultions' in Remark 4.3, 'we proof the existence' at the start of Section 5, and 'Cauchy-Schwartz' in the proof of Lemma 3.3 should be 'Cauchy-Schwarz'.","section":"Throughout"},{"comment":"The description of subdividing the time interval and restarting iterations on each segment would benefit from a precise statement of how the initial guess for each new time segment is obtained and how the contraction constants propagate across segments.","section":"Remark 4.2"}],"recommendation":"reject","confidential_remarks":"I concur with the skeptical assessment: the counterexample to the claimed uniform lower bound after (4.14) is decisive for the proof as written, and the unquantified Taylor expansions in (4.15) and (4.20) compound the problem. The paper may be salvageable with a different regularization possessing genuine strong monotonicity or with a substantially different argument, but that would go beyond a routine revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper takes on a real gap: no well-posedness analysis exists for rate-and-state friction coupled to the self-gravitating elastic system, and the iterative coupling / multi-rate time stepping perspective is genuinely useful. The variational setup, the function spaces, and the viscous regularization are all clearly presented, and the authors are honest that they only treat the regularized problem. That is the good part.\n\nBut the central contraction theorem, Theorem 4.1, has a load-bearing flaw. After (4.14) the authors assert that there is a uniform constant C_epsilon > 0 with C_epsilon ||i|-|j|| <= |sqrt(|i|^2+eps^2)-sqrt(|j|^2+eps^2)|. That is false: take j=0 and i=t e with t going to 0; the right side behaves like t^2/(2eps), which is much smaller than t, so no positive constant can work for all t. The paper's own bound (3.11) shows the regularized slip rate is smaller than the actual slip rate, so the difference can be far smaller. This invalidates the negative dissipative term -C_epsilon C_{F,s}||eps_s^k||^2 that controls the state-error coupling. The proof also uses unquantified first-order expansions, denoted by an approximate equality sign, in (4.15) and (4.20), without bounding remainders. Those are sloppy, but the false C_epsilon bound is fatal on its own.\n\nThe viscosity term could possibly provide the missing control of the slip-rate trace, since it gives -gamma||dot-eps_u||^2_H, and the trace inequality bounds the slip-rate difference. But the proof as written does not exploit that; it relies entirely on the false C_epsilon. So the contraction claim, and therefore Theorems 5.2 and 6.1, are not supported by the presented mathematics.\n\nThis paper is for mathematical seismologists and numerical analysts working on rupture simulation. It deserves a serious referee, because the gap is real and the approach is promising, but as it stands it should not be accepted. The fix may be possible—replace the C_epsilon argument with a direct use of the viscosity to control the slip-rate trace, and make the Taylor expansions rigorous—but that is substantial work. Send it to peer review, but expect major revision.","headline":"A serious and needed analysis, but the contraction proof rests on a false estimate, so the main theorems don't yet hold as written.","tokens_in":28329,"tokens_out":6069,"would_cite":false,"duration_ms":59314,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q74","74H20","74H25","86A15"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that adding a small viscosity term to the elastic-gravitational wave system makes the rate-and-state friction rupture problem well-posed, and that the natural iterative coupling scheme converges as a contraction.","keywords":["dynamic rupture","rate and state friction","seismic wave generation","self-gravitating planet","iterative coupling","viscous regularization","well-posedness","weak solution"],"falsifier":"Evaluate the ratio $R=(\\sqrt{|i|^2+\\varepsilon^2}-\\sqrt{|j|^2+\\varepsilon^2})/(|i|-|j|)$ for slip-rate vectors $i,j$ with $|i|\\ne|j|$ and seek its infimum over admissible values, for example $|i|\\to 0$ with $|j|$ vanishing; if $R$ is not bounded away from zero for fixed $\\varepsilon$, the uniform constant $C_\\varepsilon$ asserted in the estimate after (4.14) does not exist and the proof of Theorem 4.1 needs a different argument.","tokens_in":27162,"feed_emoji":"🌍","tokens_out":8468,"duration_ms":80295,"temperature":0.7,"pith_summary":"The paper addresses a gap in the mathematical theory of earthquakes: the system that couples seismic waves in a self-gravitating, prestressed elastic planet to a rate-and-state friction law (fault strength depending on slip rate and an evolving state variable) had no rigorous weak-formulation well-posedness result in three dimensions. It proposes an iterative scheme that alternates between solving the elastic-gravitational wave equation and updating the friction state variable, and proves that the iteration is a contraction in natural energy norms once a small Kelvin-Voigt viscosity is added and the slip rate is regularized. The main results are Theorem 4.1 for continuous time and Theorem 6.1 for backward-Euler time stepping, with explicit conditions linking the viscosity coefficient $\\gamma$, the time horizon, and the friction and geometry constants. If correct, the paper supplies the missing existence proof for the regularized coupled system and gives numerical analysts concrete criteria for choosing artificial viscosity in rupture simulations.","feed_headline":"With viscosity, fault-wave equations are well-posed","feed_subtitle":"A contraction proof for the coupled friction-wave system opens the way to multi-rate earthquake simulations.","key_machinery":"The load-bearing construction is the regularized slip rate $\\Psi_\\varepsilon(v)=\\sqrt{|v|^2+\\varepsilon^2}-\\varepsilon$, which smooths the non-differentiable $|v|$ term in the friction law, together with the Kelvin-Voigt viscosity term $\\gamma(\\dot u,w)_H$ added to the weak elastic-gravitational form. The iterative scheme decouples the wave equation (4.1) from the state ODE (4.2); the proof subtracts equations for two successive iterates, bounds the error using the coercivity of the prestressed elastic form $a'_3$, trace inequalities for the fault surface, and Gronwall's lemma, and chooses $\\gamma$ large enough (condition (4.4)) to absorb the friction and normal-stress coupling terms. It is this viscosity that gives the $H$-norm dissipation needed to turn the nonlinear iteration into a contraction.","core_discovery":"The central claim is that the viscosity-regularized coupled system has a unique weak solution that can be computed by a splitting iteration. In Theorem 4.1, the paper shows that, under the coefficient conditions (4.4), the map from one iterate $(u^{k-1},\\psi^{k-1})$ to the next is a contraction in $V_1 \\times V_2$, with a contraction factor $\\lambda < 1$ and constants $\\kappa_1,\\dots,\\kappa_4$ in the energy norm (4.5). Theorem 5.2 then passes to the limit and shows the fixed point solves the original weak Problem 3.1, establishing existence and uniqueness of the regularized solution. Theorem 6.1 extends the contraction to the implicitly discretized backward-Euler scheme, with the viscosity coefficient proportional to the time step. The friction law is not solved directly; it enters as a pointwise nonlinear boundary condition through the regularized functional $\\mathcal{F}_\\varepsilon$, and the viscosity term $\\gamma(\\dot u,w)_H$ is the term that makes the error estimate close.","pith_inferences":["Because the contraction constant and the coefficient conditions depend on $\\varepsilon$ through the constant $C_\\varepsilon$ (which tends to 1 as $\\varepsilon \\to 0$), the proof does not automatically extend to the unregularized friction law; a separate argument would be needed to pass $\\varepsilon \\to 0$ uniformly.","The conditions suggest that for planar faults with weak normal-stress coupling, numerical dissipation alone may satisfy the viscosity threshold, while nonplanar faults or materials with strong impedance contrast require explicit added viscosity; this could be tested systematically by measuring iteration contraction rates in a numerical rupture code.","The same splitting strategy, with the state ODE on the fault and the wave field in the volume, could be adapted to rupture problems with thermal pressurization or poroelastic stress perturbations, since those enter only through the time-regular forcing $T_\\delta$.","If the comparability of regularized and unregularized slip-rate differences fails near zero slip, the physically relevant fixed point may still exist, but its rate of convergence and the required viscosity could differ from the theorem's predictions; a numerical check on the empirical contraction factor as slip rates approach zero would reveal whether this happens in practice."],"forward_implications":["Earthquake rupture simulations on self-gravitating earth models can use multi-rate time stepping, with fine steps for the friction state ODE and coarser steps for seismic waves, while the iteration between them converges linearly.","The explicit conditions (4.4) turn the choice of artificial viscosity from a numerical tuning parameter into a criterion: $\\gamma$ must dominate products of the friction sensitivity constants and the fault-surface trace constants, with larger values needed for less regular fault geometry.","The same iterative-coupling analysis applies to the alternative Lagrangian formulation of the prestressed wave equation, so the well-posedness result is not tied to one version of the stress tensor.","With backward-Euler time discretization, the contraction survives if the viscosity coefficient scales with the time step, and the required ratio $\\gamma/\\delta t$ is again controlled by friction sensitivity and fault geometry.","The state-variable equation converges in $L^\\infty([0,T];L^2(\\Sigma_f))$ and the displacement, velocity, and gradient converge strongly in their natural spaces, giving a convergence certificate for the computed rupture process."],"supporting_citations":[{"why":"Supplies the weak formulation, coercivity, and well-posedness of the elastic-gravitational system that the rupture problem couples to.","marker":"[13]"},{"why":"Provides the iterative-coupling existence method that the proof of Theorem 5.2 follows.","marker":"[29]"},{"why":"Supplies the slip-rate regularization $\\Psi_\\varepsilon$ used to define the friction functional.","marker":"[22]"},{"why":"Formulates the rate-and-state friction laws and the Lipschitz and positivity assumptions used in (2.20) and (2.24).","marker":"[38]"},{"why":"Gives the equation of motion in a prestressed self-gravitating earth that Section 2 starts from.","marker":"[6]"},{"why":"Companion implementation of the multi-rate iterative scheme that motivates the analysis and illustrates the numerical setting.","marker":"[52]"},{"why":"Documents numerical shocks in rupture simulations, motivating the artificial-viscosity regularization.","marker":"[12]"}],"fun_headline_variants":["Viscosity keeps fault-wave equations well-posed","Contraction proof for coupled friction-wave systems","Splitting iteration converges with viscosity added","How viscosity tames earthquake rupture equations","Well-posed rupture dynamics via viscosity regularization"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the regularized slip-rate difference $\\sqrt{|i|^2+\\varepsilon^2}-\\sqrt{|j|^2+\\varepsilon^2}$ stays above a uniform positive constant times the unregularized difference $\\bigl||i|-|j|\\bigr|$ for all iterates; the entire contraction estimate collapses if that ratio can be arbitrarily small.","fun_headline_variants_meta":{"raw":{"variants":["Viscosity keeps fault-wave equations well-posed","Contraction proof for coupled friction-wave systems","Splitting iteration converges with viscosity added","How viscosity tames earthquake rupture equations","Well-posed rupture dynamics via viscosity regularization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000575,"raw_usage":{"total_tokens":2671,"prompt_tokens":860,"completion_tokens":1811,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":476,"completion_tokens_details":{"reasoning_tokens":1744}},"tokens_in":476,"tokens_out":1811,"duration_ms":12997,"temperature":1.0,"reasoning_tokens":1744,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:35:37.374295+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the ratio $R=(\\sqrt{|i|^2+\\varepsilon^2}-\\sqrt{|j|^2+\\varepsilon^2})/(|i|-|j|)$ for slip-rate vectors $i,j$ with $|i|\\ne|j|$ and seek its infimum over admissible values, for example $|i|\\to 0$ with $|j|$ vanishing; if $R$ is not bounded away from zero for fixed $\\varepsilon$, the uniform constant $C_\\varepsilon$ asserted in the estimate after (4.14) does not exist and the proof of Theorem 4.1 needs a different argument.","supporting_citations":[{"cited_title":"On the system of elastic-gravitational equations describing the oscillations of the earth","cited_arxiv_id":"1511.03200","evidence_quote":"Supplies the weak formulation, coercivity, and well-posedness of the elastic-gravitational system that the rupture problem couples to."},{"cited_title":"3, 407–428","cited_arxiv_id":null,"evidence_quote":"Provides the iterative-coupling existence method that the proof of Theorem 5.2 follows."},{"cited_title":"3, 375–390","cited_arxiv_id":null,"evidence_quote":"Supplies the slip-rate regularization $\\Psi_\\varepsilon$ used to define the friction functional."},{"cited_title":"9, 1865–1898","cited_arxiv_id":null,"evidence_quote":"Formulates the rate-and-state friction laws and the Lipschitz and positivity assumptions used in (2.20) and (2.24)."},{"cited_title":"Variational formulation of the earth's elastic-gravitational deformations under low regularity conditions","cited_arxiv_id":"1702.04741","evidence_quote":"Gives the equation of motion in a prestressed self-gravitating earth that Section 2 starts from."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Companion implementation of the multi-rate iterative scheme that motivates the analysis and illustrates the numerical setting."},{"cited_title":"B12, 1–23","cited_arxiv_id":null,"evidence_quote":"Documents numerical shocks in rupture simulations, motivating the artificial-viscosity regularization."}],"review_version":1}