[GPlates-discuss] Estimating uncertainty in plate positions and motion paths

Brandon, Mark mark.brandon at yale.edu
Tue Oct 8 13:43:26 AEDT 2019


There are several approaches for estimating uncertainties for plate reconstructions. The covariance approach discussed previously goes back to the beginning of plate reconstructions. Keep in mind that this approach is a linearized approximation based on a truncated Taylor series. A comprehensive analysis is provided by Chang, Stock and Molnar…. “ON THE STATISTICAL PROPERTIES OF ESTIMATED ROTATIONS” 1987, and “The Rotation Group in Plate-Tectonics and the Representation of Uncertainties of Plate Reconstructions” 1990. This approach assumes that you have estimated covariance matrices for the time series of interest, but you may not have that information.

An alternative approach is to use bootstrap resampling, which allows one to see how the best estimate for a reconstruction, either a single point in time or at multiple points on a path, are perturbed by the residuals (aka misfit) relative to the best-fit solution. The advantage of the bootstrap approach is that it avoids the truncated Taylor series approximation used for the estimated covariance approach of Chang and others.

Note that these methods only consider stochastic errors and are only as good as the data that is available. The “errors” discussed previously include a wider set of issues including biases and departures from assumptions (non-rigid plates), which may or may not be important. Keep in mind that scientists commonly tend to overestimate uncertainties. In fact, humans have little skill with probability. One example, of many, is the popularity of lotto tickets…. Best to rely on statistics rather than intuition for estimating uncertainties.

One more point: There are many fields, with climate analysis being the most visible example, where an ensemble of model results are used to provide a mean solution and an associated variation around the mean. Many people are skeptical of this approach, but there are some good arguments in support.

For both of the issues discussed above, search for “bootstrap resampling”, and “ensemble average”. Remember that, in the 21st century, it all about finding the right search terms….
Best,
Mark


From: GPlates-discuss <gplates-discuss-bounces at mailman.sydney.edu.au> on behalf of Ben Mather <ben.mather at sydney.edu.au>
Reply-To: GPlates general discussion mailing list <gplates-discuss at mailman.sydney.edu.au>
Date: Monday, October 7, 2019 at 6:15 PM
To: Sabin Zahirovic <sabin.zahirovic at sydney.edu.au>, "gplates-discuss at mailman.sydney.edu.au" <gplates-discuss at mailman.sydney.edu.au>, "mdarin at uoregon.edu" <mdarin at uoregon.edu>
Subject: Re: [GPlates-discuss] Estimating uncertainty in plate positions and motion paths

Hi all,

Great discussion! But just to clear up any confusion: the parameter uncertainty problem is less complicated than the inference problem some are alluding to. If you want the total uncertainty of the model given the data, then a likelihood function would need to be constructed that links a reconstruction to a set of observations (palaeomagnetic stripes, fossil records, etc.) and compares the misfit. From there you would sample the posterior probability, using a MCMC algorithm, and build up an ensemble of realisations until you can extract the posterior uncertainty. A frequentist approach may be preferable in some linear cases e.g. anomaly picks - thanks for the reading material, Douwe. Alas, as Sabin pointed out, all of this is outside the scope of what can be achieved short term.

Rather than finding the posterior uncertainty, you can compare a suite of reconstructions using all the different rotation files available and quantify their collective uncertainty through time. The uncertainty on Arabia-Eurasia convergence will not reflect the total uncertainty of the model given the data, but it will reflect the uncertainty between rotation files. Admittedly, this is far from perfect, given the small set of reconstructions we have available and the weighting bias I mentioned earlier, but it's a start...

I hope that gives a bit more insight into what I had in mind. Hopefully this will provide a better framework to assess and compare models - and, as Bruce suggests - illustrate where and when certain models are thought to be better constrained than other parts.

Cheers,
Ben

--
Dr. Ben Mather | Computational Geophysicist
Room 418, Madsen Building F09
School of Geoscience, Faculty of Science
The University of Sydney, NSW 2006
m +61 422 470 117
w benmather.info<https://protect-au.mimecast.com/s/xB71Cr8DLRtG0Y8BC7oIv0?domain=nam05.safelinks.protection.outlook.com>
t https://protect-au.mimecast.com/s/je-uCvl0PoC3j57kHX9e4v?domain=twitter.com<https://protect-au.mimecast.com/s/fOIsCwVLQmiWzmGwi9uZp2?domain=nam05.safelinks.protection.outlook.com>


