In situ U-Pb dating method for calcite

Abstract

The present disclosure discloses an in situ UPb dating method for calcite, including: cutting a calcite sample to prepare an epoxy resin sample target; placing the sample in a laser ablation sample chamber, and adjusting a position of the sample in an optical axis direction; conducting line scanning ablation on the sample target, and measuring ion signal intensity data of .sup.43Ca, .sup.88Sr, .sup.139La, and .sup.238U; conducting two-dimensional (2D) element imaging to obtain a 2D element content distribution map; according to the 2D element content distribution map, determining a high-U analysis target area, conducting point ablation on the high-U target area, and measuring ion signal intensity data of .sup.206Pb, .sup.207Pb, and .sup.238U; and after the element signal data is obtained, calculating .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb fractionation coefficients, correcting ratios of an unknown sample, constructing a Tera-Wasserbug diagram, and calculating age data and an initial Pb isotope (.sup.207Pb/.sup.206Pb) composition of the calcite sample.

Claims

1. An in situ UPb dating method for calcite, comprising the following steps: a) cutting a calcite sample or a rock sample containing calcite to prepare an epoxy resin sample target with a diameter of 1 inch and a thickness of 5 mm to adapt to a size of a laser ablation sample chamber; b) placing the sample in the laser ablation sample chamber, and adjusting a position of the sample in an optical axis direction to make a laser beam focus well; c) conducting line scanning ablation on calcite in the sample target using the laser beam, loading ablated aerosol into a quadrupole inductively coupled plasma mass spectrometry (ICP-MS) plasma source for ionization using a helium gas as a carrier gas, and measuring ion signal intensity data of .sup.43Ca, .sup.88Sr, .sup.139La, and .sup.238U; d) conducting two-dimensional (2D) element imaging using Iolite 4 software, processing data offline, conducting instrument blank deduction for a signal, and conducting instrument drift correction to calculate signal intensities of different elements; and calculating an instrument response coefficient according to a signal intensity and content of each element of NIST SRM 612, calculating a content of an unknown sample according to the instrument response coefficient, and imaging a 2D element content distribution map; e) according to the 2D element content distribution map, determining a high-U analysis target area, conducting point ablation on the high-U target area using the laser beam, loading ablated aerosol into a double-focusing sector field ICP-MS plasma source for ionization using a carrier gas, and measuring ion signal intensity data of .sup.206Pb, .sup.207Pb, and .sup.238U; f) in a measurement process, after every ten unknown samples are tested, repeatedly testing two calibration reference materials NIST SRM 614 for .sup.207Pb/.sup.206Pb correction, three calibration reference materials WC-1 for .sup.238U/.sup.206Pb correction, and two quality monitoring reference materials Duff Brown Tank to ensure same measurement conditions for the reference materials and the unknown samples; g) after the element signal data is obtained, obtaining corresponding fractionation coefficients according to measured values and standard values of .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb of the reference materials, and correcting .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb ratios of the unknown sample; and h) constructing a Tera-Wasserbug diagram through the corrected .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb ratios, and calculating age data and an initial Pb isotope (.sup.207Pb/.sup.206Pb) composition of the calcite sample.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0017] FIG. 1 is a 2D element content distribution map in an example of the present disclosure.

[0018] FIG. 2 is an age data and initial Pb isotope (.sup.207Pb/.sup.206Pb) composition diagram of a calcite sample in the example of the present disclosure.

[0019] FIG. 3 shows element fractionation effect result diagrams of NIST SRM 614, ARM-3, and WC-1 in the example of the present disclosure.

[0020] FIG. 4 is a sensitization effect diagram of auxiliary nitrogen addition in different cone groups on the instrument in the example of the present disclosure.

[0021] FIG. 5A-C show analysis result diagrams of JT in the example of the present disclosure, where FIG. 5A shows S+H cone group+Ar plasma, FIG. 5B shows Jet+H cone group+Ar plasma, and FIG. 5C shows Jet+X cone group+Ar-N.sub.2 plasma.

[0022] FIG. 6A-C show comparison diagrams of results of a two-step method (FIG. 6B) and a one-step method (FIG. 6C) in the example of the present disclosure.

[0023] FIG. 7A-D show analysis result diagrams of Duff Brown Tank in the example of the present disclosure, where FIGS. 7A, 7B, 7C, and 7D respectively represent results of four independent tests at different times, U contents of the four tests are also listed correspondingly, FIGS. 7A, 7B, and 7C use an 85-micron laser beam spot, and FIG. 7D uses a 50-micron laser beam spot.

