Analogue modelling (geology)

Analogue modelling is a laboratory experimental method using uncomplicated physical models (such as a sandbox) with certain simple scales of time and length to model geological scenarios and simulate geodynamic evolutions.[1][2]

Pure shear sandbox model of thrust fault formation

There are numerous limitations affecting the direct study of the Earth. Firstly, the timescales of geodynamic processes are exceptionally long (millions of years), and most of the processes started long before human records.[1][3] Secondly, the length scales of geodynamic processes are enormous (thousands of kilometres), and most of them happen at depth within the Earth.[1][3] Thus, scientists began making proportional small-scale simulations of features in the natural world to test geological ideas. Analogue models can directly show the whole structural pattern in 3D and cross-section. They are helpful in understanding the internal structures and the progressive development of Earth's deforming regions.[1]

Analogue modelling has been widely used for geodynamic analysis and to illustrate the development of different geological phenomena. Models can explore small-scale processes, such as folding and faulting, or large-scale processes, such as tectonic movement and interior Earth structures.[1][4]

History

Lateral compression machine made by James Hall to model geological folding. This machine is still present in the Royal Society of Edinburgh. The materials squeezed in the box are blankets or layers of clay.[2]

Analogue modelling has a development history of over 200 years.[1]

It has been used since at least 1812, when James Hall squeezed layers of clay to produce folds similar to those that he had studied at an outcrop.[2] This idea of modelling led to many other small-scale studies, such as fault-propagation fold,[5] thrust fault,[6] and folds[7] in the late 19th century. These studies were all qualitative.[1]

King Hubbert came up with the theory of scaling in 1937, meaning that the study of analogue modelling became quantitative.[8] The quantitative approach was further developed by many scientists later.[1] As the field of geodynamic study expanded, analogue modelling increased, especially for large-scale geological processes. Examples include proto-subduction[9] subduction[10][11] in plate tectonics, collision,[12] diapirism,[13] and rifting.[14][1][4]

Components

Analogue model of caldera formation using flour to represent the upper part of the crust and a balloon to represent the inflating magma chamber

Scaling

In 1937 King Hubbert described the key principles for scaling analogue models. He defined three types of similarity between models and the natural world: geometric, kinematic and dynamic.[8][15]

Geometric similarity

To be geometrically similar, lengths in the model and natural example must be proportional and angles must be equal.[15] When the length of a natural prototype (p) is ${\displaystyle l_{n}^{p}}$ (n=1, 2, 3...) and the angle is ${\displaystyle \alpha _{n}^{p}}$. Correspondingly, the length in the model (m) is ${\displaystyle l_{n}^{m}}$ and the angle is ${\displaystyle \alpha _{n}^{m}}$. They need to conform to the following formulas:[1]

${\displaystyle {\frac {l_{1}^{m}}{l_{1}^{p}}}={\frac {l_{2}^{m}}{l_{2}^{p}}}={\frac {l_{3}^{m}}{l_{3}^{p}}}={\frac {l_{n}^{m}}{l_{n}^{p}}}}$ & ${\displaystyle \alpha _{n}^{m}=\alpha _{n}^{p}}$

For example, 1 centimetre in the model represents 1 kilometre in nature.

Kinematic similarity

To be kinematically similar, they must be geometrically similar and the time needed for changes to occur must be proportional.[9] When the required time for changing is ${\displaystyle t_{n}}$:[1]

${\displaystyle {\frac {t_{1}^{m}}{t_{1}^{p}}}={\frac {t_{2}^{m}}{t_{2}^{p}}}={\frac {t_{3}^{m}}{t_{3}^{p}}}={\frac {t_{n}^{m}}{t_{n}^{p}}}}$

For example, 1 second in the model represents 1 thousand years in the nature.

As is known: ${\displaystyle v={\frac {l}{t}}}$, the velocities (${\displaystyle v}$) can be scaled by the following equation:[1]

${\displaystyle v^{p}=v^{m}{\frac {l^{p}t^{m}}{l^{m}t^{p}}}}$

Dynamic similarity

