On the criticality of return flows in viscous accretionary wedges and its implications for deep-crustal exhumation in subduction zones

Patsa, Ayan; Mandal, Nibir

doi:10.5194/egusphere-2025-2909

Preprints

https://doi.org/10.5194/egusphere-2025-2909

Preprints

12 Aug 2025

| 12 Aug 2025

On the criticality of return flows in viscous accretionary wedges and its implications for deep-crustal exhumation in subduction zones

Ayan Patsa and Nibir Mandal

Abstract. In subduction zones, the accretionary wedges play a vital role in mediating the burial processes of incoming oceanic sediments and eventually their return pathways to the surface. A direction of the previous tectonic models invoked the standard corner flow theory, assuming a slab-parallel shear and a rigid, fixed overriding plate, to elucidate the crustal recycling processes in tectonic wedges. To deal with more complex subduction-collisional settings, where they have deformable overriding plates, and associate a horizontal slab migration (advance or rollback) component during subduction, we develop a generalized corner flow model to revisit the problem of return flow mechanics, providing a criticality analysis of the return flows as a function of the geometric, kinematic, and rheological conditions in accretionary wedges. A new set of analytical solutions is presented to evaluate the limiting conditions in which a wedge can set in significant return flows, leading to focused exhumation of the deep-crustal materials. The theoretical results suggest that, for moderate wedge-taper angles (~30^o), the viscosity ratios (µ_r) between the overriding plate and the wedge ≥ ~10³ provide favourable tectonic settings for the return flow kinematics in accretionary wedges. Decrease in µ_r, or addition of slab roll back weakens the return flows, whereas slab advance greatly strengthens the return flows. The analytical solutions are also utilized to demonstrate reversals in the shear-sense patterns across the wedge. We expand this study by reproducing some of the theoretical flow patterns in laboratory experiments. It is shown from the theoretical model that the total pressure in the accretionary wedge dynamics becomes close to the lithostatic value when the rheological setting has low-viscosity (10¹⁹ Pa s) wedge materials.

Received: 18 Jun 2025 – Discussion started: 12 Aug 2025

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.

Download & links

Preprint (PDF, 69556 KB)

Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
Preprint (69556 KB)

Supplement (2538 KB)

Download & links

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Journal article(s) based on this preprint

25 Mar 2026

On the criticality of return flows in viscous accretionary wedges and its implications for deep-crustal exhumation in subduction zones

Ayan Patsa and Nibir Mandal

Solid Earth, 17, 573–600, https://doi.org/10.5194/se-17-573-2026,https://doi.org/10.5194/se-17-573-2026, 2026

Short summary

Ayan Patsa and Nibir Mandal

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2025-2909', Anonymous Referee #1, 09 Sep 2025

