scholarly journals Correlation dispersion as a measure to better estimate uncertainty of remotely sensed glacier displacements

2021 ◽  
Author(s):  
Bas Altena ◽  
Andreas Kääb ◽  
Bert Wouters

Abstract. In recent years a vast amount of glacier surface velocity data from satellite imagery has emerged based on correlation between repeat images. Thereby, much emphasis has been put on fast processing of large data volumes. The metadata of such measurements are often highly simplified when the measurement precision is lumped into a single number for the whole dataset, although the error budget of image matching is in reality not isotropic and constant over the whole velocity field. The spread of the correlation peak of individual image offset measurements is dependent on the image content and the non-uniform flow of the ice. Precise dispersion estimates for each individual velocity measurement can be important for inversion of, for instance, rheology, ice thickness and/or bedrock friction. Errors in the velocity data can propagate into derived results in a complex and exaggerating way, making the outcomes very sensitive to velocity noise and errors. Here, we present a computationally fast method to estimate the matching precision of individual displacement measurements from repeat imaging data, focussing on satellite data. The approach is based upon Gaussian fitting directly on the correlation peak and is formulated as a linear least squares estimation, making its implementation into current pipelines straightforward. The methodology is demonstrated for Sermeq Kujalleq, Greenland, a glacier with regions of strong shear flow and with clearly oriented crevasses, and Malaspina Glacier, Alaska. Directionality within an image seems to be dominant factor influencing the correlation dispersion. In our cases these are crevasses and moraine bands, while a relation to differential flow, such as shear, is less pronounced.

2004 ◽  
Vol 39 ◽  
pp. 175-180 ◽  
Author(s):  
Veijo Allan Pohjola ◽  
Jim Hedfors ◽  
Per Holmlund

AbstractHow well can we estimate the incoming ice flux by calculating the ice flux through a well-defined cross-section? We test this by comparing calculated ice flux out from the small glacier Bonnevie-Svendsenbreen with the measured accumulation rate integrated over the well-defined catchment area in the Sivorgfjella plateau, Dronning Maud Land, Antarctica (74˚45’ S, 11˚10’ W). The ice flux is calculated using ice-dynamical properties from an ice temperature model and the distribution of forces calculated using a force-budget model. The input we use includes velocity data of the glacier surface, combined with ice-thickness measurements. The result is an accumulation rate on the Sivorgfjella plateau of 0.50 ± 0.05 mw.e. a–1. We find that this is similar to the accumulation rate recorded by ground-penetrating radar work in the area. We therefore find the balance-flow method, in combination with the force-budget technique and ice temperature modeling, to be a useful tool for studies of mass fluxes in a catchment area. The most important source of uncertainty in these calculations is the quality and the spatial distribution of the ice surface velocity data. The high accumulation rate shows the effect of orographic enhancement on accumulation in montane areas in Antarctica.


BMJ Open ◽  
2021 ◽  
Vol 11 (7) ◽  
pp. e047349
Author(s):  
Ewoud ter Avest ◽  
Barbara C van Munster ◽  
Raymond J van Wijk ◽  
Sanne Tent ◽  
Sanne Ter Horst ◽  
...  

PurposeResearch in acute care faces many challenges, including enrolment challenges, legal limitations in data sharing, limited funding and lack of singular ownership of the domain of acute care. To overcome these challenges, the Center of Acute Care of the University Medical Center Groningen in the Netherlands, has established a de novo data, image and biobank named ‘Acutelines’.ParticipantsClinical data, imaging data and biomaterials (ie, blood, urine, faeces, hair) are collected from patients presenting to the emergency department (ED) with a broad range of acute disease presentations. A deferred consent procedure (by proxy) is in place to allow collecting data and biomaterials prior to obtaining written consent. The digital infrastructure used ensures automated capturing of all bed-side monitoring data (ie, vital parameters, electrophysiological waveforms) and securely importing data from other sources, such as the electronic health records of the hospital, ambulance and general practitioner, municipal registration and pharmacy. Data are collected from all included participants during the first 72 hours of their hospitalisation, while follow-up data are collected at 3 months, 1 year, 2 years and 5 years after their ED visit.Findings to dateEnrolment of the first participant occurred on 1 September 2020. During the first month, 653 participants were screened for eligibility, of which 180 were approached as potential participants. In total, 151 (84%) provided consent for participation of which 89 participants fulfilled criteria for collection of biomaterials.Future plansThe main aim of Acutelines is to facilitate research in acute medicine by providing the framework for novel studies and issuing data, images and biomaterials for future research. The protocol will be extended by connecting with central registries to obtain long-term follow-up data, for which we already request permission from the participant.Trial registration numberNCT04615065.


2019 ◽  
Vol 75 (1) ◽  
pp. 14-24 ◽  
Author(s):  
Joseph A. M. Paddison

