Relating Dike Geometry and Injection Rate in Analogue Flux-Driven Experiments

Dikes feed most eruptions, so understanding their mechanism of propagation is fundamental for volcanic hazard assessment. The variation in geometry of a propagating dike as a function of the injection rate remains poorly studied. Here we use experiments injecting water into gelatin to investigate the variation of the thickness, width and length of a flux-driven dike connected to its source as a function of the injection time and intruded volume. Results show that the thickness of vertically propagating dikes is proportional to the injection rate and remains constant as long as the latter is constant. Neither buoyancy nor injected volume influence the thickness. The along-strike width of the dike is, however, proportional to the injected volume. These results, consistent with the inferred behavior of several dikes observed during emplacement, open new opportunities to better understand how dikes propagate and also to forecast how emplacing dikes may propagate once their geometric features are detected in real-time through monitoring data.


INTRODUCTION
The dominant mechanism of magma transfer in the upper crust is via magma-filled cracks, including steeply dipping dikes. Given that dikes feed most eruptions, achieving a deeper knowledge of dike propagation is important in eruption forecasting (e.g. Rubin, 1995;Rivalta et al., 2015).
Several studies have addressed the general features related to the propagation of magma-filled cracks (Pollard, 1987;Lister and Kerr, 1991;Rubin, 1995;Rivalta et al., 2015) or more specific features investigated through field data (e.g. Gudmundsson, 2002;Daniels et al., 2012), numerical modeling (e.g. Pollard and Muller, 1976;Rubin and Pollard, 1987;Traversa et al., 2010) and laboratory experiments (Takada, 1990;Menand and Tait, 2002;Taisne et al., 2011;Kavanagh et al., 2018;Urbani et al., 2017;Urbani et al., 2018). However, any variation in the geometry of a propagating dike as a function of the injection rate and of the pressures acting on it remains poorly studied. Dike propagation is indeed a complicated phenomenon with several processes storing and dissipating energy (Lister and Kerr, 1991;Menand and Tait, 2002;Rivalta et al., 2015). Although there are several external factors potentially controlling dike propagation, such as tectonic stresses or topographic loads, the main pressures acting on a propagating dike fed by its source are the excess or elastic pressure P e , the source pressure ΔP r , the buoyancy pressure P b , the viscous pressure drop P v , and the fracture pressure P f (Pollard, 1987;Lister and Kerr, 1991;Rubin, 1995;Menand and Tait, 2002;Canon-Tapia and Merle, 2006;Taisne et al., 2011;Kavanagh et al., 2013;Rivalta et al., 2015). Several studies show that these pressures influence many features of a propagating dike, including its shape and geometry. As for the shape, when a propagating dike reaches a critical length (parallel to its direction of propagation), related to the above-mentioned pressures, it develops an inflated head, followed by a thinner tail. As for the geometry, these pressures may affect the length, width and thickness (i.e. the opening) of the dike (see Rubin, 1995;Rivalta et al., 2015 for details). In particular, the geometric parameters of a dike have been investigated using field data (e.g., Kavanagh and Sparks, 2011;Daniels et al., 2012;Geshi et al., 2020), analogue models (e.g., Takada, 1990;Menand and Tait, 2002;Taisne et al., 2011;Kavanagh et al., 2018) and numerical models (e.g., Traversa et al., 2010). However, some questions remain poorly studied and understood. A first one concerns the relationship between the thickness and the injection rate during the propagation of a dike connected to its source (flux-driven dike), which remains elusive. For example, some studies assume that a dike maintains a near constant thickness during its propagation (e.g. Aloisi et al., 2006). Although analogue models and the inversion of deformation data seem to support this assumption (Heimpel and Olson, 1994;Taisne et al., 2011;Grandin et al., 2010;Segall et al., 2001;Morita et al., 2006;Aloisi et al., 2006), this aspect has never been tested in thorough targeted analyses or laboratory experiments. A second question concerns the possible effect of the injected volume and injection rate on the width ( Figure 1E) of a propagating dike. Laboratory models suggest that a buoyancy-driven dike of finite volume not connected to the source quickly reaches a constant width, proportional to its volume (Taisne et al., 2011). Conversely, laboratory experiments modeling a flux-driven dike fed from its source show an increase in width with time and with volume (Kavanagh et al., 2018); nevertheless, the possibility that different injection rates determine different dike widths remains unclear.
Achieving a better understanding of the 3D geometry variations of a flux-driven propagating dike may have important implications in volcanology. For example, prompt detection of the geometric features of a propagating dike through monitoring data may allow estimating the associated energy, providing promising opportunities to expect how far a dike may propagate, ultimately enhancing eruptive hazard assessment (e.g. Bonaccorso et al., 2017).
In this study, we have investigated the variation of the geometric parameters (length, width and thickness) of a continuously fed flux-driven dike as a function of the injection rate. For this goal, we have carried out analogue experiments by injecting water into gelatin.

