Publications | Facundo Sapienza

2024

A Virtual Solar Wind Monitor for Mars with Uncertainty Quantification using Gaussian Processes

AR Azari, E Abrahams, F Sapienza, and 8 more authors

arXiv preprint arXiv:2402.01932, 2024

@article{azari2024virtual,
  title = {A Virtual Solar Wind Monitor for Mars with Uncertainty Quantification using Gaussian Processes},
  author = {Azari, AR and Abrahams, E and Sapienza, F and Halekas, J and Biersteker, J and Mitchell, DL and P{\'e}rez, F and Marquette, M and Rutala, MJ and Bowers, CF and others},
  journal = {arXiv preprint arXiv:2402.01932},
  year = {2024},
}

2023

Magnetic Field Draping in Induced Magnetospheres: Evidence From the MAVEN Mission to Mars

A. R. Azari, E. Abrahams, F. Sapienza, and 8 more authors

Journal of Geophysical Research: Space Physics, 2023

Abs Bib HTML PDF

Abstract The Mars Atmosphere and Volatile EvolutioN (MAVEN) mission has been orbiting Mars since 2014 and now has over 10,000 orbits which we use to characterize Mars’ dynamic space environment. Through global field line tracing with MAVEN magnetic field data we find an altitude dependent draping morphology that differs from expectations of induced magnetospheres in the vertical ( Mars Sun-state, MSO) direction. We quantify this difference from the classical picture of induced magnetospheres with a Bayesian multiple linear regression model to predict the draped field as a function of the upstream interplanetary magnetic field (IMF), remanent crustal fields, and a previously underestimated induced effect. From our model we conclude that unexpected twists in high altitude dayside draping (>800 km) are a result of the IMF component in the MSO direction. We propose that this is a natural outcome of current theories of induced magnetospheres but has been underestimated due to approximations of the IMF as solely directed. We additionally estimate that distortions in low altitude (<800 km) dayside draping along are directly related to remanent crustal fields. We show dayside draping traces down tail and previously reported inner magnetotail twists are likely caused by the crustal field of Mars, while the outer tail morphology is governed by an induced response to the IMF direction. We conclude with an updated understanding of induced magnetospheres which details dayside draping for multiple directions of the incoming IMF and discuss the repercussions of this draping for magnetotail morphology.
@article{azari2023magnetic, author = {Azari, A. R. and Abrahams, E. and Sapienza, F. and Mitchell, D. L. and Biersteker, J. and Xu, S. and Bowers, C. and Pérez, F. and DiBraccio, G. A. and Dong, Y. and Curry, S.}, title = {Magnetic Field Draping in Induced Magnetospheres: Evidence From the MAVEN Mission to Mars}, journal = {Journal of Geophysical Research: Space Physics}, volume = {128}, number = {11}, pages = {e2023JA031546}, keywords = {induced magnetospheres, hybrid magnetospheres, planetary science, Bayesian methods, data visualization, Gaussian process regression}, doi = {https://doi.org/10.1029/2023JA031546}, url = {https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/2023JA031546}, eprint = {https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023JA031546}, year = {2023}, }
Quantitative Analysis of Paleomagnetic Sampling Strategies

F. Sapienza, L. C. Gallo, Y. Zhang, and 3 more authors

Journal of Geophysical Research: Solid Earth, 2023

Abs Bib HTML PDF Code

Abstract Sampling strategies used in paleomagnetic studies play a crucial role in dictating the accuracy of our estimates of properties of the ancient geomagnetic field. However, there has been little quantitative analysis of optimal paleomagnetic sampling strategies and the community has instead defaulted to traditional practices that vary between laboratories. In this paper, we quantitatively evaluate the accuracy of alternative paleomagnetic sampling strategies through numerical experiments and an associated analytical framework. Our findings demonstrate a strong correspondence between the accuracy of an estimated paleopole position and the number of sites or independent readings of the time-varying paleomagnetic field, whereas larger numbers of in-site samples have a dwindling effect. This remains true even when a large proportion of the sample directions are spurious. This approach can be readily achieved in sedimentary sequences by distributing samples stratigraphically, considering each sample as an individual site. However, where the number of potential independent sites is inherently limited the collection of additional in-site samples can improve the accuracy of the paleopole estimate (although with diminishing returns with increasing samples per site). Where an estimate of the magnitude of paleosecular variation is sought, multiple in-site samples should be taken, but the optimal number is dependent on the expected fraction of outliers. The use of filters based on angular distance helps the accuracy of paleopole estimation, but leads to inaccurate estimates of paleosecular variation. We provide both analytical formulas and a series of interactive Jupyter notebooks allowing optimal sampling strategies to be developed from user-informed expectations.
@article{sapienza2023quantitative, author = {Sapienza, F. and Gallo, L. C. and Zhang, Y. and Vaes, B. and Domeier, M. and Swanson-Hysell, N. L.}, title = {Quantitative Analysis of Paleomagnetic Sampling Strategies}, journal = {Journal of Geophysical Research: Solid Earth}, volume = {128}, number = {11}, pages = {e2023JB027211}, keywords = {paleomagnetism, paleopole estimation, secular variation, error quantification}, doi = {https://doi.org/10.1029/2023JB027211}, url = {https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/2023JB027211}, eprint = {https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB027211}, year = {2023}, }
Embracing Uncertainty to Resolve Polar Wander: A Case Study of Cenozoic North America

L. C. Gallo, M. Domeier, F. Sapienza, and 11 more authors

Geophysical Research Letters, 2023

Abs Bib HTML PDF

Abstract Our understanding of Earth’s paleogeography relies heavily on paleomagnetic apparent polar wander paths (APWPs), which represent the time-dependent position of Earth’s spin axis relative to a given block of lithosphere. However, conventional approaches to APWP construction have significant limitations. First, the paleomagnetic record contains substantial noise that is not integrated into APWPs. Second, parametric assumptions are adopted to represent spatial and temporal uncertainties even where the underlying data do not conform to the assumed distributions. The consequences of these limitations remain largely unknown. Here, we address these challenges with a bottom-up Monte Carlo uncertainty propagation scheme that operates on site-level paleomagnetic data. To demonstrate our methodology, we present an extensive compilation of site-level Cenozoic paleomagnetic data from North America, which we use to generate a high-resolution APWP. Our results demonstrate that even in the presence of substantial noise, polar wandering can be assessed with unprecedented temporal and spatial resolution.
@article{gallo2023embracing, author = {Gallo, L. C. and Domeier, M. and Sapienza, F. and Swanson-Hysell, N. L. and Vaes, B. and Zhang, Y. and Arnould, M. and Eyster, A. and Gürer, D. and Király, Á. and Robert, B. and Rolf, T. and Shephard, G. and van der Boon, A.}, title = {Embracing Uncertainty to Resolve Polar Wander: A Case Study of Cenozoic North America}, journal = {Geophysical Research Letters}, volume = {50}, number = {11}, pages = {e2023GL103436}, keywords = {paleomagnetism, polar wander, uncertainty quantification}, doi = {https://doi.org/10.1029/2023GL103436}, eprint = {https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023GL103436}, year = {2023}, }
Universal Differential Equations for glacier ice flow modelling

J. Bolibar, F. Sapienza, F. Maussion, and 3 more authors

Geoscientific Model Development, 2023

Abs Bib HTML PDF Code

Geoscientific models are facing increasing challenges to exploit growing datasets coming from remote sensing. Universal Differential Equations (UDEs), aided by differentiable programming, provide a new scientific modelling paradigm enabling both complex functional inversions to potentially discover new physical laws and data assimilation from heterogeneous and sparse observations. We demonstrate an application of UDEs as a proof of concept to learn the creep component of ice flow, i.e. a nonlinear diffusivity differential equation, of a glacier evolution model. By combining a mechanistic model based on a 2D Shallow Ice Approximation Partial Differential Equation with an embedded neural network, i.e. a UDE, we can learn parts of an equation as nonlinear functions that then can be translated into mathematical expressions. We implemented this modelling framework as ODINN.jl, a package in the Julia programming language, providing high performance, source-to-source automatic differentiation (AD) and seamless integration with tools and global datasets from the Open Global Glacier Model in Python. We demonstrate this concept for 17 different glaciers around the world, for which we successfully recover a prescribed artificial law describing ice creep variability by solving ca 500,000 Ordinary Differential Equations in parallel. Furthermore, we investigate which are the best tools in the scientific machine learning ecosystem in Julia to differentiate and optimize large nonlinear diffusivity UDEs. This study represents a proof of concept for a new modelling framework aiming at discovering empirical laws for large-scale glacier processes, such as the variability of ice creep and basal sliding for ice flow, and new hybrid surface mass balance models.
@article{bolibar2023universal, title = {Universal Differential Equations for glacier ice flow modelling}, author = {Bolibar, J. and Sapienza, F. and Maussion, F. and Lguensat, R. and Wouters, B. and P\'erez, F.}, journal = {Geoscientific Model Development}, volume = {16}, pages = {6671--6687}, year = {2023}, number = {22}, doi = {10.5194/gmd-16-6671-2023}, }
Choosing the parameter of the Fermat distance: navigating geometry and noise

Frédéric Chazal, Laure Ferraris, Pablo Groisman, and 3 more authors

2023

Abs Bib PDF

The Fermat distance has been recently established as a useful tool for machine learning tasks when a natural distance is not directly available to the practitioner or to improve the results given by Euclidean distances by exploding the geometrical and statistical properties of the dataset. This distance depends on a parameter α that greatly impacts the performance of subsequent tasks. Ideally, the value of α should be large enough to navigate the geometric intricacies inherent to the problem. At the same, it should remain restrained enough to sidestep any deleterious ramifications stemming from noise during the process of distance estimation. We study both theoretically and through simulations how to select this parameter.
@misc{chazal2023choosing, title = {Choosing the parameter of the Fermat distance: navigating geometry and noise}, author = {Chazal, Frédéric and Ferraris, Laure and Groisman, Pablo and Jonckheere, Matthieu and Pascal, Frédéric and Sapienza, Facundo}, year = {2023}, eprint = {2311.18663}, archiveprefix = {arXiv}, primaryclass = {stat.ML}, }

2022

Heat engines with single-shot deterministic work extraction

Federico Cerisola, Facundo Sapienza, and Augusto J Roncaglia

Physical Review E, 2022

Abs Bib HTML

We introduce heat engines working in the nanoregime that allow one to extract a finite amount of deterministic work. Using the resource theory approach to themodynamics, we show that the efficiency of these cycles is strictly smaller than Carnot’s, and we associate this difference with a fundamental irreversibility that is present in single-shot transformations. When fluctuations in the extracted work are allowed there is a trade-off between their size and the efficiency. As the size of fluctuations increases so does the efficiency and optimal efficiency is attained for unbounded fluctuations, while a certain amount of deterministic work is drawn from the cycle. Finally, we show that when the working medium is composed of many particles, by creating an amount of correlations between the subsystems that scale logarithmically with their number, Carnot’s efficiency can also be approached in the asymptotic limit along with deterministic work extraction.
@article{cerisola2022heat, title = {Heat engines with single-shot deterministic work extraction}, author = {Cerisola, Federico and Sapienza, Facundo and Roncaglia, Augusto J}, journal = {Physical Review E}, volume = {106}, number = {3}, pages = {034135}, year = {2022}, publisher = {APS}, }
An optimization method for paleomagnetic Euler pole analysis

Leandro C Gallo, Facundo Sapienza, and Mathew Domeier

Computers & Geosciences, 2022

Abs Bib HTML

Owing to the axial symmetry of the Earth’s magnetic field, paleomagnetic data only directly record the latitudinal and azimuthal positions of crustal blocks in the past, and paleolongitude cannot be constrained. An ability to overcome this obstacle is thus of fundamental importance to paleogeographic reconstruction. Paleomagnetic Euler pole (PEP) analysis presents a unique means to recover such information, but prior implementations of the PEP method have incorporated subjective decisions into its execution, undercutting its fidelity and rigor. Here we present an optimization approach to PEP analysis that addresses some of these deficiencies—namely the objective identification of change-points and small-circle arcs that together approximate an apparent polar wander path. We elaborate on our novel methodology and conduct some experiments with synthetic data to demonstrate its performance. We furthermore present implementations of our methods both as adaptable, stand-alone scripts in Python and as a streamlined interactive workflow that can be operated through a web browser.
@article{gallo2022optimization, title = {An optimization method for paleomagnetic Euler pole analysis}, author = {Gallo, Leandro C and Sapienza, Facundo and Domeier, Mathew}, journal = {Computers \& Geosciences}, volume = {166}, pages = {105150}, year = {2022}, publisher = {Elsevier}, }
Efficient adjustment sets in causal graphical models with hidden variables

Ezequiel Smucler, Facundo Sapienza, and Andrea Rotnitzky

Biometrika, 2022

Abs Bib HTML Code

We study the selection of adjustment sets for estimating the interventional mean under a point exposure dynamic treatment regime, that is, a treatment rule that depends on the subject’s covariates. We assume a nonparametric causal graphical model with, possibly, hidden variables and at least one adjustment set comprised of observable variables. We provide the definition of a valid adjustment set for a point exposure dynamic treatment regime, which generalizes the existing definition for a static intervention. We show that there exists an adjustment set, referred to as optimal minimal, that yields the nonparametric estimator of the interventional mean with the smallest asymptotic variance among those that are based on observable minimal adjustment sets. An observable minimal adjustment set is a valid adjustment set such that all its variables are observable and the removal of any of its variables destroys its validity. We provide similar optimality results for the class of observable minimum adjustment sets, that is, valid observable adjustment sets of minimum cardinality among the observable adjustment sets. Moreover, we show that if either no variables are hidden or if all the observable variables are ancestors of either treatment, outcome or the variables that are used to decide treatment, a globally optimal adjustment set exists. We provide polynomial-time algorithms to compute the globally optimal, optimal minimal and optimal minimum adjustment sets. Because static interventions can be viewed as a special case of dynamic regimes, all our results also apply for static interventions.
@article{smucler2022efficient, title = {Efficient adjustment sets in causal graphical models with hidden variables}, author = {Smucler, Ezequiel and Sapienza, Facundo and Rotnitzky, Andrea}, journal = {Biometrika}, volume = {109}, number = {1}, pages = {49--65}, year = {2022}, publisher = {Oxford University Press}, }
Nonhomogeneous Euclidean first-passage percolation and distance learning

Pablo Groisman, Matthieu Jonckheere, and Facundo Sapienza

Bernoulli, 2022

Abs Bib HTML Code

Consider an i.i.d. sample from an unknown density function supported on an unknown manifold embedded in a high dimensional Euclidean space. We tackle the problem of learning a distance between points, able to capture both the geometry of the manifold and the underlying density. We define such a sample distance and prove the convergence, as the sample size goes to infinity, to a macroscopic one that we call Fermat distance as it minimizes a path functional, resembling Fermat principle in optics. The proof boils down to the study of geodesics in Euclidean first-passage percolation for nonhomogeneous Poisson point processes.
@article{groisman2022nonhomogeneous, title = {Nonhomogeneous Euclidean first-passage percolation and distance learning}, author = {Groisman, Pablo and Jonckheere, Matthieu and Sapienza, Facundo}, journal = {Bernoulli}, volume = {28}, number = {1}, pages = {255--276}, year = {2022}, publisher = {Bernoulli Society for Mathematical Statistics and Probability}, }

2019

Correlations as a resource in quantum thermodynamics

Facundo Sapienza, Federico Cerisola, and Augusto J Roncaglia

Nature communications, 2019

Abs Bib HTML PDF

The presence of correlations in physical systems can be a valuable resource for many quantum information tasks. They are also relevant in thermodynamic transformations, and their creation is usually associated to some energetic cost. In this work, we study the role of correlations in the thermodynamic process of state formation in the single-shot regime, and find that correlations can also be viewed as a resource. First, we show that the energetic cost of creating multiple copies of a given state can be reduced by allowing correlations in the final state. We obtain the minimum cost for every finite number of subsystems, and then we show that this feature is not restricted to the case of copies. More generally, we demonstrate that in the asymptotic limit, by allowing a logarithmic amount of correlations, we can recover standard results where the free energy quantifies this minimum cost.
@article{sapienza2019correlations, title = {Correlations as a resource in quantum thermodynamics}, author = {Sapienza, Facundo and Cerisola, Federico and Roncaglia, Augusto J}, journal = {Nature communications}, volume = {10}, number = {1}, pages = {2492}, year = {2019}, publisher = {Nature Publishing Group UK London}, }

2018

Anhysteretic remanent magnetization: model of grain size distribution of spherical magnetite grains

Carlos A Vasquez, Facundo Sapienza, Agustı́n Somacal, and 1 more author

Studia Geophysica et Geodaetica, 2018

Abs Bib

A phenomenological model based on a linear relationship between the magnetic coercivity field and the reciprocal of the grain diameter is applied to explain the anhysteretic remanent magnetization (ARM) imparted to artificial samples with different concentrations of a very well characterized magnetite powder. By analyses of scanning electron microscopy images, the spherically shaped single domain synthetic magnetite is found to follow a lognormal grain size distribution with 86 nm of mean diameter. The proposed model, fitted to ARM measurements up to a peak alternating field of 100 mT, yields a very good agreement. The coercivity behaviour predicted by micromagnetism theory disagrees with the experimental results of this work. A likely explanation for the discrepancy is that the magnetite particles, which consist of a mixture of grains in coherent rotation and curling modes, produce similar observations as domain processes.
@article{vasquez2018anhysteretic, title = {Anhysteretic remanent magnetization: model of grain size distribution of spherical magnetite grains}, author = {Vasquez, Carlos A and Sapienza, Facundo and Somacal, Agust{\'\i}n and Fazzito, Sabrina Y}, journal = {Studia Geophysica et Geodaetica}, volume = {62}, pages = {339--351}, year = {2018}, publisher = {Springer}, }
Weighted Geodesic Distance Following Fermat’s Principle

Facundo Sapienza, Pablo Groisman, and Matthieu Jonckheere

2018

Abs Bib PDF

We propose a density-based estimator for weighted geodesic distances suitable for data lying on a manifold of lower dimension than ambient space and sampled from a possibly nonuniform distribution. After discussing its properties and implementation, we evaluate its performance as a tool for clustering tasks. A discussion on the consistency of the estimator is also given.
@article{sapienza2018weighted, title = {Weighted Geodesic Distance Following Fermat's Principle}, author = {Sapienza, Facundo and Groisman, Pablo and Jonckheere, Matthieu}, year = {2018}, }
[thesis] Distancia de Fermat y geodésicas en percolación euclidea: teoría y aplicaciones en Machine Learning

Facundo Sapienza

Aug 2018

PDF

2017

[thesis] Teoría de recursos de la termodinámica cuántica : correlaciones en los procesos de formación de recursos

Facundo Sapienza

Dec 2017

PDF