List of the largest meteorites of the world.
The Bacubirito meteorite, although largely forgotten by the scientific community after its excavation in 1902, remains an incredible artifact and inspired generations of Mexican scientists. It is the fifth largest meteorite in the world and the longest with a length of 4.16 m. Using the Monte Carlo method, an innovative geometrical model and scanner, we have been the first to calculate the precise dimensions and mass (21 tons) of Bacubirito meteorite’s complex structure. Moreover, we are advocating that it be added to the list of world heritage sites due to its scientific, cultural, and historical significance in Mexico and the world.
- iron meteorites
- anomalous meteorite
- tridimensional model
Meteorites have provided tremendous amounts of information about the planetary system and the cores of planets or star while fascinating societies for centuries. These pieces of asteroids or fragmented planets are composed of rock, iron, or a mixture of both. Iron meteorites were part of the nucleus of planetoids or asteroids, rock from the surface and the rock-iron from the intermediate zone. Mexico’s large surface area contains a considerable number of valuable meteorites, and in the last 100 years, scientific research has evolved to study them. For instance, the Mexican Allende meteorite was used to calculate the age of our solar system .
The Bacubirito meteorite is a famous Mexican meteorite found in a small town named Camichín in the mountain range of Sinaloa in Northwestern Mexico, see Figure 1. The name comes from the closest town Bacubirito, and it has been the source of scientific research in Mexico for the last century. This enormous meteorite has an outstanding length of 4.16 m and a mass close to 20 tons.
We hope that the research presented here will also support our case for the meteorite’s inclusion in the list of Heritage Sites in the United Nations Educational, Scientific, and Cultural Organization (UNESCO) and demonstrate its value to Mexican scientific research and the world’s knowledge of meteorites. From our point of view, it meets the following criteria :
Further arguments for its significance will appear in this chapter, and we will also discuss its impact on Mexican science in addition to the innovative method utilized to determine its mass.
2. Mexican scientific interest in meteorites
The magnitude of the surface of Mexico (1,964,375 km2) allowed the discovery of a great number of meteorites. We can calculate the probability of meteorite impacts  on the Mexican territory with the cumulative distribution
The early Mexican Scientific Community advocated the study and the catalog of meteorites and other geological features, and they successfully petitioned the government to protect and categorize them as a national heritage. For this reason, any of these objects appear in museums or public institutions. The field of geology and the mining industry promoted scientific development in Mexico in the nineteenth century. Scientists trained in these disciplines began to dedicate themselves to the study of the earth and established the procedures to protect, study, and catalog the meteorites.
Two mining engineers , Antonio Del Castillo (1820–1895) and Mariano Santiago de Jesús de la Bárcena y Ramos (Mariano Barcenas for future reference) (1842–1899) featured prominently in the initial study of meteorites in Mexico and guaranteed legal protection for these objects while advancing research in the field. Castillo, an outstanding academic, created schools and institutes about geology and geography (
Mariano Barcenas was one of the most internationally well-known Mexican scientists and one of the most prominent students of Del Castillo. He combined an outstanding academic career with a successful public life as an entrepreneur and politician. In the congress of the Academy of Natural Sciences in Philadelphia in 1876 , he made the Bacubirito meteorite internationally known. Barcenas studied this meteorite by means of a sample that the governor of the state, Eustaquio Buelna-Perez from 1871 to 1875. Afterward, he obtained the composition of the meteorite and gave a mass estimation, though not precisely, since no other information was available without excavation. Despite the interest around meteorite, it remained buried for many more years, mainly due to the difficulty in arriving to the Camichin.
Professor Henry Augustus Ward (1834–1906), from the Rochester Academy of Science of New York and meteorite collector, excavated the meteorite in 1902 in Mexico [6, 7]. He describes in detail the hard work to arrive to Bacubirito and how the 20 individuals were required to unearth this enormous celestial piece, see Figure 2. He also reported the object’s dimensions: 13 feet, 1-inch length (3.96 m), 6 feet, 2 inches wide (1.8 m), and 5 feet, 4 inches in depth (1.52 m) and mass of 50 tons. He estimated the object’s mass assuming a cubic shape (although the meteorite is highly irregularly shaped). Despite of the affection that Camichin’s population had for the meteorite, it was moved to Constitution Civic Center in Culiacan (see Figure 3), the capital of the state of Sinaloa in 1959, a 245.5-km route, in order to facilitate the object’s study.
Professor Vagn Fabritius Buchwald (1929–current) from the Department of Mechanical Engineering in the Technical University of Denmark visited the meteorite in 1978. He evaluated the physical properties of Bacubirito meteorite and determined its volume and mass. He described the width, which varied between 100 and 185 cm and a depth between 30 and 75 cm. The intricate shape and ear-like curve made it difficult to calculate the exact volume. However, he used a photographical method to arrive at a volume of 2.80 m3 ± 10%, corresponding to a weight of 22 tons: a good approximation with the tools available at the time. In 1993, the meteorite was relocated to its current position in Science Center of Sinaloa in Culiacan (Figure 4).
A new generation of scientist in Mexico is analyzing and promoting Bacubirito meteorite [8, 9]. They have obtained new measurements that place it as the longest in the world (4.1 m). They have also promoted it in the State Congress of Sinaloa in 2018 and declared that the meteorite should be considered a cultural and historical patrimony of the state. This was achieved with the support of Dr. Victor Antonio Corrales-Burgueño (1954–today), a congressman and former Head Dean of the Autonomous University of Sinaloa (UAS). Thus, with this legislation, we hope to have the meteorite declared to Human Heritage, so that it will be included in this esteemed list and take place as an “exceptional and universal value.”
One of the strongest arguments to promote this initiative is the studies that we have achieve to measure the weight of this metal giant. As we have seen in its description, its size has made this difficult to establish with accuracy. Notwithstanding, advances in techniques and methods have enabled us to establish the final weight of this specimen. In the following section, we explain this in detail.
3. Geometrical model of Bacubirito to assess its mass
As stated formerly, the meteorite has suffered classification changes in size over time and is considered anomalous in density, which may be due to a lack of detailed investigations. Back in 1975, Buchwald indicated that only a few studies existed on the meteorite, and little has been done 40 years later. There have only been rough estimates of properties, such as its mass, dimensions, and densities, as noted from the large mass variations reported. Moreover, for most of the estimations, uncertainties have not been provided. Accurate measures of geometrical parameters, the mass, and statistics of the regmaglypts were determined from the elaboration of a three-dimensional model of Bacubirito.
The importance of studying this meteorite can be understood in terms of the role it plays on the entry dynamics of the object. If we consider the theoretical prediction for speed given by Regan and Anandakrishnan , two important variables appear namely the drag coefficient Cd and cross-sectional area A that are determined by its geometry. Another important dynamic variable that can be investigated is the mass using the ablation model in Revelle , and from the two parameters just mentioned, an estimate of the object mass can be obtained before entry. Finally, given that the meteorite shape is nonspherical and assuming that it was not greatly modified during entrance, we can use our geometrical model and study the drag force for an asymmetrical object as in Leith . However, these studies are beyond the scope of the present investigation, and for this research, we will limit ourselves with the model acquisition and parameters previously cited.
The Bacubirito meteorite has a rather complex shape (see Figure 5) and large weight (Tables 1 and 2), and thus, determining the surface shape of Bacubirito can be used to infer important information about its passage through the atmosphere and learn more about the meteorite itself. The geometric model can also provide reliable classification for meteorites, and comparisons between this and other recovered pieces contain additional details . This model allowed us to calculate some of its important features: geometric parameters, mass, and statistics of regmaglypts. First, we generated the model, and we defined and calculated the geometric features. An analysis of reported densities to obtain a precise mass estimate was performed. Then, we reported the size and depth of the regmaglypts, and finally, we obtained a precise geometrical model of the meteorite.
Under these restrictions, we have chosen to use a portable scanner to determine the model. This instrument allows for a vastly accurate estimation of the relative positions of surface points on a rigid object and is highly suited from a practical perspective—the outdoor location of the meteorite and its resting on a high position are benefits in this case. This device has been used in civil engineering applications, where detailed, precise and three-dimensional representations are required. Being a highly characterized instrument, it allows us to reduce the uncertainty of our results.
The scanner used was the Leica Nova MS50. This equipment radiates and collects a laser beam, which can either directly interact with the object to measure it or by means of a prism. The direct operation mode is mostly employed since it offers a higher resolution and precision to acquire the 3-d points. According to an extensive characterization from the manufacturer (see Leica Nova MS50 Datasheet), it has an angular accuracy of =1″ in horizontal and vertical angular measurements, a linear precision of
Before proceeding to computing the geometric parameters, we propose a Monte Carlo simulation to obtain even more precise positions and evaluate the uncertainty of our results from the single series of measurements. A classical reference to the Monte Carlo method and its applications in many areas is found in Rubinstein and Kroese . With this goal in mind, we first develop a simple probabilistic model and later explain its practical implementation. We start assuming that the uncertainties in our measured positions are random, namely, that our calibrated instrument presents negligible systematic deviations as supported by the small parts per million on its linear accuracy. Furthermore, both distance and angle measurements present systematic-free errors that follow a Gaussian distribution. Now, the propagation error in position due to the angular uncertainty is negligible (≪1 mm) for our short separation distances and angular precision, and we conclude that this will follow a Gaussian distribution. From this discussion, the deviation of a measured position with respect to the instrument can be written as , with a unit vector (see Figure 6), where the parameter k is the random variable in the Monte Carlo simulation with a Gaussian distribution centered at zero and . This information is valid for both the characterization provided by the manufacturer (see datasheet) and the measurement technique in which the instrument was positioned in different locations with uncertainties lower than 4 mm.
The application of the method consists of taking the
Finally, we start determining Bacubirito’s dimensions. We define its length as the separation between the two farthest-apart points and fix an axis onto these points. We next define the width as the distance between the two farthest-away points that are both perpendicular to the axis and lie on a plane containing the axis. Similarly, the thickness is defined as the largest distance between the two points that define a line orthogonal to that plane and that pass through the axis. A sketch is shown in Figure 7.
In the next step of experiment, we estimate the meteorite volume by combining the Monte Carlo simulation with Gauss’ Theorem. A Monte Carlo experiment is first employed. Gauss’ Theorem is then used to estimate the volume enclosed by the simulated surface defined by a Delaunay triangulation. Thus, we obtained the associated volumes for each simulation run.
We continued with the analysis of the meteorite’s densities reported so far and its mass determination. The percentages of its chemical elements are shown in Table 3 (see Ref. ). Those values were measured in different samples from Bacubirito. On average, the measurements should be an approximation of the full meteorite’s concentrations and densities. We obtained the average density of 7.7250 ± 0.2061 g/cm3 from the mentioned table and with the aid of the volume calculated the mass.
The diameter and depth of the regmaglypts were the next determined quantities using the model. We defined a regmaglypt’s diameter as that of the best-fit circle for the corresponding pit mouths. The depth is the difference between the bottom and the top. In order to characterize the distribution of regmaglypts, three zones around the meteorite were selected for presenting contrasting structures as shown in Table 1. The diversity in structures was assessed statistically using the Kolmogorov-Smirnov test for a thorough and rigorous treatment of the technique. To this end, we retrieved diameters and depths for each region from the model and took each as a random variable. The null hypothesis that any two of the regions originate from a same distribution for either variable was then evaluated.
3.1 Results and discussion
Renders of our tridimensional model are given in Figure 8 and can be observed in a high level of detail. A precise retrieval of any geometrical parameters and its derivatives is also possible. We focus in this section on the dimensions, mass, and statistics of the regmaglypts.
The basic defining parameters of Bacubirito are the dimensions and their uncertainties depicted in the last column of Table 4. Contrasting differences are observed with respect to one of the first studies performed by Ward and the latest one by Buchwald. While estimates of geometric parameters of Bacubirito have been reported previously, different studies exhibited noticeable changes. We made the formal definition for consistency among studies and determine them precisely from the model. Therefore, we can confirm the old claim that Bacubirito is the longest meteorite in the world, with a total length of 4.130 ± 0.005 m (see Table 4).
The previous estimate of volume was based on photographs by Buchwald  resulting in 2.8 ± 10% m3. This was very good approximation considering the prevalent instrumentation and data analysis techniques of the time. The volume calculated from a number of simulations through the Monte Carlo method is:
This represents a novel rigorous calculation of the volume of the meteorite. The remaining two parameters give the reader the idea of the aspect ratio of the body and its importance among meteorites in the world.
We emphasize that the density of points considered here (a total of 1,812,875) reflects an upper bound for an instrument with uncertainties in the millimeter range. An increase in the number of points would result in separation distances along the meteorite surface in the range of 𝜇m, which is well beyond the maximum precision of the device.
The mass is calculated from averaging the densities in Table 3 and our volume estimation.
We obtained a mass uncertainty of around 3%, and this was consistent with the estimations made by Buchwald and Sanchez-Rubio et al. . Furthermore, we can note that these densities represent the main uncertainty, and hence, we recommend improving their estimation.
Finally, an important feature of iron meteorites is their regmaglypts (Figure 9). Those meteorite traits provide important information related to its fall. For example, according to Lin and Qun , larger pits correspond to zones where the surface vector is more aligned with the meteorite velocity vector. The statistics of the regmaglypts are depicted in Table 5, which demonstrates different areas on the meteorite that exhibit contrasting average diameters and depths. To confirm those differences, a Kolmogorov-Smirnov test shows a
|Buchwald||1975||22 ± 10%|
|Dimensions||Ward [m]||Buchwald [m]||This work [m]|
|Width||1.88||1.00–1.85||2.053 ± 0.005|
|Thickness||1.63||0.30–0.75||0.911 ± 0.005|
|Length||3.99||4.25||4.130 ± 0.005|
|Region||Diameter [cm]||Depth [cm]|
|Front||8.0 ± 2.3||1.0 ± 0.3|
|Rear||12.6 ± 2.6||1.4 ± 0.3|
|Below||10.9 ± 2.6||2.2 ± 0.4|
These data obtained from Bacubirito distinguish it from other meteorites in the world in length and size and for previous generations of scholars presented a significant challenge.
In this work, we obtained a geometrical model that determined the mass of the Bacubirito meteorite and its other geometrical features. Our results are listed as follows:
The mass has a value of 19.43 ± 0.51 tons. The results are both the geometric dimensions and mass. These were achieved through state-of-the-art analysis techniques and equipment, considerable improvement on previous estimations.
The volume of the meteorite is 2.5151 ± 0.0005 m3. We utilized the Monte Carlo method to simulate a series of measurements for a tridimensional rigid object in order to reduce time and expenses related to the repetition of in-field measurements.
The maximum length of the Bacubirito is 4.130 ± 0.005 m—the longest known in the world, and we also obtained estimations of the meteorite’s main dimensions.
The depth and width of the regmaglypts of three regions were retrieved. The regions were found to exhibit different structures. The zone labeled as below is most likely to have been exposed directly to the atmosphere in the ablation altitudes.
The study also supports our proposal to include the Bacubirito meteorite on the list of the World Heritage Sites for its tremendous impact on culture, science, history, and geological studies. In addition to the impressive physical characteristics, mass, and length of the meteorite, the early Mexican scientific community developed through its research of the object and developed a series of legislative initiatives to conserve celestial objects. Moreover, local residents of Bacubirito and Sinaloa believe that the object represents national heritage and is a source of local pride.
Amelin Y, Krot A. Pb isotopic age of the Allende chondrules. Meteoritics & Planetary Science. 2007; 42(7-8):1321-1335
Leon-Portilla M. The Aztec Account of the Conquest of Mexico. Beacon Press; 2011
Jambon A. Bronze age iron: Meteoritic or not? A chemical strategy. Journal of Archaeological Science. 2017; 88:47-53
Morelos-Rodriguez L. La geología mexicana en el siglo XIX. Plaza y Valdés; 2012. 356 p
Mariano B. On certain Mexican meteorites. Proceedings of the Academy of Natural Sciences of Philadelphia. 1876; 28:122-126
Ward HA. On Bacubirito, the great meteorite of Sinaloa, Mexico. Science. 1902; 16(401):395-398
Morelos-Rodriguez L. 125 años de las Meteoritas en el Palacio de Minería. Facultad de Ingeniería; 2018
Terán-Bobadilla E, Abundis-Patiño JH, Añorve-Solano C, Moraila-Valenzuela CR, OrtegaGutierrez F. Bacubirito: The longest meteorite in the world. Astronomy & Geophysics. 2018; 59(2):2.30-22.1
Terán-Bobadilla E, Abundis-Patiño JH, Añorve C, Moraila CR, Ortega-Gutiérrez F, Aragón-Calvo MA. On a novel geometric analysis of the Bacubirito meteorite. Earth, Moon, and Planets. 2017; 120(2):101-111
Regan FJ, Anandakrishnan SM. Dynamics of Atmospheric Re-entry. American Institute of Aeronautics and Astronautics; 1993
Revelle DO. A quasi-simple ablation model for large meteorite entry: Theory vs observations. Journal of Atmospheric and Terrestrial Physics. 1979; 41(5):453-473
Leith D. Drag on nonspherical objects. Aerosol Science and Technology. 1987; 6(2):153-161
Buchwald VF. Handbook of Iron Meteorites. Berkeley: University of California Press; 1975
Dirk K, Rubinstein R. Simulation and the Monte Carlo Method. John Wiley & Sons; 2016
Sánchez-Rubio GR-KR, Lozano-Santa Cruz R, Medina-Malagón S, Flores-Gutiérrez D. In: UNAM, editor. Las meteoritas de México. Instituto de Geología; 2001
Lin TC, Qun P. On the formation of regmaglypts on meteorites. Fluid Dynamics Research. 1987; 1(3):191-199