METHODOLOGY Experimental Set-Up
In our experiments, we injected water (magma analogue) into a pigskin gelatin (type A gelatin, purchased from Italgelatine S.p.A.), selected as the crust analogue material (Kavanagh et al., 2013;Brizzi et al., 2016). The pigskin gelatin solid was prepared at 2.5 wt% by dissolving the appropriate amount of gelatin powder into water at 80°C. We also dissolved 10 wt% of NaCl, so that the gelatin rigidity decreased and the density increased (Brizzi et al., 2016), providing an adequate scaling (see Section Scaling). We poured the solution into a 30 cm diameter cylindrical Plexiglas ® container and filled it up to a height of 30 cm. We cooled the solution in a refrigerator at 8°C for 19 h.
We carried out six experiments injecting water (20°C) at a constant flux from the base of the tank via a tapered-injector using a peristaltic pump. The orientation of the tapered needle allowed us to control the initial orientation of the crack. We used a different injection rate (Q) for each experiment. Q, here used as a proxy of the source pressure (see Section Pressures for details), is the only tested parameter, to understand the variations in the dike geometry during propagation. We focused on the injection rate also as this is an easy and reliable parameter to control and measure with a peristaltic pump.
We monitored the experiments from the side views using two digital cameras, one placed in front of the dike, along its strike, and the other at 90°to the side, acquiring images every 5 s.

Pressures
Following Menand and Tait (2002; and references therein), we consider five pressures to quantitatively characterize our experiments, where a dike is constantly fed. These are: 1) The elastic (or excess) pressure P e required to contrast the stress perpendicular to a crack and to keep it open. This is expressed by a quasi-static 2D equation (Pollard, 1987;Rubin, 1995) where h is the thickness of the dike (i.e. its tensile opening), L is the length of the dike in the direction of propagation (see Figure 1E), G is the shear modulus and ν the Poisson's coefficient.
2) The source pressure ΔP r , which is the pressure of the source feeding the dike. Menand and Tait (2002) proposed the following experimental equation, valid for constantly fed dikes (Menand and Tait, 2002): Where E is the Young modulus of the host rock, Δρ is the density contrast between the host rock and the magma, g is the acceleration due to gravity, Q the volumetric flux of magma injected into the dike, ν the Poisson's ratio and u is the average velocity of magma inside the dike (i.e. the dike velocity). This formulation has the advantage to relate ΔP r to the injection rate (Q) and other parameters (u, Δρ), which can be experimentally easily controlled and measured.
3) The buoyancy (or hydrostatic) pressure P b , due to the density contrast (Δρ) between the host rock and the magma, defined as (Kavanagh et al., 2013): 4) The viscous pressure drop (P v ), which is the pressure drop due to the viscous flow of the magma within the dike, and is usually negligible in the head of the dike and dominant in the tail (Tait and Taisne, 2013): where η is the magma viscosity. 5) The fracture pressure (P f ), which is the pressure required to propagate the dike tip and is usually the main resistive pressure acting on the head of the dike (Lister and Kerr, 1991;Menand and Tait, 2002;Rivalta et al., 2015); this is expressed as: where K c is the critical fracture toughness. During dike propagation, ΔP r and P b are the main driving pressures acting on the dike, with their sum defined fluid overpressure by Menand and Tait (2002).
The buoyancy pressure (Eq. 3) is the same in all experiments, as we did not change Δρ, whereas the source pressure, directly related to the injection rate Q (Eq. 2), changes with Q in each experiment ( Table 1).