[0024] FIG. 8A-B show analysis result diagrams of JT in the example of the present disclosure, where FIG. 8A shows results of JT calcite, FIG. 8B shows results of JT calcite surrounding rock, U contents of two tests are also listed accordingly, and a 110-micron laser beam spot is used.

[0025] FIG. 9A-B show analysis results of ASH-15 in the example of the present disclosure, where FIG. 9A and FIG. 9B respectively represent results of two independent tests at different times, U contents of the two tests are also listed correspondingly, and a 110-micron laser beam spot is used.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0026] To make the objectives, technical solutions, and advantages of the present disclosure clearer, the following further describes the present disclosure in detail according to the accompanying drawings and examples.

[0027] An in situ uranium-lead dating method for calcite included the following steps.

[0028] 1) The specific implementation of the technical method of the present disclosure was introduced with four calcite age reference materials: WC-1 (254?6.4 Ma), Duff Brown Tank (64.04?0.67 Ma), JT (13.797?0.031 Ma), and ASH-15 (2.965?0.011 Ma). The U contents of the four calcite standard samples were widely distributed (0.5-20 ?g g.sup.?1), which fully demonstrated the universality of this method for various calcites. These known ages were used as validation criteria in the method to verify the accuracy and precision of the method.

[0029] 2) First, the above calcite samples were cast into sample targets (with a diameter of 1 inch and a thickness of about 5 mm) with epoxy resin, which were lightly polished to expose the calcite sections, and polished, cleaned, and dried for standby.

[0030] 3) The calcite sample targets were placed into a laser ablation instrument chamber, and the sample chamber was purged with a helium gas to fill it with the helium gas.

[0031] 4) An NIST SRM 612 glass reference material was used to optimize the conditions of the quadrupole ICP-MS instrument, making .sup.206Pb and .sup.238U signals optimal, while ensuring that the oxide yield (ThO.sup.+/Th.sup.+) was less than 1.0%, the secondary ion yield (Ca.sup.2+/Ca.sup.+) was less than 2.0%, and the Th.sup.+/U.sup.+ signal ratio was between 0.95 and 1.05.

[0032] 5) In the line scanning mode, the laser beam spot was adjusted to be circular, with a diameter of 50 ?m, a laser energy density of 2.0 J/cm.sup.2, an ablation frequency of 15 Hz, and a line scanning speed of 40 ?m/s, and area scanning was conducted. The area scanning area was generally 5.0 mm*5.0 mm.

[0033] 6) The acquisition procedure was: 10 seconds of blank before instrument analysis, several seconds of laser ablation line scanning of data (depending on the length of the line), and 10 seconds of blank after instrument analysis. After each line scanning of 10 unknown samples, one NIST SRM 612 glass reference material was repeatedly tested (for instrument drift correction and element content calculation) to ensure same measurement conditions for the reference materials and the unknown samples.

[0034] 7) In this example, during 2D element imaging, the ions to be measured and the mass numbers to be characterized were .sup.24Mg, .sup.29Si, .sup.43Ca, .sup.55Mn, .sup.57Fe, .sup.69Ga, .sup.85Rb, .sup.88Sr, .sup.90Zr, .sup.137Ba, .sup.139La, .sup.140Ce, .sup.141Pr, .sup.146Nd, .sup.147Sm, .sup.153Eu, .sup.158Gd, .sup.159Tb, .sup.163Dy, .sup.165Ho, .sup.166Er, .sup.169Tm, .sup.173Yb, .sup.175Lu, .sup.206Pb, .sup.207Pb, .sup.208Pb, .sup.232Th, and .sup.238U, and the measurement integration time was 6 milliseconds.

[0035] 8) The Iolite software (version 4.0) was used to process the data offline. First, the instrument blank deduction was conducted for the signal, and the instrument drift correction was conducted to calculate signal intensities of different elements. According to the signal intensity and content of each element of NIST SRM 612, an instrument response coefficient is calculated, as shown in the following formula (1).

[00001]$\begin{array}{cc}{k}_{\mathrm{el}}=\frac{I_{\mathrm{el}}^{\mathrm{Standard}\mathrm{sample}612}}{C_{\mathrm{el}}^{\mathrm{Standard}\mathrm{sample}612}},& \left(1\right)\end{array}$