On 2019-10-07 21:43:42+11:00 GPlates-discuss wrote:
Hi Mike and team,

I agree with Bruce, using alternative models is not likely going to capture the uncertainty present in the Arabia-Eurasia convergence. If one is just looking at the Arabia-Eurasia convergence rate during a time where Atlantic seafloor spreading is occurring, then one probably wouldn’t need to take into account the absolute reference frame, but just evaluate the uncertainties across the relative plate motion hierarchy (Arabia – Africa – North America – Eurasia). However, as Bruce pointed out, other constraints might be useful, especially if intraplate deformation is significant. The reality at the moment is that this is not easy to do, and would generally need to be done on a case-by-case basis, by trawling back through the rotation file and checking whether uncertainty ellipses are published for the poles used, and so on. Would be cool to have this in GPlates one day! In terms of Ben’s suggestion, it might be interesting to see what different models propose/predict for the Arabia-Eurasia convergence. There’s a bunch of models to download from here: https://protect-au.mimecast.com/s/ytcwCxnMRvt34V1WHwlbpS?domain=earthbyte.org<https://protect-au.mimecast.com/s/4lx6CyoNVrcRzwrBhNpBdr?domain=nam05.safelinks.protection.outlook.com> (EarthByte models, Chris Scotese model, etc.), and I think Douwe van Hinsbergen and team have some GPlates-related files here: https://protect-au.mimecast.com/s/ZsXCCzvOWKiApZMjfwp2ch?domain=geologist.nl<https://protect-au.mimecast.com/s/vN8FCANZvPiz3XNDu2swsj?domain=nam05.safelinks.protection.outlook.com>

Cheers,
Sabin

--
DR SABIN ZAHIROVIC | Postdoctoral Research Associate
School of Geosciences | Faculty of Science

THE UNIVERSITY OF SYDNEY

Rm, 404, Madsen Building F09 | The University of Sydney | NSW | 2006
M +61 416 775 589 P +61 2 9351 2467
E sabin.zahirovic at sydney.edu.au<mailto:sabin.zahirovic at sydney.edu.au> | W https://protect-au.mimecast.com/s/n18GCBNZwLiO9r7LUr4UPv?domain=earthbyte.org<https://protect-au.mimecast.com/s/RnjlCD1jy9tkKz58H8rDGX?domain=nam05.safelinks.protection.outlook.com> | R http://www.researchgate.net/profile/Sabin_Zahirovic<https://protect-au.mimecast.com/s/k2g4CE8kz9tNJ23ZC4GW0j?domain=nam05.safelinks.protection.outlook.com>
F https://protect-au.mimecast.com/s/62XgCGvmB5iQpR1zCXebRe?domain=facebook.com<https://protect-au.mimecast.com/s/rxwoCJyp0qhDmvqjCmGEaA?domain=nam05.safelinks.protection.outlook.com> | T https://protect-au.mimecast.com/s/8AcYCK1qJZtGvN2RCykC9-?domain=twitter.com<https://protect-au.mimecast.com/s/hEtlCL7rK8tZprRzHGJ7eg?domain=nam05.safelinks.protection.outlook.com>

CRICOS 00026A
This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments.


From: GPlates-discuss <gplates-discuss-bounces at mailman.sydney.edu.au> on behalf of "Eglington, Bruce" <bruce.eglington at usask.ca>
Reply-To: GPlates general discussion mailing list <gplates-discuss at mailman.sydney.edu.au>
Date: Monday, 7 October 2019 at 8:21 am
To: GPlates general discussion mailing list <gplates-discuss at mailman.sydney.edu.au>, "mdarin at uoregon.edu" <mdarin at uoregon.edu>
Subject: Re: [GPlates-discuss] Estimating uncertainty in plate positions and motion paths

Hi
 The issue is actually much more difficult to conceptualise because different parts of different models are variably uncertain and there is no reason to assume that the average of available models is any closer to reality. I doubt that motion paths are going to resolve better models from worse, other than to identify plate motions that are too fast, etc. In order to better constrain models, we will needs to have far more comprehensive data sets that can be ‘draped’ over the kinematic blocks to assess matches and mismatches. One can already start some of this with data like palaeomagnetics, geochronology, etc and, to a limited extent geochemistry and geology. In most cases, though, we need data sets with much better temporal and spatial resolution than are currently available. There are lots of ways in which machine learning and AI could help, if only we had better organised, structured and relevant data.