Scaling
Model parameters have to be geometrically, kinematically and dynamically scaled, in order to ensure similarities between natural prototypes and experimental results (Merle 2015). Kavanagh et al. 2013 showed that experiments of water injection into gelatin can be appropriately scaled to model dike propagation in nature. The characteristic length scale of the experiments is the buoyancy length L b , which represents the length over which magma buoyancy driving ascent balances the resistance from rock fracture, defined as (Kavanagh et al., 2013): Where K c is the fracture toughness of the host material, Δρ is the density difference between the intruded fluid and the host   Table 2. Q is the injection rate used in each experiment. ΔP r , P b , P e are the pressures defined by Eq. 1-3, calculated in the models at a length L L b . P 1 * is the pressure ratio defined in Section Scaling. For natural dikes we used L∼5,000 (corresponding to L b , see scaling) and h 1.5 m. ΔP r + P b are the fluid overpressures (Menand and Tait, 2002).
Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 665865 material and g is the acceleration due to gravity (9.8 m/s 2 in both model and nature). Under the applied experimental conditions (see Section Experimental Set-Up), the mechanical properties of the gelatin and of the injected water, as well as the values that we assume representative of the natural dikes, are reported in Table 2. Therefore, the length scale L * b (the asterisk denotes the model to nature ratio) is: imposing that 1 cm in the experiments corresponds to ∼200 m in nature. From the equations proposed by Kavanagh et al. 2013, we calculated also the time scale T* and the velocity scale U*. T* 7.8*10 -3 ; U* 6.2*10 -3 . Finally, for the pressure scale we propose both the P e * ( Table 1) and an adimensional ratio P 1 * between the fluid overpressures (ΔP r + P b ; Menand and Tait, 2002) and the resistive pressure P * 1 (ΔP * r + P * b )/(P * v + P * f ) As some of these pressures are related to the length L of the dike (see Eqs. 1-5), we calculated them for the experiments at a length equal to L b (23.8 cm). In Table 1 we report the P * 1 value for each experiment. The values approach 1, implying an appropriate scaling.

Method Limitations
In our experiments we made some assumptions. We neglected solidification effects (Chanceaux and Menand, 2016) and used a single-phase system (water) as analogue of magma, which is usually a three-phase system (liquid + gas + crystals). In addition, the scaling highlights how our experiments simulate the propagation of low viscosity dikes. Therefore, our results are mainly applicable to mafic dikes, widespread in different tectonic settings. Moreover, during each experiment we assumed a constant injection rate, which in nature can change in time (Rivalta, 2010;Rivalta et al., 2015). As we wanted to test the role of the injection rate on the geometry of the propagating dike, we did not change any other possible parameter affecting the dike geometry, such as the viscosity of magma and the mechanical properties of the host rock, which may also vary in nature (Lister and Kerr, 1991;Rivalta et al., 2005;Rivalta and Dahm, 2006;Rivalta et al., 2015).

Image Analysis
We used the Matlab ® image toolbox functions to measure the parameters of the dike such as the length (L), the width (W), the area and the dip (Figure 1). For this, we first subtracted each image with the dike (e. g. Figures 1A,B) to the image without the dike (the master image acquired at time t 0) in order to obtain an image in which the pixel value is zero (black) in the area without the dike and a positive value (gray) in the area with the dike (Figures 1C,D). Then, this image is converted into a binary image and we measured the parameters of the dike using the Matlab ® function regionprops ( Figures 1E,F). The measured length (L, Figure 1E) was then corrected for the dip of the dike (α).

