{"id":"620bfcc6-be59-4adc-b6a9-a425a4500c9a","arxiv_id":"2606.19065","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":1,"one_line_summary":"A nonparametric covariance estimator for fragmented partially observed 2D functional data is proposed via matrix completion with Laplacian regularization, without stationarity assumptions.","lead":"The paper introduces a Laplacian-regularized matrix completion approach to estimate covariance from 2D functional data observed in a fragmented missingness regime. This targets improved analysis of ground displacement measurements from interferometry for geophysical monitoring.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Laplacian regularization may bias recovery when true covariance lacks the implied smoothness on the domain graph","rationale":"The reader's weakest_assumption correctly isolates the critical unverified step. The full text presumably supplies simulation details, but the load-bearing risk remains whether those details include unfavorable covariance structures; the proposed check directly probes that gap.","tokens_in":1668,"tokens_out":328,"duration_ms":16072,"concrete_test":"Re-run the simulation protocol on a 2D grid using a Matérn covariance with smoothness parameter ν=0.5 (rough sample paths) under the same fragmented missingness pattern; if the reported estimation error increases by more than 30% relative to the smoother cases shown in the paper, the recovery assumption does not hold generally.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The method casts covariance estimation as matrix completion regularized by the graph Laplacian on the 2D domain. For the estimator to recover the true spatial dependence without stationarity or isotropy assumptions, the penalty must correctly impute entries under the fragmented missingness regime. This implicitly requires that the true covariance operator is well-approximated by low-frequency eigenfunctions of the Laplacian; otherwise the regularizer introduces systematic bias. The abstract states that simulations achieve low error across a range of structures, yet provides no indication that the tested covariances include cases deliberately misaligned with the Laplacian spectrum (e.g., rough or localized dependence). Without such stress tests or accompanying theory, the nonparametric claim rests on an unverified compatibility between the regularizer and the data-generating process.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper introduces the fragmented regime of systematic partial observation for 2D functional data arising from InSAR measurements. It formulates covariance estimation as a matrix-completion problem regularized by the graph Laplacian on the spatial domain, claiming the resulting estimator is nonparametric and free of stationarity or isotropy assumptions. Simulations are reported to yield consistently low error across covariance structures, and the method is applied to Phlegraean Fields displacement data to recover spatial dependence patterns.","tokens_in":1844,"tokens_out":364,"duration_ms":15832,"significance":"If the Laplacian regularizer recovers the true second-order structure without systematic bias under the fragmented regime, the approach would supply a practical, assumption-light tool for covariance estimation in geophysical remote-sensing applications where large contiguous blocks of data are missing.","major_comments":[{"comment":"Simulation study: the claim of low error 'across a range of covariance structures' does not include deliberate stress tests against covariances whose eigenstructure is misaligned with the low-frequency modes of the domain Laplacian (e.g., rough or spatially localized dependence). Without such cases the nonparametric claim remains unverified and the bias concern raised by the regularizer is unaddressed.","section":"Simulation study"},{"comment":"Estimator definition: the regularization parameter appears as a free hyper-parameter whose selection procedure is not shown to preserve the nonparametric character of the estimator; any data-driven choice must be demonstrated not to re-introduce implicit parametric assumptions.","section":"Methods / estimator definition"}],"minor_comments":[{"comment":"The precise mathematical definition of the 'fragmented regime' should be stated explicitly with a diagram or equation before the estimator is introduced.","section":"Introduction / setup"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their detailed and constructive report. We address each major comment below and outline the revisions we will make to the manuscript.","responses":[{"response":"We agree that the existing simulations, while spanning multiple covariance structures, do not contain explicit stress tests for eigenstructures deliberately misaligned with the low-frequency modes of the Laplacian (e.g., rough or spatially localized dependence). To verify the nonparametric claim and directly address potential bias introduced by the regularizer, we will add new simulation scenarios with such covariance structures in the revised manuscript. These additions will include quantitative error comparisons under the fragmented observation regime.","revision_made":"yes","referee_comment":"[Simulation study] Simulation study: the claim of low error 'across a range of covariance structures' does not include deliberate stress tests against covariances whose eigenstructure is misaligned with the low-frequency modes of the domain Laplacian (e.g., rough or spatially localized dependence). Without such cases the nonparametric claim remains unverified and the bias concern raised by the regularizer is unaddressed."},{"response":"The regularization parameter is treated as a hyper-parameter whose value is chosen via cross-validation on the observed entries. We acknowledge that the manuscript does not explicitly demonstrate that this data-driven procedure preserves the nonparametric character of the estimator. In the revision we will expand the methods section with a clearer description of the selection algorithm together with additional empirical checks (and, if feasible, a brief theoretical argument) showing that the procedure does not implicitly re-introduce stationarity or other parametric assumptions.","revision_made":"yes","referee_comment":"[Methods / estimator definition] Estimator definition: the regularization parameter appears as a free hyper-parameter whose selection procedure is not shown to preserve the nonparametric character of the estimator; any data-driven choice must be demonstrated not to re-introduce implicit parametric assumptions."}],"tokens_in":1281,"tokens_out":404,"duration_ms":19019,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main takeaway is a practical method for covariance estimation when observations are missing in large systematic blocks across replicates, framed as matrix completion with a graph Laplacian penalty on the 2D domain.\n\nWhat stands out as new is the explicit definition of the fragmented regime together with the regularized completion setup for nonparametric covariance in spatial functional data. The work does well by targeting a real problem in interferometric SAR data, where standard approaches fail due to the structured missingness. Simulations report consistently low error across several covariance structures, and the Phlegraean Fields application recovers spatial dependence patterns that align with known geophysical features.\n\nThe soft spot is the implicit smoothness assumption built into the Laplacian penalty. The estimator recovers the true structure only if that structure is well approximated by the low-frequency eigenfunctions of the domain graph; rougher or localized dependence would likely produce systematic bias. The abstract gives no sign that the simulations included deliberate mismatches with the Laplacian spectrum, so the robustness claim rests on untested compatibility between the regularizer and the data-generating process.\n\nThis is for statisticians or applied researchers working with spatial functional data under structured missingness, especially in remote sensing or environmental monitoring. A reader facing similar incomplete spatial curves would find usable ideas and a concrete example.\n\nIt deserves a serious referee because the problem is well-motivated, the method is implementable, and the application provides external grounding even if additional stress tests would help.","headline":"The paper defines a fragmented missingness regime for 2D functional data and estimates covariance via Laplacian-regularized matrix completion without stationarity assumptions.","tokens_in":2320,"tokens_out":361,"would_cite":false,"duration_ms":18229,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Covariance from partially observed 2D interferometric data is recovered by Laplacian-regularized matrix completion without stationarity assumptions.","keywords":["covariance estimation","matrix completion","Laplacian regularization","functional data analysis","interferometric data","fragmented regime","spatial dependence","nonparametric estimation"],"falsifier":"Apply the estimator to synthetic data generated from a known covariance that cannot be well approximated by Laplacian regularization and check whether the recovered matrix deviates substantially from the generating covariance in Frobenius norm.","tokens_in":2589,"feed_emoji":"","tokens_out":593,"duration_ms":11631,"temperature":0.7,"pith_summary":"The paper treats ground displacement measurements from the Small BAseline Subset technique as partially observed functional data on a two-dimensional domain. It defines a fragmented regime where large portions of the curves are systematically missing across replicates. The authors formulate covariance estimation as a matrix completion problem with Laplacian regularization to recover spatial dependence. The resulting nonparametric estimator avoids any stationarity or isotropy requirements. Simulations across varied covariance structures show consistently low error, and the method applied to Phlegraean Fields data recovers interpretable dependence patterns useful for risk monitoring.","feed_headline":"Laplacian matrix completion recovers covariance from fragmented 2D data","feed_subtitle":"The nonparametric estimator avoids stationarity assumptions and yields low error on simulated and real interferometric displacement measurem","key_machinery":"Laplacian-regularized matrix completion applied to the covariance matrix of fragmented functional data on a 2D domain.","core_discovery":"The paper establishes a nonparametric covariance estimator for functional data observed under the fragmented regime of partial observation. The estimator is obtained by solving a matrix completion problem whose objective includes a Laplacian regularization term; this term is shown to recover the true spatial dependence structure of the data-generating process on the two-dimensional domain.","pith_inferences":["The same matrix-completion framing could be tested on other remote-sensing modalities that exhibit block-wise missingness.","Performance comparisons against kriging or other spatial covariance estimators would clarify when the Laplacian term adds value.","If the fragmented regime accurately describes additional data sources, the estimator could be reused directly in those domains."],"forward_implications":["The estimator produces low error for a wide range of covariance structures in simulation studies.","Application to real interferometric data yields spatial dependence patterns that align with geophysical expectations.","The approach supports downstream analysis of geophysical processes without requiring stationarity or isotropy.","The method extends covariance estimation tools to functional data with systematic large-scale missingness."],"fun_headline_variants":["Laplacian-regularized matrix completion for 2D covariance estimation","Covariance recovery from fragmented observations using Laplacian regularization","Nonparametric covariance estimator with Laplacian matrix completion","Fragmented 2D covariance estimated via Laplacian matrix completion"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The true spatial dependence structure is recoverable by Laplacian regularization when the missingness follows the fragmented regime of systematic partial observation.","fun_headline_variants_meta":{"raw":{"variants":["Laplacian-regularized matrix completion for 2D covariance estimation","Covariance recovery from fragmented observations using Laplacian regularization","Nonparametric covariance estimator with Laplacian matrix completion","Fragmented 2D covariance estimated via Laplacian matrix completion"]},"model":"grok-4.3","cost_usd":0.003199,"raw_usage":{"total_tokens":1614,"prompt_tokens":615,"num_sources_used":0,"completion_tokens":62,"cost_in_usd_ticks":31990500,"prompt_tokens_details":{"text_tokens":615,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":937,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":615,"tokens_out":62,"duration_ms":7167,"temperature":1.0,"reasoning_tokens":937,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-26T19:54:20.758576+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Apply the estimator to synthetic data generated from a known covariance that cannot be well approximated by Laplacian regularization and check whether the recovered matrix deviates substantially from the generating covariance in Frobenius norm.","supporting_citations":[],"review_version":1}