Patsa & Mandal present a paper to study return flow in accretionary wedges, a topic that has been widely studied in the recent decades and is still relevant today. This process has been modelled with theoretical solutions, numerical and analog models to explain field-based observations on pressure temperature conditions on metamorphic rocks along accretionary prisms. The authors provide a more generalized theoretical solution to include a non-parallel component on the slab to reproduce slab advance/rollback. In addition, they conducted analog models to enhance these results. Finally, an overview with natural observations is made.
The integration of analytical solutions with analog modelling represents a novel and valuable approach, with the potential to be further strengthened by incorporating existing numerical results. However, the current manuscript does not fully emphasize this novelty: the introductory section underplays the contribution, and the analytical and analog results are presented somewhat independently, without sufficient cross-comparison.
To improve the manuscript, I recommend (i) revising the introduction to more clearly articulate the novelty and significance of the combined approach, (ii) providing greater coverage and integration of the analog modelling, and (iii) addressing several technical issues within the analytical solution. If these major concerns are resolved, the manuscript would meet the standards for publication in Solid Earth and I look forward to see the revised version.
Major Comment 1
The use of a non-dimensional parameter to quantify the strength or weakness of return flow is both effective and intuitive. This approach allows the balance of incoming and outgoing material along the top boundary to be measured, corresponding to sediment influx and exhumed units in an accretionary wedge. In the simplest case—rigid walls with parallel subduction—the value must equal 1, as dictated by wedge geometry and mass conservation.
However, when a non-parallel component is introduced to the slab, the theoretical models also incorporate material flux along the slab boundary. This introduces bias in the calculated rates of burial and exhumation, since part of the return flow ratio (FR) is influenced by this artificial slab-sourced material. In the case of slab advance, FR appears anomalously high due to two factors: (1) the top-boundary influx of sediments, which the authors correctly identify, and (2) the additional, unaddressed influx of material along the slab boundary.
For the rollback scenario, the formulation permits material to exit through the slab boundary, which is physically unrealistic. With small rollback velocities, wedge geometry would still force material to return to the surface; however, in the present formulation, this is instead channelled out through the slab, leading to FR = 0. Furthermore, at high rollback velocities, the assumption of a downward-tapered wedge may break down, opening the system to the upper mantle. To remain consistent with the wedge geometry assumptions, I suggest restricting the analysis to small perturbations of the non-parallel component.
Finally, the comparison of slab advance, rollback, and normal subduction must be conducted under consistent assumptions regarding material influx. One possible way forward would be to fix the total incoming volume and instead vary return velocities, though I acknowledge that implementing this within the analytical framework may be non-trivial.
Major Comment 2:
I strongly recommend undertaking a systematic comparison between the analytical solution and the analog modelling, as this would greatly strengthen the validation of the theoretical framework. At present, the two sets of results are described independently, which makes the paper look unbalanced. Establishing a one-to-one correspondence between the analytical predictions and the analog experiments for identical setups would provide a more rigorous test of the model and highlight the novelty of the combined approach.
That said, some modifications and additional discussion will be required to enable such a comparison. As noted in lines 415–416, the analog experiments allow material to exit through the bottom boundary, a feature not incorporated into the analytical solution. This discrepancy must be explicitly acknowledged and its implications discussed, since it directly affects the comparability of the results. One option is to consider modified boundary conditions in the analytical framework, or alternatively to constrain the analog results so that they are evaluated under conditions more consistent with the theoretical assumptions.
In addition, I recommend expanding the analog modelling section. At present, the manuscript dedicates considerably more space to the analytical solution, leaving the analog results underdeveloped. A more balanced treatment would not only give greater weight to the experiments but also allow for meaningful side-by-side comparisons. Importantly, you could draw inspiration from the approach of Moulas et al. (2021), who validated their analytical solution against numerical models with a similar setup. Extending your study to include a three-way comparison of analytical, analog, and numerical would significantly increase the robustness and originality of the manuscript.
To facilitate these improvements, some reorganization of the manuscript structure is advised. For example, moving Section 4 earlier in the text, immediately after the presentation of the analytical results would allow for more direct comparisons between the different methods. This restructuring would make the narrative more cohesive and highlight the integrative character of the study, which is currently one of its main strengths but not fully emphasized.
Line to Line comments:
Line 20: Delete “the” or simply state “facilitate subduction.”
Line 21: References are missing for geophysical observations; see Abers (2005).
Lines 25–28: I recommend mentioning the P–T–t path for consistency with the rest of the sentence, and introducing the concept of recycling here.
Line 33: Retain only geochronological, since this sentence refers exclusively to exhumation rates and not geochemical constraints.
Line 39: Verify the reference “?, for review.”
Line 45: A reference is required. The corner flow model also accounts for both prograde and retrograde metamorphism; please mention this. Additionally, note that the model supports the possibility of sediments reaching mantle depths where partial melting may occur.
Lines 48–49: Add one or two sentences on the thermal regime of subduction zones, as this strongly influences eclogite formation.
Lines 56–57: Clarify the rationale of this sentence, or consider removing it.
Line 59: Quantify exhumation rates, providing values from numerical models and natural estimates. Restrict the discussion to the specific tectonic setting under study (wedge geometry)
Line 61: At present, the introduction does not clearly define the scientific gap. While the questions posed are valid, they appear abruptly. The gap would be clearer if you outlined: (i) the discrepancies between modelled and observed exhumation rates, (ii) the wide variability in return-flow models, and (iii) the influence of parameters such as channel width, rheology, and boundary conditions. Emphasize the mismatch between observations and models.
Lines 62–63: Clarify whether questions (1) and (2) are essentially identical.
Line 64: Revise to “theoretical and analog study.”
Line 65: Specify the depth range of both the theoretical and analog models. Although this is mentioned later, it should also appear here. State the main assumptions explicitly—for example, that the accretionary wedge is closed and material cannot enter the mantle—since this represents a special-case scenario.
Line 70: Integrate content from later lines: prior work has already tested different boundary conditions in numerical models (e.g., Gerya et al., 2002). The novelty here lies in the analytical treatment of non-parallel slab boundary conditions, which allows replication of slab rollback and advance. When combined with analog modelling, this provides a unique contribution.
Line 76: Clarify the phrase “oblique to the slab.” If it refers to the trench, rephrase as “a non-parallel component of slab velocity.”
Lines 86–88: This is the first mention of analog experiments, which are central to the manuscript. Introduce them earlier in the introduction and highlight the novelty of combining analytical solutions with analog modelling.
Line 96: Since some material may be dragged down, this is an important limitation—discuss explicitly.
Line 107: Even if non-linearity is not first-order, note that complexity may also arise from contrasting lithologies within the accretionary wedge.
Line 117: Clarify whether this component is oblique to the trench or simply non-parallel.
Lines 120–121: If the trench-oblique component is merely a reduction of subduction velocity, avoid presenting it as a trench-oblique term, since this implies a 3D model.
Lines 127–128: Replace “trench-perpendicular vertical plane” with “non-parallel component of slab velocity.”
Line 130: In Fig. 2 it seems only the oceanic plate’s fixed wall is constrained. Confirm whether the top boundary is also fixed.
Line 181: This derivation follows Moulas et al. (2021). Add a phrase such as “Following the approach of Moulas et al. (2021)...” You may shorten this section and direct readers to that reference until the non-parallel extension is introduced.
Line 214: Same as line 181.
Lines 229–230: Indicate that these models assume either a rigid overriding plate or very strong subduction channels, which generate extremely high overpressures.
Line 245: Clarify whether this is the discretization used to evaluate equations. If so, specify resolution and grid type.
Lines 261–263: Add references or case studies linking models to natural observations.
Line 267: In this simple case, the flow ratio (FR) should equal 1 due to mass conservation and wedge geometry. Clarify whether this depends on discretization.
Line 271: Is this because you are adding material through the slab or because of the "squeezing" of the wedge?
Line 286: Is this difference with the rigid case because of the viscosity ratio only? what if mu_r is even higher (i.e., 10^5 or 10^7), do you reach the rigid wall solution? Also, see later comments to define high/low FR.
Lines 290–291: Figures 4–5 show material still returning to the surface but further from the trench. Confirm whether FR is calculated only for the accretionary wedge (if so, specify at line 245).
Lines 292–294: Consider adding a figure similar to FR vs. obliquity (Fig. 7) to illustrate this result, and extend the same approach to other variables.
Line 323: Since deformation is not described, either remove the vorticity figures or move them to supplementary material with an explanation.
Line 344 and Fig. 7: FR = 0.5 is presented as a threshold between significant and negligible return flow. Explain how this value was determined, or move section 3.2 to the discussion. Kerswell et al. (2023) may provide guidance.
Line 351: Revise “wedge” to “downward-tapered wedge.”
Lines 356–357: Move this sentence to the discussion section.
Lines 392–395: Provide scaling for analog experiments. Do they correspond with analytical models or plate tectonic velocities? Indicate scaling parameters (e.g., Schellart & Strak, 2016) to demonstrate consistency with natural systems and analytical calculations.
Lines 411, 413: Replace “30%” with “0.3U” and “one-sixth” with “U/6” to align with line 415 (“0.6U”).
Line 413: Clarify whether oblique shortening with slab advance and oblique extension correspond to slab rollback. If so, use consistent terminology.
Line 419: Remove “grossly.” Add a one-to-one comparison with the analytical solution.
Lines 422–423: Gravity influences the analog model if the bottom boundary is open. Add a brief discussion of this effect.
Line 425: If slab rollback or oblique extension precludes return flow, clarify how comparisons were made, since the analog and analytical models differ in bottom boundary conditions.
Line 429: Change “crusts” to “crust.”
Lines 449–450: Clarify how comparisons between theoretical models and natural examples are made. Estimate taper angles and subduction dynamics (advance vs. rollback) for each exhumation case, then compare calculated FR values with reviews (e.g., Agard et al., 2009) or case studies (e.g., Franciscan Complex; Ring, 2008).
Line 453: Are there documented cases of absent accretionary wedges in modern rollback settings?
Line 454: Revise: note that HP units are present in Chile (Willner, 2005), with localized pressures of 2–2.5 GPa (González-Jiménez et al., 2017).
Lines 503–504: Mention the thermal regime, as it controls the brittle–ductile transition depth and influences viscosity.
Line 510: Replace “subduction” with “subduction zones.”
Line 519: Expand to “Multiple structural fabrics and fluid-assisted deformation (e.g., Muñoz-Montecinos & Behr, 2023).”
Lines 530–545: This section would benefit greatly from plotting the geological cases in Figure 8, enabling direct comparison between tectonic settings, model outputs, and natural data. Analog model results could also be added for completeness.
Line 590: Include rock strength in the brittle regime. In favorable conditions, rocks can sustain tens to hundreds of MPa before failure (Platt, 2019).
Line 625: Remove mu_r here, as it denotes viscosity ratio.
References in this review:
Abers (2005): https://doi.org/10.1016/j.pepi.2004.10.002