When the models and the natural world are geometrically and kinematically similar, dynamic similarity additionally requires that the various forces acting on a point in the model are proportional to those at a corresponding point in nature.[15] When the forces (${\displaystyle F_{n}}$) acting on the system are ${\displaystyle F_{g}}$ (gravity), ${\displaystyle F_{v}}$ (viscous force), and ${\displaystyle F_{f}}$ (friction):[15]

${\displaystyle {\frac {F_{g}^{m}}{F_{g}^{p}}}={\frac {F_{v}^{m}}{F_{v}^{p}}}={\frac {F_{f}^{m}}{F_{f}^{p}}}={\frac {F_{n}^{m}}{F_{n}^{p}}}}$

However, since the forces acting in the nature are unmeasurable, it is impossible to scale the forces and stresses directly. Scientists have been using different formulas to convert forces into the parameters that can be measured. Cauchy momentum equation is usually used for showing the relationship between forces and densities (${\displaystyle \rho }$ is density):[1]

${\displaystyle {\frac {F^{m}}{F^{p}}}={\frac {\rho ^{m}(l^{m})^{3}}{\rho ^{p}(l^{p})^{3}}}}$(Generating from Cauchy momentum equation[16])

Stokes' law is usually used for showing the relationship between forces and density contrasts (${\displaystyle \Delta \rho }$ is density constant):[1]