[0036] where k is the instrument response coefficient, el is an element, I is the signal intensity, and C is the element content.

[0037] A content of an unknown sample is calculated by formula (2) according to the instrument response coefficient, and a 2D element content distribution map is imaged.

[00002]$\begin{array}{cc}{C}_{\mathrm{el}}^{\mathrm{Unknown}\mathrm{sample}}=\frac{I_{\mathrm{el}}^{\mathrm{Unknown}\mathrm{sample}}}{k_{\mathrm{el}}}& \left(2\right)\end{array}$

2D element content map, as shown in FIG. 1.

[0038] 9) According to the 2D element content distribution map, a high-U analysis target area was determined, point ablation was conducted on the high-U target area using the laser beam (a laser beam spot diameter of 50-110 ?m, a laser energy density of 2.0 J/cm.sup.2, and an ablation frequency of 10 Hz), ablated aerosol was loaded into a double-focusing sector field ICP-MS plasma source for ionization using a carrier gas, and ion signal intensity data of .sup.206Pb, .sup.207Pb, and .sup.238U were measured.

[0039] 10) The acquisition procedure was: 5 seconds of blank before instrument analysis, 30 seconds of laser ablation line scanning of data, and 10 seconds of blank after instrument analysis. After each operation on 10 unknown samples, two NIST SRM 614 were repeatedly tested (for .sup.207Pb/.sup.206Pb correction), three WC-1 were repeatedly tested (for .sup.238U/.sup.206Pb correction), and two Duff Brow Tank were repeatedly tested (for data quality monitoring), so as to ensure same measurement conditions for the reference materials and the unknown samples.

[0040] 11) In the present disclosure, during UPb dating experiment, the ions to be measured and the mass numbers to be characterized were .sup.202Hg, .sup.204Pb, .sup.206Pb, .sup.207Pb, .sup.208Pb, .sup.232Th, .sup.235U, and .sup.238U, and the measurement integration times were 2 milliseconds, 2 milliseconds, 15 milliseconds, 25 milliseconds, 2 milliseconds, 2 milliseconds, 2 milliseconds, and 10 milliseconds.

[0041] 12) The Iolite software (version 3.7) was used to process the data offline. First, the instrument blank deduction was conducted for the signal, and the instrument drift correction was conducted. The .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb ratios were calculated. The element fractionation curve was simulated according to NIST SRM 614 and used to correct the element fractionation of unknown samples. According to NIST SRM 614 and WC-1, instrument fractionation coefficients of .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb are calculated, as shown in the following formulas (3) and (4).

[00003]$\begin{array}{cc}{k}_{7/6}=\frac{\frac{{207}_{\mathrm{Pb}}}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Standard}\mathrm{sample}614\\ \mathrm{Measured}\mathrm{value}\end{array}}{\frac{{207}_{\mathrm{Pb}}}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Standard}\mathrm{sample}614\\ \mathrm{Recommended}\mathrm{value}\end{array}},& \left(3\right)\end{array}$

[0042] where k is the instrument fractionation coefficient, and 7/6 is the abbreviation for the .sup.207Pb/.sup.206Pb ratio.

[00004]$\begin{array}{cc}{k}_{8/6}=\frac{\frac{{238}_U}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Standard}\mathrm{sample}\mathrm{WC}-1\\ \mathrm{Measured}\mathrm{value}\end{array}}{\frac{{238}_U}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Standard}\mathrm{sample}\mathrm{WC}-1\\ \mathrm{Recommended}\mathrm{value}\end{array}},& \left(4\right)\end{array}$ [0043] where 8/6 is the abbreviation for the .sup.238U/.sup.206Pb ratio.

[0044] The .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb ratios of unknown samples are corrected according to the fractionation coefficient, as shown in the following formulas:

[00005]$\begin{array}{cc}\frac{{207}_{\mathrm{Pb}}}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Unknown}\mathrm{sample}\\ \mathrm{Corrected}\mathrm{value}\end{array}=\frac{\frac{{207}_{\mathrm{Pb}}}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Unknown}\mathrm{sample}\\ \mathrm{Measured}\mathrm{value}\end{array}}{k_{7/6}},\mathrm{and}& \left(5\right)\\ \frac{{238}_U}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Unknown}\mathrm{sample}\\ \mathrm{Corrected}\mathrm{value}\end{array}=\frac{\frac{{238}_U}{{206}_{\mathrm{Pb}}}\begin{array}{c}\mathrm{Unknown}\mathrm{sample}\\ \mathrm{Measured}\mathrm{value}\end{array}}{k_{8/6}}.& \left(6\right)\end{array}$