Agard et al. (2009): https://doi.org/10.1016/j.earscirev.2008.11.002

Gerya et al. (2002): https://doi.org/10.1029/2002TC001406

González-Jiménez et al. (2017): https://doi.org/10.1127/ejm/2017/0029-2668

Kerswell et al. (2023): https://doi.org/10.1029/2022GC010834

Moulas et al. (2021): https://doi.org/10.1093/gji/ggab246

Muñoz-Montecinos & Behr (2023): https://doi.org/10.1029/2023GL104244

Platt (2019): https://doi.org/10.5194/se-10-357-2019

Ring (2008): https://doi.org/10.1130/2008.2445

Schellart & Strak (2016): https://doi.org/10.1016/j.jog.2016.03.009

Willner (2005): https://doi.org/10.1093/petrology/egi035

Citation: https://doi.org/10.5194/egusphere-2025-2909-RC1
- AC1: 'Reply on RC1', Nibir Mandal, 17 Oct 2025
  
  We thank the reviewer for all the insightful comments and suggestions. In the revised version, we carefully address the issues raised by the reviewer and incorporate all the suggestions. Please find the response letter in the attached file.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2909-AC1
RC2:
'Comment on egusphere-2025-2909', Stefan Markus Schmalholz, 11 Sep 2025

In the manuscript “On the criticality of return flows …”, the authors investigate scenarios of rock burial and exhumation during subduction and orogeny. They present analytical solutions for corner flow involving two adjacent deformable wedges: one representing the overriding plate and the other the accretionary wedge (or, in some cases, the subduction channel), where burial and exhumation may occur. The presented analytical model extends the model of Moulas et al. (GJI, 2021). While Moulas et al. considered only subduction velocities parallel to the wedge base, the present study also examines non-parallel velocities. The authors systematically analyze conditions under which significant return flow arises. They further compare the analytical velocity fields with velocity fields from analogue laboratory experiments for similar configurations and show first-order agreement between experimental and analytical velocity fields. Research on burial and exhumation in subduction zones remains highly relevant, as the controlling mechanisms are still debated. Comparing analytical predictions with analogue experiments is also very relevant. However, I have major concerns about the applicability of the presented analytical solution to natural subduction and orogenic burial-exhumation cycles.
Major comments
1) A key assumption of the analytical solution is that wedge geometry remains constant throughout the burial-exhumation cycle. Both the analytical model and analogue experiments yield essentially instantaneous velocity fields. However, the boundary conditions in the analytical model imply that wedge geometry must change over time, as observed in the experiments. If the velocity at the wedge base is not parallel to its base, the wedge will be squeezed or extended. Assuming average burial and exhumation rates of 10 mm/yr, a 50 km burial followed by 50 km exhumation would last ~10 Myr. In the scenarios presented (slab advance and rollback), non-parallel velocity components of just a few mm/yr would displace the wedge base by several tens of kilometers over this timescale. Such large geometric changes would strongly alter the internal corner flow field. Moreover, the issue of corner “stability” was examined by Moulas et al. (2021), but their findings are not fully considered here. For example, Fig. 4 includes results for viscosity ratios ≤1000, yet at such values the internal wedge boundary would deform significantly during the burial-exhumation cycle, effectively destroying the corner geometry. For these reasons, I am not convinced that the presented velocity fields for slab advance and rollback can reliably be integrated to predict a complete burial-exhumation path.
2) Analytical corner flow models can be useful for certain geodynamic scenarios. However, regarding (U)HP rock exhumation in subduction zones, most 2D thermo-mechanical simulations do not generate wedges with forced corner flow when the subduction zone develops self-consistently (i.e., without a pre-imposed weak zone or wedge). In such models, (U)HP exhumation is typically driven by buoyancy or plate divergence/extension, rather than by forced return flow. Thus, the presented corner flow model may be applicable to burial and exhumation at crustal depths, but for rocks buried deeper than ~35 km, exhumation is more likely controlled by buoyancy. Another limitation is the assumption of constant linear viscosity across a wedge spanning the entire crust or even deeper. In reality, significant temperature variations produce large variations in effective rock viscosity. Also, deeper ductile regions may localize strain into shear zones, allowing the subducting plate to slide beneath the wedge without initiating a distributed corner flow. Hence, the authors should clearer discuss the range of applicability of their model. In particular, I find the reference to both an accretionary wedge and a subduction channel in Fig. 1a problematic. An application to shallow accretionary wedges seems far more realistic than to subduction channels in sub-crustal depths.
Minor comments
The introduction is a bit confusing because the authors do not use a consistent terminology. They mention accretionary wedges in places but seem to refer to subduction channels. In Fig. 1a the corner flow seems to be representative for both accretionary wedges and subduction channels. Buoyancy may be ignored in a shallow accretionary wedge but is likely important for subduction channels at sub-crustal depths.
Line 79-90: In section 2.1. the authors state that they generalize the corner-flow model of Moulas et al. (2021), which is based on coupling two corner flow solutions to consider the deformation of the overriding plate. Also, dynamic pressure fields, velocity fields and a velocity profile showing burial and return velocities have all been shown and discussed in Moulas et al. (2021). Also, Fig. 2 is very similar to figures 2 and 3 in Moulas et al. (2021). Maybe the authors could clearer state in the Introduction that their model is a modification of the model of Moulas et al. (2021). In lines 79-90 the study of Moulas et al. (2021) is not mentioned.
Line 103: Buoyancy is excluded but may play the dominant role for rock exhumation from sub-crustal depths in many orogenies.
Line 236-237: I would argue that even more studies consider the positive buoyancy of subducted rocks as a more important mechanisms for the exhumation of (U)HP rocks that have been subducted to depth larger than the average crustal thickness.
Line 270: If I understood correctly, a value of F>1 implies that more material is exhumed by return flow than is subducted. How can this be applied to natural scenarios? It means that more material is exhumed than buried. How is this in agreement with mass balance and where is the additional material that is exhumed generated? F>1 also implies that the wedge disappears after some time. Likely I missed something, but the authors should better explain the meaning, implications and applications for models with F>1.
Line 274: The lower boundary of the wedge in Fig. 3 indicates the top of the subducting plate (see also Fig. 2a). When the velocity arrows are not parallel to this wedge boundary and point towards the subducting plate, then this implies that material from the wedge flows into the subducting plate, because the lower wedge boundary does not move with the velocity. This makes no sense to me. Keeping the wedge geometry constant but applying boundary velocities that would change the wedge geometry is not consistent.
Line 315-316: The mentioned interface migration has been studied by Moulas et al. (2021) with a 2D dimensionless “regime diagram“ distinguishing a “stable corner“ (with negligible interface migration) and an “unstable corner“ (with significant interface migration).
Section 4.1: Please provide a scaling analysis for the experiment so that it is clear how the laboratory experiments can be scaled to the natural situation.
Best regards,
Stefan Schmalholz