Estimate of the Average Thickness from the Volume/Area Ratio
With the method in section Image Analysis we calculated the area (A) of the dike for each image, which we corrected for the dip of the dike (α). As for each image we also know the associate injected volume (V i ), we measured the average thickness of the dike from the V i /A ratio.
However, this method has a limitation. As soon as the dike length (L, parallel to the direction of propagation) overcomes a critical value, which depends on the driving and resistive pressures, the dike develops an inflated head followed by a thinner tail (see Rubin, 1995;Rivalta et al., 2015 for details). Once the experimental dike shows a well-developed tail in the second part of each experiment (see Figures 2E,G and Result Section), our method progressively led to an underestimate of the average thickness (Supplementary Figure S2). Indeed, the estimated area of the dike also included the area of the tail, but the injected volume mainly lied within the head ( Figures  2E,G). Therefore, to better estimate the average thickness of the dike with a tail, we propose a further method based on the ratio V head /A head , where A head and V head are respectively the area and the volume and of the dike head. We estimated the area and the volume of the head from the black and gray images ( Figure 1C) using the following approach. In the images from the first half of the experiment, with no or a poorly developed tail, the intensity value of each gray pixel (I) associated with the dike showed just a limited variability. As in these images the average thickness (h) estimated with the method V i /A was almost constant (see result section), we associated that thickness to the average value of the gray intensity (I a ). Knowing the area of the pixel (A pixel ), corrected for the dike dip, we calculated the average volume (V pixel_mean ) associated with the average value of the gray intensity using the formula V pixel_mean A pixel *h. Then, we associated the gray intensity (I) of each pixel to a volume V x using equation V x (I*V pixel_mean )/I a . To test the validity of this method and its assumptions, we recalculated the total volume of the dike by summing the V x estimated for each pixel and we compared this value with the real injected volume. As these values were almost equal within a reasonable error (<5%), this method proved reliable. At this point, we applied the calibrated gray scale to calculate the volume within the head of the dike in the images where the dike showed a head and a tail. The area of the head was selected in the summit region in which the gray intensity of the pixels was similar to that shown before the formation of the tail. From the V head /A head ratio we calculated the average thickness of the head.

Models Description
We briefly describe the first experiment (Dike 1), whose general evolution is representative of all the experiments. When we switch the peristaltic pump on, a sub-vertical water-filled crack forms and propagates. The propagation of our analogue dike can be schematically divided into two main stages, similar to those identified by previous studies with analogous set-up (Menand and Tait, 2002;Kavanagh et al., 2018). At the beginning (first phase), the crack propagates mainly radially in a sub-vertical plane, with sub-circular shape and the length almost equal to the width (Figures 2A,B, 3). In this phase, the crack lacks any preferential direction of propagation (see Figures 2A,B, 3). Then, in the second phase, the crack starts propagating mainly vertically and the length (L) increases faster than the width (W), with the crack assuming an elliptical/inverted teardrop geometry ( Figures 2C-H, 3). After few tens of seconds from the onset of the second phase, a thicker head and a thinner tail are clearly identified from the side and front view images ( Figures 2E-H).
The length (L) of the dike at which the transition from the first to the second phase occurs is proportional to the injection rate of the experiment (Figure 3) and ranges from ∼10 cm in the experiment at the lowest injection rate to ∼16 cm in the experiment at the highest rate ( Figure 3). In addition, the front and side view images suggest that in the higher Q experiments the amount of magma left over in the tail increases (see images in the data repository Galetto, 2021). Figure 3 reports the temporal variation of the average thickness of the crack in each experiment. For the first and second phase, we used the method based on the Volume/Area ratio. At the onset of injection (beginning of phase one), when the dike nucleates, the thickness quickly increases with time. Then, the thickness remains approximately constant until the end of the experiment. After few tens of seconds from the onset of the second phase, when the images of the experiments show a thicker head and a thinner tail, we start measuring the thickness with the method V head /A head based on the color intensity to avoid the possibility that the tail overestimates the area, leading to an underestimate in the thickness (which would linearly decrease with time). With the V head /A head method, we find that the thickness of the head remains almost constant until the end of the experiment. However, there is a small drop in the average thickness at the transition between the two methods, when the tail becomes identifiable; this drop is negligible in the three experiments at lower Q, whereas it is slightly more evident in the three experiments at higher Q.

