Supplementary Code S1_MagmaticMigrationRateCalculator

Code Usage Instructions

Data Preparation

Please prepare a CSV file (e.g., magma_data.csv) that includes the following columns:

Age (Ma): The age of the magmatic rocks in millions of years (e.g., 160, 155, ..., 100).

Longitude (°E): The longitude in degrees east (e.g., 108.2, 108.8, ..., 117.5).

Latitude (°N): The latitude in degrees north.

Phase: The phase label (1, 2, 3 correspond to 160–140 Ma, 140–125 Ma, 125–100 Ma, respectively).

Place this CSV file in the current working directory of MATLAB.

Sample data format:

Age

Longitude

Latitude

Phase

160

108.2

33.5

1

155

108.8

33.7

1

...

...

...

...

100

117.5

34.8

3

Parameter Adjustment

outputFileName: Customize the name of the output result file.

Running the Code

Execute the script in MATLAB, and it will automatically generate:

Graphical Output: A scatter plot with regression lines, labeled with the average positions of each phase.

Textual Output: A text file saving the phase-specific migration rates and regression analysis results.

Output example

Figure presents the spatiotemporal migration of magmatism, displaying an age-longitude plot featuring magmatic samples color-coded by phase and a regression line.

Interpreting the Results

Phase-Specific Rates: The migration rates for each phase are based on the weighted average longitude.

Overall Rate: The statistically significant trend rate derived from linear regression analysis.

Detailed Steps

Prepare Your Data: Ensure your data is formatted correctly in a CSV file with the necessary columns as specified above.

Customize Parameters: Open the MATLAB script and modify the outputFileName variable to your desired output file name if necessary.

Execute the Script: Run the script in MATLAB. This will perform all necessary calculations and generate the outputs.

Review Outputs:

The graphical output will display the migration of magmatic activity over time, with different colors representing different phases.

The textual output will contain detailed migration rates for each phase and the overall trend rate, including statistical significance.

Interpret Results: Use the provided outputs to interpret the migration patterns of magmatic activity across the specified phases.

This script is designed to provide a comprehensive analysis of magmatic migration rates, offering both visual and statistical insights into the data. Ensure MATLAB is properly configured to execute the script without errors. If issues arise, check the data file path and format, as well as MATLAB's current working directory settings.