{"id":"05231f55-47e5-48ee-8a93-a8dfae2be025","arxiv_id":"2506.05758","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The proposed 'adjacent correlation map' is a restatement of the gradient structure tensor and produces spatial maps of a local correlation angle and degree.","lead":"This paper makes maps that show, at each location, how two measured quantities change together, by combining the local gradients with a polarization-style trick. The authors apply the maps to rainfall and temperature, to simulated magnetized turbulence, and to Turing patterns, and argue the maps separate distinct physical regions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The scalar p in Eq. (8) is a coherence index of the gradient structure tensor, not a correlation degree; it can approach 1 for independent fields, so the central map does not establish local correlations or regime separation.","rationale":"The reader's stated weakest assumption is the unspecified gradient scale and lack of robustness tests. I agree that scale is a substantive concern, but the more load-bearing issue is that the central scalar is not measuring correlation at all. If p can be high for independent fields, no choice of smoothing scale or window size can repair the claim. The reader did identify the Eq. (8) typo and the eigenvalue misinterpretation in the rationale, but did not elevate this to the primary weakness; hence partial agreement. This is an internal mathematical mismatch rather than a dispute about novelty or consensus, so it directly invalidates the central assertion as written. The GitHub code availability is helpful, but it does not resolve the interpretation; the proposed synthetic experiment would settle whether p really behaves as a correlation degree. Because the central claim fails on this more basic ground, the reader's REJECT verdict remains appropriate.","tokens_in":6648,"tokens_out":6931,"duration_ms":63400,"concrete_test":"Generate a 256x256 synthetic image pair: p1(x,y) = x and p2(x,y) = 0.01 * eta(x,y), where eta is independent Gaussian noise. Apply Eqs. (5)-(8) with a fixed stencil or window, for example a 9x9 averaging of the per-pixel Stokes parameters, and compute the median p over interior pixels. If the median p is close to 1 while the Pearson correlation between p1 and p2 is near 0, then p is not a correlation degree and the central claim fails. As a second check, rerun with p2 = 2*x: a coherence index should again give p near 1, but the two cases can be distinguished by also plotting the joint distribution of the gradient vectors, which is collinear in the second case and nearly isotropic in the first.","verdict_should_be":"UNCHANGED","load_bearing_attack":"To support the central claim, the map must measure the nature and strength of the local correlation between p1 and p2. The quantity called 'correlation degree' in Eqs. (7)-(8) does not do this. After correcting the obvious missing square in Eq. (8), p = sqrt((Q/I)^2 + (U/I)^2) = (lambda_plus - lambda_minus)/(lambda_plus + lambda_minus), where lambda_plus and lambda_minus are the eigenvalues of the summed outer product of the two gradient components. This is the fractional anisotropy, or coherence, of the gradient distribution, not a correlation coefficient between the two fields. A concrete counterexample: take p1 = x on a uniform grid and p2 = epsilon * eta(x,y), where eta is independent noise and epsilon is small. At every interior pixel the gradient vector is approximately (1, 0) in the (p1, p2) plane, so p is near 1, while the Pearson correlation between p1 and p2 is near 0. Thus the map can label regions as 'perfectly correlated' where the two physical quantities are statistically independent. The claimed connection to the correlation matrix in Eq. (11) is also inconsistent: p equals (lambda_plus - lambda_minus)/(lambda_plus + lambda_minus), not lambda_plus/(lambda_plus + lambda_minus) as stated. The regime interpretations in Sections 4.1-4.3 are therefore assignments to a quantity that has not been shown to be a correlation measure.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an 'adjacent correlation map' that uses Stokes parameters constructed from local gradients of two scalar fields to produce spatially resolved maps of the 'type' and 'degree' of correlation between those fields. The method is applied to climate data (temperature vs. precipitation), a simulation of compressive MHD turbulence (density vs. magnetic field), and a Turing pattern (activator vs. inhibitor). The paper claims that the resulting maps separate a domain into subregions corresponding to physical regimes or climate zones, and it introduces a distinction between 'correlated' (Type 1) and 'stiff' (Type 2) regularities. A connection to Buckingham's Pi theorem is also sketched.","tokens_in":6951,"tokens_out":5091,"duration_ms":50113,"significance":"If the central claim were correct, the method would provide a simple and widely applicable tool for visualizing local correlations in spatially inhomogeneous data, with potential uses in climate science, astrophysics, and pattern formation. The Stokes-parameter construction is elegant, and the author provides code on GitHub. However, the paper's central quantity, the 'correlation degree' p, is not a measure of correlation between the two fields; it is the fractional anisotropy (coherence) of the gradient structure tensor. This fundamentally undermines the interpretation of every map in the paper. The mathematical errors in Eq. (8) and in the stated relation to eigenvalues of a correlation matrix reinforce the concern. The method may still be useful as a coherence visualizer, but that is not what the paper claims or demonstrates.","major_comments":[{"comment":"Equation (8) is missing a square on the U/I term: the formula as written, p = ((Q/I)^2 + (U/I))^{1/2}, is dimensionally inconsistent and algebraically wrong. The correct expression is p = sqrt((Q/I)^2 + (U/I)^2). With this correction, p equals (lambda_plus - lambda_minus)/(lambda_plus + lambda_minus), where lambda_plus and lambda_minus are the eigenvalues of the averaged outer product of the gradient vectors, i.e., the fractional anisotropy of the gradient distribution. The text's interpretation that p = 0 means 'non-correlation' and p = 1 means 'perfect correlation' (Section 3.1, after Eq. (10)) is not supported by this algebraic form.","section":"3.1, Eq. (8)"},{"comment":"The 'correlation degree' p does not measure the degree of correlation between the two fields p1 and p2. Consider p1 = x on a uniform grid and p2 = epsilon * eta(x,y), where eta is independent noise and epsilon is small. At every interior point the gradient vectors in the (p1,p2) plane are approximately (1,0), so Q/I is near 1, U/I is near 0, and p is near 1, while the Pearson correlation between p1 and p2 is near 0. The proposed map would label statistically independent fields as 'perfectly correlated'. This directly contradicts the abstract and Section 6, which claim the method provides 'a spatially resolved view of the nature and strength of correlations'.","section":"3.1, Eq. (8) and Figs. 2-5"},{"comment":"The claimed relation to the correlation matrix M_ij = <(p_i - \\bar p_i)(p_j - \\bar p_j)> is not correct. The Stokes parameters are constructed from spatial gradients of the fields, not from the fields themselves, so the eigenvalues of the summed outer product of gradients do not coincide with the eigenvalues of the correlation matrix. The statement that the correlation degree 'reflects l_max/(l_max + l_min)' is also algebraically wrong; the correct expression after fixing Eq. (8) is (l_max - l_min)/(l_max + l_min). This incorrect theoretical justification is load-bearing because it is used to argue that p is a correlation measure.","section":"3.1, Eq. (11)"},{"comment":"The scale at which the gradients in Eq. (2) are evaluated is never specified. The maps in Figs. 2-5 depend on this scale (e.g., nearest-neighbor finite differences versus a smoothed derivative), yet no robustness tests are reported. Without such tests, the claimed separation into 'physical regimes' in Sections 4.1-4.3 could be an artifact of the gradient stencil rather than a property of the data. This is particularly problematic for a paper whose central claim is that the method is a systematic approach to regularity discovery.","section":"3.1 and Section 4"}],"minor_comments":[{"comment":"The title 'Regularizes from spatially-inhomogeneous systems' appears to be a typo for 'Regularities from spatially-inhomogeneous systems'.","section":"Section 2 title"},{"comment":"The phrase 'correction degree' is used repeatedly (e.g., Figures 2 and 4 captions, lower panels) and should read 'correlation degree'.","section":"Figures 2-5 captions and text"},{"comment":"In the description of Fig. 1, 'perception' should be 'precipitation'.","section":"Section 2, text"},{"comment":"The reference to '(Li 2015 submitted)' appears to be a typo; the companion paper is cited elsewhere as '(Li 2025 submitted)'.","section":"Section 4.2"},{"comment":"The title 'Relation to the Buckinghum's Pi theorem' misspells Buckingham, and the text contains 'visa versa' instead of 'vice versa'.","section":"Section 5.1"},{"comment":"The paper relies heavily on an unpublished companion paper (Li 2025 submitted) for the concept of locally-correlated variations and for the phase-space interpretation of the Stokes parameters; this is not ideal for a self-contained submission and the relevant definitions should be repeated here or the companion paper made available.","section":"General"},{"comment":"The notation E_{p1} and E_{p2} for the pseudo-Stokes parameters is confusing, as it suggests electric field components; a different symbol or explicit definition would improve clarity.","section":"Eq. (7)"}],"recommendation":"reject","confidential_remarks":"The paper's main methodological claim fails because p is a coherence index of gradient directions, not a correlation coefficient between fields; the counterexample with p1=x and p2=small independent noise shows that high p can occur for statistically independent fields. The algebraic errors in Eq. (8) and in the eigenvalue relation would be correctable, but the interpretational gap is fundamental and would require rewriting the central claim. Additionally, the lack of any robustness analysis regarding the gradient scale, and the heavy dependence on an unpublished companion paper, make it difficult to assess the novelty and reliability of the results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe paper is a clean and well-written re-derivation of the gradient structure tensor, renamed as \"adjacent correlation map.\" The Stokes-parameter packaging is neat, and the maps are visually appealing. But the central claim—that the resulting \"correlation degree\" p measures the strength of correlation between two fields—does not hold up. p is actually the coherence (fractional anisotropy) of the local gradient vectors, not a statistical correlation. A simple counterexample: if p1 = x and p2 is small independent noise, p is near 1 everywhere while the Pearson correlation is near 0. So the map cannot be interpreted as a spatially resolved correlation coefficient.\n\nThe paper also contains two concrete errors: Eq. (8) is missing the square on (U/I), and the stated relation to eigenvalues lmax/(lmax+lmin) should be (lmax-lmin)/(lmax+lmin). Those are easy to fix but underline the conceptual confusion.\n\nWhat is genuinely useful: the method provides a simple way to visualize regions where the local gradient of one quantity relative to another has a consistent orientation. The three applications—climate, MHD turbulence, reaction-diffusion—show that the maps can segment data into visually coherent patches. If reframed as an application of the classical structure tensor (e.g., in image processing) and the \"correlation degree\" were renamed something like \"gradient coherence,\" the examples would make a decent short application note. The code is on GitHub, which is good for reproducibility.\n\nThe main soft spot is the novelty claim. The gradient structure tensor has been used for decades in computer vision for edge and corner detection, and the Stokes parameters are just a basis for the same second-moment matrix. The paper does not cite any of that literature, and instead refers to an unpublished companion paper. That is a significant gap in engagement.\n\nThe physical-regime interpretations in Sections 4.1–4.3 feel post hoc: the boundaries depend on an unspecified gradient stencil scale, and no robustness checks are given. That said, the core algebra, once corrected, is straightforward.\n\nWho is this for? A reader with a spatially-varying dataset who wants a quick look at local gradient orientations might find the code handy. But as a contribution to physics or data visualization, it needs substantial revision. I would not cite it as a new method.\n\nRecommendation: Send to peer review? Yes—the paper is clearly written, technically testable, and the flaws are the kind referees should catch. It will likely be rejected as submitted, but it deserves proper review rather than a desk reject.","headline":"A well-written re-derivation of the gradient structure tensor as an 'adjacent correlation map,' but the key quantity is coherence, not correlation, and the novelty claim needs major reframing.","tokens_in":7448,"tokens_out":4877,"would_cite":false,"duration_ms":44468,"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":"A pixel-by-pixel map of how two quantities correlate, built from spin-2 Stokes parameters, splits data into physically distinct regions.","keywords":["adjacent correlation map","Stokes parameters","spatially resolved correlations","locally-correlated variations","correlation angle","correlation degree","regularity discovery","data visualization"],"falsifier":"Take the North America temperature-precipitation map and recompute the correlation-angle map with the gradient stencil spacing doubled, then halved: if the boundaries between negative-, positive-, and weakly-correlated regions shift by more than the stated climate-zone uncertainty, the claim that the map identifies physical regimes is not robust.","tokens_in":6438,"feed_emoji":"🗺️","tokens_out":5388,"duration_ms":55037,"temperature":0.7,"pith_summary":"Global correlation plots wash out structure when data are spatially inhomogeneous, because opposite local correlations cancel. This paper proposes a new representation, the adjacent correlation map, that computes local gradient pairs of two quantities at each location and treats them as spin-2 vectors, combining them through Stokes parameters. The resulting maps show the correlation angle and correlation degree at every point, so a dataset separates into coherent patches that share the same local correlation. Applied to weather data, MHD turbulence simulations, and Turing patterns, the maps highlight regions that correspond to climate zones, magnetically versus kinematically dominated regimes, and self-regulating regions. The central claim is that this procedure is a systematic, broadly applicable tool for discovering spatially resolved regularities that global correlation analysis misses.","feed_headline":"A pixel-by-pixel map reveals where two fields correlate","feed_subtitle":"Treating local gradient pairs as spin-2 vectors splits climate, turbulence, and Turing patterns into coherent regions.","key_machinery":"The machinery is the adjacent correlation vector: at each location the gradients of the two quantities form a vector $(\\partial p_1/\\partial x_i, \\partial p_2/\\partial x_i)$ for each spatial direction $i$. Because reversing the order of the two quantities changes the sense of the vector, the vector behaves like a spin-2 quantity under 180-degree rotation, so it is combined with Stokes parameters $I_i = G_{1,i}^2 + G_{2,i}^2$, $Q_i = G_{1,i}^2 - G_{2,i}^2$, and $U_i = 2G_{1,i}G_{2,i}$, summed over directions and neighbors. The normalized $Q/I$ and $U/I$ encode the correlation angle and the correlation degree; the paper notes this is equivalent to extracting eigenvalues and eigenvectors of the local correlation matrix, with the correlation degree equal to $l_{\\max}/(l_{\\max}+l_{\\min})$. This construction is what turns local correlation information into a spatial map rather than a single global number.","core_discovery":"The paper's central claim is that locally correlated variations carry the information that global correlations destroy, and that this information becomes visible when adjacent gradient pairs are summed as spin-2 objects. For two measured fields $p_1$ and $p_2$, each location contributes gradient vectors $(G_1,G_2)$; the Stokes parameters $I$, $Q$, and $U$ are accumulated over spatial directions, and the normalized maps $Q/I$ and $U/I$ define a correlation angle $\\theta = \\tfrac{1}{2}\\arctan(U/Q)$ and a correlation degree $p = \\sqrt{(Q/I)^2+(U/I)^2}$. Perfectly correlated fields give $p \\approx 1$, uncorrelated fields $p \\approx 0$, and the angle tells how the two fields vary relative to each other. The paper demonstrates on temperature-precipitation data, MHD turbulence, and reaction-diffusion patterns that these maps reveal large coherent regions with a uniform correlation type, allowing the system to be divided into physically distinct subregions.","pith_inferences":["The gradient scale used to compute the derivatives is a free parameter; a natural next test is to measure how patch boundaries shift as that scale changes, since the paper does not report such sensitivity.","For more than two quantities, the author's suggested route through the local correlation matrix could yield maps of local principal-component structure, a multivariate extension not demonstrated in the paper.","Because the Stokes sum weights directions equally, anisotropic data might be better served by a directional variant that records the dominant gradient orientation, which the paper does not explore."],"forward_implications":["A dataset that appears uncorrelated globally can be broken into subregions, each with a well-defined local correlation type, so regime boundaries become visible without manually drawing boxes.","The correlation degree provides a spatially resolved measure of how strongly two fields regulate each other, with $p$ near 1 indicating mutual dependence and $p$ near 0 indicating local independence.","Type 1 (correlated) and Type 2 (stiff) regularities, distinguished by whether variations follow a relation $\\delta p_1 = k\\,\\delta p_2$ or are dominated by a single field, give a language for describing emergent simplicity versus dynamical detachment.","Because the method works for any two measured fields on a grid, it extends naturally to climate science, plasma physics, image analysis, and any other field with spatial or temporal data."],"supporting_citations":[{"why":"Supplies the climate-zone classification (semi-arid versus humid continental) that the weather-data application is compared against.","marker":"[1]"},{"why":"Provides the compressive MHD turbulence simulation data used to demonstrate the separation into magnetically and kinematically dominated regions.","marker":"[2]"},{"why":"Is the turbulence simulation catalogue from which the MHD data are drawn, grounding the second application.","marker":"[4]"},{"why":"Defines the reaction-diffusion system that generates the Turing pattern used as the third application.","marker":"[5]"}],"fun_headline_variants":["Spin-2 mapping splits data into coherent correlation regions","Adjacent gradients map correlations, expose hidden structure","Stokes parameters turn gradient pairs into correlation maps","Correlation maps reveal coherent regimes in complex data","Pixel-wise correlation maps unmask hidden regularities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The output map depends on the spatial scale over which the gradients in Equations (2)–(6) are computed and summed, and the paper neither fixes that scale nor tests whether the detected regions are robust to it.","fun_headline_variants_meta":{"raw":{"variants":["Spin-2 mapping splits data into coherent correlation regions","Adjacent gradients map correlations, expose hidden structure","Stokes parameters turn gradient pairs into correlation maps","Correlation maps reveal coherent regimes in complex data","Pixel-wise correlation maps unmask hidden regularities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000611,"raw_usage":{"total_tokens":2853,"prompt_tokens":967,"completion_tokens":1886,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":1814}},"tokens_in":583,"tokens_out":1886,"duration_ms":13990,"temperature":1.0,"reasoning_tokens":1814,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:13:13.882237+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the North America temperature-precipitation map and recompute the correlation-angle map with the gradient stencil spacing doubled, then halved: if the boundaries between negative-, positive-, and weakly-correlated regions shift by more than the stated climate-zone uncertainty, the claim that the map identifies physical regimes is not robust.","supporting_citations":[{"cited_title":"The Two States of Star Forming Clouds","cited_arxiv_id":"1202.2594","evidence_quote":"Is the turbulence simulation catalogue from which the MHD data are drawn, grounding the second application."}],"review_version":1}