Citation: https://doi.org/10.5194/egusphere-2025-2909-RC2
- AC2: 'Reply on RC2', Nibir Mandal, 17 Oct 2025
  
  We are grateful to the reviewer for thoroughly reviewing the manuscript, and raising critical questions. In the revised version, we carefully address these issues. The comments have greatly helped us improve the quality of our theoretical and experimental results presentations. Please find the response letter in the attached file.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2909-AC2

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2025-2909', Anonymous Referee #1, 09 Sep 2025

Patsa & Mandal present a paper to study return flow in accretionary wedges, a topic that has been widely studied in the recent decades and is still relevant today. This process has been modelled with theoretical solutions, numerical and analog models to explain field-based observations on pressure temperature conditions on metamorphic rocks along accretionary prisms. The authors provide a more generalized theoretical solution to include a non-parallel component on the slab to reproduce slab advance/rollback. In addition, they conducted analog models to enhance these results. Finally, an overview with natural observations is made.
The integration of analytical solutions with analog modelling represents a novel and valuable approach, with the potential to be further strengthened by incorporating existing numerical results. However, the current manuscript does not fully emphasize this novelty: the introductory section underplays the contribution, and the analytical and analog results are presented somewhat independently, without sufficient cross-comparison.
To improve the manuscript, I recommend (i) revising the introduction to more clearly articulate the novelty and significance of the combined approach, (ii) providing greater coverage and integration of the analog modelling, and (iii) addressing several technical issues within the analytical solution. If these major concerns are resolved, the manuscript would meet the standards for publication in Solid Earth and I look forward to see the revised version.
Major Comment 1
The use of a non-dimensional parameter to quantify the strength or weakness of return flow is both effective and intuitive. This approach allows the balance of incoming and outgoing material along the top boundary to be measured, corresponding to sediment influx and exhumed units in an accretionary wedge. In the simplest case—rigid walls with parallel subduction—the value must equal 1, as dictated by wedge geometry and mass conservation.
However, when a non-parallel component is introduced to the slab, the theoretical models also incorporate material flux along the slab boundary. This introduces bias in the calculated rates of burial and exhumation, since part of the return flow ratio (FR) is influenced by this artificial slab-sourced material. In the case of slab advance, FR appears anomalously high due to two factors: (1) the top-boundary influx of sediments, which the authors correctly identify, and (2) the additional, unaddressed influx of material along the slab boundary.
For the rollback scenario, the formulation permits material to exit through the slab boundary, which is physically unrealistic. With small rollback velocities, wedge geometry would still force material to return to the surface; however, in the present formulation, this is instead channelled out through the slab, leading to FR = 0. Furthermore, at high rollback velocities, the assumption of a downward-tapered wedge may break down, opening the system to the upper mantle. To remain consistent with the wedge geometry assumptions, I suggest restricting the analysis to small perturbations of the non-parallel component.
Finally, the comparison of slab advance, rollback, and normal subduction must be conducted under consistent assumptions regarding material influx. One possible way forward would be to fix the total incoming volume and instead vary return velocities, though I acknowledge that implementing this within the analytical framework may be non-trivial.
Major Comment 2:
I strongly recommend undertaking a systematic comparison between the analytical solution and the analog modelling, as this would greatly strengthen the validation of the theoretical framework. At present, the two sets of results are described independently, which makes the paper look unbalanced. Establishing a one-to-one correspondence between the analytical predictions and the analog experiments for identical setups would provide a more rigorous test of the model and highlight the novelty of the combined approach.
That said, some modifications and additional discussion will be required to enable such a comparison. As noted in lines 415–416, the analog experiments allow material to exit through the bottom boundary, a feature not incorporated into the analytical solution. This discrepancy must be explicitly acknowledged and its implications discussed, since it directly affects the comparability of the results. One option is to consider modified boundary conditions in the analytical framework, or alternatively to constrain the analog results so that they are evaluated under conditions more consistent with the theoretical assumptions.
In addition, I recommend expanding the analog modelling section. At present, the manuscript dedicates considerably more space to the analytical solution, leaving the analog results underdeveloped. A more balanced treatment would not only give greater weight to the experiments but also allow for meaningful side-by-side comparisons. Importantly, you could draw inspiration from the approach of Moulas et al. (2021), who validated their analytical solution against numerical models with a similar setup. Extending your study to include a three-way comparison of analytical, analog, and numerical would significantly increase the robustness and originality of the manuscript.
To facilitate these improvements, some reorganization of the manuscript structure is advised. For example, moving Section 4 earlier in the text, immediately after the presentation of the analytical results would allow for more direct comparisons between the different methods. This restructuring would make the narrative more cohesive and highlight the integrative character of the study, which is currently one of its main strengths but not fully emphasized.
Line to Line comments:
Line 20: Delete “the” or simply state “facilitate subduction.”
Line 21: References are missing for geophysical observations; see Abers (2005).
Lines 25–28: I recommend mentioning the P–T–t path for consistency with the rest of the sentence, and introducing the concept of recycling here.
Line 33: Retain only geochronological, since this sentence refers exclusively to exhumation rates and not geochemical constraints.
Line 39: Verify the reference “?, for review.”
Line 45: A reference is required. The corner flow model also accounts for both prograde and retrograde metamorphism; please mention this. Additionally, note that the model supports the possibility of sediments reaching mantle depths where partial melting may occur.
Lines 48–49: Add one or two sentences on the thermal regime of subduction zones, as this strongly influences eclogite formation.
Lines 56–57: Clarify the rationale of this sentence, or consider removing it.
Line 59: Quantify exhumation rates, providing values from numerical models and natural estimates. Restrict the discussion to the specific tectonic setting under study (wedge geometry)
Line 61: At present, the introduction does not clearly define the scientific gap. While the questions posed are valid, they appear abruptly. The gap would be clearer if you outlined: (i) the discrepancies between modelled and observed exhumation rates, (ii) the wide variability in return-flow models, and (iii) the influence of parameters such as channel width, rheology, and boundary conditions. Emphasize the mismatch between observations and models.
Lines 62–63: Clarify whether questions (1) and (2) are essentially identical.
Line 64: Revise to “theoretical and analog study.”
Line 65: Specify the depth range of both the theoretical and analog models. Although this is mentioned later, it should also appear here. State the main assumptions explicitly—for example, that the accretionary wedge is closed and material cannot enter the mantle—since this represents a special-case scenario.
Line 70: Integrate content from later lines: prior work has already tested different boundary conditions in numerical models (e.g., Gerya et al., 2002). The novelty here lies in the analytical treatment of non-parallel slab boundary conditions, which allows replication of slab rollback and advance. When combined with analog modelling, this provides a unique contribution.
Line 76: Clarify the phrase “oblique to the slab.” If it refers to the trench, rephrase as “a non-parallel component of slab velocity.”
Lines 86–88: This is the first mention of analog experiments, which are central to the manuscript. Introduce them earlier in the introduction and highlight the novelty of combining analytical solutions with analog modelling.
Line 96: Since some material may be dragged down, this is an important limitation—discuss explicitly.
Line 107: Even if non-linearity is not first-order, note that complexity may also arise from contrasting lithologies within the accretionary wedge.
Line 117: Clarify whether this component is oblique to the trench or simply non-parallel.
Lines 120–121: If the trench-oblique component is merely a reduction of subduction velocity, avoid presenting it as a trench-oblique term, since this implies a 3D model.
Lines 127–128: Replace “trench-perpendicular vertical plane” with “non-parallel component of slab velocity.”
Line 130: In Fig. 2 it seems only the oceanic plate’s fixed wall is constrained. Confirm whether the top boundary is also fixed.
Line 181: This derivation follows Moulas et al. (2021). Add a phrase such as “Following the approach of Moulas et al. (2021)...” You may shorten this section and direct readers to that reference until the non-parallel extension is introduced.
Line 214: Same as line 181.
Lines 229–230: Indicate that these models assume either a rigid overriding plate or very strong subduction channels, which generate extremely high overpressures.
Line 245: Clarify whether this is the discretization used to evaluate equations. If so, specify resolution and grid type.
Lines 261–263: Add references or case studies linking models to natural observations.
Line 267: In this simple case, the flow ratio (FR) should equal 1 due to mass conservation and wedge geometry. Clarify whether this depends on discretization.
Line 271: Is this because you are adding material through the slab or because of the "squeezing" of the wedge?
Line 286: Is this difference with the rigid case because of the viscosity ratio only? what if mu_r is even higher (i.e., 10^5 or 10^7), do you reach the rigid wall solution? Also, see later comments to define high/low FR.
Lines 290–291: Figures 4–5 show material still returning to the surface but further from the trench. Confirm whether FR is calculated only for the accretionary wedge (if so, specify at line 245).
Lines 292–294: Consider adding a figure similar to FR vs. obliquity (Fig. 7) to illustrate this result, and extend the same approach to other variables.
Line 323: Since deformation is not described, either remove the vorticity figures or move them to supplementary material with an explanation.
Line 344 and Fig. 7: FR = 0.5 is presented as a threshold between significant and negligible return flow. Explain how this value was determined, or move section 3.2 to the discussion. Kerswell et al. (2023) may provide guidance.
Line 351: Revise “wedge” to “downward-tapered wedge.”
Lines 356–357: Move this sentence to the discussion section.
Lines 392–395: Provide scaling for analog experiments. Do they correspond with analytical models or plate tectonic velocities? Indicate scaling parameters (e.g., Schellart & Strak, 2016) to demonstrate consistency with natural systems and analytical calculations.
Lines 411, 413: Replace “30%” with “0.3U” and “one-sixth” with “U/6” to align with line 415 (“0.6U”).
Line 413: Clarify whether oblique shortening with slab advance and oblique extension correspond to slab rollback. If so, use consistent terminology.
Line 419: Remove “grossly.” Add a one-to-one comparison with the analytical solution.
Lines 422–423: Gravity influences the analog model if the bottom boundary is open. Add a brief discussion of this effect.
Line 425: If slab rollback or oblique extension precludes return flow, clarify how comparisons were made, since the analog and analytical models differ in bottom boundary conditions.
Line 429: Change “crusts” to “crust.”
Lines 449–450: Clarify how comparisons between theoretical models and natural examples are made. Estimate taper angles and subduction dynamics (advance vs. rollback) for each exhumation case, then compare calculated FR values with reviews (e.g., Agard et al., 2009) or case studies (e.g., Franciscan Complex; Ring, 2008).
Line 453: Are there documented cases of absent accretionary wedges in modern rollback settings?
Line 454: Revise: note that HP units are present in Chile (Willner, 2005), with localized pressures of 2–2.5 GPa (González-Jiménez et al., 2017).
Lines 503–504: Mention the thermal regime, as it controls the brittle–ductile transition depth and influences viscosity.
Line 510: Replace “subduction” with “subduction zones.”
Line 519: Expand to “Multiple structural fabrics and fluid-assisted deformation (e.g., Muñoz-Montecinos & Behr, 2023).”
Lines 530–545: This section would benefit greatly from plotting the geological cases in Figure 8, enabling direct comparison between tectonic settings, model outputs, and natural data. Analog model results could also be added for completeness.
Line 590: Include rock strength in the brittle regime. In favorable conditions, rocks can sustain tens to hundreds of MPa before failure (Platt, 2019).
Line 625: Remove mu_r here, as it denotes viscosity ratio.
References in this review:
Abers (2005): https://doi.org/10.1016/j.pepi.2004.10.002