[0045] A Tera-Wasserbug diagram is constructed according to the corrected .sup.207Pb/.sup.206Pb and .sup.238U/.sup.206Pb ratios, and age data and an initial Pb isotope (.sup.207Pb/.sup.206Pb) composition of the calcite sample are calculated, as shown in FIG. 2.

[0046] The element fractionation effect of LA-ICP-MS was one of the main reasons that affected the accuracy and precision of data results. Therefore, the present disclosure took NIST SRM 614, ARM-3, and WC-1 as examples to explore the element fractionation effect. The results are shown in FIG. 3. It can be seen from the figure that the fractionation effects of the three different matrix elements are similar, so the fractionation model of NIST STM 614 can be directly used for the element fractionation correction of calcite.

[0047] The low U content of calcite required high instrument sensitivity. Therefore, the present disclosure explored the sensitization effect of auxiliary nitrogen addition in different cone groups on the instrument. Under six modes (S+H cone group+Ar plasma, S+H cone group+AR-N.sub.2 plasma, Jet+H cone group+Ar plasma, Jet+H cone group+AR-N.sub.2 plasma, Jet+X cone group+Ar plasma, and Jet+X cone group+AR-N.sub.2 plasma), instrument parameters were adjusted, making .sup.206Pb and .sup.238U signals optimal, while ensuring that the oxide yield (ThO.sup.+/Th.sup.+) was less than 1.0%, the secondary ion yield (Ca.sup.2+/Ca.sup.+) was less than 2.0%, and the Th.sup.+/U.sup.+ signal ratio was between 0.95 and 1.05. The results are shown in FIG. 4. It can be seen from the figure that the signal intensity is optimal under the condition of Jet+X cone group+Ar-N.sub.2 plasma instrument, and the instrument sensitivity is increased by 3-5 times.

[0048] In order to further verify the influence of different instrument sensitivity modes on the experimental results, the present disclosure took JT calcite as the research object, and conducted UPb dating experiments under three modes. The experimental results are shown in FIG. 5A-C below. It can be seen from the figure that JT has the optimal result under the Jet+X cone group+nitrogen addition mode. In addition, the figure shows that only this mode can give accurate age data. To sum up, the experiments showed that the instrument had the optimal performance under the Jet+X-cone group+Ar-N.sub.2 plasma mode, and could accurately determine the UPb age of calcite. With the optimized instrument sensitivity, JT could be dated with only 110 ?m, which was significantly superior to the 150-200 ?m previously reported.

[0049] In order to further verify the advantages of the two-step dating method for calcite provided by the present disclosure, namely (1) 2D element imaging, and (2) selection of an area with a high U content for UPb dating, the results of the two-step method and the one-step method are compared, as shown in FIG. 6. It can be seen from the figure that the success rate of the two-step method is obviously superior to that of the one-step method.

[0050] In order to further verify the reliability of the method, three calcite samples are analyzed, as shown in FIG. 7A-D, FIG. 8A-B, and FIG. 9A-B. It can be seen from the figure that the age results of the three calcite reference materials are consistent with the recommended values within the error range. The analysis accuracy is 2-3%.

[0051] In combination with the above data, the present disclosure provides an in situ UPb dating method for calcite, which can provide more favorable technical support for the application of calcite geochronology.

[0052] The data in the above examples were all completed on LA-ICP-MS of Photo Machine Analyst G2 excimer laser series thermoelectric Element XR ICP-MS. The cited examples were only used to explain the present disclosure, and not to limit the present disclosure. Those skilled in the art can obtain the same results on similar LA-ICP-MS instruments according to the method.

[0053] It will be appreciated by those of ordinary skill in the art that the examples herein are used for helping a reader to understand the implementation method of the present disclosure, and it should be understood that the protection scope of the present disclosure is not limited to such special statement and examples. Those of ordinary skill in the art may make other various specific modifications and combinations according to the technical teachings disclosed in the present disclosure without departing from the essence of the present disclosure, and such modifications and combinations still fall within the protection scope of the present disclosure.