${\displaystyle {\frac {F^{m}}{F^{p}}}={\frac {\Delta \rho ^{m}(l^{m})^{3}}{\Delta \rho ^{p}(l^{p})^{3}}}}$(Generating from Stokes' law[17])

(While the gravitational acceleration ${\displaystyle g^{m}=g^{p}}$)

Since the densities and density contrasts are proportional to forces and stresses, it is easy to scale densities or density contrasts instead of scaling forces and stresses.[1]

However, these two equations can lead to different topography scales.[1]

The simple analogue modelling of the fold-thrust belt system taken in the Nanjing University physical laboratory. The white and blue materials are quartz sands.

Experimental apparatus

Different geodynamic processes are simulated by different experimental apparatus.

For example, lateral compression machines are commonly used in simulating deformations involving lithospheric shortening, such as folding,[2] thrust faulting, collision, and subduction. Longitudinal compression machines are usually used for fracturing.[18] There is a large variety of devices based on the different sources of forces applied to the material. Some devices have multiple forcing systems because nature is not homogeneous.[1]

Systems

For experimental systems, the energy can be supplied externally (at the boundary) and internally (buoyancy forces). If the deformation is only caused by internal forces, it is a closed system. Conversely, if the deformations are caused by external forces or a combination of internal and external forces, it is an open system.[1]

For the open system, the extrusion or stretching forces are imposed externally. However, the buoyancy forces can be generated both externally or internally. The materials and thermal energy can be added to or remove from the system. For the closed system, there is no energy and materials added to the system. Thus, all the deformations are caused by internal buoyancy forces. Only buoyancy-driven deformation can be simulated in a closed system.[1]

Gravity field

The simple analogue modelling of a subduction zone. The materials this model uses are sand mixture and silicone putty for the continental crust (left in layered brown) and oceanic crust (right in layered brown), and glucose syrup for the asthenosphere (greenish-blue liquid in the glass tank). There is a heater in the tank for heating the liquid.[2][19][20]

Because the major research object of analogue modelling is Earth, the gravity field that most experiments utilize is ordinarily the Earth's field of gravity. However, many models are carried out using simulated gravity, for example, by use of a centrifuge. These technologies are usually used in studying development of gravity-controlled structures, such as dome formation,[21] and diapirism.[1]

Materials

Analogue modelling uses various materials, such as sand, clay, silicone and paraffin wax.[2] Different materials came into use for quantitative analysis of analogue modelling experiments, compared to qualitative analysis.[22] Before Hubbert's scaling theory, scientists used natural materials (e.g. clays, soil, and sand) for analogue modelling.[1] For large-scale simulation, analogue modelling should have geometric, kinematic, and dynamic similarity with nature. If the model has these similarities, results from simulation will be more accurate.[8] All these different materials represent the natural features of Earth (such as crust, mantle, and river).[22] Selection of analogue materials is difficult, because of the largely rheology-dependent deformation and inconstant rheology influenced by the thermal gradient in nature. The rheological characteristic of internal layering was developed by the study of seismology and geochemistry.[1]

To simulate layers with different properties, different materials are chosen:

Materials for analogue modelling[1]

Categories		Examples	Simulation
Granular materials (various in density, shape, and size)		Quartz sand, glass microbeads, feldspar powder	Brittle upper crust[8]
Low-viscous materials		Water, sugar solution, honey	Asthenosphere, Sub-lithospheric mantle
		Corn syrup, glucose syrup	Sinking slabs[23]
High linear viscous materials		Syrup, silicone putty	Ductile lithosphere
Visco-elastic materials		Amorphous polymers, biopolymers, bitumen
Non-linear viscous materials	Plastic materials	Plasticine
	Visco-plastic materials	Wax, paraffin
	Visco-elasto-plastic materials	Gelatin

Advantages

There are many useful properties of analogue modelling:

1. Analogue models can directly show whole geodynamic processes from start to finish.[1]
2. Geodynamic processes can be stopped at any time for investigation, and allow the study of 3D structures.[24]
3. The scales of the model can be controlled in a practicable range for the laboratory.[1]
4. The simulation can show different results of geodynamic processes by altering the parameters, and the influence of each parameter is clarified.[24]
5. The results of analogue modelling can be directly used for interpreting nature if the accuracy of the model is high.[1]
6. Analogue modelling can provide new ways of thinking about geological problems.[24]

Disadvantages

Because analogue modelling involves the simplification of geodynamic processes, it also has several disadvantages and limitations:[15]

1. The study of natural rock properties still needs more research. The more accurate the input data, the more accurate the analogue modelling.[15]
2. There are many more factors in nature that affect the geodynamic processes (such as isostatic compensation and erosion), and these are most likely heterogeneous systems. Thus they are challenging for simulations (some factors are not even known).
3. The variation of natural rocks is greater than in simulated materials; therefore it is difficult to fully model the real situation.[15]
4. Analogue modelling cannot simulate chemical reactions.[15]
5. There are systematic errors in the apparatus, and random errors due to human factors.[1]

Applications

Simple analogue modelling of the growth and erosion of an orogenic wedge. This simulation is done in a glass tank, with layered different granular materials that represent the crust.

Analogue modelling can be used to simulate different geodynamic processes and geological phenomena, such as small-scale problems – folding, fracturing, boudinage and shear zone, and large-scale problems – subduction, collision, diapirism, and mantle convection.[1][4] The following are some examples of applications of analogue modelling.

Compressional tectonics

The first analogue model was built by James Hall for simulating folds. He used a lateral compression machine for the simulation, and this machine is still shown in the Royal Society of Edinburgh.[2] The final result from the model is quite close to observation of the Berwickshire coast.[2] Although the model he used is simpler than current ones, the idea remains in use.

The use of more complex compression machines substantially increases the number of simulations of compressional tectonics, including subduction, collision, lithospheric shortening, fracture formation, thrust and accretionary wedge. If the simulation only focuses on the upper crustal, the model is always built in the glass box (or two lateral glass walls) with a piston and/or wedges to supply forces to layers of granular materials (normally called sandbox). Depending on the different natural features, erosion (removal of top materials at a certain angle), décollement (inserted layers with low cohesion, normally glass microbeads), and any other parameters can put into the model, producing various results.[25]

Simulations of mantle influences vary. Because of the different physical and chemical properties between the asthenosphere and lithosphere, viscous materials and a heater (for mantle convection) are also used.[2]

The simple analogue modelling of the extension tectonics which showing the formation of normal fault and salt dome (diapirism). This model is built in a glass box. The darker greyish layer is silicone which represents salt, and brownish layers are dry quartz sands which represent the brittle sedimentary rocks.[13]

Extensional tectonics

Compression machines can also be used in reverse for simulating extensional tectonics, such as lithospheric extension, the formation of rifts, normal faulting, boudinage and diapirs. These models can also be built in a glass box which is similar to the above, but instead of thrust force, tensile force is applied.[13]

Strike-slip tectonics

The simplified analogue modelling setting of shear deformation. This model is built on two separate horizontal plates. The brownish layers are dry sand, wet clay, and viscous materials, such as silicone or polydimethylsiloxane.[26]

Strike-slip tectonics differ from the dominantly vertical crust movements associated with shortening and extension, being dominantly horizontal in character (in relative terms sinistral or dextral). This kind of horizontal movement will create a shear zone and several types of fractures and faults. A typical model used for strike-slip tectonics has two (or more) horizontal basal plates moving in opposite directions (or only move one of the plates, other are fixed). The visual results are shown from bird's-eye view. Scientists used CT-analysis to collect the cross-section images for the observation of the most influenced area during the simulation.[26]

See also

• Geologic modelling
• Numerical modeling (geology)
• Earth analog

References

1. Schellart, Wouter P.; Strak, Vincent (2016). "A review of analogue modelling of geodynamic processes: Approaches, scaling, materials and quantification, with an application to subduction experiments". Journal of Geodynamics. 100: 7–32. Bibcode:2016JGeo..100....7S. doi:10.1016/j.jog.2016.03.009. ISSN 0264-3707.
2. Ranalli, Giorgio (2001). "Experimental tectonics: from Sir James Hall to the present". Journal of Geodynamics. 32 (1–2): 65–76. Bibcode:2001JGeo...32...65R. doi:10.1016/s0264-3707(01)00023-0. ISSN 0264-3707.
3. Schreurs, Guido; Buiter, Susanne J. H. (Susanne Janita Henriët) (2006). Analogue and numerical modelling of crustal-scale processes. Geological Society. ISBN 978-1862391918. OCLC 191801955.
4. Strak, Vincent; Schellart, Wouter P. (2016). "Introduction to the special issue celebrating 200 years of geodynamic modelling". Journal of Geodynamics. 100: 1–6. Bibcode:2016JGeo..100....1S. doi:10.1016/j.jog.2016.08.003. ISSN 0264-3707.
5. Hall, Sir James. "Geological Studies in the Pays-D'Enhaut Vaudois". Cite journal requires |journal= (help)
6. Cadell, Henry M. (1889). "VII.—Experimental Researches in Mountain Building". Earth and Environmental Science Transactions of the Royal Society of Edinburgh. 35 (1): 337–357. doi:10.1017/S0080456800017658. ISSN 2053-5945.
7. Bailey Willis (1894). The Mechanics of Appalachian Structure. Harvard University. Govt. print. off.
8. HUBBERT, M. K. (1937-10-01). "Theory of scale models as applied to the study of geologic structures". Geological Society of America Bulletin. 48 (10): 1459–1520. Bibcode:1937GSAB...48.1459H. doi:10.1130/gsab-48-1459. ISSN 0016-7606.
9. Ph.H., Kuenen (1937). The negative isostatic anomalies in the East Indies (with Experiments). OCLC 945425263.
10. JACOBY, WOLFGANG R. (1973). "Model Experiment of Plate Movements". Nature Physical Science. 242 (122): 130–134. Bibcode:1973NPhS..242..130J. doi:10.1038/physci242130a0. ISSN 0300-8746.
11. Kincaid, Chris; Olson, Peter (1987-12-10). "An experimental study of subduction and slab migration". Journal of Geophysical Research: Solid Earth. 92 (B13): 13832–13840. Bibcode:1987JGR....9213832K. doi:10.1029/jb092ib13p13832. ISSN 0148-0227.
12. Tapponnier, P.; Peltzer, G.; Le Dain, A. Y.; Armijo, R.; Cobbold, P. (1982). "Propagating extrusion tectonics in Asia: New insights from simple experiments with plasticine". Geology. 10 (12): 611. Bibcode:1982Geo....10..611T. doi:10.1130/0091-7613(1982)10<611:petian>2.0.co;2. ISSN 0091-7613.
13. Vendeville, B.C.; Jackson, M.P.A. (1992-01-01). "The Rise and Fall of Diapirs During Thin-Skinned Extension". Report Investigation. doi:10.23867/ri0209d. ISSN 2475-367X.
14. Brune, James N.; Ellis, Michael A. (1997-05-01). "Structural features in a brittle–ductile wax model of continental extension". Nature. 387 (6628): 67–70. Bibcode:1997Natur.387...67B. doi:10.1038/387067a0. ISSN 0028-0836.
15. Koyi, H. (2007-12-18). "Analogue modelling: From a qualitative to a quantitative technique — A historical outline". Journal of Petroleum Geology. 20 (2): 223–238. Bibcode:1997JPetG..20..223K. doi:10.1111/j.1747-5457.1997.tb00774.x. ISSN 0141-6421.
16. Davy, Ph.; Cobbold, P.R. (1991-03-10). "Experiments on shortening of a 4-layer model of the continental lithosphere". Tectonophysics. 188 (1–2): 1–25. Bibcode:1991Tectp.188....1D. doi:10.1016/0040-1951(91)90311-f. ISSN 0040-1951.
17. JACOBY, WOLFGANG R. (1973). "Model Experiment of Plate Movements". Nature Physical Science. 242 (122): 130–134. Bibcode:1973NPhS..242..130J. doi:10.1038/physci242130a0. ISSN 0300-8746.
18. Mead, Warren J. (1920). "Notes on the Mechanics of Geologic Structures". The Journal of Geology. 28 (6): 505–523. Bibcode:1920JG.....28..505M. doi:10.1086/622731. JSTOR 30063760.
19. Shemenda, Alexander I. (1994). Subduction. Modern Approaches in Geophysics. 11. doi:10.1007/978-94-011-0952-9. ISBN 978-94-010-4411-0. ISSN 0924-6096.
20. Rossetti, Federico; Ranalli, Giorgio; Faccenna, Claudio (1999). "Rheological properties of paraffin as an analogue material for viscous crustal deformation". Journal of Structural Geology. 21 (4): 413–417. Bibcode:1999JSG....21..413R. doi:10.1016/s0191-8141(99)00040-1. ISSN 0191-8141.
21. Ramberg, H. (2010-01-26). "Model Experimentation of the Effect of Gravity on Tectonic Processes". Geophysical Journal of the Royal Astronomical Society. 14 (1–4): 307–329. doi:10.1111/j.1365-246x.1967.tb06247.x. ISSN 0016-8009.
22. Klinkmüller, M.; Schreurs, G.; Rosenau, M.; Kemnitz, H. (2016-08-02). "Properties of granular analogue model materials: A community wide survey". Tectonophysics. 684: 23–38. Bibcode:2016Tectp.684...23K. doi:10.1016/j.tecto.2016.01.017. ISSN 0040-1951.
23. Griffiths, Ross W.; Hackney, Ronald I.; van der Hilst, Rob D. (1995). "A laboratory investigation of effects of trench migration on the descent of subducted slabs". Earth and Planetary Science Letters. 133 (1–2): 1–17. Bibcode:1995E&PSL.133....1G. doi:10.1016/0012-821x(95)00027-a. hdl:1874/7889. ISSN 0012-821X.
24. Gelder, Inge. "Analogue Modelling".
25. Konstantinovskaia, Elena; Malavieille, Jacques (2005-02-26). "Erosion and exhumation in accretionary orogens: Experimental and geological approaches". Geochemistry, Geophysics, Geosystems. 6 (2): Q02006. Bibcode:2005GGG.....6.2006K. doi:10.1029/2004gc000794. ISSN 1525-2027.
26. Dooley, Tim P.; Schreurs, Guido (2012-10-29). "Analogue modelling of intraplate strike-slip tectonics: A review and new experimental results". Tectonophysics. 574-575: 1–71. Bibcode:2012Tectp.574....1D. doi:10.1016/j.tecto.2012.05.030. ISSN 0040-1951.