age-gridding gaps

[michaelchin]

Here is an example of the gaps in the age-gridding process using the 1 Ga model.

[michaelchin]

Investigation is in process. Set to "medium" priority. Let me know if "high" priority is needed.

[michaelchin]

Customer mentioned "the 1 Ga model". Need to confirm with customer which model exactly is mentioned here.

[michaelchin]

I had sent an email to customer to clarify the meaning of "the 1 Ga model". Waiting for customer's response.

[michaelchin]

while waiting for customer's response, I tested with merdith2021 model which goes back to 1Ga as well. There are no seed points in the area which the customer had reported. See the screenshot below. Need to investigate why the missing seed points.

Note: the seed points come from two sources, ocean seed points at the initial time and MOR points at each time step.

[michaelchin]

I have put the intermediate files at https://repo.gplates.org/webdav/mchin/test-agegrid-output-20240613/.

The model being used is merdith2021.

sample_points/: sea floor seed points in both .gpmlz and .pkl format SEAFLOOR_AGE/ : age grid netcdf files agegrids.mp4 : the animation from 1000Ma to 0Ma images/ : .png images from 1000Ma to 0Ma

You may use macOS "Finder" to copy the files (login as Guest).

Alternatively, you may use CyberDuck with "webdav" protocol to download the files.

You may open the .gpmlz file in GPlates and colour the points by "FeatureAge".

[michaelchin]

Investigation is still in progress. Need to find out why there are no "seed points" in certain areas.

[jcannon-gplates]

Thanks for putting that all together.

I see what you mean now! The seed point GPML files are very large and slow (as you mentioned in the last meeting).

Here’s a better way to do this – as done in the original isopolate - and should help a lot…

Instead of saving as individual points you can save as a multipoint coverage. Basically, each point gets one or more scalar values – such as age, spreading rate, direction, obliquity, etc. But it’s all compactly stored in the GPML file – so it’s much smaller and loads faster. In isopolate, you can see how it’s done here (where this function shows how the age, spreading rate, direction, obliquity, etc, are added). This ultimately hinges around pygplates.Feature.set_geometry() (see the bottom example in that function).

The following screenshot shows some old isopolate output (it’s actually lines instead of multipoints, but it’s the same idea, each point gets a different colour). And it’s easy to switch between the different output types (eg, age, spreading obliquity, etc) and colour each output type with a different palette. This is in grey layers (that aren’t commonly seen in GPlates).

Also, if you could set the feature time-period on each GPML to [time+0.5, time-0.5] that would be great! Because it means we can load all the GPMLs at once and animate the time slider - and only one GPML will be visible at a time. That allows us to more easily animate the seed points to see where they go over time. Currently, when you animate the slider you also have to disable some layers at later times so they don’t overlap the current time.

[michaelchin]

At around 428Ma, Laurussia plate started to move to right and away from Farallon Ocean plate. See the screenshots below. While Laurussia moving to right, the seed points in Farallon Ocean plate did not follow and there was no middle ocean ridge between the two plates. So no points to fill the void and a gap started forming. At around 410Ma, a big gap formed. I have very limited plate tectonic knowledge and I don't want to pretend. I could not understand what was going on there.

[michaelchin]

The seed points come from two sources. One is the initial oceanic points and the other is the middle ocean ridges. So either the old seed points move into the gap or the middle ocean ridges generate the new points, otherwise the gap is inevitable.

[michaelchin]

I had sent an email to customer to clarify the meaning of "the 1 Ga model". Waiting for customer's response.

Still no response from customer.

Assume "the 1 Ga model" means any models which go back to 1000Ma, such as merdith2021, unless the customer says otherwise.

[jcannon-gplates]

Thanks for making those changes! Much easier to view now.

I agree with your assessment. It appears the ocean plate is moving away from the continent. Which leaves a gap that cannot be captured by the seed points (for the reasons you stated).

I guess that’s where the interpolation ideas bounced around in recent meetings (a separate concept than seed point reconstruction) might need to come into play.

Just a minor point: in future if you add more than just seafloor ages to the sample points (eg, spreading rate/obliquity/etc) then it might be easier to view if all sample points are in a single file – only because with hundreds of file/layers (ie, one per time step) you’d have to go into each layer and select spreading rate, for example, if you wanted to switch from viewing seafloor age to spreading rate. And if that single file is too big then could always have one file every 100Myr or so, for example.

[michaelchin]

it might be easier to view if all sample points are in a single file

yes, it is a good idea.

[michaelchin]

@jcannon-gplates I tried to recall what decision had been made in the auscope meeting on 2 July, 2024, but my memory failed me again. Did we conclude anything in the meeting regarding this age gridding gaps?