Diffuse scattering is a rich source of information about disorder in crystalline materials, which can be modelled using atomistic techniques such as Monte Carlo and molecular dynamics simulations. Modern X-ray and neutron scattering instruments can rapidly measure large volumes of diffuse-scattering data. Unfortunately, current algorithms for atomistic diffuse-scattering calculations are too slow to model large data sets completely, because the fast Fourier transform (FFT) algorithm has long been considered unsuitable for such calculations [Butler & Welberry (1992). J. Appl. Cryst. 25, 391–399]. Here, a new approach is presented for ultrafast calculation of atomistic diffuse-scattering patterns. It is shown that the FFT can actually be used to perform such calculations rapidly, and that a fast method based on sampling theory can be used to reduce high-frequency noise in the calculations. These algorithms are benchmarked using realistic examples of compositional, magnetic and displacive disorder. They accelerate the calculations by a factor of at least 102, making refinement of atomistic models to large diffuse-scattering volumes practical.


2018 ◽  
Vol 64 (248) ◽  
pp. 969-976 ◽  
Author(s):  
J. W. SANDERS ◽  
K. M. CUFFEY ◽  
K. R. MACGREGOR ◽  
J. L. KAVANAUGH ◽  
C. F. DOW

ABSTRACTFollowing pioneering work in Norway, cirque glaciers have widely been viewed as rigidly rotating bodies. This model is incorrect for basin-filling cirque glaciers, as we have demonstrated at West Washmawapta Glacier, a small glacier in the Canadian Rocky Mountains. Here we report observations at the same glacier that assess whether complex temporal variations of flow also occur. For parts of three summers, we measured daily displacements of the glacier surface. In one year, four short-duration speed-up events were recorded. Three of the events occurred during the intervals of warmest weather, when melt was most rapid; the fourth event occurred immediately following heavy rain. We interpret the speed-up events as manifestations of enhanced water inputs to the glacier bed and associated slip lubrication by increased water volumes and pressures. No further speed-ups occurred in the final month of the melt season, despite warm temperatures and several rainstorms; the dominant subglacial water system likely transformed from one of poorly connected cavities to one with an efficient channel network. The seasonal evolution of hydrology and flow resembles behaviors documented at other, larger temperate glaciers and indicates that analyses of cirque erosion cannot rely on simple assumptions about ice dynamics.


1986 ◽  
Vol 32 (110) ◽  
pp. 101-119 ◽  
Author(s):  
Almut Iken ◽  
Robert A. Bindschadler

AbstractDuring the snow-melt season of 1982, basal water pressure was recorded in 11 bore holes communicating with the subglacial drainage system. In most of these holes the water levels were at approximately the same depth (around 70 m below surface). The large variations of water pressure, such as diurnal variations, were usually similar at different locations and in phase. In two instances of exceptionally high water pressure, however, systematic phase shifts were observed; a wave of high pressure travelled down-glacier with a velocity of approximately 100 m/h.The glacier-surface velocity was measured at four lines of stakes several times daily. The velocity variations correlated with variations in subglacial water pressure. The functional relationship of water pressure and velocity suggests that fluctuating bed separation was responsible for the velocity variations. The empirical functional relationship is compared to that of sliding over a perfectly lubricated sinusoidal bed. On the basis of the measured velocity-pressure relationship, this model predicts a reasonable value of bed roughness but too high a sliding velocity and unstable sliding at too low a water pressure. The main reason for this disagreement is probably the neglect of friction from debris in the sliding model.The measured water pressure was considerably higher than that predicted by the theory of steady flow through straight cylindrical channels near the glacier bed. Possible reasons are considered. The very large disagreement between measured and predicted pressure suggests that no straight cylindrical channels may have existed.


1985 ◽  
Vol 31 (108) ◽  
pp. 115-119 ◽  
Author(s):  
L. A. Rasmussen

Abstract An algorithm is developed for adjusting glacier surface-velocity vectors, given on the nodes of a square grid, so that they obey a central-difference approximation of the continuity equation. Also required on the grid nodes are the glacier thickness, the ratio of the surface-velocity to the average velocity in the column, and the difference between the mass balance and the thickness change. All these other variables are assumed to be known exactly, and only the surface-velocity field is adjusted. The result is optimum in the sense that the magnitude of the adjustment is minimized. Either the relative or the absolute adjustment can be minimized, depending on how weights are specified. No restriction is placed on the shape of the solution region, and no boundary condition is required. The algorithm is not iterative. The algorithm first forms a parallel flow field that satisfies the continuity equation, and then uses a stream function to add a divergenceless field to it. The stream function that leads to the minimum velocity adjustment is obtained as four independent, interlacing solutions covering the solution region. For each of the four, a well-conditioned, sparse-matrix system of simultaneous linear equations is solved. A compact, sub-optimum, well-behaved iterative procedure is also developed for transforming part of the velocity adjustment into an adjustment of the thickness field.


