Siderite data within the Rock Magnetic Bestiary

import pmagpy.rockmag as rmag
import pmagpy.contribution_builder as cb
import pmagpy.ipmag as ipmag
%config InlineBackend.figure_format = 'retina'

from bokeh.plotting import figure, show
from bokeh.io import output_notebook, save
output_notebook()

Import data¶

The Rock Magnetic Bestiary siderite data are available within MagIC (UndefinedundefinedNick Swanson-Hysell (2025)). Here we download these data and unpack them to get the measurements

magic_id = '20384'
dir_path = 'example_data/siderite'

result, magic_file = ipmag.download_magic_from_id(magic_id, directory=dir_path)
ipmag.unpack_magic(magic_file, dir_path, print_progress=False)
siderite_contribution = cb.Contribution(dir_path)
measurements = siderite_contribution.tables['measurements'].df

Download successful. File saved to: example_data/siderite/magic_contribution_20384.txt
1  records written to file  /Users/penokean/0000_GitHub/RockmagPy-notebooks/RMB_notebooks/example_data/siderite/contribution.txt
1  records written to file  /Users/penokean/0000_GitHub/RockmagPy-notebooks/RMB_notebooks/example_data/siderite/locations.txt
1  records written to file  /Users/penokean/0000_GitHub/RockmagPy-notebooks/RMB_notebooks/example_data/siderite/sites.txt
3  records written to file  /Users/penokean/0000_GitHub/RockmagPy-notebooks/RMB_notebooks/example_data/siderite/samples.txt
6  records written to file  /Users/penokean/0000_GitHub/RockmagPy-notebooks/RMB_notebooks/example_data/siderite/specimens.txt
1243  records written to file  /Users/penokean/0000_GitHub/RockmagPy-notebooks/RMB_notebooks/example_data/siderite/measurements.txt

-I- Using online data model

-I- Getting method codes from earthref.org

-I- Importing controlled vocabularies from https://earthref.org

Summarize experiments¶

rmag.make_experiment_df(measurements)

specimen siderite_Dalizi-YP-01¶

specimen_name = 'siderite_Dalizi-YP-01'
fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data = rmag.extract_mpms_data_dc(measurements, specimen_name)
rmag.plot_mpms_dc(fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data,
                  plot_derivative=True,interactive=True, return_figure=False, show_plot=True)

FC_ZFC_ratio = fc_data['magn_mass'].iloc[0]/zfc_data['magn_mass'].iloc[0]
print(f"FC/ZFC ratio for {specimen_name}: {round(FC_ZFC_ratio,1)}")

FC/ZFC ratio for siderite_Dalizi-YP-01: 12.1

experiment_name = 'IRM-OldBlue-LP-X:LP-X-T:LP-X-F-9811'
Dalizi_ac_data = measurements[measurements['experiment'].str.contains(experiment_name, na=False)]
rmag.plot_mpms_ac(Dalizi_ac_data,
                  interactive=True,phase='in',
                  figsize=(6,4),legend_location='upper right',)

specimen siderite_Roman-01¶

specimen_name = 'siderite_Roman-01'
fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data = rmag.extract_mpms_data_dc(measurements, specimen_name)
rmag.plot_mpms_dc(fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data,
                  plot_derivative=True,interactive=True)
FC_ZFC_ratio = fc_data['magn_mass'].iloc[0]/zfc_data['magn_mass'].iloc[0]
print(f"FC/ZFC ratio for {specimen_name}: {round(FC_ZFC_ratio,1)}")

FC/ZFC ratio for siderite_Roman-01: 17.4

experiment_name = 'siderite_Roman-01-LP-MST-DC-9810'
Roman_Ms_T = measurements[measurements['experiment'].str.contains(experiment_name, na=False)]
rmag.plot_M_T(Roman_Ms_T, interactive=True, legend_location='upper right',)

specimen siderite_unk-01¶

specimen_name = 'siderite_unk-01'
fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data = rmag.extract_mpms_data_dc(measurements, specimen_name)
rmag.plot_mpms_dc(fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data,
                  plot_derivative=True,interactive=True)
FC_ZFC_ratio = fc_data['magn_mass'].iloc[0]/zfc_data['magn_mass'].iloc[0]
print(f"FC/ZFC ratio for {specimen_name}: {round(FC_ZFC_ratio,1)}")

FC/ZFC ratio for siderite_unk-01: 94.5

Export figure for siderite chapter¶

specimen_name = 'siderite_Roman-01'
fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data = rmag.extract_mpms_data_dc(measurements, specimen_name)
fig = rmag.plot_mpms_dc(fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data,
                        return_figure=True,show_plot=False)
fig.savefig('../book/images/siderite_Roman-01_MPMS.png', dpi=300, bbox_inches='tight')

bokeh_fig = rmag.plot_mpms_dc(fc_data, zfc_data, rtsirm_cool_data, rtsirm_warm_data,
                        interactive=True,return_figure=True,show_plot=False)
save(bokeh_fig, '../book/images/siderite_Roman-01.html')

'/Users/penokean/0000_GitHub/RockmagPy-notebooks/book/images/siderite_Roman-01.html'

References¶

UndefinedundefinedNick Swanson-Hysell. (2025). Institute-for-Rock-Magnetism/RMB-carbonates: RMB-carbonates data (Dataset). Magnetics Information Consortium (MagIC). 10.7288/V4/MAGIC/20384