Power-law viscous materials for analogue experiments: New data on the rheology of highly-filled silicone polymers

The selection of appropriate analogue materials is a central consideration in the design of realistic physical models. We investigate the rheology of highly-filled silicone polymers in order to find materials with a power-law strain-rate softening rheology suitable for modelling rock deformation by dislocation creep and report the rheological properties of the materials as functions of the filler content. The mixtures exhibit strain-rate softening behaviour but with increasing amounts of filler become strain-dependent. For the strain-independent viscous materials, flow laws are presented while for strain-dependent materials the relative importance of strain and strain rate softening/hardening is reported. If the stress or strain rate is above a threshold value some highly-filled silicone polymers may be considered linear visco-elastic (strain independent) and power-law strain-rate softening. The power-law exponent can be raised from 1 to ~3 by using mixtures of high-viscosity silicone and plasticine. However, the need for high shear strain rates to obtain the power-law rheology imposes some restrictions on the usage of such materials for geodynamic modelling. Two simple shear experiments are presented that use Newtonian and power-law strain-rate softening materials. The results demonstrate how materials with power-law rheology result in better strain localization in analogue experiments.

Networking of shear zones at the brittle-to-viscous transition (Cap de Creus, NE Spain)

A crustal-scale shear zone network at the fossil brittle-to-viscous transition exposed at Cap de Creus, NE Spain evolved by coeval fracturing and viscous, mylonitic overprinting of an existing foliation. Initial fracturing led to mylonitic shearing as rock softened in ductilely deformed zones surrounding the fractures. Mylonitic shear zones widened by lateral branching of fractures from these shear zones and by synthetic rotation of the existing foliation between the fractures and shear zones. Shear zones lengthened by a combination of fracturing and mylonitic shearing in front of the shear zone tips. Shear zones interconnected along and across their shearing planes, separating rhomb-shaped lozenges of less deformed rock. Lozenges were subsequently incorporated into the mylonitic shear zones by widening in the manner described above. In this way, deformation became homogeneous on the scale of initial fracturing (metre- to decametre-scale). In contrast, the shear zone network represents localisation of strain on a decametre-length scale. The strength of the continental crust at the time of coeval fracturing and viscous shearing is inferred to have decreased with time and strain, as fracturing evolved to mylonitic shearing, and as the shear zones coalesced to form a through-going network subparallel to the shearing plane. Crustal strength must therefore be considered as strain- and scale-dependent.

Pore formation during dehydration of polycrystalline gypsum observed and quantified in a time-series synchrotron radiation based X-ray micro-tomography experiment

We conducted an in-situ X-ray micro-computed tomography heating experiment at the Advanced Photon Source (USA) to dehydrate an unconfined 2.3 mm diameter cylinder of Volterra Gypsum. We used a purpose-built X-ray transparent furnace to heat the sample to 388 K for a total of 310 min to acquire a three-dimensional time-series tomography dataset comprising nine time steps. The voxel size of 2.2 μm3 proved sufficient to pinpoint reaction initiation and the organization of drainage architecture in space and time. We observed that dehydration commences across a narrow front, which propagates from the margins to the centre of the sample in more than four hours. The advance of this front can be fitted with a square-root function, implying that the initiation of the reaction in the sample can be described as a diffusion process. Novel parallelized computer codes allow quantifying the geometry of the porosity and the drainage architecture from the very large tomographic datasets (20483 voxels) in unprecedented detail. We determined position, volume, shape and orientation of each resolvable pore and tracked these properties over the duration of the experiment. We found that the pore-size distribution follows a power law. Pores tend to be anisotropic but rarely crack-shaped and have a preferred orientation, likely controlled by a pre-existing fabric in the sample. With on-going dehydration, pores coalesce into a single interconnected pore cluster that is connected to the surface of the sample cylinder and provides an effective drainage pathway. Our observations can be summarized in a model in which gypsum is stabilized by thermal expansion stresses and locally increased pore fluid pressures until the dehydration front approaches to within about 100 μm. Then, the internal stresses are released and dehydration happens efficiently, resulting in new pore space. Pressure release, the production of pores and the advance of the front are coupled in a feedback loop.

Role of kilometer-scale weak circular heterogeneities on upper crustal deformation patterns: Evidence from scaled analogue modeling and the Sudbury Basin, Canada

The Sudbury Basin is a non-cylindrical fold basin occupying the central portion of the Sudbury Impact Structure. The impact structure lends itself excellently to explore the structural evolution of continental crust containing a circular region of long-term weakness. In a series of scaled analogue experiments various model crustal configurations were shortened horizontally at a constant rate. In mechanically weakened crust, model basins formed that mimic several first-order structural characteristics of the Sudbury Basin: (1) asymmetric, non-cylindrical folding of the Basin, (2) structures indicating concentric shortening around lateral basin termini and (3) the presence of a zone of strain concentration near the hinge zones of model basins. Geometrically and kinematically this zone corresponds to the South Range Shear Zone of the Sudbury Basin. According to our experiments, this shear zone is a direct mechanical consequence of basin formation, rather than the result of thrusting following folding. Overall, the models highlight the structurally anomalous character of the Sudbury Basin within the Paleoproterozoic Eastern Penokean Orogen. In particular, our models suggest that the Basin formed by pure shear thickening of crust, whereas transpressive deformation prevailed elsewhere in the orogen. The model basin is deformed by thickening and non-cylindrical synformal buckling, while conjugate transpressive shear zones propagated away from its lateral tips. This is consistent with pure shear deformation of a weak circular inclusion in a strong matrix. The models suggest that the Sudbury Basin formed as a consequence of long-term weakening of the upper crust by meteorite impact.

