Climate change triggers a wide range of hydrometeorological, glaciological, and geophysical processes that span across vast spatiotemporal scales. With the advances in technology and analytics, a multitude of remote sensing (RS), geodetic, and in situ instruments have been developed to effectively monitor and help comprehend Earth’s system, including its climate variability and the recent anomalies associated with global warming. A huge volume of data is generated by recording these observations, resulting in the need for novel methods to handle and interpret such big datasets. Managing this enormous amount of data extends beyond current computer storage considerations; it also encompasses the complexities of processing, modeling, and analyzing. Big datasets present unique characteristics that set them apart from smaller datasets, thereby posing challenges to traditional approaches. Moreover, computational time plays a crucial role, especially in the context of geohazard warning and response systems, which necessitate low latency requirements.

The choice of an appropriate metric is mandatory to perform deformation analysis between two point clouds (PC)—the distance has to be trustworthy and, simultaneously, robust against measurement noise, which may be correlated and heteroscedastic. The Hausdorff distance (HD) or its averaged derivation (AHD) are widely used to compute local distances between two PC and are implemented in nearly all commercial software. Unfortunately, they are affected by measurement noise, particularly when correlations are present. In this contribution, we focus on terrestrial laser scanner (TLS) observations and assess the impact of neglecting correlations on the distance computation when a mathematical approximation is performed. The results of the simulations are extended to real observations from a bridge under load. Highly accurate laser tracker (LT) measurements were available for this experiment: they allow the comparison of the HD and AHD between two raw PC or between their mathematical approximations regarding reference values. Based on these results, we determine which distance is better suited in the case of heteroscedastic and correlated TLS observations for local deformation analysis. Finally, we set up a novel bootstrap testing procedure for this distance when the PC are approximated with B-spline surfaces.

The global navigation satellite system (GNSS) daily position time series are often described as the sum of stochastic processes and geophysical signals which allow to study global and local geodynamical effects such as plate tectonics, earthquakes, or ground water variations. In this work, we propose to extend the Generalized Method of Wavelet Moments (GMWM) to estimate the parameters of linear models with correlated residuals. This statistical inferential framework is applied to GNSS daily position time-series data to jointly estimate functional (geophysical) as well as stochastic noise models. Our method is called GMWMX, with X standing for eXogenous variables: it is semi-parametric, computationally efficient and scalable. Unlike standard methods such as the widely used maximum likelihood estimator (MLE), our methodology offers statistical guarantees, such as consistency and asymptotic normality, without relying on strong parametric assumptions. At the Gaussian model, our results (theoretical and obtained in simulations) show that the estimated parameters are similar to the ones obtained with the MLE. The computational performances of our approach have important practical implications. Indeed, the estimation of the parameters of large networks of thousands of GNSS stations (some of them being recorded over several decades) quickly becomes computationally prohibitive. Compared to standard likelihood-based methods, the GMWMX has a considerably reduced algorithmic complexity of order $$\mathcal {O}\{\log (n) n\}$$O{log(n)n} for a time series of lengthn. Thus, the GMWMX appears to provide a reduction in processing time of a factor of 10–1000 compared to likelihood-based methods depending on the considered stochastic model, the length of the time series and the amount of missing data. As a consequence, the proposed method allows the estimation of large-scale problems within minutes on a standard computer. We validate the performances of our method via Monte Carlo simulations by generating GNSS daily position time series with missing observations and we consider composite stochastic noise models including processes presenting long-range dependence such as power law or Matérn processes. The advantages of our method are also illustrated using real time series from GNSS stations located in the Eastern part of the USA.

For a trustworthy least-squares (LS) solution, a good description of the stochastic properties of the measurements is indispensable. For a terrestrial laser scanner (TLS), the range variance can be described by a power law function with respect to the intensity of the reflected signal. The power and scaling factors depend on the laser scanner under consideration, and could be accurately determined by means of calibrations in 1d mode or residual analysis of LS adjustment. However, such procedures complicate significantly the use of empirical intensity models (IM). The extent to which a point-wise weighting is suitable when the derived variance covariance matrix (VCM) is further used in a LS adjustment remains moreover questionable. Thanks to closed loop simulations, where both the true geometry and stochastic model are under control, we investigate how variations of the parameters of the IM affect the results of a LS adjustment. As a case study, we consider the determination of the Cartesian coordinates of the control points (CP) from a B-splines curve. We show that a constant variance can be assessed to all the points of an object having homogeneous properties, without affecting the a posteriori variance factor or the loss of efficiency of the LS solution. The results from a real case scenario highlight that the conclusions of the simulations stay valid even for more challenging geometries. A procedure to determine the range variance is proposed to simplify the computation of the VCM.