Agard et al. (2009): https://doi.org/10.1016/j.earscirev.2008.11.002

Gerya et al. (2002): https://doi.org/10.1029/2002TC001406

González-Jiménez et al. (2017): https://doi.org/10.1127/ejm/2017/0029-2668

Kerswell et al. (2023): https://doi.org/10.1029/2022GC010834

Moulas et al. (2021): https://doi.org/10.1093/gji/ggab246

Muñoz-Montecinos & Behr (2023): https://doi.org/10.1029/2023GL104244

Platt (2019): https://doi.org/10.5194/se-10-357-2019

Ring (2008): https://doi.org/10.1130/2008.2445

Schellart & Strak (2016): https://doi.org/10.1016/j.jog.2016.03.009

Willner (2005): https://doi.org/10.1093/petrology/egi035

Citation: https://doi.org/10.5194/egusphere-2025-2909-RC1
- AC1: 'Reply on RC1', Nibir Mandal, 17 Oct 2025
  
  We thank the reviewer for all the insightful comments and suggestions. In the revised version, we carefully address the issues raised by the reviewer and incorporate all the suggestions. Please find the response letter in the attached file.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2909-AC1
RC2:
'Comment on egusphere-2025-2909', Stefan Markus Schmalholz, 11 Sep 2025

In the manuscript “On the criticality of return flows …”, the authors investigate scenarios of rock burial and exhumation during subduction and orogeny. They present analytical solutions for corner flow involving two adjacent deformable wedges: one representing the overriding plate and the other the accretionary wedge (or, in some cases, the subduction channel), where burial and exhumation may occur. The presented analytical model extends the model of Moulas et al. (GJI, 2021). While Moulas et al. considered only subduction velocities parallel to the wedge base, the present study also examines non-parallel velocities. The authors systematically analyze conditions under which significant return flow arises. They further compare the analytical velocity fields with velocity fields from analogue laboratory experiments for similar configurations and show first-order agreement between experimental and analytical velocity fields. Research on burial and exhumation in subduction zones remains highly relevant, as the controlling mechanisms are still debated. Comparing analytical predictions with analogue experiments is also very relevant. However, I have major concerns about the applicability of the presented analytical solution to natural subduction and orogenic burial-exhumation cycles.
Major comments
1) A key assumption of the analytical solution is that wedge geometry remains constant throughout the burial-exhumation cycle. Both the analytical model and analogue experiments yield essentially instantaneous velocity fields. However, the boundary conditions in the analytical model imply that wedge geometry must change over time, as observed in the experiments. If the velocity at the wedge base is not parallel to its base, the wedge will be squeezed or extended. Assuming average burial and exhumation rates of 10 mm/yr, a 50 km burial followed by 50 km exhumation would last ~10 Myr. In the scenarios presented (slab advance and rollback), non-parallel velocity components of just a few mm/yr would displace the wedge base by several tens of kilometers over this timescale. Such large geometric changes would strongly alter the internal corner flow field. Moreover, the issue of corner “stability” was examined by Moulas et al. (2021), but their findings are not fully considered here. For example, Fig. 4 includes results for viscosity ratios ≤1000, yet at such values the internal wedge boundary would deform significantly during the burial-exhumation cycle, effectively destroying the corner geometry. For these reasons, I am not convinced that the presented velocity fields for slab advance and rollback can reliably be integrated to predict a complete burial-exhumation path.
2) Analytical corner flow models can be useful for certain geodynamic scenarios. However, regarding (U)HP rock exhumation in subduction zones, most 2D thermo-mechanical simulations do not generate wedges with forced corner flow when the subduction zone develops self-consistently (i.e., without a pre-imposed weak zone or wedge). In such models, (U)HP exhumation is typically driven by buoyancy or plate divergence/extension, rather than by forced return flow. Thus, the presented corner flow model may be applicable to burial and exhumation at crustal depths, but for rocks buried deeper than ~35 km, exhumation is more likely controlled by buoyancy. Another limitation is the assumption of constant linear viscosity across a wedge spanning the entire crust or even deeper. In reality, significant temperature variations produce large variations in effective rock viscosity. Also, deeper ductile regions may localize strain into shear zones, allowing the subducting plate to slide beneath the wedge without initiating a distributed corner flow. Hence, the authors should clearer discuss the range of applicability of their model. In particular, I find the reference to both an accretionary wedge and a subduction channel in Fig. 1a problematic. An application to shallow accretionary wedges seems far more realistic than to subduction channels in sub-crustal depths.
Minor comments
The introduction is a bit confusing because the authors do not use a consistent terminology. They mention accretionary wedges in places but seem to refer to subduction channels. In Fig. 1a the corner flow seems to be representative for both accretionary wedges and subduction channels. Buoyancy may be ignored in a shallow accretionary wedge but is likely important for subduction channels at sub-crustal depths.
Line 79-90: In section 2.1. the authors state that they generalize the corner-flow model of Moulas et al. (2021), which is based on coupling two corner flow solutions to consider the deformation of the overriding plate. Also, dynamic pressure fields, velocity fields and a velocity profile showing burial and return velocities have all been shown and discussed in Moulas et al. (2021). Also, Fig. 2 is very similar to figures 2 and 3 in Moulas et al. (2021). Maybe the authors could clearer state in the Introduction that their model is a modification of the model of Moulas et al. (2021). In lines 79-90 the study of Moulas et al. (2021) is not mentioned.
Line 103: Buoyancy is excluded but may play the dominant role for rock exhumation from sub-crustal depths in many orogenies.
Line 236-237: I would argue that even more studies consider the positive buoyancy of subducted rocks as a more important mechanisms for the exhumation of (U)HP rocks that have been subducted to depth larger than the average crustal thickness.
Line 270: If I understood correctly, a value of F>1 implies that more material is exhumed by return flow than is subducted. How can this be applied to natural scenarios? It means that more material is exhumed than buried. How is this in agreement with mass balance and where is the additional material that is exhumed generated? F>1 also implies that the wedge disappears after some time. Likely I missed something, but the authors should better explain the meaning, implications and applications for models with F>1.
Line 274: The lower boundary of the wedge in Fig. 3 indicates the top of the subducting plate (see also Fig. 2a). When the velocity arrows are not parallel to this wedge boundary and point towards the subducting plate, then this implies that material from the wedge flows into the subducting plate, because the lower wedge boundary does not move with the velocity. This makes no sense to me. Keeping the wedge geometry constant but applying boundary velocities that would change the wedge geometry is not consistent.
Line 315-316: The mentioned interface migration has been studied by Moulas et al. (2021) with a 2D dimensionless “regime diagram“ distinguishing a “stable corner“ (with negligible interface migration) and an “unstable corner“ (with significant interface migration).
Section 4.1: Please provide a scaling analysis for the experiment so that it is clear how the laboratory experiments can be scaled to the natural situation.
Best regards,
Stefan Schmalholz