The causes of sinuous crustal-scale deformation patterns in hot orogens: Evidence from scaled analogue experiments and the southern Central Andes

The cause of upper-crustal segmentation into rhomb-shaped, shear zone-bound domains associated with contractional sedimentary basins in hot, wide orogens is not well understood. Here we use scaled multilayered analogue experiments to investigate the role of an orogen-parallel crustal-strength gradient on the formation of such structures. We show that the aspect ratio and size of domains, the sinuous character and abundance of transpressional shear zones vary with the integrated mechanical strength of crust. Upper-crustal deformation patterns and the degree of strain localization in the experiments are controlled by the ratio between the brittle and ductile strength in the model crust as well as gradients in tectonic and buoyancy forces. The experimental results match the first-order kinematic and structural characteristics of the southern Central Andes and provide insight on the dynamics of underlying deformation patterns in hot, wide orogens.

The analogue shear zone : from rheology to associated geometry

The geometry of ductile strain localization phenomena is related to the rheology of the deformed rocks. Both qualitative and quantitative rheological properties of natural rocks have been estimated from finite field structures such as folds and shear zones. We apply physical modelling to investigate the relationship between rheology and the temporal evolution of the width and transversal strain distribution in shear zones, both of which have been used previously as rheological proxies. Geologically relevant materials with well-characterized rheological properties (Newtonian, strain hardening, strain softening, Mohr-Coulomb) are deformed in a shear box and observed with Particle Imaging Velocimetry (PIV). It is shown that the width and strain distribution histories in model shear zones display characteristic finite responses related to material properties as predicted by previous studies. Application of the results to natural shear zones in the field is discussed. An investigation of the impact of 3D boundary conditions in the experiments demonstrates that quantitative methods for estimating rheology from finite natural structures must take these into account carefully.

Understanding treatment characteristics of constructed stormwater wetlands

Constructed wetlands are a common structural treatment measure employed to remove stormwater pollutants and forms an important part of the Water Sensitive Urban Design (WSUD) treatment suite. In a constructed wetland, a range of processes such as settling, filtration, adsorption, and biological uptake play a role in stormwater treatment. Occurrence and effectiveness of these processes are variable and influenced by hydraulic, chemical and biological factors. The influence of hydraulic factors on treatment processes are of particular concern. This paper presents outcomes of a comprehensive study undertaken to define the treatment performance of a constructed wetland highlighting the influence of hydraulic factors. The study included field monitoring of a well established constructed wetland for quantity and quality factors, development of a conceptual hydraulic model to simulate water movement within the wetland and multivariate analysis of quantity and quality data to investigate correlations and to define linkages between treatment performance and influential hydraulic factors. Total Suspended Solids (TSS), Total Nitrogen (TN) and Total Phosphorus (TP) concentrations formed the primary pollutant parameters investigated in the data analysis. The outcomes of the analysis revealed significant reduction in event mean concentrations of all three pollutants species. Treatment performance of the wetland was significantly different for storm events above and below the prescribed design event. For events below design event, TSS and TN load reduction was comparatively high and strongly influenced by high retention time. For events above design event, TP load reduction was comparatively high and was found to be influenced by the characteristics of TP wash-off from catchment surfaces.

WSUD application audit at Gold Coast City

In recent years, Water Sensitive Urban Design (WSUD) has been strongly promoted in South East Queensland to mitigate quantity and quality issues in relation to stormwater. Gold Coast City Council has implemented WSUD devices widely for stormwater management for a number of years and is planning to continue this practice into the future. According to the planning policy of Gold Coast City Council, the adoption of WSUD practices is now mandatory for any new development within the city. As a result, Council is expected to be in possession of tens of millions of dollars of these assets in the future and will be responsible for their maintenance and long-term management. Any shortcoming in the implementation of best practice can potentially result in substantial liability for the Council in the future. However, there has been limited evaluation of WSUD systems in relation to their performance, long-term maintenance, and current knowledge gaps. It was considered that periodical audits of WSUD applications on the Gold Coast is vital to ensure that Council’s WSUD policies are continually improved to new learning and best practice is implemented and risk to Council is mitigated. After a series of stakeholder interviews within Council to understand current practical issues (weaknesses and strengths) in relation to the implementation of WSUD on the Gold Coast, a field audit comprising of condition assessment of eleven WSUD systems within four suburbs was undertaken to identify weaknesses and strengths in WSUD implementation on the Gold Coast. The outcomes of this study are presented in this paper.

Performance characterisation of a constructed wetland

Performance of a constructed wetland is commonly reported as variable due to the site specific nature of influential factors. This paper discusses outcomes from an in-depth study which characterised treatment performance of a wetland based on the variation in runoff regime. The study included a comprehensive field monitoring of a well established constructed wetland in Gold Coast, Australia. Samples collected at the inlet and outlet was tested for Total Suspended Solids (TSS), Total Nitrogen (TN) and Total Phosphorus (TP). Pollutant concentrations in the outflow were found to be consistent irrespective of the variation in inflow water quality. The analysis revealed two different treatment characteristics for events with different rainfall depths. TSS and TN load reduction is strongly influenced by hydraulic retention time where performance is higher for rainfall events below the design event. For small events, treatment performance is higher at the beginning of the event and gradually decreased during the course of the event. For large events, the treatment performance is comparatively poor at the beginning and improved during the course of the event. The analysis also confirmed the variable treatment trends for different pollutant types.