2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Michele Allegra ◽  
Elena Facco ◽  
Francesco Denti ◽  
Alessandro Laio ◽  
Antonietta Mira

Abstract One of the founding paradigms of machine learning is that a small number of variables is often sufficient to describe high-dimensional data. The minimum number of variables required is called the intrinsic dimension (ID) of the data. Contrary to common intuition, there are cases where the ID varies within the same data set. This fact has been highlighted in technical discussions, but seldom exploited to analyze large data sets and obtain insight into their structure. Here we develop a robust approach to discriminate regions with different local IDs and segment the points accordingly. Our approach is computationally efficient and can be proficiently used even on large data sets. We find that many real-world data sets contain regions with widely heterogeneous dimensions. These regions host points differing in core properties: folded versus unfolded configurations in a protein molecular dynamics trajectory, active versus non-active regions in brain imaging data, and firms with different financial risk in company balance sheets. A simple topological feature, the local ID, is thus sufficient to achieve an unsupervised segmentation of high-dimensional data, complementary to the one given by clustering algorithms.


2018 ◽  
Vol 119 (5) ◽  
pp. 1863-1878 ◽  
Author(s):  
Vahid Rahmati ◽  
Knut Kirmse ◽  
Knut Holthoff ◽  
Stefan J. Kiebel

Calcium imaging provides an indirect observation of the underlying neural dynamics and enables the functional analysis of neuronal populations. However, the recorded fluorescence traces are temporally smeared, thus making the reconstruction of exact spiking activity challenging. Most of the established methods to tackle this issue are limited in dealing with issues such as the variability in the kinetics of fluorescence transients, fast processing of long-term data, high firing rates, and measurement noise. We propose a novel, heuristic reconstruction method to overcome these limitations. By using both synthetic and experimental data, we demonstrate the four main features of this method: 1) it accurately reconstructs both isolated spikes and within-burst spikes, and the spike count per fluorescence transient, from a given noisy fluorescence trace; 2) it performs the reconstruction of a trace extracted from 1,000,000 frames in less than 2 s; 3) it adapts to transients with different rise and decay kinetics or amplitudes, both within and across single neurons; and 4) it has only one key parameter, which we will show can be set in a nearly automatic way to an approximately optimal value. Furthermore, we demonstrate the ability of the method to effectively correct for fast and rather complex, slowly varying drifts as frequently observed in in vivo data. NEW & NOTEWORTHY Reconstruction of spiking activities from calcium imaging data remains challenging. Most of the established reconstruction methods not only have limitations in adapting to systematic variations in the data and fast processing of large amounts of data, but their results also depend on the user’s experience. To overcome these limitations, we present a novel, heuristic model-free-type method that enables an ultra-fast, accurate, near-automatic reconstruction from data recorded under a wide range of experimental conditions.


2015 ◽  
Vol 40 (2) ◽  
pp. 305-321 ◽  
Author(s):  
Lydia Sam ◽  
Anshuman Bhardwaj ◽  
Shaktiman Singh ◽  
Rajesh Kumar

Changes in ice velocity of a glacier regulate its mass balance and dynamics. The estimation of glacier flow velocity is therefore an important aspect of temporal glacier monitoring. The utilisation of conventional ground-based techniques for detecting glacier surface flow velocity in the rugged and alpine Himalayan terrain is extremely difficult. Remote sensing-based techniques can provide such observations on a regular basis for a large geographical area. Obtaining freely available high quality remote sensing data for the Himalayan regions is challenging. In the present work, we adopted a differential band composite approach, for the first time, in order to estimate glacier surface velocity for non-debris and supraglacial debris covered areas of a glacier, separately. We employed various bandwidths of the Landsat 8 data for velocity estimation using the COSI-Corr (co-registration of optically sensed images and correlation) tool. We performed the accuracy assessment with respect to field measurements for two glaciers in the Indian Himalaya. The panchromatic band worked best for non-debris parts of the glaciers while band 6 (SWIR – short wave infrared) performed best in case of debris cover. We correlated six temporal Landsat 8 scenes in order to ensure the performance of the proposed algorithm on monthly as well as yearly timescales. We identified sources of error and generated a final velocity map along with the flow lines. Over- and underestimates of the yearly glacier velocity were found to be more in the case of slow moving areas with annual displacements less than 5 m. Landsat 8 has great capabilities for such velocity estimation work for a large geographic extent because of its global coverage, improved spectral and radiometric resolutions, free availability and considerable revisit time.


Sign in / Sign up

Export Citation Format

Share Document