The best that can be said at present seems to be that each model can be tested for validity with localised data sets where sufficient detail is available, from which we can start to improve models.

None of this is to say that we should not attempt to develop ways to test model comparisons in anticipation of improved data compilations. By doing so, we will have a better idea of how to assess and compare models. It would also be useful to somehow be able to illustrate which parts of individual models are thought to be better constrained than other parts (time, polygons, motion, etc). Something like this would certainly help address the very real need to quantify model uncertainties to help potential end users.

Bruce Eglington

Get Outlook for iOS<https://protect-au.mimecast.com/s/2I_cCMwvLQTA01qnf4s6N1?domain=nam05.safelinks.protection.outlook.com>
________________________________
From: GPlates-discuss <gplates-discuss-bounces at mailman.sydney.edu.au> on behalf of Ben Mather <ben.mather at sydney.edu.au>
Sent: Saturday, October 5, 2019 4:39:49 AM
To: gplates-discuss at mailman.sydney.edu.au <gplates-discuss at mailman.sydney.edu.au>; mdarin at uoregon.edu <mdarin at uoregon.edu>
Subject: Re: [GPlates-discuss] Estimating uncertainty in plate positions and motion paths

Hi Mike,

Thanks for reaching out.

In principle, if you had enough reconstructions (rotation files) you could build up an ensemble of plate models that would enable you to interrogate uncertainty information from (standard deviation, confidence intervals, etc.), but one of the key problems is it's unclear how to weight each reconstruction. For instance, one could always argue that the latest reconstruction is the "best" because it's the result of painstaking tweaks and adjustments of all those that came before it as a product of additional geological constraints.
Nonetheless, if you see beyond this caveat, a qualitative comparison between the different models is still useful and some quantitative metrics can also be employed. Comparing the plate velocities and cumulative plate distance is relatively simple to do with motion paths, but it can give the false impression that the plate reconstructions do not agree at all. A test we are working on is to place a bunch of points equally spaced across the globe and measure the distance between these points through time for different rotation files. Essentially this is just comparing multiple motion paths across the globe. From there you would calculate (pseudo-)uncertainty maps of where the reconstructions agree and where they do not.

This workflow is particularly suitable for pyGPlates because swapping rotation files a breeze, and it's easy to import other python modules that deal with statistical analysis. Once we have something in place (that works!) I'll send you the python scripts so you can customise it to your problem. Is that OK? Let me know if you have any other suggestions and if you would like to be involved with testing.

Best regards,
Ben
--
Dr. Ben Mather | Computational Geophysicist
Room 418, Madsen Building F09
School of Geoscience, Faculty of Science
The University of Sydney, NSW 2006
m +61 422 470 117
w benmather.info<https://protect-au.mimecast.com/s/yqbNCNLwM9iGr80gCxHfSc?domain=nam05.safelinks.protection.outlook.com>
t https://protect-au.mimecast.com/s/je-uCvl0PoC3j57kHX9e4v?domain=twitter.com<https://protect-au.mimecast.com/s/bNJlCOMxNyty9jpBCXrco1?domain=nam05.safelinks.protection.outlook.com>


On 2019-10-05 05:25:49+10:00 GPlates-discuss wrote:
Hello,
I'm using GPlates to reconstruct the Arabia-Eurasia collision zone and to make a plot showing Total Arabia-Eurasia Relative Motion vs. Time. To do this, I fixed the Eurasia plate and then created a set of motion paths for points on Arabia. This worked great, as I was then able to export a text file containing the coordinates of each reference point at different time steps, which allowed me to calculate the distance that Arabia had traveled between each time step.
However, there must be some uncertainty in the rotation model that is not currently reflected in the generated motion paths. I am wondering if and how I could generate error ellipses in map view for the motion paths, which would allow me to add error bars to each point on the plots that reflect the uncertainties in the locations of each reference point through time? This doesn't seem straightforward to me, so I am reaching out to the discussion list for suggestions. ​
Thank you,
Mike
Michael H. Darin, Ph.D.
Postdoctoral Research Associate
Department of Earth Sciences
University of Oregon
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.sydney.edu.au/pipermail/gplates-discuss/attachments/20191008/ae29eb57/attachment-0001.html>


More information about the GPlates-discuss mailing list