Matuyama–Brunhes geomagnetic reversal record and associated key tephra layers in Boso Peninsula: extraction of primary magnetization of geomagnetic fields from mixed magnetic minerals of depositional, diagenesis, and weathering processes

We report paleomagnetic records of Matuyama–Brunhes geomagnetic polarity reversal and associated key tephra layers from the Early–Middle Pleistocene marine sedimentary succession in the Boso Peninsula. The outcrop is in Terasaki, Chiba, Japan and ~ 25 km northeast of the Chiba section. The sediment succession consists of a massive siltstone layer of the Kokumoto Formation, Kazusa Group. A tephra layer was identified in the middle of the outcrop with chemical composition comparable to that of the Byk-E tephra layer from the Chiba section defining the base of the Chibanian Stage. Oriented paleomagnetic samples were collected at intervals of 1–10 cm from the siltstone. To identify the primary remanent magnetization, progressive alternating field demagnetization (PAFD) and progressive thermal demagnetization (PThD) were conducted on pilot samples. Identification of primary magnetization with PAFD was not successful, especially for reversely magnetized samples. In addition, magnetization during PThD showed sharp drops around 175 °C, which decreased gradually between 175 °C and ~ 300 °C, and became unstable above ~ 350 °C. To extract the primary remanent magnetization while avoiding laboratory alteration by heating, a PThD up to 175 °C followed by PAFD was conducted. Combined analysis of remagnetization circles enables extraction of primary magnetization with improved reliability. Rock magnetic experiments were conducted during stepwise heating to understand the magnetic minerals involved and to evaluate the influence of laboratory heating. During heating, FORC-PCA revealed significant changes of magnetic minerals at 200 °C, 400 °C, 450 °C and 550 °C. Rock magnetic analyses and electron microscopy indicate that titanomagnetite/magnetite are magnetic minerals contributing to primary remanent magnetization. Greigite was also identified preserving secondary magnetizations during sub-seafloor diagenesis. The presence of feroxyhyte is suggested as secondary magnetization through the weathering of pyrite by exposure to the air after the Boso Peninsula uplift. The correlation of relative paleointensity with the Chiba section provides an age model with sedimentation rates of 30 cm/kyr and 18 cm/kr for the intervals above and below the Byk-E tephra. VGP latitudes are highly consistent with those from the Chiba section based on the age model, which assigns the main directional swing from reversed to normal polarities as 772.8 ± 0.5 ka.