Citation: https://doi.org/10.5194/egusphere-2025-2909-RC2
- AC2: 'Reply on RC2', Nibir Mandal, 17 Oct 2025
  
  We are grateful to the reviewer for thoroughly reviewing the manuscript, and raising critical questions. In the revised version, we carefully address these issues. The comments have greatly helped us improve the quality of our theoretical and experimental results presentations. Please find the response letter in the attached file.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2909-AC2

Peer review completion

AR – Author's response | RR – Referee report | ED – Editor decision | EF – Editorial file upload

AR by Nibir Mandal on behalf of the Authors (27 Nov 2025) Author's response Author's tracked changes Manuscript

ED: Referee Nomination & Report Request started (08 Dec 2025) by Taras Gerya

RR by Anonymous Referee #2 (04 Jan 2026)

Suggestions for revision or reasons for rejection

In the revised manuscript and the authors’ response, the authors have adequately addressed my previous comments. A key strength of the manuscript is the comparison between analytical solutions and results from analogue experiments for corner flow involving a deformable overriding plate and subduction velocities that are not necessarily parallel to the wedge base. However, the specific model geometry, the required strength contrast between the wedge and the overriding plate, and the omission of buoyancy and thermal evolution still make it difficult to assess the significance of the proposed forced return flow for natural rock exhumation. Nevertheless, the study contributes to an improved understanding of forced return flow under a range of corner-flow configurations.