Dike Thickness (h)
The results highlight two main aspects related to the thickness of the dike: i) this is proportional to the injection rate, here used as a proxy of the source pressure (ΔP r ) ( Figure 4A and Table 1); ii) after a short initial stage of growth, the thickness remains approximately constant. In addition, Figure 4B shows that in each experiment the thickness/length ratio (h/L) and the related P e (see Eq. 1), calculated for the same L of the dike, are also directly proportional to the injection rate. According to Pollard (1987) and Rubin and Pollard (1987), we calculated P e using the length L parallel to the direction of propagation of the dike. Some Authors (e.g. Lister and Kerr, 1991) proposed to consider the shortest dimension between L and W (which would be W in our models), and using the h/W ratio we obtain the same linear correlation with P e (see Supplementary Figure S3). Finally, the P e calculated in each experiment for equal intruded volumes (but different length) is again proportional to Q ( Figure 4C). Figure 3 shows the temporal variation of the dike width W in each experiment. During the first phase, the width increases Results highlight that the dike width is proportional to the injected volume, whereas there is no clear relationship with the FIGURE 3 | Temporal variation of the width (W), length (L) and thickness (h) in each experiment arranged (A-F) with increasing injection rates. The V/A Method is based on the ratio Volume/Area applied in the first phase and at the beginning of the second phase. Then, we apply the method V head /A head based on the ratio Volume of the head/Area of the dike head.

Dike Width (W)
Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 665865 6 variation in the injection rate ( Figure 5). Therefore, similar dike widths are reached with the same injected volumes, even at different injection rates.
Finally, the total volume injected when the dike erupts is proportional to the injection rate ( Figure 6A). Therefore, the final maximum width of the dike, when this erupts (W er ), is also somehow proportional to the injection rate ( Figure 6B), as a higher rate will generate a dike with larger volume that in turn is related to a larger width.

General Features
Our models reproduce near vertical dikes fed by a constant injection rate (Q). During the propagation, we observe different relationships between the geometry of the dike (L, W, and h) and the injection rate. During the first phase, the length and width of the dike grow about at the same rate and the thickness soon approaches a constant value. Then, during the second phase, the dike starts propagating vertically faster than laterally, while the thickness remains almost constant, although there may be a minor drop when the tail appears ( Figure 3). These behaviors are explained considering the pressures acting on the dike ( Figure 7A) (Menand and Tait, 2002). During the first phase, the resistive pressures, mostly related to the fracture pressure P f , are higher than the fluid overpressure (P b + ΔP r ) driving dike propagation ( Figure 7B) and indeed the dike grows radially, without any preferred direction. At the beginning of the second phase, the fluid overpressure approaches and then overcomes P f , mainly because of the increase of P b, that becomes the main driving component, and the dike starts  propagating mostly upward ( Figure 7B; Menand and Tait, 2002). This explains why L increases faster than W (Menand and Tait, 2002). The progressive decrease of P f and P e ( Figure 7A) is expected by the linear elastic theory, as, the longer the dike, the lower the energy P f required to propagate its tip, and the lower the elastic energy P e required to keep the fracture open (e.g. Rubin and Pollard, 1987;Rubin, 1995). Figure 7A also shows that during the propagation of a constantly fed dike the viscous pressure drop is a negligible resistive pressure at the dike head, in agreement with theoretical models (Lister and Kerr, 1991;Tait and Taisne, 2013). Finally, the small drop in the thickness of the dike head at the appearance of the tail (Figure 3) may be related to the amount of fluid left over in the tail (Rivalta et al., 2015). This is supported by the fact that this drop is negligible in the experiments at lower injection rate, where the dike tail is almost closed, while it becomes more evident in the experiments at higher rate, with thicker tail (see data in Galetto, 2021). The fact that the amount of magma left over in the tail increases with the injection rate Q, as suggested by the front and side images, is expected by fluid-mechanical models of crack propagation (see Lister and Kerr, 1991 for a detailed analysis).

