{"id":"e7f461fd-43fb-44e7-a943-eba9f9d5bb55","arxiv_id":"2501.08021","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The Sun's near-surface shear layer contains a distinct 8 Mm deep leptocline with strong rotational gradients, enhanced turbulence, and cycle-dependent variations.","lead":"This paper combines helioseismic observations and 3D simulations to describe the leptocline, a shallow 8 Mm layer of sharp rotational shear just below the Sun's surface. It also links changes in this layer's rotation gradient and seismic radius to the solar activity cycle.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Global inversions with l <= 300 cannot resolve an 8 Mm leptocline; the quantitative depth/gradient structure in Figs. 2-4 needs a synthetic-recovery test before being accepted.","rationale":"The paper's central claim—that the NSSL contains an approximately 8 Mm deep leptocline with enhanced shear and turbulent convection—is supported by three independent strands: global helioseismic gradients, local ring-diagram measurements, and 3D RHD simulations. The strand the paper itself presents as new analysis is the global inversion, and it is also the strand with the weakest resolving power near the surface. The authors honestly flag the smoothing from averaging kernels but do not quantify how much this affects the depth, amplitude, or latitude dependence shown in Figs. 2-4. Because averaging kernels in the top 20 Mm are likely several Mm wide, a synthetic-inversion recovery test is needed before the apparent 8 Mm scale and the high-latitude sign reversal can be interpreted as physical. This is not a fatal objection: local helioseismology and simulations independently indicate a shallow shear layer, so the existence claim survives. It does, however, reinforce the reader's CONDITIONAL verdict: the quantitative leptocline properties should be accepted conditionally on a resolution verification. I did not make the f-mode radius contamination the primary concern because it is ancillary to the central claim and already acknowledged in the paper. The simulation-versus-observation gradient discrepancy (-4 versus -1.5/-2.6) is a real internal inconsistency, but it changes the amplitude rather than the existence of the leptocline. On balance, the reader's conditional assessment already captures the main risk, so no verdict change is needed.","tokens_in":11426,"tokens_out":10392,"duration_ms":104181,"concrete_test":"Compute the radial averaging-kernel widths for the l <= 300 inversion in the 0-20 Mm depth range, and run a synthetic recovery test: generate noiseless rotational frequency splittings from a forward model that contains a leptocline-like profile (a sharp gradient layer about 8 Mm deep, with d log Omega / d log r near -4 at low latitudes and positive values poleward of about 60 degrees), then invert these splittings with the same pipeline used for Figs. 1-4. If the recovered gradients resemble Fig. 2, the global data are consistent with a leptocline but cannot determine its true sharpness or depth. If the synthetic leptocline is smoothed into a broad NSSL-like gradient or the high-latitude sign reversal is not reproduced, then the feature labeled 'leptocline' in Figs. 2-4 is not resolved, and the quoted depth and gradient amplitudes should not be used as quantitative constraints.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing weakness is that the paper's own global helioseismic analysis cannot resolve the feature it characterizes. The inversions use modes with angular degree up to l = 300, and Section 2 concedes that these data 'do not resolve sharp variations near the surface' because gradients are smoothed by averaging kernels. Despite this, Figs. 2-4 are interpreted as showing the leptocline's 8 Mm depth, its latitudinal gradient reversal, and its solar-cycle variations. If the radial averaging-kernel width near r/R = 0.985-1.0 is comparable to or larger than 8 Mm, the apparent gradient enhancement and its latitude/time structure could be inversion artifacts rather than a physically resolved leptocline. This does not overturn the central claim, because local ring-diagram measurements (Komm 2022; Rabello Soares et al. 2024) and 3D RHD simulations (Kitiashvili et al. 2023) independently point to a shallow near-surface shear layer. However, it means the specific quantitative claims about the leptocline's depth and gradient are not established by the analysis presented in this paper and should be treated as conditional on a resolution check.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper analyzes global helioseismic rotation inversions from SoHO/MDI and SDO/HMI for 1996–2024, supplemented by local helioseismology and 3D radiative hydrodynamic simulations, to argue that the Near-Surface Shear Layer (NSSL) contains a distinct shallow substructure, the 'leptocline', approximately 8 Mm deep. The paper reports that the radial gradient of rotation, d log Ω / d log r, is strongly negative at low latitudes in the leptocline, reverses sign at high latitudes, varies with the solar cycle in sunspot-forming regions, and that f-mode-derived seismic radius variations are strongest near 5 Mm depth, inside the leptocline. The manuscript is explicitly built on prior work by the same authors and collaborators, including the naming of the leptocline and the simulations, and it is candid about the limited spatial resolution of global helioseismic inversions.","tokens_in":11627,"tokens_out":5252,"duration_ms":55673,"significance":"If the quantitative claims hold, the leptocline would be a dynamically important shallow layer controlling angular momentum transport, the formation of sunspots and active regions, and the solar-cycle variation of the seismic radius. The paper benefits from the convergence of global helioseismology, local ring-diagram measurements, and 3D simulations, and it explicitly cautions that inversion kernels smooth sharp gradients. However, the central quantitative content—the 8 Mm depth, the gradient amplitudes, and the latitude/time structure shown in Figs. 2–4—is not self-contained: it depends on low-resolution global inversions and on prior publications. The paper is a useful synthesis and a plausible qualitative case, but its specific quantitative characterization of the leptocline needs a resolution validation before it can be accepted as established by the new analysis.","major_comments":[{"comment":"The global helioseismic inversions use only modes with angular degree up to l=300, and the paper itself states that 'these data do not resolve sharp variations near the surface' and that gradients are smoothed by averaging kernels. This limitation is load-bearing because the central quantitative claims—the 8 Mm leptocline depth, the gradient value of about -1.5 at low latitudes, and the sign reversal at high latitudes—are read directly from these smoothed profiles. Please add a synthetic-recovery test: construct a rotation profile containing a sharp 8-Mm leptocline with known gradient values, synthesize rotational splittings for the same mode set and inversion procedure, and compare the recovered gradient profiles with the input. If the averaging-kernel width near r/R = 0.985–1.0 is comparable to or larger than 8 Mm, the apparent depth, amplitude, and latitudinal structure in Figs. 2–3 are not established by this data set alone. The qualitative existence of a shallow shear layer is independently supported by local helioseismology and simulations, but the quantitative characterization presented here requires this test.","section":"§2, Figs. 2–3"},{"comment":"The time-latitude maps of d log Ω / d log r in the leptocline (Fig. 4c) are derived from the same low-resolution global inversions. The claim that the gradient is enhanced in the leptocline during activity minima as well as maxima, and the contrast between the 'below leptocline' (Fig. 4b) and 'in leptocline' (Fig. 4c) patterns, could be affected by changes in the radial averaging kernels with time or with the mode set. Please show that these patterns survive the synthetic-recovery test described above, or at least provide the radial averaging-kernel widths for the relevant depths and epochs. Without this, the apparent solar-cycle contrast may reflect a varying sensitivity rather than a physical change in the leptocline gradient.","section":"§3, Fig. 4"},{"comment":"The inference that solar-cycle variations of the seismic radius are strongest at about 5 Mm depth, inside the leptocline, relies on f-mode frequency shifts, and the paper itself notes that magnetic field and temperature effects are difficult to separate. Since Fig. 5d is used to support the leptocline's role in cyclic changes, the authors should either quantify the likely systematic error from magnetic and thermal contamination (for example, using sensitivity kernels or comparing with p-mode-based radius estimates) or explicitly soften the claim to say that the variations are 'consistent with' a maximum inside the leptocline rather than 'strongest at' 5 Mm. As written, the quantitative depth localization is not supported by the analysis presented.","section":"§4, Fig. 5"}],"minor_comments":[{"comment":"The word 'superadiabiatic' should be 'superadiabatic'.","section":"§1"},{"comment":"The phrase 'near-shear shear layer' should be 'near-surface shear layer'.","section":"§3"},{"comment":"The expansion 'Near Sub Surface Layer' should be 'Near-Surface Shear Layer' to match the acronym NSSL introduced in the abstract and used throughout the paper.","section":"§6"},{"comment":"Several references are incomplete, e.g., Kholikov and Hill (2008) lacks an article title, and Javaraiah (2003) lacks an article title; the reference list should be completed for consistency.","section":"References"},{"comment":"The captions and panel descriptions for Figures 2 and 3 appear somewhat scrambled: the text refers to a latitude dependence at three depths, but the Figure 2 caption only describes depth. Please ensure each figure panel is described by its own caption.","section":"Figure 2 and 3 captions"}],"recommendation":"major_revision","confidential_remarks":"The paper is a short proceedings contribution that relies heavily on the authors' previous publications for the naming, simulation, and seismic-radius analyses. The referee's main concern is that the new global-inversion figures are presented with quantitative depth and amplitude claims that the mode set cannot resolve; a synthetic-recovery test is feasible and should be added within the manuscript's scope. If the authors prefer to keep the proceedings format brief, they should weaken the quantitative claims and present the global-inversion results as qualitative support, deferring the quantitative characterization to higher-resolution local helioseismology."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a proceedings-style synthesis of the leptocline argument, with one genuinely new piece—a 1996-2024 averaged map of the radial gradient of rotation from public SoHO/SDO inversions—plus time-latitude diagrams connecting the gradient to the butterfly diagram. What is new is largely descriptive. The conceptual load is carried by earlier papers (Kitiashvili et al. 2023; Komm 2022; Rabello Soares et al. 2024; Kosovichev & Rozelot 2018), and the paper is transparent about that.\n\nCredit where due: the new averaged maps are a reasonable thing to have in the literature, the text explicitly warns that global inversions with l≤300 smooth the near-surface gradients and that true values may be larger, and the authors flag the mismatch between global and local helioseismic patterns rather than papering over it. The writing is clear and the citation pattern is fine.\n\nThe soft spots are real, and the main one is exactly what the stress test says. The paper spends Figures 2-4 presenting the leptocline's 8 Mm depth, its latitudinal gradient reversal, and its solar-cycle variations as if these were read directly off the inversions, despite the admission that the averaging kernels cannot resolve 8 Mm near the surface. If the kernel width there is comparable to the claimed layer thickness, the apparent features—especially the sign reversal at high latitude—could be inversion artifacts. A synthetic recovery test (invert a known leptocline-like profile with the same mode set and see what comes back) would settle this. Until then, the quantitative depth and amplitude claims are conditional. That is a load-bearing soft spot, but not a fatal one, because the local helioseismology and the 3D simulations independently point to a shallow sharp shear layer.\n\nA second, smaller issue: Section 5 says the simulated gradient of about -4 is \"in agreement with observations,\" but the global inversions show about -1.5 and the local ring-diagram values about -2.6. Those differ by a factor of 1.5-2.5. Agreement is too strong a word.\n\nWho is this for? Readers who want a compact orientation to the leptocline discussion and a first look at the long-baseline gradient maps. It is not the definitive paper on the subject. If submitted as a regular journal article, it should be refereed: the central claim is important enough, and the resolution issue is serious enough, that the quantitative conclusions need to be checked before being cited as established. Recommendation: send it to review, and ask for the resolution test and a more careful simulation-observation comparison.","headline":"A candid conference summary of the leptocline case: useful new gradient maps, honest caveats, but the quantitative depth claims need a resolution check before they are taken as established.","tokens_in":12183,"tokens_out":2652,"would_cite":false,"duration_ms":25080,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The Sun's Near-Surface Shear Layer contains a distinct, about 8-Mm-deep substructure—the leptocline—where the rotational shear sharply intensifies and the strongest solar-cycle variations of the seismic radius occur.","keywords":["Near-Surface Shear Layer","leptocline","solar rotation gradient","helioseismology","seismic radius","solar cycle","torsional oscillations"],"falsifier":"An observation that resolves the top 10 Mm of the Sun with high-degree modes (angular degree far beyond 300, as local helioseismology approaches) and finds that the logarithmic radial gradient of rotation does not steepen near the H i/He i ionization zone—or that the steepening is a fixed artifact of the inversion—would falsify the claim that the leptocline is a distinct shear layer. Concretely, if a resolved inversion shows $d\\log\\Omega/d\\log r$ remaining near $-1$ continuously from 30 Mm to the surface, with no sharp excursion around 3–8 Mm, the leptocline as a physical layer would be ruled out.","tokens_in":11230,"feed_emoji":"☀️","tokens_out":12235,"duration_ms":95703,"temperature":0.7,"pith_summary":"This paper assembles global helioseismic observations from two decades (1996–2024) to argue that the Sun's Near-Surface Shear Layer (NSSL) is not a simple monotonic shear but contains a distinct, about 8-Mm-deep substructure just below the surface, which the authors call the leptocline. In this layer the radial gradient of rotation steepens sharply, enhanced turbulent convection from the hydrogen ionization zone modifies the stratification, and the strongest solar-cycle variations of the Sun's seismic radius occur. A sympathetic reader would care because the leptocline may be the place where the dynamo-generated magnetic field is concentrated and where sunspot-forming regions first feel the rotational shear, potentially linking the observed magnetic butterfly diagram to subsurface dynamics. The paper combines global helioseismology, local ring-diagram analyses, and 3D radiative hydrodynamics simulations to make the case, and notes that the true gradients may be even stronger than inversion-resolved values.","feed_headline":"Solar rotation shear is sharpest inside an 8-Mm leptocline","feed_subtitle":"The layer drives solar-cycle changes in the Sun's radius and may seed sunspot formation.","key_machinery":"The central object is the leptocline: a shallow, sharp rotational-shear layer about 8 Mm deep at the top of the solar convection zone (roughly $0.985$–$1.0\\,R_\\odot$), named by analogy with the tachocline. It is characterized by a steep radial gradient of rotation, enhanced turbulent convection rooted in the hydrogen and helium ionization zones, and self-organized meridional flows. The argument is carried by three complementary tools: (1) global helioseismic inversions of rotation splittings from SoHO/MDI and SDO/HMI data, which provide the average gradient structure and its time-latitude evolution; (2) f-mode frequency inversions, which track the solar-cycle displacement of subsurface layers and define the seismic radius; and (3) 3D radiative hydrodynamics simulations, which reproduce the leptocline and suggest its origin in anisotropic turbulent downdrafts that overshoot from the H i/He i ionization zone into the layer below. The simulations also indicate that the true rotational gradient in the leptocline may be as steep as about $-4$, much larger than the smoothed inversion values.","core_discovery":"The paper's central claim is that the Sun's Near-Surface Shear Layer has a distinct substructure, the leptocline, which occupies the top approximately 8 Mm of the convection zone (roughly between $0.985\\,R_\\odot$ and $1.0\\,R_\\odot$). In this layer the logarithmic radial gradient of the rotation rate, $d\\log\\Omega/d\\log r$, increases sharply from about $-1$ in the deeper NSSL to about $-1.5$ at the equator and low latitudes, while at higher latitudes it becomes positive, so rotation speeds up toward the surface; the inversion results are smoothed by averaging kernels, so the actual gradients may be considerably larger. The leptocline coincides with the hydrogen and helium ionization zones, where enhanced anisotropic turbulent overshooting produces strong density fluctuations and turbulent mixing. The paper further reports that the rotational gradient in and just below the leptocline varies with the solar cycle: below the leptocline the enhanced gradient follows the magnetic butterfly diagram, whereas inside the leptocline the pattern is more complex and resembles the overlapping extended solar cycles of the torsional oscillations. Finally, inversion of f-mode frequencies shows the solar-cycle variations of the seismic radius are strongest at a depth of about 5 Mm, in the middle of the leptocline, with the Sun's average radius shrinking by a few kilometers near activity maxima.","pith_inferences":["If the leptocline is a generic feature of cool stars with near-surface convection and hydrogen ionization zones, similar sharp shear layers should appear in other solar-like stars and could be probed through high-degree asteroseismic rotational splittings.","The anti-correlation between seismic radius and activity, strongest at 5 Mm, offers a testable energy-reservoir mechanism: if the leptocline stores or releases energy with the cycle, total irradiance variations and radius changes should track the layer's density and temperature perturbations in f-mode inversions.","The discrepancy between global-inversion and ring-diagram reversal latitudes of the gradient is a concrete place to look: a resolution-aware joint inversion could show whether the high-latitude positive gradient is real or a smoothing artifact, and whether it drifts with the extended cycle.","If the leptocline's gradient is as steep as simulations indicate, the magneto-rotational instability in this layer could be strong enough to matter for near-surface dynamo models, making the leptocline a boundary condition for theories that place dynamo action close to the surface."],"forward_implications":["Solar rotation models must treat the leptocline as a distinct layer: the sharp shear there means the NSSL is not a single power-law gradient but a two-step structure near the surface.","The enhancement of the rotational gradient in sunspot-forming regions, with a butterfly-like pattern of enhanced gradient below the leptocline, implies the leptocline participates in the magnetic activity cycle, likely in the formation or emergence of active regions.","The solar-cycle seismic-radius variations peak at about 5 Mm depth, inside the leptocline, meaning the Sun's small activity-related contractions are driven by subsurface stratification changes in this layer rather than by a global radial change.","Because inversion averaging kernels smooth sharp structures, the actual leptocline gradient may be significantly steeper than the reported $-1.5$, which strengthens the layer's dynamical importance for shear-driven instabilities and dynamo action.","The convergence of global helioseismology, local ring-diagram analysis, and 3D simulations on a shallow shear-enhanced layer suggests the leptocline is a physical feature, not an inversion artifact."],"supporting_citations":[{"why":"3D radiative hydrodynamics simulations that first revealed the leptocline as a shallow substructure and identified its origin in the hydrogen ionization zone with overshooting.","marker":"Kitiashvili et al. 2023"},{"why":"Local helioseismology (ring-diagram) analysis resolving the NSSL substructure and finding gradients reaching about $-2.6$ near 3 Mm depth.","marker":"Rabello Soares et al. 2024"},{"why":"Ring-diagram measurements of the radial gradient of rotation in the NSSL and leptocline, including the quiet-Sun enhancement and the latitude-radius diagram used for comparison.","marker":"Komm 2022"},{"why":"F-mode inversion establishing solar-cycle variations of the seismic radius with depth, with the strongest variations around 5 Mm.","marker":"Kosovichev and Rozelot 2018"},{"why":"Global-mode analysis pipeline used to produce the 72-day rotation inferences from SoHO/MDI and SDO/HMI data.","marker":"Larson and Schou 2018"},{"why":"Helioseismic inversion methodology and earlier determination of differential rotation and the near-surface shear layer, including the averaging-kernel smoothing caveat.","marker":"Schou et al. 1998"},{"why":"Initial MDI medium-l results establishing the rotation structure of the solar interior used as a baseline.","marker":"Kosovichev et al. 1997"},{"why":"Coined the term 'leptocline' and provided first evidence of the shallow layer at the top of the convection zone.","marker":"Godier and Rozelot 2001"}],"fun_headline_variants":["Sun's leptocline drives cycle changes in rotation and radius","Sharp shear layer at 8 Mm reshapes solar rotation","Leptocline: key to Sun's rotation shear and activity cycle","Solar leptocline links rotation shear to radius changes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results rest on helioseismic inversions whose averaging kernels smooth the sharpest near-surface structure, so the leptocline's reported gradients and their cycle variations could be diluted versions of stronger underlying values.","fun_headline_variants_meta":{"raw":{"variants":["Sun's leptocline drives cycle changes in rotation and radius","Sharp shear layer at 8 Mm reshapes solar rotation","Leptocline: key to Sun's rotation shear and activity cycle","Solar leptocline links rotation shear to radius changes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000211,"raw_usage":{"total_tokens":1450,"prompt_tokens":1013,"completion_tokens":437,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":629,"completion_tokens_details":{"reasoning_tokens":367}},"tokens_in":629,"tokens_out":437,"duration_ms":4293,"temperature":1.0,"reasoning_tokens":367,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:28:47.726366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An observation that resolves the top 10 Mm of the Sun with high-degree modes (angular degree far beyond 300, as local helioseismology approaches) and finds that the logarithmic radial gradient of rotation does not steepen near the H i/He i ionization zone—or that the steepening is a fixed artifact of the inversion—would falsify the claim that the leptocline is a distinct shear layer. Concretely, if a resolved inversion shows $d\\log\\Omega/d\\log r$ remaining near $-1$ continuously from 30 Mm to the surface, with no sharp excursion around 3–8 Mm, the leptocline as a physical layer would be ruled out.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"3D radiative hydrodynamics simulations that first revealed the leptocline as a shallow substructure and identified its origin in the hydrogen ionization zone with overshooting."}],"review_version":1}