5.1. Mineralogical Evidence
One other issue of note is that a key piece of evidence for the purely impact based formation of craters, as follows from
Section 1.1.2, is the presence of various shocked minerals. However, the very same effect can be produced by a strong electric current, as “... intense heat and pressure of impact events are similar to those experienced during lightning strikes” [
21]. The similarity is so strong that the fulgurites are assumed to be the closest analogs to the products of an impact shock metamorphism [
22,
23], down to microscopic features such as shock lamellae in quartz [
24].
In theory, there may be differences in redox states of minerals affected by electric currents and impact shocks, however, the studies of this issue demonstrate that both may lead to reduction [
21] and in general the results are rather inconclusive. The same situation is with magnetization of minerals by impacts and electric currents: it seems that both can lead to magnetization [
25,
26,
27], as well as demagnetization [
28,
29]—mostly due to intense heat—so the difference between the two is rather hard to pinpoint. The only obvious case where the magnetization could not have been caused by an impact and is most likely produced by electric currents is the sinuous rilles discussed in
Section 1.2.
Cristobalite seems to be an interesting candidate for a “signature” mineral present in fulgurites [
30], though some studies reveal that it’s only observed in younger samples, i.e. over time it relaxes back to quartz [
31]. If there would be significant amount of evidence (for now there is not enough) that cristobalite is indeed the marker of electrically altered minerales, perhaps the fraction of cristobalite that has relaxed back into quartz might become a novel way to date an EDM produced crater. Overall, perhaps further research (which goes beyond the scope of this paper) would facilitate the search for detectable differences, so that impact and EDM hypotheses of crater production would be more distinguishable in terms of their mineralogical results.
There is some evidence for a strictly EDM mechanism of crater production, which comes from an inherently fractal nature of electric currents. For example, the microscopic structures of crystals taken from an “impact site” are “stringer”-like and shocked at the edges, whereas at the center they are more homogeneous and glassified [
32], which mirrors a typical microscopic structure of fulgurites (Ende et al. [
22]. The same (given fractality) can be observed on larger scale—for the craters themselves [
33]. In this context, formation of “stringers” could be caused by the microscopic dendritic pattern of current itself rather than defects and structures in the mineral before the impact, as is assumed by some researchers [
34]. It is curious that cristobalite (see above) is mostly formed at the edges [
35], where we’d expect the current to go.
5.3. Piezoelectricity
Piezoelectricity is another possibility for coupling of mechanical deformation of planetary surfaces (including crater production) with electromagnetic forces. E.g. it has been recently shown [
41] that piezo effects are a viable method of power generation from the ambient mechanical vibration in an urban environment.
In the context of our study, we may consider a
direct piezoelectric effect, where a deformation of the surface of a celestial body (regardless of its cause) may provoke the production of an electric field in it that would lead to a discharge and the formation of a crater. Likewise, a
reverse piezoelectric effect may play its role, where an external electric field (e.g. the field of an infalling plasma cloud, as described in
Section 3.1 and
Section 5.4 would produce a mechanical deformation of the surface. It is hard to imagine an exact arrangement where such deformation would produce a crater, but it is possible that this mechanism might be involved in the production of various accompanying structures.
Understandably, the impact events proper might produce extreme electric potentials in the crust via the direct piezoelectric effect causing discharges during the compression and/or excavation phases (see
Section 1.1.2–
Section 1.1.3. The presence of silicates on Mars and the Moon would accommodate this idea [?], silica being a crystal capable of the piezo effect. Same is true for quartz, which is extremely common on Earth.
Some studies of the Rochechouart crater, France [
42], might lend support to EDM or piezo impact origins: “These minerals are of impact, not meteoritic, origin and probably formed by a combination of shock and plasma processes with the contribution due to each mechanism yet to be resolved”. Silicon carbide was among the minerals observed in the cited study. It is also found in fulgurites [
43], comets, and meteorites (no confirmation on asteroids). Diligent study should be applied to consider any circumstantial evidence that might exist linking minerals such as silicon carbide to a possible impact/EDM hybrid process involving piezoelectricity.
For example, under the pressures of 100 GPa assumed for the impact (see
Section 1.1.2 quartz—given its favorable orientation—would experience piezoelectric fields an order of magnitude larger than its breakdown strength:
where
E is the strength of an electric field between two parallel surfaces with an area
S within the quartz,
q is the electric charge that appears on these surfaces,
C/N [?] is the piezoelectric coefficient,
F is the force acting perpendicularly to the surfaces,
P is the pressure, and
is the dielectric permittivity. As we can see, it is an order of magnitude above the breakdown strength of quartz (
V/m at 300 K—[
44].
Hence an electric breakdown under the influence of an impact is likely. The reverse is not necessarily true (i.e. a simple electric breakdown is arguably not enough to replicate the results of an impact), yet this is another curious similarity between the two processes. To the best of our knowledge, none of the previous literature has taken these effects into account neither in the EDM nor in the impact cratering scenarios, whereas, as shown above, they may prove to be important for both.
5.4. Large Scale Craters
Here we shall consider the energies required to produce a large crater on the surface of a celestial body through the mechanism outlined in
Section 3.1, as well as discuss the potential alternatives. To derive an upper boundary on the needed current we note that there are no craters larger than 10
4 km in the Solar System [?], and all the craters over a 10
3 km threshold are “basins” on the Moon and Mars. So we would use 10
3 km as a benchmark for the largest crater in the System. In order to produce such a crater, an energy of 10
26 J is required [
20]—see Formula (
8) for an exact value.
Let’s consider a plasma cloud with flat parallel boundaries infalling perpendicularly on the surface (oriented parallel to the boundaries of the cloud itself). For an order of magnitude estimation we shall consider the following assumption: the energy available for the crater production is equal to the energy of the electric current, and thus is (per mechanism described in
Section 3.1 of the same order of magnitude as the energy required to deviate the initial plasma cloud—i.e. its upper boundary is equal to the kinetic energy of the particles of the cloud itself (see (
15) and further for a more comprehensive analysis). Hence
where
W is the energy required for the crater production (here
J),
M is the mass of the deviated particles, and
v the cloud’s velocity. Obviously,
, where
m is the mass of a single particle (we shall assume they are protons, so
kg),
n is their number density (
;
N is the number of particles, and
V is the volume of the cloud involved in the cratering process). In such consideration
, where
L is the cloud thickness, and
is the surface area of the part of the cloud involved in the cratering. Clearly,
R is of the same order as the crater/celestial body size (
m). Hence
.
Thus, our estimation of values in (
10) leads to (we shall neglect the numerical coefficients and only focus on the order of magnitude factors):
or, expressing the parameters
:
If we assume that the plasma comes from the Sun, then an obvious upper constraint on
would come from an expression
, where
r is the distance from the body to the Sun,
is the mass of the Sun, and
is the fraction of this mass that is carried away by the cloud. Hence
, or (taking
m, i.e. an astronomical unit, so that
):
Considering that
kg, we arrive at
, or, in terms of
:
An obvious upper limit on
v is the speed of light:
m/s, so the lower limit on the fraction of the Sun’s mass carried by the cloud is
. Currently observed coronal mass ejections may have velocities of the order of
m/s [
45], which would correspond to a more realistic estimate of the mass fraction:
(this is the same order of magnitude as the mass of the Earth). Solar CME may have a mass from
kg up to
kg [
46], though recent advancements in observational techniques have made possible measurements of plasma clouds with masses of
kg being ejected by Sun-like stars [
47]. The upper figure of the given solar CME mass (
kg) corresponds to the mass fraction of
, hence the CME required for the production of the largest craters in the Solar System per mechanism described in
Section 3.1 needs to be
times larger in terms of the mass carried away from the Sun.
To estimate how rare such an event may be, unfortunately, we can only rely on extrapolation of the present day
in situ measurements of the solar plasma (e.g. based on the distribution characteristics presented in [
48]. And following from these considerations such event is virtually impossible (according to our calculations, already for
cm
−3—whereas the actual densities required would be closer to
cm
−3—the expected frequency of occurrence is
yr
−1, i.e. the corresponding time is of the order of the assumed age of the Solar System). Of course, due to the extremely short span of the observations (with respect to the assumed lifetime of the Solar System), not even including the well known extreme events of the last few centuries, such statistical estimation in meaningless. In principle, solar activity may have undergone significant changes during the evolution of the System, and we cannot exclude a higher (or even much higher) frequency of such events in the past.
Additionally, we shall note here that obviously these particular craters should not necessarily be the result of the cratering mechanism proposed in this paper. They might as well be the result of impacts. Our point, however, is to evaluate the possibility of even the largest craters being produced through the outlined mechanism.
Conversely, a strong modern day CME with a mass of
kg traveling at velocities of 1000 km/s would have total kinetic energy of
J, though only a tiny fraction of it would be accessible for cratering (determined by the size of the celestial body and its distance from the Sun). If we assume that the CME covers 1/4 of the heliosphere in directional sense (solid angle is equal to
), then at distances of the order of an astronomical unit and the celestial body size of
m, the available energy would be determined by the solid angle taken by the body as seen from the Sun:
, hence only
J is available for cratering. If all of this energy is being spent on the production of a single crater, the Formula (
8) then determines its size as
km. This is a relatively large crater, which would be easily noticeable e.g. on the Moon or Mars, the surfaces of which are under constant observation of artificial satellites.
So we must conclude that not all of the CME energy is being spent on the cratering, so only smaller sized craters (if any) are produced by modern day CMEs. Perhaps the key point is that in order for the crater to appear through the outlined mechanism, the current has to be strong enough to undergo significant pinching (see Formula (
6) and the following discussion), so that its thermal energy can be deposited in a localized manner. Hence only the stronger, fast and dense solar outbursts may lead to cratering. We have also omitted the consideration of other phenomena which may play a role in stronger CMEs: relativistic particle implantation being just one example.
On the other hand, a single energetic event may potentially produce multiple craters at once due to dendritic fractioning of the current into multiple branches (
Figure 2 obviously only shows a single cross-section, whereas the event possesses a cylindrical symmetry with respect to the X axis, if we disregard the intrinsic plasma magnetic field as a symmetry breaker), which may also cause a polygonal shape in individual craters (see
Section 5.5).
In [
49] it was proposed that a strong CME in the Sun’s past could result in reduction of magnetic field strength. The decreased magnetic field would fail to deflect charged particles from the CME and solar wind. “... the incoming gas bringing its strong field into the virtually insulating atmosphere would then result in very large electric fields so directed that the resulting currents would maintain those fields”. These fields would lead to “electrical breakdown” on each side of the incoming CME cloud. “Hundreds of millions of amperes” would flow in the crust. It’s possible that smaller currents acting in lower pressures might machine for duration, producing a crater much larger than the plasma current may imply. Mars does indeed have a global electric current [
50]. This current may have come under fluctuation in periods of increased solar activity [
51] that might produce machining in the crust. In modern conditions the most obvious examples of continuous surface erosion under the influence of electric currents are observed in the magnetospheres of large planets: with electric currents flowing between Saturn and Enceladus [
52] or Io and Jupiter [
53] etc.
We acknowledge that a model of a solar storm presented in broad strokes in [
49] in its essence is analogous to the scheme outlined in
Section 3.1, but without the more specific analysis given therein. More recently, the same scenario was proposed to explain the tectonic uplift on the Earth itself [
54], where the energies accessible to very strong CMEs [
55,
56] are transformed into heat, which alters the planetary surface.
The considerations given here, as well as the
Section 3.1 are obviously over-simplified and potentially ignore some higher order effects and dynamic phenomena (e.g. plasma waves). As we’re dealing with plasma, the actual dynamics of the process might be extremely complicated. For example, we have ignored the fact that the CME cloud should be electrically polarized due to the influence of the magnetic field of the Sun itself [
57]. Thus, an infalling plasma cloud would produce an electric field between the surface of the celestial body and its immediate surroundings (at the time of impact occupied by the cloud’s plasma), which may provoke an electric discharge between the two, as the charge carriers are provided by the cloud itself. This is another potential source for cratering, though a simplified analysis shows that it may only be significant if the magnetic field within the cloud would reach unrealistically high values.
E.g. for a typical CME cloud that is
m thick [
58] the available volume of plasma is
cm
3, and we shall estimate the energy available for cratering as
, where
Q is the available charge (
), and
U is the potential difference (as for the measure of
E, we still take the largest Solar System crater size estimate of
J, which follows from (
8). Hence
. In terms of field strength this corresponds to (at the given cloud thickness)
V/m, which should be equal to the product
:
Thus for the velocities of
m/s the product of magnetic field
B (in T) and number density
n (in cm
−3) should be of the order of 10 (typically it is
for a calm solar wind and is assumed to be around 10
−7 for a Carrington event CME—[
59]. So in our case it should be extremely high, though it is somewhat facilitated by the fact that particle density, solar wind velocity and magnetic field strength are usually connected and may scale upwards together [
48].
Likewise, we have omitted other potential sources of electric current energy deposition to the surface (such as the formation of plasmoids in the celestial body’s magnetotail and their subsequent migration towards the body itself or the interaction with neighboring bodies through solar plasma—e.g. in the case of the Moon, as described in [
60], yet all these higher order effects may play a significant role at higher plasma energies required by our model.
Another thing we’ve ignored is the energy of the magnetic field of the plasma cloud, which may also theoretically be available for cratering, as the full energy of the cloud is the sum of its kinetic, electric and magnetic energies:
if we assume
m/s. Only the first component was included in the initial estimation (
10). The reason for this is that the second one heavily depends on magnetic field strength of the plasma and is typically 2–3 orders of magnitude smaller (e.g. in the case of a calm solar wind, where
m
−3,
T), yet this may not be the case for more extreme situations (as has been demonstrated by [
59] for Carrington event).