New Insights on Dike Thickness and Width
Our experiments show a different behavior of the thickness and the width of the dikes with respect to the injection rate and volume change. A first important aspect is that the thickness h tends to remain approximately constant during dike propagation. The fact that after the minor drop at the appearance of the tail the thickness of the head remains almost constant also confirms that under the applied experimental conditions the thickness does not change in time. Therefore, the total injected volume, which increases with time, does not seem to influence the thickness. On the contrary, the thickness is directly proportional to the injection rate Q, here used as a proxy of the source pressure ΔP r ( Figure 4A and Table 1; Menand and Tait, 2002). FIGURE 7 | A) Variation of the pressures as a function of the length of the dike in experiment "Dike 1" (Q 0.162 ± 0.008 cm 3 /s). P e , P v , P b , P f , and ΔP r are the elastic pressure, the viscous pressure drop, the buoyancy pressure, the fracture pressure and the source pressure, respectively (see equation 2-5). (B) Variation of the fluid overpressure (P b + ΔP r ) and of the resistive pressure (P f + P v ) as a function of the dike length in "Dike 1".
Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 665865 Another important point is that the buoyancy pressure (P b ) does not affect the dike thickness. Indeed P b (see Eq. 3) increases with the length of the dike in each experiment ( Figure 7A), whereas the thickness remains constant. The experiments also highlight that a higher Q (i.e. higher ΔP r ) generates a higher h/L ratio and therefore higher P e ( Figures  4A,B), in agreement with the linear elastic theory (Pollard, 1987;Rubin, 1995).
A novelty revealed by the experiments is that, contrary to the thickness, the dike width is proportional to the injected volume, and therefore increases over time, whereas there is no clear relationship with the injection rate ( Figure 5). This new aspect shows that similar crack widths are obtained when same crack volumes are injected, even with different injection rate. Under the applied experimental conditions, the volume of the dike and the P b progressively grow. Taisne et al. 2011 found a relationship between the buoyancy force, the volume, and width of the dike (see equation 12 in Taisne et al., 2011) for buoyancy-driven dikes. Their equation derives from the balance between the fracture and buoyancy pressures, as their buoyancy-driven dikes resulted from a limited volume of magma rising only for buoyancy. By using their equation, we may estimate the contribution of the buoyancy to the dike width in our models. The difference between our measured dike width and the buoyancy width estimated following Taisne et al. 2011 should be the contribution to the width due to the constant injection (Figure 8), which is the only driving force in our models that is absent in Taisne et al. 2011. This suggests that even if the width of the dike mainly depends on the injected volume, for equal injected volumes the two widths may differ depending on whether the dike is continuously fed (flux-driven dike) or not (buoyancy-driven dike) from its source.

