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

Eglington, Bruce bruce.eglington at usask.ca
Mon Oct 7 08:20:51 AEDT 2019


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/XBRWCQnzP0tV26XJSx9Aku?domain=aka.ms>
________________________________
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/GsbACROAQotYMrGysNTW9q?domain=benmather.info>
t https://protect-au.mimecast.com/s/h2m6CVAGXPtBWl2zsyg4f1?domain=twitter.com<https://protect-au.mimecast.com/s/zc3_CWLJY7i3Mj6NunhMTW?domain=twitter.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/20191006/a45355f9/attachment-0001.html>


More information about the GPlates-discuss mailing list