I have a few minor comments that the authors may wish to consider:
Line 77: It could also be stated here that buoyancy is explicitly excluded from the model.
Lines 236–238: This sentence and the cited references are potentially misleading. In the studies by Agard et al. (2009) and Yamato et al. (2007), buoyancy plays a fundamental role in driving return flow. In contrast, the present study investigates forced return flow in the absence of buoyancy. In my view, the studies of Agard et al. (2009) and Yamato et al. (2007) do not support return flow without buoyancy, but rather argue for buoyancy as the dominant mechanism controlling return flow (see, for example, Fig. 13 in Agard et al., 2009, and Fig. 11d in Yamato et al., 2007).
Section 5.2: Regarding the Western Alps, the study by Vaughan-Hammon et al. (2022, Geochemistry, Geophysics, Geosystems) could be cited. That study models syn-convergent burial–exhumation cycles and produces pressure–temperature paths and metamorphic facies gradients that agree, to first order, with observations from exhumed rocks in the Western Alps.

Hide

RR by Anonymous Referee #1 (05 Jan 2026)

Suggestions for revision or reasons for rejection

In the revised version of this manuscript, the authors took care to solve all the comments and did a good job implementing the other reviewer and my reviews. As a result, the introduction and the presentation of the analogue models improved substantially and strengthened the analytical approach. I still have some minor comments based on the newly added sections and especially the meaning of the return flow ratio Fr. If these comments are resolved, this manuscript is suitable for publication in Solid Earth.
Minor comments:
1. As far as I understood, the return flow ratio Fr takes into account both burial/exhumation plus the effect of squeezing/unsqueezing of the accretionary wedge. So, a sentence or two should be added near Line 240 to clarify this parameter. This would avoid my initial confusion and help the readers to understand that Fr>>1 does not imply that more material is being exhumed than subducted and is due to the extrusion of the wedge. This should be reminded in parts of the discussions (e.g., in section 5.1, line 510)
2. When the authors compare exhumation rates between their models and nature, they should use an appropriate scaling for the veylocity. In Line 536, U=3 cm/yr is used to calculate exhumation rates between 4-7 mm/yr (U/4 to U/8) but oceanic plates are much faster (2 to even 4 times faster in some plate reconstructions). Therefore, I recommend using plate velocities appropriate for the tectonic setting being compared.
3. In the analogue models, section 5.1 and especially Fig 10b, how realistic are the OP-wedge viscosity ratios and the scaled wedge viscosity (5e20 Pa*s; line 387) compared to estimates through power-laws? A more in-depth comparison is missing here, the composition, thermal structure and deformation mechanisms of the overriding plate (oceanic or continental) and the accretionary wedge provide clues on the plausible range of viscosity ratios for different settings, which ultimately impact the conditions for return flows. Also, is it realistic to consider viscosity ratios < 1?

