A 200 km suspected impact crater Kotuykanskaya near Popigai, Siberia, in the light of new gravity aspects from EIGEN 6C4, and other data

We provide arguments in favour of impact origin of a 200 km suspected impact crater Kotuykanskaya near Popigai, Siberia, Russia. We use the gravity aspects (gravity disturbances, the Marussi tensor of the second derivatives of the disturbing geopotential, the gravity invariants and their specific ratio, the strike angles and the virtual deformations), all derived from the combined static gravity field model EIGEN 6C4, with the ground resolution of about 10 km and a precision of about 10 milliGals. We also use the magnetic anomalies from the model EMAG2 and emphasize the evidence of much deeper sources in the suspected area, constraining the impact origin of this structure.

complete review of the theory (with examples) is in [5][6][7] . The examples devoted to the impact craters can be found in 5,[8][9][10] and cannot be repeated here owing to space reasons and selfplagiarism; thus, an extensive Supplement (with theory and examples) is added.
The Marussi tensor Γ having five independent components delivers more complex information than the gravity anomalies Δg alone; the location of the target body is related to the radial component T zz , the orientation and shape of the causative structure is characterised through the non-diagonal components T ij of Γ 11 . The gravity invariants I 1 and I 2 act as non-linear filters augmenting the sources of substantial volumes; they help to divide larger density anomalies into individual units. Their specific ratio tells us about 2D structures. The strike angles θ inform us on how the observations (used to derive Γ) correlate with the main directions of the underground structure; these are the main direction(s) of Γ. Analogously to the tidal deformation, the virtual deformation vd describes the "tensions" (compression and dilatation) produced by the body in question; one can conceive such a deformation as being due to "erosion" caused solely by gravity (first described in 4 ).
In the literature, many examples of the local use of usual gravity aspects (mostly Δg, sometimes T ij , seldom I 1 and I 2 or θ) can be found (see, e.g 3,11-13 .). Our methodology and various results covering different parts of the world can be found in 4,[6][7][8][9][14][15][16] . The values of all the gravity aspects were obtained by using the software by B. Bucha 17  Data. Gravity. The EIGEN 6C4 (European Improved Gravity model of the Earth by New techniques) 1 is a global, static, combined, comprehensive and detailed gravity field model including gradiometry from the whole GOCE mission (Gravity field and steady-state Ocean Circulation Explorer, ESA, in orbit in 2009-2013). The important fact is that EIGEN 6C4 is a better and higher resolution gravity model in a comparison with all its predecessors. Its terrestrial data (very important component of the model from the high-resolution part) came from the older EGM 2008 gravity model and other sources (like satellite altimetry). Its worldwide resolution is 5 × 5 arcmin, which is about 10 km on the ground and its precision, expressed as a typical standard deviation in the gravity anomalies, is about 10 mGal. The quality of EIGEN 6C4 is not given only by its resolution and precision, but also by a global and nearly regular coverage (excluding small polar gaps) achieved by the five years of GOCE gradiometry (http://www.esa.int/GOCE). In this study, we solely make use of EIGEN 6C4.
Impact structures inventories. The Earth Impact Database is a database of confirmed, proved impact structures, now with nearly 200 objects, see http://www.passc.net/EarthImpactDatabase/index.html, developed and maintained by the Planetary and Space Science Centre, University of New Brunswick, Canada. Unconfirmed impact craters were summarized by Rajmon in the previous SEIS (Suspected Earth Impact Sites) catalogue, now in 18 . Kotuykanskaya objects are not listed there. Another inventory has been worked out by Russian authors, http:// labmpg.sscc.ru/impact/index1.html; see Mikheeva 19 . It contains 260 proven, 260 probable, 2190 potential, 630 questionable and 45 discredited impact structures. A few of them are located near the well-known, confirmed Popigai (I) impact crater. Popigai itself might be in fact a multiple crater 8 , which is, however, not yet confirmed (but it was accepted and applied in 20 ). Kotuykanskaya objects are listed there (using gravity anomalies from old gravity field model EGM 2008 8 ).
Magnetic anomalies. We add magnetic anomalies from the EMAG2 project, which is a global 2-arc-minute resolution grid of the anomaly of the magnetic intensity at an altitude of 4 km above mean sea level, compiled from satellite, marine, aeromagnetic and ground magnetic surveys 21 . For the area of Russia, the authors mention the resolution of magnetic anomalies to be around 5 km. While data are provided in 2 arcmin resolution, the resolution of magnetic signatures is dictated by the distance from the source that exceeds the altitude of the survey [22][23][24][25] . Wherever available, the original shipborne and airborne data replaced the precompiled oceanic magnetic grids and interpolation between sparse track lines in the oceans was improved by directional gridding and extrapolation, based on an oceanic crustal age model 21 .
Surface topography. Topography model ETOPO 1 26 is a global 1-arcmin (cell size) relief model of the Earth surface that integrates land topography and ocean bathymetry from huge number of satellite measurements. Its precision is about 10-15 m in heights (but not everywhere). Topography from ETOPO 1 model for our area of interest is plotted in Fig. 1a. It is not in a conflict with our other results (see below) but, in a contrast to Popigai itself, the existence of a hole after the impact is not shown clearly. It depends on subsequent effects of geologic forces acting on the spot, here probably mainly water erosion, and on resolution of the topography model in the given area. The surface topography for Kotuykanskaya II is available also from a local Russian remote sensing aircraft material 19 but the data themselves are not accessible. Topography from ETOPO1 for the area is in Fig. 1.
Geological map. A general geologic map of the world is available from the Geological Survey of Canada at approximately 1:35,000,000 scale. We have the area of our interest in Fig. 2.
In situ data. Such geological data as shocked quartz and other impact related grains, shatter cones, coesit, stishovite, diamonds (see the case of Popigai) and other signs of the impact origin are not available, at least in open

Results and discussion
Gravity signal expected for impact craters. What gravity signal is expected (according to our experience) for the proven impact crater? It is expected that Δg and T zz would be negative inside the crater (in the "hole" on open air or covered by less-dense newer sediments), changing positive and negative values for the respective rims around the crater and in space between them. The gravity invariants have extreme values inside and around the crater itself and they are significantly concentrated to "spots" (blotches?), creating isolated extremes in the rim(s) along the crater's ring. The strike angles θ usually show a certain degree of orientation (they are "combed") into one prevailing direction inside, around and near the crater, which is not true in an arbitrary location -the  Results for Kotuykanskaya. First, we refer to the gravity anomalies derived from an older gravity model EGM 2008 8 , as amended by circles around the impact and suspected impact structures in 19 . Then, in the same area, with a focus to Kotuykanskaya II, we present a choice of the gravity aspects, computed from EIGEN 6C4 (Fig. 3a-c, 4a-d).
Finally, we add magnetic anomalies (Fig. 5) from the EMAG2 21 . Both structures, the Popigai proven impact crater as well as the Kotuykanskaya II putative impact crater, exhibit a decline in magnitude (more strongly for Kotuykanskaya II) in the magnetic anomalies when moving towards the centre of the proposed crater (the question is their interpretation, see Discussion). www.nature.com/scientificreports www.nature.com/scientificreports/

Discussion
Gravity field. Now we discuss our "gravity observations" for the Popigai and Kotuykanskaya area, first for both crater together (Fig. 3a-c), and then in zooms for Kotuykanskaya II (Fig. 4a-d). The gravity disturbances Δg (Fig. 3a) and the radial component T zz (Fig. 3b) are negative and circular inside the crater or putative crater and Kotuykanskaya II clearly shows a large central peak, namely in the second derivative (which is expected, because the second derivative is more sensitive to shallower density variations). The craters Popigai II, III…, predicted in 8 , are lined up in NW-SE direction from the most intensive Popigai (I) on NW, well visible in Fig. 3a,b, their right upper corner.
Kotuykanskaya has in its southern part a crater inside the main crater also with a central peak well shown mainly by T zz . The rim around Popigai is fragmented (by erosion?). The rim (or rims?) around Kotuykanskaya are very fragmented and partly missing on the eastern side (the same is true for the surface topography, Fig. 1). The vd show dilatation at the rims and compression inside the craters (Fig. 3c).
The zooms for Kotuykanskaya ( Fig. 4a-b) repeat Δg and T zz with more details about the inner crater, central peaks and fragmented ring(s). They show the smaller crater inside the crater also with a central peak. The I 2 values in Fig. 4c concentrate power into isolated peaks in the ring, completely missing on east side. The vd in Fig. 4d exhibit, as also expected from previous experience with other craters, compression inside the crater and dilatation around it in the ring. The central peaks are denoted in all Fig. 4. These are not artefacts but real features; the resolution is still sufficient.
All these observations support the hypothesis about the impact origin of Kotuykanskaya II and Kotuykanskaya I. We accept of course that the gravity-based results are not unique. For this reason we discuss also "magnetic www.nature.com/scientificreports www.nature.com/scientificreports/ results" for the Popigai and Kotuykanskaya area (Fig. 5) trying to strengthen or to challenge our results with the gravity aspects.
While the general expression of magnetic anomalies on Earth is driven by induced magnetization due to main dipolar geomagnetic field areas containing an even moderate amount of titanohematite and/or hematite, reflecting the oxidation of the crust, are dominated by remanent magnetization 25,28,[31][32][33][34] . When the crust is impacted by meteorite whose mass is such that the atmosphere has negligible effect on the slowing down the impactor, the speed of the impact has tens of km/s. The resulting impact and explosion create a shock pressure exposure of the crust that decays exponentially from the impact epicentre 35,36 . Such pressure exposure necessary demagnetizes the crust. As the shock wave compresses the magnetic carriers, it randomizes the overall magnetic energy (on shock wave passing, the carrier is behaving like superparamagnetic grain when magnetic moments are fluctuating chaotically 24,31 ). Different minerals have different sensitivities to the impact demagnetization 36 . In general, the multidomain magnetite is most sensitive, followed by single domain magnetite, multidomain pyrrhotite, and multidomain hematite. The most resistant is exsolved titanohematite, that can be abundant in the Earth's crust 27,37 . The crustal rocks in this area are exposed to the main component of the dipolar field. As a result, there are two components of magnetization that form magnetic anomalies. Positive anomalies are in the direction of the main dipolar field. First look at Popigai in Fig. 5. There are three clusters of positive (red) anomalies around the crater. These are due to iron rich bodies with inherent induced magnetization. Note that Popigai is centred on the edge of the large positive anomaly to the north, cluster of positive anomalies to the south west and small positive anomaly to the south east. The space in between these three anomalies is filled with negative anomaly centred over the crater area. This interruption of the induced magnetic signature we interpret is due to demagnetization of the impact. The demagnetization penetrates down at least the crater's radius.
The Kotuykanskaya structure shows that around the edges of the presumed crater are well defined induced magnetic anomalies (in red). Since the crater is 200 km wide the demagnetization penetrates at least to the same depth of 200 km. Such a process is connected with plasma generation that is capable to penetrate and permanently magnetize the deep previously demagnetized material resulting in four to five blurred magnetic anomalies similar like in the case of the Chicxulub impact crater 38 .
The magnetic map shows details of magnetic field over the suspected impact structures. Note the contrast of sharpness of magnetic anomalies around the Kotuykanskaya II structure and inside the structure. This difference separates the magnetic sources that are close to the surface that shows the sharp magnetic boundaries from sources that are deep-seated and whose boundaries are not as sharp. Positive anomalies (red colour) show www.nature.com/scientificreports www.nature.com/scientificreports/ induced magnetization, locations where there is likely a large amount of magnetite, the main carrier of induced magnetization 27,34 . Geologic units containing these sources are distributed not only in the north from the crater but also southward and eastward from the suspected crater. The centre of the crater where we see the diffused magnetic boundaries have been demagnetized to a depth of several tens of kilometres and left magnetization only of magnetic sources that are deep-seated. This location correlates well with the location of the suspected crater and leaves an evidence of impact demagnetization process 35,36 .
The magnetic map shows similar signs of demagnetization over the Popigai structure. The structure is bounded by four distinct positive anomalies, the largest from the north and two moderate southwest but also one from the southeast. The large negative field detected over the centre is consistent with impact demagnetization similarly like Kotuykanskaya and supports the impact origin.
Assuming geology from the Geological Survey of Canada (back to Fig. 2), the position of the proposed crater indicates the Mesozoic age, while the age of Popigai is of the Cenozoic age. An impact crater that happen to Mesozoic sediment creates a hole in it and exposes the underlying Mesozoic geology. There is an evidence that the impact created an actual topographic relief 39 . The black lines are Cenozoic geologic contacts, and they interestingly round around the Kotuykanskaya structure. This would be an indication that the Cenozoic sedimentary terrain were filling the dimple made by the cratering process in Mesozoic and some of these Cenozoic sedimentary units were subsequently eroding away leaving the exposures of Cenozoic geological contacts. The structure can be a Mesozoic event, from about 251 to 66 million years ago. While we are aware that the age estimates are not precise, the geology map in Fig. 2 and the relationship of the units offers the most simple and straightforward interpretation of the Kotuykanskaya II crater age.

Conclusion
While there is no direct evidence of shock in terms of shatter cones, planar deformation features, high pressure minerals, etc., we provide evidence in terms of new gravity and magnetic data and new methodology that at least Kotuykanskaya II, the suspected impact structure near Popigai, Siberia, is likely created by impact process. Kotuykanskaya would be one of the largest impact craters on the Earth -with a diameter about 200 km, so further studies are recommended.

Data availability
The input data from the models global models EIGEN 6C4, EMAG v2, and ETOPO 1 are publicly available (see the list of reference below). The outputs (the files to plot our figures) are available upon request from the authors.