|Wosnitza, E. 2002. Data Analysis in Thermomechanical Analogue Modelling. Schellart, W. P. and Passchier, C. 2002. Analogue modelling of large-scale tectonic processes. Journal of the Virtual Explorer.|
Data Analysis in Thermomechanical Analogue Modelling
The experiment using sand as brittle upper crust implies the necessity for correct scaling of density. A heavy orogen building up causes an instability, sinking down to the model base. In a natural scenario, the density inversion would lead to subsidence of the dense material, down into the mantle. However, such a scenario is extremely unrealistic. Nevertheless, the rise of diapirs could be modelled using this combination of materials.
the relaxation phase of the experiments, the stresses in the upper mantle
changed from a horizontally compressive to an extensional regime. This
is in accordance with the predictions by Dewey (1988) or England &
Thompson (1984), who, however, suggested time scales of up to 120Ma for the extension and erosion of
the orogenes. The relaxation phase also suggests that the paraffin waxes
used as analogue materials do not show completely Newtonian behaviour.
This is in accordance with measurements
A scenario containing a large thrust can be found in the urals, at the boundary between Europe and Asia. The collision of the Eastern European and the Kazachian plate started around 310Ma ago in the Middle Carboniferous (Zonenshain et al., 1994). During the collision, the shortening along the Main Uralian Front (MUF) was between 12.5km and 16.7km (Perez-Estaun et al., 1997), and convergence rates between 1 cm/a and 3 cm/a can be assumed (Zonenshain et al., 1994). The interpretation of the "Urals Seismic Experiment and Integrated Studies" (URSEIS) profile given by Steer et al. (1998) shows various structures, the MUF and the Moho being the most prominent ones (Figure 11). According to Steer et al. (1998), the Moho is a structural detachment, as is the case in the analogue experiments. The dip of the MUF corresponds well to the dip of the thrust in the analogue experiments presented, the scaling of length being correct within 30%. Since the Uralian Orogenesis ended about 210Ma ago (Echtler & Hetzel, 1997) and shows no significant remobilisation (Seward et al., 1997), the Urals had more time to equilibrate than the corresponding relaxation time of the models. Therefore, the topography of the Uralian Moho is less intensive than in the analogue model adjacent to the MUF.
The Alexandrovka Reflection Sequence (ARS) in the URSEIS profile at a depth of 80-100 km coincides with the base of the experiments. It can be interpreted as a continuation of the MUF, as a rheological boundary due to material or phase changes, or as a local fault zone within the mantle (Steer et al., 1998). The analogue experiments presented do not show evidence for penetration of the thrust zone into the mantle. If the ARS is a continuation of the MUF, it has not been formed synorogenically. Either the Alexandrovka Reflection Sequence is an older structure just as the MUF, or one of the alternative interpretations has to be used.
Although the length scales of the Higher Himalaya are larger than the lengths modelled by the analogue experiments presented, the geometry of the models is similar to that of the Himalayas according to the interpretation of the "International Deep Profiling of Tibet and the Himalaya" (INDEPTH) seismic profile (Figure 12). The convergence rate between the Indian and the Asian plates has been estimated by Henry et al. (1997) from various sources to be in the range 1.5-2.0 cm/a, which is in the range modelled by the analogue models presented. Since in the wedge experiments the deformation data obtained from the photos did not represent the deformation in the whole model, the data should be interpreted with care. Nevertheless, the wedge of weak material has been extruded between the two colliding blocks (Figure 13).
Newer numerical models of Beaumont et al. (2001) and analysis of topographical data by Clark & Royden (2000) show evidence for a channel of weak material instead of a wedge. To verify this concept by using analogue models, a way has to be found to model larger lengths either by choosing different materials, or by modifying the parameters of the apparatus.
The main points from the data analysis presented in the previous chapter apply to both types of orogenes. The steepness of the natural roots suggests convergence rates rather on the higher side of the ranges quoted above: 3 cm/a for the Urals, 2 cm/a for the Himalayas.
Dextral shear has been observed in the experiments at both ends of the thrust or the weak zone. The tip of the footwall block is not accessible in nature, but in the area of the upper crust where the thrust cuts the surface indicators for dextral simple shear should occur.
of the bulk stress data as well as the distribution of differential stresses
stresses in the orogenes, at least for some Ma after the deformation.
In the Urals, these remanent stresses might have decayed by now, but in
The method of thermomechanical analogue modelling has proven the ability to produce, analyse and present data for continental collision on lithospheric scales. It has been shown that the set-up used is able to reproduce and visualize the deformation of isotherms. The machine is ready to use and the experimental techniques as well as the procedures for data analysis were tested and proved to give reproducable results.
From the work with the apparatus for thermomechanical analogue modelling and from the first experiments, several suggestions arise for future machines and applications.
One important point is to learn more about the materials involved. For the analogue materials, measurement of the stress exponent n is necessary, as well more accurate estimates of the other flow properties. On the other hand, the data for natural rocks also lacks certain points. The data quoted in the literature is gained mostly from samples of some cubic centimetres, while it is applied to layers some ten kilometres thick and at strain rates some orders of magnitude less than those used to obtain the data.
Maybe in future experiments it might also be possible to introduce the important surface process erosion and sedimentation, without disturbance of the thermal equilibrium. The use of microwaves could be a way to include local heat production into the model.
Presently (December 2001), the apparatus is used for a Masters Thesis (Diplomarbeit of Philip Schiwek) to examine the onset of subduction at the (still) passive continental margin of Galicia. For these experiments, the scaling of density is extremely important, and the use of thin layers of lead powder embedded into the wax layers seems to be promising.
The machine presented is able to produce extensional regimes as well as compressive scenarios. The rise of diapirs or rifting processes can be modelled, maybe using sand as a brittle overburden. Additionally, the collisional experiments could be continued, examining the extensional phase of the collapse
In future thermomechanical experiments, the bottom temperature sensor should be fixed to the model base directly. All obstacles in the experimental domain must be avoided. The small floating PT100 temperature probes are permissible. Other sensors could be inserted into the models, such as piezoelectric stress sensors. Thus, the stress analysis presented in this work could be verified. Anyway, a second strain gauge for the second piston is necessary. Due to different material fluxes behind the pistons, the bulk stresses imposed to the plates are not equal. By also measuring the material level on the rear side of the pistons, the pure loading component could be identified.
A less steep temperature gradient in the model should be able to provide larger viscosity contrasts at the model Moho. Nevertheless, it is necessary to have the base of the model molten to decrease friction at the bottom of the model and to allow material flux behind the pistons for isostatic equilibration.
An important factor in the experiments is the analysis of the data obtained. As errors induced by the process of digitization of the marker positions can influence the results, an automated procedure for the strain analysis would largely improve both the quality and the quantity of data. It would also be easier to take into account more marker particles at more time steps.
For further experiments modelling the Higher Himalayas, the scaling of the experiments must be recalculated. Using a different set of gears, the motor could provide higher deformation rates, resulting in a smaller scaling factor for length, so that larger natural scenarios could be modelled. However, the experimental techniques for these models must also be refined to reduce boundary effects within the weak wedge or the newly suggested channel.
Probably one of the most important applications for thermomechanical analogue modelling lies in its interplay with numerical models. Although simulating large internal deformations in combination with discrete detachment zones still seems to pose numerical problems, the increasing power of computers keeps up with the refinement of data analysis in analogue experiments. The disadvantage of laboratory models to be limited to physical material properties does not exist in numerical models. On the other hand, analogue models include the physics implicitly, whereas numerical models can only simulate the processes they are programmed to simulate. Therefore, analogue models could be used to calibrate numerical models, which can then extend the results to nature. Only if a numerical model is able to reproduce a simple analogue model close enough, the computer code is justified to model geological scenarios.
Analogue modelling needs to be advanced further and further to provide an additional instrument for the Geologist. In this work, a major step forward has been done.
I want to express my gratitude towards Djordje Grujic, David Tanner and Jan Behrmann, who contibuted to this paper with many discussions. Beate Söhngen did the invaluable work of digitizing the data. Grant BE1041 by the Deutsche Forschungsgemeinschaft financed large parts of the project.
J., Gardien, V., Allemand, P. & Grandjean, P. (2000): Analog modelling
of the deformation of anisotropic partially molten rocks using paraffin
wax. Geophysical Research Abstracts 2: SE24/p27.
C., Jamieson, R. A., Nguyen, M. H. & Lee, B. (2001): Himalayan tectonics
explained by extrusion of a low-viscosity crustal channel coupled to focused
surface denudation. Nature 414(6565): 738 - 742.
P. D., Jessell, M. W. & Passchier, C. W. (1993): The analysis of progressive
deformation in rock analogues. Journal of Structural Geology 15(3-5):
403 - 411.
J. N. & Ellis, M. A. (1997): Structural features in a brittle-ductile
wax model of continental extension. Nature 387(6628): 67 - 70.
H. M. (1889): Experimental researches in mountain building. Transactions
of the royal society Edinburgh 35: 337 - 357.
N. L. & Tsenn, M. C. (1987): Flow properties of continental lithosphere.
Tectonophysics 136: 27 - 63.
A. I., Burg, J. P. & Mattauer, M. (2000): Evolutionary model of the
Himalaya-Tibet system; geopoem based on new modelling, geological and
geophysical data. Earth and Planetary Science Letters 174(3 - 4): 397-409.
A. I., Mattauer, M., Malavieille, J. & Bokun, A. N. (1995): A mechanism
for syn-collisional rock exhumation and associated normal faulting: results
from physical modelling. Earth and Planetary Science Letters 132: 225-232.
M. K. & Royden, L. H. (2000): Topographic ooze; building the eastern
margin of Tibet by lower crustal flow. Geology 28(8): 703 - 706.
P. R. & Jackson, M. P. A. (1992): Gum rosin (colophony): a suitable
material for thermomechanical modelling of the lithosphere. Tectonophysics
210: 255 - 271.
P. R., Rossello, E. & Vendeville, B. (1989): Some experiments on interacting
sedimentation and deformation above salt horizons. Bulletin de la Societe
Geologique de France 8: 453 - 460.
P. & Cobbold, P. R. (1991): Experiments on shortening of a 4-layer
model of the continental lithosphere. Tectonophysics 188: 1 - 25.
E. R., Heasler, H. P., Buelow, K. L., Baker, K. H. & Hallin, J. S.
(1988): Signicance of past and recent heat-flow and radioactivity studies
in the Southern Rocky Mountains region. Geological Society of America
Bulletin 100: 1851 - 1885.
J. F. (1988): Extensional collapse of orogenes. Tectonics 7(6): 1123 -
Z. (1984): Some formulae for calculating parameters of the strain ellipse.
Tectonophysics 110: 167 - 175.
H. P. & Hetzel, R. (1997): Main Uralian thrust and main Uralian normal
fault; nonextensional Palaeozoic high-P rock exhumation, oblique collision,
and normal faulting in the Southern Urals. Tera Nova 9(4): 158 - 162.
T. (1994): Deviatoric stressitis; a virus infecting the earth science
community. EOS, Transactions 75(18): 209 - 210.
P. C. (1987): Diffuse continental deformation; length scales, rates and
metamorphic evolution. Philosophical Transactions of the Royal Society
of London, Series A: Mathematical and Physical Sciences 321(1557): 3 -
P. C. & Holland, T. J. B. (1979): Archimedes and the Tauern eclogites;
the role of buoyancy in the preservation of exotic eclogite blocks. Earth
and Planetary Science Letters 44(2): 287 - 294.
P. C. & Molnar, P. (1997): Active deformation of Asia: from kinematics
to dynamics. Science 278(5338): 647 - 650.
P. C. & Thompson, A. B. (1984): Pressure-temperature-time paths of
regional metamorphism; I, Heat transfer during the evolution of regions
of thickened continental crust. Journal of Petrology 25(4): 894 - 928.
A., Vasseur, G. & Daignieres, M. (1981): An inverse problem in microtectonics
for the determination of stress tensors from fault striation analysis.
Journal of Structural Geology 3(1): 51 - 65.
C. M. R. (1990): The Solid Earth. Cambridge University Press, Cambridge,
D., Casey, M., Davidson, C., Hollister, L. S., Kuendig, R., Pavlis, T.
L. & Schmid, S. M. (1996): Ductile extrusion of the Higher Himalayan
Crystalline in Bhutan; evidence from quartz microfabrics. Tectonophysics
260: 21 - 43.
M., Braun, J., Brown, M., Kukowski, N., Paterson, M. S., Schmid, S. M.,
Stockhert, B., Stuwe, K., Thompson, A. B. & Wosnitza, E. M. (2001):
Rheology and geodynamic modelling: the next step forward. International
Journal of Earth Sciences 90: 149 - 156.
P., Le Pichon, X. & Goffe, B. (1997): Kinematic, thermal and petrological
model of the Himalayas: constraints related to metamorphism within the
underthrust Indian crust and topographic elevation. Tectonophysics 273(1-2):
31 - 56.
K., Koide, H., Inami, K., Iwamura, S. & Mitsui, S. (1972): Mechanical
properties of Japanese Tertiary sedimentary rocks under high confining
pressures. Geological Survey of Japan Report 244.
M. K. (1937): Theory of scale models as applied to the study of geologic
structures. Geological Society of America Bulletin 48: 1459 - 1520.
M. K. & Rubey, W. W. (1959): Mechanics of fluid-filled porous solids
and its application to overthrust faulting. Geological Society of America
Bulletin 70(2): 115 - 166.
R. A. & Beaumont, C. (1988): Orogeny and metamorphism; a model for
deformation and pressure-temperature-time paths with applications to the
Central and Southern Appalachians. Tectonics 7: 417 - 445.
R. A., Beaumont, C., Hamilton, J. & Fullsack, P. (1996): Tectonic
assembly of inverted metamorphic sequences. Geology 24(9): 839 - 854.
S. H. & Kronenberg, A. K. (1987): Rheology of the lithosphere; selected
topics. Reviews of Geophysics 25(6): 1219 - 1244.
H.-J. (1979): Lineare Algebra. de Gruyter, Berlin, New York, 9 edition,
H. A. (2001): Modeling the influence of sinking anhydrite blocks on salt
diapirs targeted for hazardous waste disposal. Geology 29(5): 387 -390.
H. A. & Skelton, A. (2001): Centrifuge modelling of the evolution
of low-angle detachment faults from high-angle normal faults. Journal
of Structural Geology 23(8): 1179 - 1185.
S. E., Schnuerle, P. & Malavieille, J. (1994): Coulomb theory applied
to accretionary and nonaccretionary wedges: Possible causes for tectonic
erosion and/or frontal accretion. Journal of Geophysical Research 99(B6):
12 033 - 12 055.
N. S. (1988): The rheology of paraffin wax and its usefulness as an analogue
for rocks. Bulletin of the Geological Institutions of the University of
Uppsala N.S.14: 181 - 193.
N. S. (1991): The analysis of progressive deformation from an inscribed
grid. Journal of Structural Geology 13(7): 859 - 864.
N. S. (1995): Nonlithostatic pressure during sediment subduction and the
development and exhumation of high-pressure metamorphic rocks. Journal
of Geophysical Research 100(B1): 571 - 583.
J. W., Hussain, A. & Yeats, R. S. (1993): The Main Boundary Thrust
and propagation of deformation into the foreland fold-and-thrust beld
in northern Pakistan near the Indus River. In: Himalayan Tectonics, Treloar,
P. J. & Searle, M. P., eds., volume 74 of Special Publications, pp.
581 - 588. Geological Society, London.
W. D. (1990): Kinematics, stress, deformation and material behavior. Journal
of Structural Geology 12(8): 953 - 971.
W. D., Hobbs, B. E., Lister, G. S. & Williams, P. F. (1980): Vorticity
and non- coaxiality in progressive deformations. Journal of Structural
Geology 2(3): 371 - 378.
R. P. (1993): Pressure, temperature and time constraints on metamorphism
across the Main Central Thrust zone and High Himalayan Slab in the Garhwal
Himalaya. In: Himalayan Tectonics, Treloar, P. J. & Searle, M. P.,
eds., volume 74 of Special Publications, pp. 485 - 509.Geological
H. C. & Richter, F. M. (1982): Convection experiments in fluids with
highly temperature-dependent viscosity and the thermal evolution of the
planets. Physics of the Earth and Planetary Interiors 29: 320 - 329.
K. D., Zhao, W., Brown, L. D., Kuo, J. T., Che, J., Liu, X., Klemperer,
S. L., Makovsky, Y., Meissner, R., Mechie, J., Kind, R., Wenzel, F., Ni,
J. F., Nabelek, J., Chen, L., Tan, H., Wie, W., Jones, A. G., Booker,
J. R., Unsworth, M. J., Kidd, W. S. F., Hauck, M. L., Alsdorf, D., Ross,
A., Cogan, M. J., Wu, C., Sandvol, E. A. & Edwards, M. A. (1996):
Partially molten middle crust beneath southern Tibet; synthesis of Project
INDEPTH results. Science 274(5293): 1684 - 1688.
D. W. & Brune, J. N. (1972): Ridge transform fault spreading pattern
in freezing wax. Science 178(4058): 301 - 304.
D. W. & Brune, J. N. (1975): An explanation for the orthogonality
of ocean ridges and transform faults. Journal of Geophysical Research
80(17): 2575 - 2585.
A. & Hobbs, B. E. (1989): The strength of the continental crust, detachment
zones and the development of plastic instabilities. Tectonophysics 158:
269 - 289.
C. W. & Trouw, R. A. J. (1996): Microtectonics. Springer Verlag, Berlin,
Heidelberg, New York, London, Paris, Tokyo, Hong Kong, 289 pp.
A., Alvarez-Marron, J., Brown, D., Puchkov, V., Gorozhanina, Y. &
Baryshev, V. (1997): Along-strike structural variations in the foreland
thrust and fold belt of the southern Urals. Tectonophysics 276(1-4): 265
ner, O. A. & Ramsay, J. G. (1982): Constraints on geological strain
rates; arguments from finite strain states of naturally deformed rocks.
Journal of Geophysical Research 87(B1): 311 - 321.
H. (1967): Gravity, Deformation and the Earth's Crust. Academic Press,
London, New York, 214 pp.
J. G. & Huber, M. I. (1983): The techniques of modern structural geology.
Volume 1: Strain analysis. Academic Press, London, New York, 307 pp.
G. (1987): Rheology of the Earth. Allen & Unwin, Boston, London, Sydney,
Wellington, 366 pp.
G. (2001): Experimental tectonics: from Sir James Hall to the present.
Journal of Geodynamics 32: 65 - 76.
F., Faccenna, C., Ranalli, G. & Storti, F. (2000): Convergence rate-dependent
growth of experimental viscous orogenic wedges. Earth and Planetary Science
Letters 178(3 - 4): 367 - 372.
F., Ranalli, G. & Faccenna, C. (1999): Rheological properties of paraffin
as an analogue material for viscous crustal deformation. Journal of Structural
Geology 21(4): 413 - 417.
M. & Koyi, H. A. (2001): Modeling the role of erosion in diapir development
in contractional settings. In: Tectonic Modeling: A Volume in Honor of
H. A. & Mancktelow, N. S., eds., Geological Society of America Memoir
193, pp. 111 - 122. Geological Society of America, Boulder, Colorado.
R. (1988): Algorithms. Addison Wesley, Reading, Mass., 657 pp.
D., Perez-Estaun, A. & Puchkov, V. (1997): Preliminaryfi ssion-track
results from the southern Urals | Sterlitamak to Magnitogorsk. Tectonophysics
276(1-4): 281 - 290.
A. I. & Grocholsky, A. L. (1994): Physical modeling of slow sea floor
spreading. Journal of Geophysical Research 99(B5): 9137 - 9153.
D. N., Knapp, J. H., Brown, L. D., Echtler, H. P., Brown, L. D. &
Berzin, R. (1998): Deep structure of the continental lithosphere in an
unextended orogen; an explosive-source seismic reflection profile in the
Urals (Urals Seismic Experiment and Integrated Studies (URSEIS 1995)).
Tectonics 17(2): 143 - 157.
A. B. & England, P. C. (1984): Pressure-temperature-time paths of
regional metamorphism; II, Their inference and interpretation using mineral
assemblages in metamorphic rocks. Journal of Petrology 25(4): 929 - 955.
R. J. & Moores, E. M. (1992): Structural Geology. Freeman, New York,
B. & Cobbold, P. R. (1988): How normal faulting and sedimentation
interact to produce listric fault profiles and stratigraphic wedges. Journal
of Structural Geology 10(7): 649 - 659.
R. & Schmeling, H. (1986): Scaling of Newtonian and non-Newtonian
fluid dynamics without inertia for quantitative modelling of rock flow
due to gravity (including the concept of rheological similarity). Physics
of the Earth and Planetary Interiors 43: 316 - 330.
S., Beaumont, C. & Fullsack, P. (1993): Mechanical model for the tectonics
of doubly vergent compressional orogens. Geology 21(4): 371 - 374.
S. (1991): The Mathematica Book. Cambridge University Press, Cambridge,
E. M., Grujic, D., Hofmann, R. & Behrmann, J. H. (2001): New apparatus
for thermomechanical analogue modeling. In: TectonicModeling:
A Volume in Honor of Hans Ramberg, Koyi, H. A. & Mancktelow, N. S.,
eds., Geological Society of America
Zonenshain, L. P., Korienevsky, V. G., Kzamin, V. G., Pechersky, D. M., Khain, V. V. & Matveenkov, V. V. (1994): Plate tectonic model of the South Urals development. Tectonophysics 109(1/2): 95 - 135.