Comparison to Natural Dikes and Implications
According to our results, a constantly fed dike increases its width over time, but not its thickness, if the injection rate remains constant. On the contrary, if the injection rate increases, the thickness also increases. Thus, during dike propagation any volume change does not necessarily imply a variation in all the geometric parameters of the dike.
A first main novelty of our experiment regards the dike thickness. Our results confirm that thicker dikes are related to higher source pressures, and therefore they will release more energy. This agrees with Bonaccorso et al. (2017), who found a simple but effective equation that relates the squared value of the dike thickness with the expected released mechanical energy and, in turn, to the recorded seismic moment produced by the earthquakes associated with dike propagation. These Authors showed that the available mechanical elastic strain energy has to be entirely released seismically, and this can be preliminarily estimated by the dike thickness. This means that the dike thickness is the most important and unequivocal parameter related to the energy to be released and, therefore, from the comparison between the expected energy and the recorded released seismicity, it may possible to evaluate in real-time the state of the propagating dike (Bonaccorso et al., 2017). This general approach was retrospectively applied to two recent dike intrusions feeding flank eruptions at Etna volcano on October 2002 (Bonaccorso et al., 2017) and December 2018 (Bonaccorso and Giampiccolo, 2020). In this rationale, a crucial aspect concerns the hypothesis that the thickness of the dike remains approximately constant during its propagation, which is a feature that has never been proven so far and that is now confirmed by our models (if injection rates do not change significantly). Furthermore, it is also interesting to note that the modeling of the continuous ground deformation data (borehole tilt and GPS) modeled over time of the 2002 (Aloisi et al., 2006) and 2018 eruptions (Aloisi et al., 2020) suggested a near constant thickness during the dike propagation, as expected by our experiments.
Our results show that a constantly fed dike should propagate varying its width and volume, but not the thickness. Conversely, any increase in the thickness would imply an increase in the injection rate. Nevertheless, if the width of a propagating dike does not change significantly, the dike would not be continuously fed from its source and would be propagating only through buoyancy (Taisne et al., 2011). Some examples of well-studied propagating dikes in nature support these findings. For example, at Cerro Azul (Galápagos) in 2008 a propagating radial dike underwent a sudden increase in the injection rate and volume, also increasing its thickness and width (Galetto et al., 2020). Similarly, the increase in width and thickness of the 2014 Bárðarbunga dike (Iceland) mainly occurred at the onset of its propagation, when the increase in volume change and injection rate were higher (Sigmundsson et al., 2015;Heimisson and Segall, 2020).
Finally, our results may allow reinterpretation of previous eruptive events. In fact, the relationship between the thickness FIGURE 8 | Variation of the dike width (W) with the intruded volume (V) in our models (flux-driven dike; blue crosses) and as predicted by Taisne et al., (2011) for buoyancy driven dikes (orange triangles).
Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 665865 and injection rate in flux-driven dikes implies that the greater the thickness the higher the injection rate of the dike and, as the injection rate is a proxy of the source pressure (Menand and Tait, 2002), the higher the source pressure. This relationship is particularly useful, as direct estimation of the source pressure feeding a dike is difficult (if not impossible) and even the estimation of injection rate may be difficult if geodetic data are not detected in a continuous mode. For example, the 2005-2009 sequence of rifting at the Manda-Hararo Rift (Afar, Ethiopia) was triggered by the intrusion of at least thirteen dikes. Although the 3D geometry has been estimated for all these dikes (Hamling et al., 2009;Grandin et al., 2010), their injection rates remain unknown. Using our approach, it may be possible to establish a hierarchy of the dikes fed by the highest and lowest injection rates (and source pressures), simply based on their thickness. These possible in-depth studies are beyond the scope of our work, although we believe they may be useful to further explore the potential for applying the results of our experiments.
Our results highlight the importance of studying the 3D geometry variation during dike propagation. Future studies should investigate whether any other parameter (e.g. viscosity, Young's modulus and cooling) may also affect the geometry of a propagating dike.

CONCLUSION
Our models show that in experimental flux-driven dikes: 1) Dike thickness is related to the injection rate and therefore to the source pressure ΔP r , whereas it does not seem to be related to buoyancy pressure. 2) Dike thickness tends to remain constant if injection rate does not change. 3) Higher injection rates (i.e. higher ΔPr) generate higher h/L ratios and therefore higher elastic pressures ( Figures 4B,C), as expected from the linear elastic theory. 4) Similar crack widths are obtained when same volumes are injected, even with different injection rates.
These results increase our understanding of the 3D geometric variation during the propagation of a dike continuously fed from its source and are consistent with the behavior inferred from observations of several historically emplaced dikes. In addition, our results offer new opportunities to better understand how dikes propagated during previous events and also how propagation will proceed for dikes with geometry delineated by real-time monitoring during an emplacement event.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Images of the analogue experiments are available in the data repository https://doi.org/ 10.17605/OSF.IO/Z42PH (Galetto, 2021). Further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
AB conceived this work and traced its scheme. FG carried out the experiments and processed the data. VA supervised the experiments. All the authors contributed with original ideas to the drafting of the article and shared the discussion of the results.

ACKNOWLEDGMENTS
We thank Fabio Corbi, Stefano Urbani and Eleonora Rivalta for their help with the experimental set-up and the scaling. FG is grateful to Matteo Galetto for assistance in developing the method based on the color intensity and for developing the associated Matlab ® code. FG was supported by the INGV project Premiale-MIUR acronym "Transienti". We thank S. Conway for improving the quality of the English. We thank the Editor J. White and two anonymous Reviewers whose comments significantly improved the quality of this manuscript.