Introduction
The Matuyama-Brunhes (M-B) boundary, the Earth's latest geomagnetic field reversal event, is an important calibration point on the geological timescale, providing a clear marker in the Pleistocene, which has been the focus of many paleomagnetic studies. Studies have been conducted on the M-B boundary using sediments (Clement and Kent 1991;Oda et al. 2000;Channell and Kleiven 2000;Channell et al. 2010) and volcanic rocks (Mochizuki et al. 2011;Ricci et al. 2020), both of which have advantages and complement each other. During the polarity transition, the Earth's geomagnetic field intensity dropped significantly (Valet et al. 2005;Valet and Fournier 2016). In addition, the reduction in geomagnetic field intensity has been recorded as increased production of cosmogenic radionuclides in the upper atmosphere, including 10 Be in marine sediments (Suganuma et al. 2010;Valet et al. 2014) and in an Antarctic ice core (Raisbeck et al. 2006;Dreyfus et al. 2008).
The geomagnetic field intensity data as well as directional records during polarity transitions contain essential information about the Earth's core and its boundary conditions, which are expected to lead to geodynamic models that can explain geomagnetic reversals (e.g., Merrill and McFadden 1999;Leonhardt and Fabian 2007;Nakagawa 2020;Tassin et al. 2021). In particular, preferred paths of transitional virtual geomagnetic poles (VGPs) are considered to hold invaluable information on the conditions at the core-mantle boundary, which have been made evident during the transitions (Clement 1991;Laj et al. 1991;Hoffman and Mochizuki 2012;Hoffman et al. 2020). Importantly, the reversal process is rapid but dynamic compared with the normal process of geomagnetism. Thus, retrieving continuous paleomagnetic records of adequate temporal resolution from sediments could allow us to accumulate the information necessary to understand the reversal process.
In the Boso Peninsula, thick Pleistocene marine sediments, the Kazusa Group and the overlying Shimosa Group were deposited (Additional file 1: Figure S1). The Kazusa Group is well exposed and contains a continuous stratigraphic succession with well-preserved marine microfossils, pollen, paleomagnetic reversal events, and many tephra layers Haneda et al. 2020;Nishida et al. 2016;Okada et al. 2017;Suganuma et al. 2018Suganuma et al. , 2021. The M-B boundary is widely distributed in the area (e.g., Okada and Niitsuma 1989;Kazaoka et al. 2015;Hyodo et al. 2016). The high rate of sedimentation with numerous tephra layers allowed high-resolution magnetostratigraphy studies on the Boso Peninsula.
In 2020, the Executive Committee of the International Union of Geological Sciences ratified the Global Boundary Stratotype Section and Point (GSSP), defining the base of the Chibanian Stage in the Chiba section, Japan (Additional file 1: Figure S1; Suganuma et al. 2021). The sediments captured both terrestrial and marine environmental signals, as well as the last geomagnetic field reversal. The Chiba section reveals a tightly defined Matuyama-Brunhes paleomagnetic polarity boundary, geomagnetic field paleointensity proxies, and numerous tephra beds, allowing the establishment of a robust and precise chronostratigraphic framework.
The M-B boundary was identified in the middle part of the Kokumoto Formation of the Kazusa Group in the central Boso Peninsula by Nakagawa et al. (1969) and Niitsuma (1971Niitsuma ( , 1976. A series of paleomagnetic studies identified the M-B boundary approximately 2 m below the Ontake-Byakubi-E (Byk-E) tephra bed using alternating field (AF) demagnetization (Niitsuma 1971;Okada and Niitsuma 1989;Tsunakawa et al. 1995Tsunakawa et al. , 1999. In contrast, recent paleomagnetic records using thermal demagnetization by Suganuma et al. (2015), Okada et al. (2017), and Haneda et al. (2020) for the Chiba section indicate that the M-B boundary is located slightly above the Byk-E tephra bed, which provides magnetostratigraphic bases for the proposal and the final ratification of the GSSP.
The age of 780 ka for the M-B boundary, which has been frequently cited in the literature, was derived from astronomically tuned benthic and planktonic oxygen isotope records from the eastern equatorial Pacific (Shackleton et al. 1990). This astronomical age of the M-B boundary is supported by the 40 Ar/ 39 Ar ages of 775.6 ± 1.9 ka obtained from lavas recording the reversal in Hawaii (Coe et al. 2004;Singer et al. 2005), which has been updated to 781-783 ka based on the revised age of the Fish Canyon Tuff sanidine standards (Kuiper et al. 2008;Renne et al. 2011). Furthermore, it has been shown that the lock-in of the geomagnetic signal occurs below the sediment-water interface in marine sediments (e.g., Roberts et al. 2013;Suganuma et al. 2011), which means a delayed magnetization acquisition yielding older ages for geomagnetic events than the actual depositions. Assuming a constant lock-in depth for a delayed magnetization acquisition, higher rates of sedimentation should minimize the age offset for geomagnetic field events (deMenocal et al. 1990;Suganuma et al. 2010). This is evidenced by the fact that the youngest astrochronological M-B boundary ages of 772-773 ka were reported for high sedimentation rate records (Channell et al. 2010;Valet et al. 2014). These M-B boundary ages are consistent with the records of cosmogenic nuclides in marine sediments (Suganuma et al. 2010;Simon et al. 2016Simon et al. , 2018 and an Antarctic ice core (Raisbeck et al. 2006;Dreyfus et al. 2008), which are free from magnetization lock-in. Suganuma et al. (2015) reported a U-Pb zircon age of 772.7 ± 7.2 ka for the Byk-E tephra bed. The Byk-E tephra is 1.1 m below the directional midpoint in the Chiba section, where the GSSP is positioned at the base. Astronomical calibration using oxygen isotope records further pinned down the age of the Byk-E tephra to 774.1 ka, which is immediately below the top of Marine Isotope Substage 19c (Suganuma et al. 2021). On the other hand, the M-B boundary in the Chiba section has an astronomically estimated age of 772.9 ka (Suganuma et al. 2018), which is consistent with astronomically tuned paleomagnetic records (Channell et al. 2010;Channell 2017;Valet et al. 2019), and cosmogenic nuclide records (Raisbeck et al. 2006;Suganuma et al. 2010;Simon et al. 2018;Valet et al. 2019).
After the application of simple thermal demagnetization by Suganuma et al. (2015), improvements were made by Okada et al. (2017) to extract primary remanent magnetizations from weakly magnetized sediments influenced by diagenesis using thermal demagnetization (ThD) up to 300 °C, followed by AF demagnetization (AFD). Haneda et al. (2020) followed the procedure of Okada et al. (2017) and provided high-quality paleomagnetic data, which finalized paleomagnetic evidence for the GSSP. Here, we report paleomagnetic records with geomagnetic reversals corresponding to the M-B boundary from silty clay sediments recovered from an outcrop around Terasaki in the Boso Peninsula (Nanayama et al. 2016). Three tephra layers have been identified in the studied outcrop, and the middle one could be correlated with the Byk-E tephra bed, a stratigraphic marker defined as the base of Chibanian (Suganuma et al. 2021). We also demonstrate the applicability of ThD followed by AFD on mixed polarity intervals with unstable magnetization affected by diagenesis, which is comparable to the records for the GSSP reported by Haneda et al. (2020). Further, we show the details of the magnetic minerals involved and evaluate the influence of laboratory heating on magnetic minerals. Finally, we investigated the secondary magnetization carried by iron sulfide minerals associated with sub-seafloor diagenesis and iron oxides generated by iron sulfide oxidation during weathering after uplifting of Boso Peninsula sediments.

Geological background
Thick marine sediments were deposited around the Japanese Islands in response to the subduction of the Pacific Plate beneath the Philippine Sea and North American plates during the Pleistocene (Additional file 1: Figure S1). In the Boso Peninsula, deep-and shallowwater marine succession of ~ 3,000 m thick, the Kazusa Group, was deposited in the early and middle Pleistocene (Additional file 2: Figure S2a). Based on calcareous nannofossils (Sato et al. 1988), planktonic foraminifera (Oda 1977), diatoms (Cherepanova et al. 2002), magnetostratigraphy (Niitsuma 1976), and oxygen isotope stratigraphy (Okada and Niitsuma 1989;Pickering et al. 1999;Tsuji et al. 2005), depositional ages were estimated as ca. 2.4 to 0.5 Ma for the Kazusa Group (Ito, 1992;Ito et al. 2016). In addition, a number of tephra beds provide detailed stratigraphic correlations and the compilation of different types of age data (Machida et al. 1980;Satoguchi and Nagahashi, 2012). These tephras from the Kazusa Group were dated using zircon fission-track dating (Tokuhashi et al. 1983;Kasuya, 1990;Watanabe and Danhara, 1996;Suzuki et al. 1998), U-Pb dating using SIMS , and LA-ICP-MS (Ito et al. 2017). The Byk-E tephra bed, originating from the Older Ontake volcano in Central Japan (Takeshita et al. 2016), was defined as the base of the Chibanian Stage/Age (Suganuma et al. 2021). The Kazusa Group is subdivided into 12 formations (Additional file 2: Figure S2a): the Kurotaki, Katsuura, Namihana, Ohara, Kiwada, Otadai, Umegase, Kokumoto, Kakinokidai, Chonan, Kasamori, and Kongochi Formations in stratigraphic ascending order .

Samples
Paleomagnetic samples were taken from an outcrop located within the area of the Geological Map of Mobara (Nanayama et al. 2016), which is dominated by siltstones exposed in Terasaki Shinden-Nishi along a roadside in the Mobara district, Chiba Prefecture, Japan (Lat. = 35.381058°N, Lon. = 140.311016°E; Additional file 1: Figure S1). The outcrop is ~ 4 m in height, facing west with a tilt angle of ~ 52° (Additional file 3: Figure S3). The lithology of the outcrop is mainly massive siltstone of the Kokumoto Formation (Additional file 2: Figure  S2b), Kazusa Group. The outcrop contains three tephra layers A, B and C. The very fine sand to silt-sized tephra B layer (0-3 mm thick) is located in the middle of the outcrop, which is presumed to be the Byk-E layer and should be slightly below the M-B boundary (Suganuma et al. , 2018Okada et al. 2017;Simon et al. 2019;Haneda et al. 2020). In addition, the pumiceous tephra A and C layers (both are 2-3 cm thick) are located 47 cm above and 60 cm below of tephra B, respectively. The dips of thin sand layers within the thick siltstone are less than ten degrees. Thus, we consider the bedding plane of the strata for the studied outcrop is subhorizontal; no tilt correction on the paleomagnetic results. To identify the paleomagnetic polarity for magnetostratigraphy, we collected oriented paleomagnetic drill cores (1-in. diameter) at intervals of 1-10 cm from 7 blocks (A-G, Additional file 3: Figure S3(c)) in the outcrop (Additional file 4: Figure S4; Additional file 6: Table S1a). Tephra layer B is in the middle between paleomagnetic samples 60 (1.2 cm below of tephra B) and 61 (1.2 cm above of tephra B). Paleomagnetic samples 76 (49.3 cm above of tephra B) and 30 (57.9 cm below of tephra B) are located several centimeters from the base of tephra layers A and C, respectively. Tephra layers, particularly for tephra B, were also collected for chemical analyses.

Tephra analyses
Chemical analyses were performed to confirm that tephra B is correspondent to the Byk-E layer. Since no glass shards were recognized after washing, hornblende grains were picked up from tephra B and analyzed for chemical compositions and compared with those for Byk-E from the Chiba section. Volcanic glass shards taken from a tephra layer, corresponding to the Byk-E, ~ 7 km west of the Terasaki section (M44 in the map of Nanayama et al. 2016; Additional file 1: Figure S1b) were also analyzed and compared with those from the Byk-E in the Chiba section. For chemical analyses, the following two types of scanning electron microscope (SEM) were used (Table 1)

Scanning electron microscope analyses
Highly polished thin sections of samples 77A and 83A were prepared for scanning electron microscope analyses. Observations and mineral identification in the thin sections were carried out using a SEM (HITACHI SU3500) equipped with EDS (Xmax 80 in Oxford Instruments) and electron back-scattered diffraction (EBSD; HKL NordlysNano in Oxford Instruments) at GSJ-Lab. These measurements were performed using the Aztec software (Oxford Instruments). EBSD measurements were conducted under an accelerating voltage of 15 kV, working distance of 18-22 mm, specimen tilting to 70°, and low-vacuum mode. All index data represent points with a mean angular deviation of ≤ 1°. The EDS analysis was also conducted under an acceleration voltage of 15 kV, a working distance of 10 mm, and the high-vacuum mode.

Paleomagnetic and rock magnetic measurements
Two or three specimens (1.1 or 2.2 cm in length) were obtained from each drill core and subjected to progressive demagnetization experiments. To identify the primary remanent magnetization and for subsequent stratigraphic correlation, the following processes were conducted at GSJ-Lab, AIST. Low-field magnetic susceptibility was measured for all specimens using a Kappabridge susceptibility meter (KLY-4S; AGICO). Subsequently, the natural remanent magnetization (NRM) was measured using a three-axis cryogenic magnetometer (SRM-760R; 2G Enterprises) in a magnetically shielded room. For each drill core or for each sampling horizon from block G, a specimen was selected to perform progressive alternating field demagnetization (PAFD) at 1-10 mT increments up to 80 mT using an AF demagnetizer in line with the magnetometer. Several pilot specimens for progressive thermal demagnetization (PThD) were selected and then heated in air at 25-50 °C increments up to 500 °C using a thermal demagnetizer (TDS-1; Natsuhara-Giken). For the specimens other than the pilot ones, PThD in vacuum up to 175 °C followed by AF demagnetization up to 80 mT was conducted. For this process, low-temperature demagnetization with liquid nitrogen (e.g., Ozima et al. 1964) was conducted before each AF demagnetization step to minimize the influence of multidomain magnetization. When three specimens from each core were available, ThD in air up to 300 °C and subsequent PAFD up to 80 mT were conducted. For specimens from block G collected during the second sampling, PThD in air up to 175 °C followed by PAFD up to 80 mT was conducted. We also applied other complex combined demagnetization methods on specimens from block G for the purpose of minimizing laboratory alteration by heating while maximizing the extraction of primary remanent magnetization. The details on the methods, and motivations, advantages and disadvantages are shown in Additional file 6: Table S1b. Paleomagnetic data were processed using Paleomagnetism.org 2.0 (Koymans et al. 2016(Koymans et al. , 2020 including features on combined analyses of remagnetization circles (McFadden and McElhinny 1988).
Several rock magnetic experiments were conducted on the studied specimens to identify their magnetic mineral components. Thermal demagnetization of the threecomponent IRM (isothermal remanent magnetization) was performed on selected specimens based on Lowrie (1990). For each specimen, IRM was imparted with a magnetic field of 2.5 T in the Z-axis using a pulse magnetizer (Model 660; 2G Enterprises) at GSJ-Lab, AIST. A secondary IRM was put at 0.4 T in the Y-axis, which was followed by a field of 0.12 T applied in the X-axis. The procedure makes each specimen acquire magnetization for high (0.4-2.5 T), medium (0.12-0.4 T), and low (0-0.12 T) coercivity components in the Z, Y, and X axes, respectively. The specimens were thermally demagnetized stepwise to monitor the unblocking temperatures of different coercivities and deduce possible magnetic mineral components in the specimens.
Chips of selected specimens were weighted and their magnetic hysteresis parameters (M s , M rs , B c , and B cr ) were measured using an alternating gradient magnetometer (PMC MicroMag 2900 AGM; Lake Shore Cryotronics Inc.) or vibrating sample magnetometer (Model 8604 VSM; Lake Shore Cryotronics Inc.) at GSJ, AIST. The ratio of saturation remanent magnetization to saturation magnetization (M rs /M s ) was plotted versus the ratio of coercivity of remanence to coercivity (B c /B cr ) based on the method proposed by Day et al. (1977). In addition, measurements of the first-order reversal curve (FORC) were performed, which provide enhanced mineral Table 1 Chemical analysis of a tephra layer (tephra B) from Terasaki section (this study), Byk-E from Yoro river (the Chiba section) and Byk-E from M44 section (Nanayama et al. 2016)  and domain state discrimination (Roberts et al. 2014). FORC results were processed using FORCinel (Harrison and Feinberg 2008), and principal component analysis (FORC-PCA) was conducted (Lascu et al. 2015;Harrison et al. 2018). Low-temperature magnetic properties were measured on a sample chip using a magnetic property measurement system (MPMS-5XL; Quantum Design Inc.) at GSJ. The magnetic moment was monitored during warming in a 4-mT field following a zero-field cooling, which provides information on the magnetic transition temperatures of magnetic minerals (e.g., Verwey transition for magnetite (Verwey 1939) and Morin transition for hematite (Morin 1950)). (2) by Borchardt et al. (1971). The lower bound of the SC value for the acceptance of correlation is 0.92 (Froggatt 1992). The SC values calculated from the analytical results of tephra B and Byk-E were 0.96, suggesting that these two tephras are correlated. The volcanic glasses for Byk-E of M44 and the Chiba section provided an SC value of 0.97, which confirms the same origin.  Fig. 1f and g). The morphology of titanomagnetite grains is angular (Fig. 1a). Based on their morphology, these grains are considered to be of detrital or volcanic origin and were transported from the Japanese Islands. Lower minerals with bright reflections and framboidal textures ( Framboidal sulfide minerals are ubiquitous and are typically observed as dense fillings in the chambers of microfossil shells, such as foraminifera in sample 83A (Fig. 2a). Figure 2b is a close-up image of an area in Fig. 2a. Sulfide minerals were composed of coarse (~ 1 µm) and fine (~ 0.3 µm) grained mineral assemblages and their morphologies suggest that they are originated from diagenesis. Figure 2c and d shows the EDS spectra of point A (coarse grain) and B (fine grain) with characteristic peaks of Fe and S, respectively. The peaks for S are lower in Fig. 2d than in c, indicating a lower atomic ratio of S in the grain for Point B than that for Point A.

Scanning electron microscope analyses
The microstructure of these sulfide minerals closely resembles to those reported by Roberts et al. (2011) in their Fig. 2e with coarse grains (~ 1 µm; interpreted as pyrite) and fine grains (~ 0.3 µm; interpreted as greigite). EBSD analysis of coarse grains shows Kikuchi patterns indicative of pyrite, supporting this interpretation. Although EBSD for fine grains was not successful, the EDS of a fine grain showed a subdued peak of sulfur ( Fig. 2d) compared with that of a coarse grain ( Fig. 2c) relative to the peak of iron, which also suggests that the fine-grained framboids are greigite. Although the presence of fine-grained greigite is suggested in the range of SD (~ 0.3 µm; e.g., Roberts et al. 2011), possible presence of MD greigite cannot be excluded.

Paleomagnetism
PAFD was conducted on the specimens from each drill core sample. Typical behaviors against PAFD are shown in Fig. 3b  175 °C, the magnetization decreases gradually up to 300-350 °C (Fig. 3a) and becomes unstable above ~ 350 °C.
Considering the significant drop of NRM by heating up to 175 °C without instability, we conducted PAFD following PThD up to 175 °C in vacuum ( Fig. 3d and f ). Although not all specimens exhibit ideal behavior decaying linearly toward the origin, some samples show satisfactory results indicative of primary magnetization (e.g., Fig. 3d and f ).
In addition, some specimens were conducted with PAFD following PThD up to 300 °C in air; Fig. 3c shows one of the specimens exhibiting successful results with negative inclination decaying linearly to the origin. Figure 4 shows the plots of the volume magnetic susceptibility and paleomagnetic results versus depth, where a horizontal broken line at zero indicates the position of the Byk-E tephra. The horizontal lines above and below the Byk-E tephra are tephra A and tephra C, respectively. The two samples slightly above tephra A (sample 76) and tephra C (sample 30) show peaks in susceptibility ( Fig. 4a) and NRM intensity (blue solid circles, Fig. 4b), suggesting that these samples contain volcanic materials.
The NRM intensity before demagnetization were significantly reduced after AFD at 30 mT (blue open diamonds, Fig. 4b). PAFD at 30 mT for samples 24, 27 and 30 with stratigraphic levels between −70 and −60 cm show declinations around 180°, and negative or low positive inclinations ( Fig. 4c and d). Susceptibility and NRM after AFD at 30 mT for this interval ( Fig. 4a and b) show much higher values than the other stratigraphic intervals. This suggests that these samples acquired more stable primary magnetization during a reversed polarity interval that could be resolved using only AFD, whose behaviors against demagnetizations are different from the other intervals. On the other hand, the paleomagnetic directions obtained by the linear regression fitting for the PAFD experiments after heating in air or in vacuum suggest that the samples with stratigraphic levels from −90 to +20 cm have negative inclination for most of the interval (Fig. 4f ). Since thermal demagnetizations were conducted on specimen B for samples 1 ~ 10 as pilot studies, which resulted in instability above 350 °C, no specimens were left for the PAFD experiments following heating. Also, demagnetization behavior against heating between 200 °C and 350 °C are not reliable enough, based on the interpretations using combined demagnetization methods and remagnetization circles on other samples. To avoid misinterpretation, we decided not to use these samples from stratigraphic levels between −120 and −100 cm for the interpretation of paleomagnetic directions ( Fig. 4e and f; Additional file 7: Table S2). The paleomagnetic directions of the linear regression fitting for the PAFD experiments after heating in air or in vacuum for intervals above + 20 cm has mostly positive inclinations except samples 76 (~ 49 cm), 87(~ 119 cm), and 93 (164 cm) ( Fig. 4f; Tables S2). Despite this, the results of linear regression fitting on PAFD experiments after partial heating to 175 °C (in air and vacuum) and 300 °C (in air) are not always satisfactory, with large MAD values (> 15°), and transitional behavior of paleomagnetic directions. Further investigation will be conducted in the following subsection 'combined analyses of remagnetization circles' . Figure 5 shows the results of the low-temperature magnetic property measurements of a sediment sample with the MPMS. Magnetization during warming in the 4 mT field shows a clear inflection point at ~ 117 K. This could be considered as a magnetic transition temperature representative of magnetite (Verwey transition;Verwey 1939). The transition temperature close to that of stoichiometric magnetite (~ 125 K) may suggest the presence of nearly pure magnetite (or low-Ti magnetite based on EDS analysis, Fig. 1d) grains in the sample. A subdued signal of the transition may indicate a partial oxidation (e.g., Ozdemir et al. 1993). Figure 6 is a typical example of the results of the thermal demagnetization experiments of the three-axis IRM. The low coercivity component (0-0.12 T) shows a significant decrease in intensity at approximately 480-560 °C. The low coercivity component also shows a broad decrease around 200-400 °C. The medium coercivity component (0.12-0.4 T) shows a significant decrease of IRM at temperatures around 200-300 °C. The high coercivity component (0.4-2.5 T) is not significant and shows a gradual decrease during heating up to 600 °C. Magnetic susceptibility shows an increase from 400 to 440 °C, which may indicate the production of magnetic (See figure on next page.) Fig. 3 Typical Zijderveld diagrams for sample 2 (a), sample 27 (b, c, and d) and sample 77 (e and f). Samples 2, 27, and 77 correspond to stratigraphic levels of −116 cm, −63 cm, and +57 cm relative to the Byk-E tephra layer, respectively. a PThD for specimen 2B up to 500 °C. Inset is a plot showing magnetization versus temperature. b PAFD for specimen 27A-1. c PThD for specimen 27B up to 300 °C (left) followed by PAFD (right). d PThD for specimen 27A-2 up to 175 °C (left), and subsequent PAFD following LTD at each step (right). e PAFD for specimen 77A-1. f PThD for specimen 77A-2 up to 175 °C (left), and subsequent PAFD following LTD at each step (right). Sample 27 (b, c, and d) and sample 77 (e and f) are considered to have reversed and normal polarity of primary remanent magnetizations, respectively. Stratigraphic distance of paleomagnetic samples from Byk-E are after correction of slope of the outcrop. For each specimen, stratigraphic distance from Byk-E (Additional file 6: Table S1a) is shown after correction of slope of the outcrop in the parentheses Oda et al. Earth, Planets and Space (2022)  minerals due to laboratory heating. Magnetic susceptibility is reduced by further heating to temperatures between 440 and 560 °C. The increase and subsequent decrease of magnetic susceptibility with a peak around 440-480 °C during heating in air might be originated from production of superparamagnetic magnetite by oxidation of sulfide minerals and further oxidation to hematite. Thermal alteration of sediment samples during laboratory heating in air will be discussed in detail in subsection 'laboratory heating experiments' in the "Discussion" section.

Magnetic hysteresis and FORC
The hysteresis parameters of the measured samples are listed in Table 2. Figure 7 shows a Day plot for the ratios of the hysteresis parameters (Day et al. 1977). The values of the measured samples could not be explained by a mixture of SD and MD magnetite (Fig. 7a). On the other hand, these values fall within the region of mixtures between single domain (SD) and superparamagnetic (SP) magnetite with particle sizes between 5 and 10 nm (red curves in Fig. 7b; Dunlop et al. 2002). The measured values also fall within the area of greigite-bearing marine sediments with diagenesis (colored area in Fig. 7b; Roberts et al. 2011). Alternatively, the measured values were compared with those calculated using micromagnetic simulations for greigite ( Fig. 7c; Valdez-Grijalva et al. 2020). It is suggested that the measured values are also consistent with those for non-interacting greigite of 80-90 nm grain size.  The results of the first-order reversal curve (FORC) measurements of the selected samples are shown in Fig. 8. Figure 8e and j shows the results for the samples taken from the sediments containing tephra C (sample 30-4) and A (sample 76-1), respectively. The FORC diagrams are typical of marine sediments with central ridge along horizontal axis and a moderate vertical spread close to vertical axis (e.g., Roberts et al. 2018b). The FORC diagrams suggest that the interactions between the magnetic particles are moderate, but not strong. To determine the magnetic mineral components in the sediments of the outcrop, we conducted FORC PCA ( Fig. 9; Harrison et al. 2018). A 94% variance is explained by two principal components, PC1 and PC2 (Fig. 9a). Figure 9b shows a plot of PC2 versus PC1. Based on the plots, the samples can be explained by three end members, EM1, EM2, and EM3. FORC diagrams for EM1, EM2, and EM3 are shown in Fig. 9c, d, and e, respectively. The horizontal profiles for EM1, EM2, and EM3 are also shown in Fig. 9f, g, and h, respectively. EM1 is very close to sample No. 12 ( Fig. 8l; 89-2), which has a relatively long central ridge up to ~ 250 mT. The negative region is obvious around B c = 0 and the vertical spread along B u is small, indicating a lower contribution of multidomain grains. EM2 is very close to sample No. 10 ( Fig. 8j; 76-1) corresponding to tephra A, which has a short and subdued central ridge up to ~ 80 mT. The vertical spread along B u is large, suggesting a higher contribution of multidomain titanomagnetite and/or magnetite. EM3 is close to sample No. 5 ( Fig. 8e; 30-4) and No. 6 ( Fig. 8f; 39-1), which has a central ridge up to ~ 150 mT. The vertical spread along B u may indicate a moderate contribution of multidomain titanomagnetite and/or magnetite.

Magnetic minerals
Based on electron microscope observations (Fig. 1a, c, d f, and g), we identified titanomagnetite as ferrimagnetic iron oxides in the sediments. Moreover, the observation of the Verwey transition (Fig. 5) confirms  Lowrie (1990). Bottom diagram shows the magnetization for high (red circles; 0.4-2.5 T), medium (pink triangles; 0.12-0.4 T) and low (blue rectangles; 0-0.12 T) coercivity components after each heating step in air. Middle diagram shows the magnetization removed at each heating step for high (red), medium (pink) and low (blue) coercivity components. Top diagram shows the magnetic susceptibility measured after each heating step the existence of low-Ti magnetite in the studied sediments. Thermal demagnetization experiments of threeaxis IRM (Fig. 6) for low coercivity component are also consistent with this finding, suggesting the presence of a magnetic mineral component with unblocking temperatures of ~ 560 °C. The observation of titanomagnetite grains might be consistent with the fact that the low coercivity component shows a broad decrease around 200-400 °C. FORC diagrams (Fig. 8) and FORC PCA analysis ( Fig. 9) show that there is a central ridge in the typical coercivity range of several tens of mT and 100 mT, suggesting the presence of non-interacting SD magnetite/ titanomagnetite as well as MD magnetite/titanomagnetite (EM2; Fig. 9d). In addition, there might be possible minor contribution of SD and MD greigite grains for EM2 and EM3 (see Fig. 2 and discussions in the following paragraph for SD greigite grains). The contribution of non-interacting SD magnetite may possibly originate from detrital grains and/or fossil magnetotactic bacteria. The FORC diagrams for EM2 (Fig. 9d) and EM3 (Fig. 9e) with significant signatures in the region of negative B u around B c = 30 mT may also suggest that magnetite contributes to the vortex state (pseudo-single domain).
Thermal demagnetization experiments of the threeaxis IRM for the medium coercivity component show a significant decrease around 200-300 °C, which is consistent with the observation of greigite in the sediments. FORC PCA analysis (Fig. 9) revealed that EM1 has a relatively long central ridge up to ~ 250 mT, which might be an indication of the presence of greigite as SD particles. Greigite might be present as vortex state particles, which is evidenced by a considerably broadened distribution around B c = 30-120 mT. Although the interpretation of the Day diagram is not straightforward (e.g., Roberts et al. 2018a), the distribution on the diagram is consistent with the area of greigite bearing marine sediments (Fig. 7b).
A striking feature of the thin sections is the strong oxidation observed in the outer part of framboidal pyrite (Fig. 1a, b, c and e). The elemental distribution pattern obtained by EDS suggests that iron remained in the pyrite grains, whereas sulfur in the pyrite grains in the outer part of the framboid was almost completely replaced with oxygen. The replacement of sulfur by oxygen depends on the grain size, the size of the framboid, and the surrounding conditions in the sediment. It has been reported that oxidative weathering of pyrite to an iron (oxyhydr)oxide (either ferrihydrite (See figure on next page.) Fig. 7 Day plots of the measured and compiled hysteresis parameters. a Day plot (Day et al. 1977) of the measured hysteresis parameters (black circles). Numbers correspond to those in Table 2. Nos 5 and 7 are tephra C and A, respectively. SD, PSD and MD shown in blue are domain states corresponding to magnetite. Blue curves and numbers are two SD-MD mixing curves and mixing ratios (Dunlop et al. 2002). b Day plot of the parameters for the measured samples (black circles) together with that for greigite bearing marine sediments around New Zealand . Red curves are mixing curves of SD with 5 nm and 10 nm SP particles and with mixing ratios (Dunlop et al. 2002).  Oda et al. Earth, Planets and Space (2022) (Gu et al. 2020). Feroxyhyte is a planar antiferromagnet with the net sublattice moments aligned parallel or antiparallel to c-axis (Koch et al. 1995). Each particle acquires a net moment due to small number of layers along the c-direction, and the presence of surface steps causing the formation of ferrimagnetic domains with an odd number of layers. The Curie temperature, saturation magnetization, and IRM after saturation for a synthetic feroxyhyte are reported as 182 °C, 14 Am 2 /kg, and 6.7 Am 2 /kg, respectively (Koch et al. 1995). Considering the observation of pyrite oxidation by electron microscopy and the evidence that there is substantial loss of magnetization by thermal demagnetization at 175 °C (e.g., Fig. 3a), feroxyhyte could be one of the principal magnetic minerals in the studied sediments, which might have acquired secondary magnetization after the sediments changed to the oxidative condition associated with the uplifting of the Boso Peninsula. The absence of a concentric distribution of FORC diagrams representative of interacting SD greigite is considered to result from partial oxidation of the outer shell, which will be explained in the next paragraph. This interpretation is also consistent with the distribution on the Day diagram for a theoretical curve for non-interacting greigite particles (open circles; Fig. 7c).
Although the contribution of the total magnetization is minor, thermal demagnetization experiments of three-axis IRM for high coercivity component (Fig. 6) show a gradual decrease during heating up to 600 °C. This may be attributed to the presence of hematite with distributed unblocking temperatures. There is no evidence of unblocking temperatures indicative of goethite (α-FeOOH) lower than its Curie temperature of ~ 120 °C (e.g., Ozdemir and Dunlop 1996).

Laboratory heating experiments
To understand the coercivity distributions of magnetic minerals contributing to magnetizations in the sediments and laboratory alteration during heating, a series of heating experiments were conducted. The first half of the experiment comprises ARM measurements during sample heating in air according to the protocol that resembles the microcoercivity unblocking temperature diagram proposed by Sato et al. (2019). First, ARM heating experiments were performed on selected samples in air. The ARM is acquired before each heating step, which is AF demagnetized and plotted in Fig. 10a. The coercivity spectrum of magnetization remaining after each heating step was calculated by taking the derivative of Fig. 10a in terms of coercivity (Fig. 10b). The coercivity spectrum indicates that the peak coercivity is centered around 25 mT, which was reduced to < 80% by heating up to ~ 200 °C. This is consistent with the reduction in the magnetization intensity observed during PThD experiments (Fig. 3a). The magnetization centered at approximately 25 mT was further reduced to < 70% by heating up to ~ 300 °C, which is also consistent with the PThD experiments. Coercivity after heating above 300 °C shows a considerably broad spectrum extending from 5 to 75 mT, which was gradually reduced to zero above ~ 550 °C.
After AFD following each heating step, ARM is acquired, followed by AFD to identify the thermal alteration during heating (Fig. 10c), and the corresponding coercivity spectrum was calculated (Fig. 10d). The diagram shows that the coercivity of heating products (distribution in Fig. 10b should be subtracted from that in Fig. 10d to observe the actual thermal alteration effect on magnetization) is prominent in the range of 5 mT and 35 mT. The heating product steadily increased with increasing temperature up to ~ 400 °C. Heating above ~ 400 °C to ~ 500 °C produced a significant amount of heating induced magnetic materials, which was slightly reduced by heating above ~ 500 °C up to ~ 550 °C. It should be noted that the spectrum in Fig. 10b is not purely representative of the coercivity distributions of the natural state before heating because ARM was acquired repeatedly after each heating step, introducing a minor amount of magnetization of the heating products. However, this effect could be ignored assuming that the magnetization of the heating product was not significant relative to the total magnetization and that the unblocking temperatures of the heating products were mostly lower than the corresponding temperature during the experiments.
The second experiment was composed of a set of measurements using VSM during the stepwise heating of a sediment sample in air. After each heating step, the hysteresis parameters were measured, followed by FORC measurements. Figure 11a shows plots of hysteresis measurements versus temperature. B cr gradually decreases from ~ 200 °C to ~ 400 °C, and then decreases significantly at 450-500 °C followed by a linear increase at 600 °C. B c shows a gradual decrease from ~ 25 °C to 400 °C and a sudden increase at 450 °C, followed by a slight decrease from 500 °C to 600 °C. M rs decreases slightly from 200 °C to 400 °C, increases suddenly at 450 °C, and then decreases from 500 °C to (See figure on next page.) Fig. 9 Results of FORC PCA (Harrison et al. 2018) on 12 samples. a Variance explained by principal components (solid circles) is shown together with cumulative variance (blue columns). 94% of variance is explained with PC1 and PC2. b FORC for each sample is plotted on PC2 versus PC1 diagram together with end members EM1, EM2 and EM3. FORC diagrams are shown for c EM1, d EM2 and e EM3. Color scales for d and e are the same as c. Horizontal profiles on FORC diagrams along B u = 0 are shown for f EM1, g EM2 and h EM3 Oda et al. Earth, Planets and Space (2022)  The data point at 550 °C moves to the right (B cr /B c = ~ 4.5) along the same trend as that between 400 °C and 450 °C, then further moves to the right (B cr /B c = ~ 5.3) at 600 °C. Figure 12 shows a series of FORC diagrams during stepwise heating in air. The diagrams show a gradual shrinkage of the central ridge from 25 °C to 400 °C, followed by a clear stepwise change between 400 °C and 450 °C. The diagram further changed from 450 to 600 °C. The serial change of the FORC diagram during heating was better captured by FORC-PCA analysis (Fig. 13). Up to 97% of the total variance of FORCs could be explained by two principal components, PC1 and PC2 (Fig. 13a). The FORC diagrams are explained by three end members, EM1, EM2, and EM3 (Fig. 13b). The FORC starts from around EM1 (25 °C in Fig. 13b), then moves gradually to around EM2 with increasing temperatures up to 400 °C. By heating to 450 °C, the FORC suddenly moved to EM3, which was also recognized by direct observation of the two FORC diagrams. The FORC stays around EM3 for a temperature of 500 °C, and then moves stepwise back to EM2 up to a temperature of 600 °C.
FORC diagrams corresponding to EM1, EM2, and EM3 are shown in Fig. 13c-e, respectively. The horizontal profiles of the FORC diagrams along B u = 0 for EM1, EM2, Fig. 10 Coercivity spectrums of ARM monitored during heating of sample 102-1Z-1-B. a The first ARM is acquired before each heating step, which is AF demagnetized and the magnetization is plotted with color for the corresponding temperature and field steps. b Coercivity spectrum of magnetization after each heating step is obtained by dividing the magnetization difference of successive AFD steps by the AFD step difference in a. c The second ARM is acquired and PAFD demagnetized up to 80 mT after removal of the first ARM by AFD up to 80 mT following each heating step, which is useful in identifying thermal alteration. d Coercivity spectrum of magnetization in c. ARM acquisition was conducted with DC field of 50 μT and AC field of 80 mT and EM3 are also shown in Fig. 13f-h, respectively. EM1, which corresponds to the FORC diagram before heating, is characterized by an extended central ridge up to B c ~ 300 mT ( Fig. 13c and f ). The higher coercivity of the central ridge could be attributed to the presence of SD greigite. EM2 (representative of FORC after heating to 400 °C) is characterized by a subdued central ridge extending to B c ~ 200 mT ( Fig. 13d and g), which may suggest the decomposition of greigite by heating. The negative region for B c ~ 0 mT and B u < ~ −5mT typical for the FORCs of EM2 and EM3 (Fig. 13d and e), which is clearly recognized for FORCs at temperatures above 350 °C (Fig. 12), might be associated with SD magnetite and greigite (Roberts et al. 2014). A slight increase in the central ridge of EM2 around 400 °C, just close to the origin (B c ~ 0 mT; Fig. 13g) may suggest the production of superparamagnetic particles by heating, possibly due to the decomposition of feroxyhyte and/or greigite. Removal of SD greigite and/or production of SP particles is also consistent with a slight shift from 25 °C to 400 °C on the Day diagram to the lower right (Fig. 11b).
EM3 corresponding to temperature of 450-500 °C is characterized by subdued central ridge extending to B c ~ 150 mT ( Fig. 13e and h). Another feature is a bump appearing on the central ridge at approximately 30 mT suggesting the production of magnetic particles (Fig. 13h). In addition, a characteristic positive region on the FORC diagram around B c ~ 40 mT and B u ~ −50 mT appeared, which was not obvious for EM1 and EM2. This might be due to the production of magnetic particles by heating, some of which could interact with each other. The production of magnetic particles could be best visualized by the previous ARM diagram for 400-550 °C and ~ 10-30 mT (Fig. 10d). The shift to the PSD region in the Day diagram at 450-500 °C is also consistent with the interpretation that (possibly interacting) SD particles are produced by heating. Finally, the FORC moves from EM3 back to EM2 (Fig. 13b), which may suggest that the particles produced by heating to 450-500 °C were decomposed by heating to 600 °C. The removal of magnetic particles produced by heating at 450-500 °C could also be observed on the ARM coercivity diagram (Fig. 10d) and the Day diagram (Fig. 11b). The magnetic particles produced by heating at 450-500 °C could be interpreted as magnetite and the removal of the signal by heating further to 600 °C could be considered to be due to the oxidation of magnetite to hematite, which has much less magnetization than magnetite.
Microcoercivity unblocking temperature diagram experiments using AF demagnetization of ARM acquired Fig. 11 Change of hysteresis parameters during heating. a Hysteresis parameters during heating in air of sample 102-1Z-3. b Day plot for the hysteresis parameters during heating before each temperature step (Sato et al. 2019) combined with thermal alteration monitoring using ARM acquired after each temperature step is a powerful tool for understanding the coercivity spectrum at each unblocking temperature due to newly produced magnetic minerals by heating. It should be noted that this method is valid, assuming that the unblocking temperature of newly produced magnetic minerals is usually less than the heating temperature that produces the secondary magnetic minerals. On the other hand, FORC and FORC PCA of stepwise heating of sediment samples provide detailed information on the changes in the magnetic mineral and domain state during heating in the laboratory. Van Velzen and Zijderveld (1992) proposed monitoring the coercivity spectrum of the IRM during stepwise thermal demagnetization. Alternatively, Torii et al. (1996) proposed a method to apply SIRM before and after each laboratory heating step to monitor unblocking and alteration temperatures, which is similar to the concept of our study of monitoring ARM acquired before and after each heating step. Overall, the combination of two types of laboratory heating experiments using ARM and FORC  Fig. 12. a Variance explained by principal components (solid circles) shown together with cumulative variance (blue columns). 97% of variance is explained with PC1 and PC2 and three end members are recognized. b FORC for each heating step is plotted on PC2 versus PC1 diagram together with end members EM1, EM2 and EM3. FORC diagrams are shown for c EM1, d EM2 and e EM3. Color scales for d and e are the same as c. Horizontal profiles on FORC diagrams along B u = 0 are shown for f EM1, g EM2 and h EM3. EM1 represents the component before heating. During heating, it gradually moves to EM2 at 400 °C. After 400 °C it moves to EM3 around 450-500 °C then back to EM2 at 600 °C Oda et al. Earth, Planets and Space (2022)  is quite effective in diagnosing unblocking temperatures and thermal alteration critical temperatures for sediments as a complex mixture of magnetic minerals of various origins.

Combined analyses of remagnetization circles
Paleomagnetic directions based on linear regression of the PAFD experiments after heating up to 175 °C and 300 °C are shown in Fig. 4e and f (Additional file 7: Table S2). To maximize the success rate and reliability of paleomagnetic directions, combined analyses of remagnetization circles (McFadden and McElhinny, 1988) were performed on the paleomagnetic results of multiple specimens of each sample drill core (examples are shown in Additional file 5: Figure S5). If there are two components of magnetization with coercivity distributions that do not overlap with each other, we could successfully separate the components by PAFD and linear regression fitting on the corresponding component. On the other hand, if the coercivity spectra of the two components overlap with each other, the demagnetization vector on an equal area projection should fall in a great circle (e.g., Kirschvink 1980). If secondary magnetizations recorded by two specimens have slightly different directions and overlap with the primary magnetization, then the great circles for the directions during progressive demagnetization intersect in the primary direction. Assuming experimental errors or components other than the two components, the great circles may intersect at multiple points. McFadden and McElhinny (1988) formulated to calculate the maximum likelihood estimate of the primary magnetization direction utilizing remagnetization circles. Figure 14b through d shows the paleomagnetic results based on combined analyses of remagnetization circles by McFadden and McElhinny (1988) using Paleomagnetism. org 2.0 (Koymans et al. 2020). Additional file 7: Table S2 shows the details of lines or great circles of individual specimens fitted with linear regression and the combined analyses of remagnetization circles for each sample. The number of specimens used for the combined analyses of each sample was between two and four. For 22 samples, no reliable line could be obtained by fitting to the directional data of progressive demagnetization (reliability category is '1'). If α 95 of the mean direction is less than 15° and t95 is less than 50°, reliability category is '3' . The reliability category of other samples is '2' . VGP latitudes suggest that the directional midpoint of the M-B boundary could be placed at the stratigraphic level between 25 and 55 cm, approximately centered at 40 cm (Fig. 14d). Although samples 91 (195.5 cm) and 93 (208.5 cm) in Fig. 14c show negative inclinations, the comparison of the VGP latitudes with those for Haneda et al. (2020) shows that they are consistent with each other (see Fig. 15e and discussion in the following section).

Correlation with the Chiba section and age model
The paleomagnetic results of the combined analysis of this study were compared with those of the Chiba section (Haneda et al. 2020). In order to facilitate the correlation on various aspects, the relative paleointensity and rock magnetic parameters, including NRM 30-50 /ARM 30-50 , k LF , k ARM , S −0.1T , S −0.3T , used by Haneda et al. (2020) were also measured in this study (Fig. 14e-h; Tables S3 and S4). Haneda et al. (2020) estimated relative paleointensity using NRM and ARM after thermal demagnetization at 300 °C in air. They used a DC field of 0.3 Gauss and an AC field of 80 mT to impart ARM; therefore, we used the same protocol for direct comparison. Then, relative paleointensity was estimated using NRM 30-50 /ARM 30-50 = (NRM 30 -NRM 50 )/(ARM 30 -ARM 50 ) ( Fig. 14e; Additional File 8: Table S3), where ARM x is ARM after AFD at X mT and NRM x is NRM after AFD at X mT, respectively. They also used k LF and k ARM to calculate k ARM /k LF as a proxy of magnetic grain size (Fig. 14f, g; Additional File 9: Table S4), where k LF is the volume magnetic susceptibility and k ARM is the ARM susceptibility. For the calculation of k ARM , ARM was imparted at DC field of 0.5 Gauss and AC field of 80 mT. S ratio (Fig. 14h; Additional File 9: Table S4) was measured for two field values of −0.1 T and −0.3 T. First, SIRM was acquired with a pulse magnetizer (2G Enterprises Model 660 at GSJ-Lab) and IRMs at 0.1 T and 0.3 T fields were acquired in the opposite direction of the SIRM. The S-ratio for −0.1 T (S −0.1 T ) and −0.3 T (S −0.3 T ) were calculated as.
In Fig. 15, paleomagnetic and rock magnetic parameters are plotted versus age (black and blue solid circles) in comparison with Haneda et al. (2020) (purple open circles). Age of Byk-E is fixed at 774.1 ka based on Suganuma et al. (2018). The age model according to Haneda et al. (2020) was estimated based on the correlation between oxygen isotope stratigraphy and a sea-level proxy curve obtained from ODP Site 1123, excluding sand layers. The sedimentation rate of the Chiba section excluding sand layers for the interval spanning M-B polarity boundary is 89 cm/kyr (Suganuma et al. 2018), which is applicable for ages younger than 776 ka (1.65 m below Byk-E). On the other hand, the sedimentation rate before 776 ka is 44 cm/kyr (Suganuma et al. 2018). The relative paleointensity variations (Fig. 15d) were used primarily for the correlation with those for Haneda et al. (2020), which is partly assisted by rock magnetic parameters. The best correlation was obtained by assuming constant sedimentation rates of 30 cm/kyr and 18 cm/kyr for the intervals above and below the Byk-E tephra layer, respectively. Based on this age model, the ages of tephra A and tephra C are 772.6 ka and 777.5 ka, respectively.
The similarity of VGP latitude variations between this study and that of Haneda et al. (2020) is striking. This could be a confirmation of the reliability of the paleomagnetic directions recorded in the sediments of the studied area and the Chiba section, and hence justifies the applicability of combined analyses of remagnetization circles. k LF and k ARM are considerably different, whereas k ARM /k LF is comparable. S ratio (S −0.1 T and S −0.3 T ) are quite similar to each other for the interval centered around Byk-E. The discrepancy in the concentration  Table S2) plotted versus level above or below the Byk-E key tephra layer (corrected for the slope of the outcrop) together with various rock magnetic parameters. From left to right, a NRM intensity, b declination, c inclination, d VGP latitude, e relative paleointensity, f k LF (solid circles) and k ARM (solid triangles) g k ARM / k LF , h S ratio (S −0.1 T and S −0.3 T shown by solid circles and solid triangles, respectively). Black symbols are the results for the 1 st sampling and the blue ones are those for the 2nd sampling (Block G in Additional file 3: Figures S3(c) and S4(e)). In b-d, reliability categories are expressed by the size of the symbol (the largest is the most reliable category '3'). Horizontal pink lines are tephra layers, which correspond to tephra A, Byk-E tephra, and tephra C from top to bottom, respectively (See figure on next page.) Fig. 15 Summary of paleomagnetic results and rock magnetic results plotted versus age compared with those by Haneda et al. (2020). From top to bottom, a S ratio (S −0.1 T and S −0.3 T shown by solid circles and solid triangles, respectively), b k ARM /k LF , c k LF (solid circles) and k ARM (solid triangles), d relative paleointensity (NRM 30-50 /ARM 30-50 ), and e VGP latitude. Details on measurements and calculations of relative paleointensity is explained in the text. Black symbols indicate the results for the 1st sampling and the blue ones those for the 2nd sampling (Block G in Additional file 3: Figure  S3(c) and Additional file 4: Figure S4(e)). Large and small symbols (black and blue circles) correspond to the highest (3) and second highest (2) reliabilities (Additional file 7: Table S2), respectively. Purple open circles are paleomagnetic results from Haneda et al. (2020) and small symbols are paleomagnetic data with MAD > 15°. Vertical dashed line is the M-B boundary. Vertical pink lines are tephra layers, which correspond to tephra A, Byk-E tephra, and tephra C from younger to older ages, respectively. Age model for this study is primarily based on the age of Byk-E tephra (774.1 ka; Suganuma et al. 2018) as a fixed point and correlation of relative paleointensity variations with those from Haneda et al. (2020). The best correlation is obtained by assuming sedimentation rates of 30 cm/kyr and 18 cm/kyr for the intervals above and below Byk-E tephra layer, respectively parameters (k LF and k ARM ) in spite of the similar values of grain size parameter (k ARM /k LF ) may suggest that the studied area had received input of magnetic minerals with similar grain sizes but in smaller amounts. The age for the directional swing from reversed to normal polarities (stratigraphic level of 40 ± 15 cm above the Byk-E tephra) could be estimated as 772.8 ± 0.5 ka using the estimated sedimentation rate of 30 cm/kyr for the  Oda et al. Earth, Planets and Space (2022) 74:80 interval above the Byk-E tephra, which is consistent with the age 772.9 ka reported for the Chiba section (Suganuma et al. 2018). The paleomagnetic results obtained in this study show remarkable similarity with those reported by Haneda et al. (2020), which provide an opportunity for one-to-one correlation with a simple model of constant sedimentation rates (Additional file 10).

Conclusions
A high-resolution paleomagnetic record spanning the Matsuyama-Brunhes polarity reversal boundary was obtained from silty clay sediments from an outcrop at Terasaki in the Boso Peninsula, Japan. The following are concluding remarks on identified tephra layers, rock magnetic signatures, and paleomagnetic results: (1) A tephra layer was identified in the middle of the studied outcrop, which was assigned as the Byk-E tephra (774.1 ka) based on chemical analysis. In addition, two pumiceous tephra layers were recognized 47 cm above and 60 cm below the Byk-E tephra.
(2) Rock magnetic and SEM-EDS analyses revealed that the primary carriers of remanent magnetization were magnetite and titanomagnetite. In addition, greigite was recognized as a diagenetic product in the studied sediment. Furthermore, feroxyhytes may exist as weathering products of pyrite.
(3) Paleomagnetic analyses of the PAFD experiments were not successful in identifying the primary remanent magnetization. PThD experiments suffer from alterations by heating in the laboratory above 175 °C. PAFD experiments following PThD up to 175 °C increased the success rate of extraction of primary magnetizations. To maximize the success rate and reliability of weak magnetization during the polarity transition, combined analyses of remagnetization circles were performed on multiple specimens from each drill core sample. The main directional swing from reversed to normal polarities was recognized at 40 ± 15 cm above the Byk-E tephra. (4) During stepwise heating in the laboratory, the ARM spectra acquired before and after heating were monitored. A significant reduction in magnetization occurs by heating to 200 °C with a peak coercivity of ~ 25 mT for ARM acquired before heating.
Heating to 400 °C resulted in further reduction of magnetization. ARM acquired after heating to 400-500 °C showed a significant increase in the coercivity range of 535 mT, suggesting the production of magnetic minerals by laboratory heating. Further-more, monitoring of FORC diagrams and FORC PCA analyses revealed a gradual shrinkage of the central ridge extending to B c ~ 300 mT by heating from 25 °C to 400 °C, suggesting the decomposition of SD greigite. Heating to 450 °C introduces a bump centered around 30 mT, which may suggest the formation of aggregates of SD magnetic minerals with interactions. (5) The paleomagnetic directions of primary magnetizations could be compared with those from the Chiba section (Haneda et al. 2020). The correlation with the Chiba section based on relative paleointensity provides an age model with sedimentation rates of 30 cm/kyr and 18 cm/kr for the intervals above and below the Byk-E tephra. The age for the main directional swing from reversed to normal polarities could be estimated as 772.8 ± 0.5 ka, which is consistent with the age 772.9 ka reported for the Chiba section (Suganuma et al. 2018).
Additional file 1: Figure S1. (a) Map of Japanese Islands. Red open rectangle is the area shown in Figure S1b. (b) Simplified geological map of Mobara area, Boso Peninsula. The Kiwada, Otadai, Umegase, Kokumoto, Kakinokidai, Chonan, Kasamori and Kongochi Formations are shown in stratigraphic ascending order from southeast to northwest. Red solid rectangles are the localities, where Byk-E tephra layer is identified, including the Terasaki Shinden-Nishi site for this study. The M44, Chonan, and Chiba sections were studied by Nanayama et al. (2016), Okada and Nitsuma (1989), and Suganuma et al. (2021), respectively.
Additional file 2: Figure S2. (a) Schematic columnar section of Kazusa Group. Key tephras, magnetostratigraphy, biostratigraphy, marine isotope stages (MIS) stratigraphy, sequence stratigraphy, ages of tephras based on radiometric datings and correlations are shown. Base of Kazusa Group is suggested as 2.3 Ma based on tephrastratigraphy (Tamura et al., 2019). (b) Schematic lithology column of the Kokumoto Formation. Key tephras, magnetostratigraphy, and biostratigraphy are also shown. Red bar shows massive mudstone used for this study. Magnetozones in the magnetostratigraphy column in a are based on Niitsuma (1976). Geomagnetic