Line to line comments:
Line 20: delete ‘the’ in ‘the subduction’
Line 28: add ‘the’ to ‘the Earth’s…’
Line 28: ‘in course of the convergence movements’ this sentence could be replaced with ‘during active subduction’
Line 50: remove ‘the’ from ‘the accretionary wedge’
Line 84: ‘Subduction’ to ‘subduction’
Lines 87-88: Here I recommend just describing it as a non-parallel component on the slab velocity, just like it was said in lines 122-123.
Lines 116-117: Rewrite this sentence, as rocks deform but flow in a solid-state manner rather than being a fluid, and the viscosity is an approximation to replicate such deformation.
Lines 239-240: In this sentence I recommend including that the return flux Fr is the ratio of material being returned and the effect of squeezing the wedge to the material being buried plus the ‘unsqueeze’ of the wedge
Line 363: typo ‘configuration’
Line 379: remove ‘The’
Line 471: here I recommend including what each parameter mean: slab advance/rollback, high viscosity contrasts, etc.
Line 483: I would change the word ‘explains’ to ‘predicts’ as there are not natural examples mentioned in the text
Line 490: ‘characteristics’ to ‘processes’
Line 501: Only the viscosity of the overriding plate or also the wedge?
Lines 521-523: Here I would reconsider mentioning subduction erosion as a process since it is mostly linked to the removal of the base of the overriding plate and not the same as sediment subduction, which I think is what the authors mean in this sentence. In addition, subduction erosion is facilitated in sediment starved margins and therefore not applicable in the context of accretionary wedges.
Line 576: typo ‘horizontal’

Hide

ED: Publish subject to minor revisions (review by editor) (09 Jan 2026) by Taras Gerya

AR by Nibir Mandal on behalf of the Authors (16 Jan 2026) Author's response Author's tracked changes Manuscript

ED: Publish as is (25 Feb 2026) by Taras Gerya

ED: Publish as is (07 Mar 2026) by Florian Fusseis (Executive editor)

AR by Nibir Mandal on behalf of the Authors (14 Mar 2026) Manuscript

Journal article(s) based on this preprint

25 Mar 2026

On the criticality of return flows in viscous accretionary wedges and its implications for deep-crustal exhumation in subduction zones

Ayan Patsa and Nibir Mandal

Solid Earth, 17, 573–600, https://doi.org/10.5194/se-17-573-2026,https://doi.org/10.5194/se-17-573-2026, 2026

Short summary

Ayan Patsa and Nibir Mandal

Supplement

https://doi.org/10.5194/egusphere-2025-2909-supplement

Ayan Patsa and Nibir Mandal

Viewed

Total article views: 1,968 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	Supplement	BibTeX	EndNote
1,708	230	30	1,968	88	95	99

HTML: 1,708
PDF: 230
XML: 30
Total: 1,968
Supplement: 88
BibTeX: 95
EndNote: 99

Views and downloads (calculated since 12 Aug 2025)

Month	HTML	PDF	XML	Total
Aug 2025	295	21	2	318
Sep 2025	1,051	11	3	1,065
Oct 2025	78	13	10	101
Nov 2025	44	15	2	61
Dec 2025	55	41	4	100
Jan 2026	88	25	3	116
Feb 2026	42	29	3	74
Mar 2026	55	75	3	133
Apr 2026	0

Cumulative views and downloads (calculated since 12 Aug 2025)

Month	HTML	PDF	XML	Total
Aug 2025	295	21	2	318
Sep 2025	1,051	11	3	1,065
Oct 2025	78	13	10	101
Nov 2025	44	15	2	61
Dec 2025	55	41	4	100
Jan 2026	88	25	3	116
Feb 2026	42	29	3	74
Mar 2026	55	75	3	133
Apr 2026	0

Viewed (geographical distribution)

Total article views: 1,945 (including HTML, PDF, and XML) Thereof 1,945 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 09 Apr 2026

Download

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Preprint (69556 KB)
Metadata XML

Short summary

Accretionary wedges are the prime locations of exhumed high-pressure (HP) and low-temperature (LT) metamorphic rocks. Previous tectonic models invoked the corner flow theory with a premise of slab-parallel motion to explain the upward return flow of buried metasediments. In this study, we develop a generalized corner flow model with additional kinematic and rheological factors and evaluate the limiting conditions in which a wedge can set in significant return flows.


Total:	0
HTML:	0
PDF:	0
XML:	0