Multiscale coupling and multiphysics approaches in earth sciences : applications

Geoscientists are confronted with the challenge of assessing nonlinear phenomena that result from multiphysics coupling across multiple scales from the quantum level to the scale of the earth and from femtoseconds to the 4.5 Ga of history of our planet. We neglect in this review electromagnetic modelling of the processes in the Earth’s core, and focus on four types of couplings that underpin fundamental instabilities in the Earth. These are thermal (T), hydraulic (H), mechanical (M) and chemical (C) processes which are driven and controlled by the transfer of heat to the Earth’s surface. Instabilities appear as faults, folds, compaction bands, shear/fault zones, plate boundaries and convective patterns. Convective patterns emerge from buoyancy overcoming viscous drag at a critical Rayleigh number. All other processes emerge from non-conservative thermodynamic forces with a critical critical dissipative source term, which can be characterised by the modified Gruntfest number Gr. These dissipative processes reach a quasi-steady state when, at maximum dissipation, THMC diffusion (Fourier, Darcy, Biot, Fick) balance the source term. The emerging steady state dissipative patterns are defined by the respective diffusion length scales. These length scales provide a fundamental thermodynamic yardstick for measuring instabilities in the Earth. The implementation of a fully coupled THMC multiscale theoretical framework into an applied workflow is still in its early stages. This is largely owing to the four fundamentally different lengths of the THMC diffusion yardsticks spanning micro-metre to tens of kilometres compounded by the additional necessity to consider microstructure information in the formulation of enriched continua for THMC feedback simulations (i.e., micro-structure enriched continuum formulation). Another challenge is to consider the important factor time which implies that the geomaterial often is very far away from initial yield and flowing on a time scale that cannot be accessed in the laboratory. This leads to the requirement of adopting a thermodynamic framework in conjunction with flow theories of plasticity. This framework allows, unlike consistency plasticity, the description of both solid mechanical and fluid dynamic instabilities. In the applications we show the similarity of THMC feedback patterns across scales such as brittle and ductile folds and faults. A particular interesting case is discussed in detail, where out of the fluid dynamic solution, ductile compaction bands appear which are akin and can be confused with their brittle siblings. The main difference is that they require the factor time and also a much lower driving forces to emerge. These low stress solutions cannot be obtained on short laboratory time scales and they are therefore much more likely to appear in nature than in the laboratory. We finish with a multiscale description of a seminal structure in the Swiss Alps, the Glarus thrust, which puzzled geologists for more than 100 years. Along the Glarus thrust, a km-scale package of rocks (nappe) has been pushed 40 km over its footwall as a solid rock body. The thrust itself is a m-wide ductile shear zone, while in turn the centre of the thrust shows a mm-cm wide central slip zone experiencing periodic extreme deformation akin to a stick-slip event. The m-wide creeping zone is consistent with the THM feedback length scale of solid mechanics, while the ultralocalised central slip zones is most likely a fluid dynamic instability.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Wet tropics reef rescue impact report 2009

The Wet Tropics region has a unique water asset and is also considered a priority region for the improvement of water quality entering the Great Barrier Reef due to a combination of high rainfall, intensive agricultural use, urban areas and the proximity of valuable reef assets to the coast. Agricultural activities are one of many identified threats to water quality and water flows in the Wet Tropics in terms of sediment and pollutant-related water quality decline. Information describing the current state of agricultural management practices across the region is patchy at best. Based on the best available information on agricultural management practices in the Wet Tropics in 2008, it is clear that opportunities exist to improve nutrient, sediment and pesticide management practice to reduce the impact on the water asset and the Great Barrier Reef. Based on current understandings of practices and the relationship between practices and reef water quality, the greatest opportunities for improved water quality are as follows: · nutrients – correct rate and the placement of fertilisers; · pesticides – improve weed control planning, herbicide rates and calibration practice; and · soil and sediment – implement new farming system practices. The 2008-09 Reef Rescue program sought to accelerate the rate of adoption of improved management practices and through Terrain invested $6.8M in the 2008-09 year for: · landholder water quality improvement incentive payments; · cross regional catchment repair of wetlands and riparian lands in areas of high sediment or nutrient loss; and · partnerships in the region to lever resources and support for on-ground practice change. The program delivered $3,021,999 in onground incentives to landholders in the Wet Tropics to improve farm practices from D or C level to B or A level. The landholder Water Quality Incentives Grants program received 300 individual applications for funding and funded 143 individual landholders to implement practice change across 36,098 ha of farm land. It is estimated that the Reef Rescue program facilitated practice change across 21% of the cane industry, and 20% of the banana industry. The program levered an additional $2,441,166 in landholder cash contributions and a further $907,653 in non-cash in-kind contributions bringing the total project value of the landholder grants program in the Wet Tropics to $6,370,819. Most funded projects targeted multiple water quality objectives with a focus on nutrient and sediment reduction. Of the 143 projects funded, 115 projects addressed nutrient management either as the primary focus or in combination with strategies that targeted other water quality objectives. Overall, 82 projects addressed two or more water quality targets. Forty-five percent of incentive funds were allocated to new farming system practices (direct drill legumes, zonal tillage equipment, permanent beds, min till planting equipment, GPS units, laser levelling), followed by 24% allocated to subsurface fertiliser applicators (subsurface application of fertiliser using a stool splitter or beside the stool, at the correct Six Easy Steps rate). As a result, Terrain estimates that the incentive grants achieved considerable reductions in nitrogen, phosphorus, sediment and pesticide loads. The program supported nutrient management training of 167 growers managing farms covering over 20% of the area harvested in 2008, and 18 industry advisors and resellers. This resulted in 115 growers (155 farms) developing nutrient management plans. The program also supported Integrated Weed Management training of 80 growers managing farms covering 8% of the area harvested in 2008, and 6 industry advisors and resellers. This report, which draws on the best available Reef Rescue Management Monitoring, Evaluation, Reporting, and Improvement (MERI) information to evaluate program performance and impact on water quality outcomes, is the first in a series of annual reports that will assess and evaluate the impact of the Reef Rescue program on agricultural practices and water quality outcomes. The assessment is predominantly focused on the cane industry because of data availability. In the next stage, efforts will expand to: · improve practice data for the banana and grazing industry; · gain a better understanding of the water quality trends and the factors influencing them in the Wet Tropics; in particular work will focus on linking the results of the Paddock to Reef monitoring program and practice change data to assess program impact; · enhance estimations of the impact of practice change on pollutant loads from agricultural land use; · gain a better understanding of the extent of ancillary practice (change not directly funded) resulting from Reef Rescue training/ education/communication programs; and · provide a better understanding of the economic cost of practice change across the Wet Tropics region. From an ecological perspective, water quality trends and the factors that may be contributing to change, require further investigation. There is a critical need to work towards an enhanced understanding of the link between catchment land management practice change and reef water quality, so that reduced nutrient, sediment, and pesticide discharge to the Great Barrier Reef can be quantified. This will also assist with future prioritisation of grants money to agricultural industries, catchments and sub catchments. From a social perspective, the program has delivered significant water quality benefits from landholder education and training. It is believed that these activities are giving landholders the information and tools to implement further lasting change in their production systems and in doing so, creating a change in attitude that is supportive and inclusive of Natural Resource Management (NRM). The program in the Wet Tropics has also considerably strengthened institutional partnerships for NRM, particularly between NRM and industry and extension organisations. As a result of the Reef Rescue program, all institutions are actively working together to collectively improve water quality. The Reef Rescue program is improving water quality entering the Great Barrier Reef Lagoon by catalysing substantial activity in the Wet Tropics region to improve land management practices and reduce the water quality impact of agricultural landscapes. The solid institutional partnerships between the regional body, industry, catchment and government organisations have been fundamental to the successful delivery of the landholder grant and catchment rehabilitation programs. Landholders have generally had a positive perception and reaction to the program, its intent, and the practical, focused nature of grant-based support. Demand in the program was extremely high in 2008-09 and is expected to increase in 2009-2010.

CFD simulations of flow and heat transfer through the porous interface of a metal foam heat exchanger

This paper offers numerical modelling of a waste heat recovery system. A thin layer of metal foam is attached to a cold plate to absorb heat from hot gases leaving the system. The heat transferred from the exhaust gas is then transferred to a cold liquid flowing in a secondary loop. Two different foam PPI (Pores Per Inch) values are examined over a range of fluid velocities. Numerical results are then compared to both experimental data and theoretical results available in the literature. Challenges in getting the simulation results to match those of the experiments are addressed and discussed in detail. In particular, interface boundary conditions specified between a porous layer and a fluid layer are investigated. While physically one expects much lower fluid velocity in the pores compared to that of free flow, capturing this sharp gradient at the interface can add to the difficulties of numerical simulation. The existing models in the literature are modified by considering the pressure gradient inside and outside the foam. Comparisons against the numerical modelling are presented. Finally, based on experimentally-validated numerical results, thermo-hydraulic performance of foam heat exchangers as waste heat recovery units is discussed with the main goal of reducing the excess pressure drop and maximising the amount of heat that can be recovered from the hot gas stream.

Megaprojects managements in Ecuador : challenges and opportunities

Megaprojects are described as large, complex and expensive construction projects. Recent studies have shown that megaprojects often result in cost overruns, time extensions and undesired outcomes. Regardless, megaprojects are common, particularly in developing countries, as they are a trigger for social and economic development (Li et al., 2010). Since 2007, the Government of Ecuador has begun an unprecedent investment in infrastructure. Through the National Water Secretary, the government has 16 projects in agenda accounting for over $ 3 billion, with 6 projects currently under construction. These projects are considered flagship infrastructure in the endeavour to enhance the country´s productivity.The Bulubulu-Naranjal-Cañar project, a $406 million multi-purpose hydraulic project for irrigation and flood control, consists of over 1,000 activities and was proposed to be completed by 2015. This novel project for Ecuador, presented as a case study, represents a challenge for project management and financing. The purpose of this preliminary study is to provide an insight to megaproject management in Ecuador, and propose improvements to megaproject management through optimization of stochastic project schedules.

The capture and retention evaluation of a stormwater gross pollutant trap design

Gross pollutant traps (GPT) are designed to capture and retain visible street waste, such as anthropogenic litter and organic matter. Blocked screens, low/high downstream tidal waters and flows operating above/below the intended design limits can hamper the operations of a stormwater GPT. Under these adverse operational conditions, a recently developed GPT was evaluated. Capture and retention experiments were conducted on a 50% scale model with partially and fully blocked screens, placed inside a hydraulic flume. Flows were established through the model via an upstream channel-inlet configuration. Floatable, partially buoyant, neutrally buoyant and sinkable spheres were released into the GPT and monitored at the outlet. These experiments were repeated with a pipe-inlet configured GPT. The key findings from the experiments were of practical significance to the design, operation and maintenance of GPTs. These involved an optimum range of screen blockages and a potentially improved inlet design for efficient gross pollutant capture/retention operations. For example, the outlet data showed that the capture and retention efficiency deteriorated rapidly when the screens were fully blocked. The low pressure drop across the retaining screens and the reduced inlet flow velocities were either insufficient to mobilise the gross pollutants, or the GPT became congested.

The Cressbrook Creek alluvial aquifer system, Southeast Queensland : hydrochemistry and isotopes to determine hydrological processes and response to floods

This study developed an understanding of hydrological processes within the Cressbrook Creek catchment of the upper Brisbane River, in particular for the alluvial aquifers. Those aquifers within the lower catchment are used for intensive irrigation, and have been impacted by long-term drought followed by flooding. The study utilised water chemistry, isotopic characters and hydraulic measurements to determine factors such as recharge, links between creeks and groundwater, and variations in water quality. The catchment-wide study will enable improved management of the local water resources.

Performance characterisation of a stormwater treatment bioretention basin

Treatment performance of bioretention basins closely depends on hydrologic and hydraulic factors such as rainfall characteristics and inflow and outflow discharges. An in-depth understanding of the influence of these factors on water quality treatment performance can provide important guidance for effective bioretention basin design. In this paper, hydraulic and hydrologic factors impacting pollutant removal by a bioretention basin were assessed under field conditions. Outcomes of the study confirmed that the antecedent dry period plays an important role in influencing treatment performance. A relatively long antecedent dry period reduces nitrite and ammonium concentrations while increasing the nitrate concentration, which confirms that nitrification occurs within the bioretention basin. Additionally, pollutant leaching influences bioretention basin treatment performance, reducing the nutrients removal efficiency, which was lower for high rainfall events. These outcomes will contribute to a greater understanding of the treatment performance of bioretention basins, assisting in the design, operation and maintenance of these systems.

Development and application of Palygorskite porous ceramsite in a biological aerated filter (BAF)

Novel filter Palygorskite porous ceramsite (PC) was prepared using Palygorskite clay, poreforming material sawdust, and sodium silicate with a mass ratio of 10:2:1 after sintering at 700°C for 180 min. PC was characterized with X-ray diffraction, X-ray fluorescence, scanning electron microscopy, elemental, and porosimetry. PC had a total porosity of 67% and specific surface area of 61 m2/g. In order to assess the usefulness of PC as a medium for biological aerated filters (BAF), PC and (commercially available ceramsite) CAC were used to treat wastewater city in two laboratory-scale upflow BAFs. The results showed that the reactor containing PC was more efficient than the reactor containing CAC in terms of total organic carbon (TOC), ammonia nitrogen (NH3-N), and the removal of total nitrogen (TN) and phosphorus (P). This system was found to be more efficient at water temperatures ranging from 20 to 26°C, an air–water (A/W) ratio of 3:1, dissolved oxygen concentration >4.00 mg/L, and hydraulic retention time (HRT) ranging from 0.5 to 7 h. The interconnected porous structure produced for PC was suitable for microbial growth, and primarily protozoan and metazoan organisms were found in the biofilm. Microorganism growth also showed that, under the same submerged culture conditions, the biological mass in PC was significantly higher than in CAC (34.1 and 2.2 mg TN/g, respectively). In this way, PC media can be considered suitable for the use as a medium in novel biological aerated filters for the simultaneous removal of nitrogen and phosphorus.

Axial capacity of barrette piles embedded in gravel layer

In this paper, the axial performance of two heavily instrumented barrette piles, with and without grouting, socket into gravel layer in Taipei are evaluated based on the results of pile load tests. Both piles are 44 m long with the same dimension of 0.8 by 2.7 m, installed by hydraulic long bucket. One of the piles with toe grouting was socket 6 m into gravel layer and the other pile without toe grouting was socket 3 m into gravel layer. The load versus displacement relationships at pile head, the t-z curves of upper soil layers and of bottom gravel layer, and the tip resistance versus displacement relationships are important concerns and are presented in the paper. The t-z curves interpreted from the measured data along depth are also simulated by the hyperbolic model.

Three-dimensional geological modelling of the Galilee and central Eromanga basins, Australia: New insights into aquifer/aquitard geometry and potential influence of faults on inter-connectivity

Study region The Galilee and Eromanga basins are located in central Queensland, Australia. Both basins are components of the Great Artesian Basin which host some of the most significant groundwater resources in Australia. Study focus This study evaluates the influence of regional faults on groundwater flow in an aquifer/aquitard interbedded succession that form one of the largest Artesian Basins in the world. In order to assess the significance of regional faults as potential barriers or conduits to groundwater flow, vertical displacements of the major aquifers and aquitards were studied at each major fault and the general hydraulic relationship of units that are juxtaposed by the faults were considered. A three-dimensional (3D) geological model of the Galilee and Eromanga basins was developed based on integration of well log data, seismic surfaces, surface geology and elevation data. Geological structures were mapped in detail and major faults were characterised. New hydrological insights for the region Major faults that have been described in previous studies have been confirmed within the 3D geological model domain and a preliminary assessment of their hydraulic significance has been conducted. Previously unknown faults such as the Thomson River Fault (herein named) have also been identified in this study.

Predicting flow and transport in highly heterogeneous alluvial aquifers

Successful prediction of groundwater flow and solute transport through highly heterogeneous aquifers has remained elusive due to the limitations of methods to characterize hydraulic conductivity (K) and generate realistic stochastic fields from such data. As a result, many studies have suggested that the classical advective-dispersive equation (ADE) cannot reproduce such transport behavior. Here we demonstrate that when high-resolution K data are used with a fractal stochastic method that produces K fields with adequate connectivity, the classical ADE can accurately predict solute transport at the macrodispersion experiment site in Mississippi. This development provides great promise to accurately predict contaminant plume migration, design more effective remediation schemes, and reduce environmental risks. Key Points Non-Gaussian transport behavior at the MADE site is unraveledADE can reproduce tracer transport in heterogeneous aquifers with no calibrationNew fractal method generates heterogeneous K fields with adequate connectivity

Perturbed Lignification Impacts Tree Growth in Hybrid Poplar-A Function of Sink Strength, Vascular Integrity, and Photosynthetic Assimilation

The effects of reductions in cell wall lignin content, manifested by RNA interference suppression of coumaroyl 3'-hydroxylase, on plant growth, water transport, gas exchange, and photosynthesis were evaluated in hybrid poplar trees (Populus alba 3 grandidentata). The growth characteristics of the reduced lignin trees were significantly impaired, resulting in smaller stems and reduced root biomass when compared to wild-type trees, as well as altered leaf morphology and architecture. The severe inhibition of cell wall lignification produced trees with a collapsed xylem phenotype, resulting in compromised vascular integrity, and displayed reduced hydraulic conductivity and a greater susceptibility to wall failure and cavitation. In the reduced lignin trees, photosynthetic carbon assimilation and stomatal conductance were also greatly reduced, however, shoot xylem pressure potential and carbon isotope discrimination were higher and water-use efficiency was lower, inconsistent with water stress. Reductions in assimilation rate could not be ascribed to increased stomatal limitation. Starch and soluble sugars analysis of leaves revealed that photosynthate was accumulating to high levels, suggesting that the trees with substantially reduced cell wall lignin were not carbon limited and that reductions in sink strength were, instead, limiting photosynthesis.

Sectional analysis of stormwater treatment performance of a constructed wetland

Constructed wetlands are among the most common Water Sensitive Urban Design (WSUD) measures for stormwater treatment. These systems have been extensively studied to understand their performance and influential treatment processes. Unfortunately, most past studies have been undertaken considering a wetland system as a lumped system with a primary focus on the reduction of the event mean concentration (EMC) values of specific pollutant species or total pollutant load removal. This research study adopted an innovative approach by partitioning the inflow runoff hydrograph and then investigating treatment performance in each partition and their relationships with a range of hydraulic factors. The study outcomes confirmed that influenced by rainfall characteristics, the constructed wetland displays different treatment characteristics for the initial and later sectors of the runoff hydrograph. The treatment of small rainfall events (<15 mm) is comparatively better at the beginning of runoff events while the trends in pollutant load reductions for large rainfall events (>15 mm) are generally lower at the beginning and gradually increase towards the end of rainfall events. This highlights the importance of ensuring that the inflow into a constructed wetland has low turbulence in order to achieve consistent treatment performance for both, small and large rainfall events.

Microscale consolidation analysis of relaxation behavior of single living chondrocytes subjected to varying strain-rates

Besides the elastic stiffness, the relaxation behavior of single living cells is also of interest of various researchers when studying cell mechanics. It is hypothesized that the relaxation response of the cells is governed by both intrinsic viscoelasticity of the solid phase and fluid-solid interactions mechanisms. There are a number of mechanical models have been developed to investigate the relaxation behavior of single cells. However, there is lack of model enable to accurately capture both of the mechanisms. Therefore, in this study, the porohyperelastic (PHE) model, which is an extension of the consolidation theory, combined with inverse Finite Element Analysis (FEA) technique was used at the first time to investigate the relaxation response of living chondrocytes. This model was also utilized to study the dependence of relaxation behavior of the cells on strain-rates. The stress-relaxation experiments under the various strain-rates were conducted with the Atomic Force Microscopy (AFM). The results have demonstrated that the PHE model could effectively capture the stress-relaxation behavior of the living chondrocytes, especially at intermediate to high strain-rates. Although this model gave some errors at lower strain-rates, its performance was acceptable. Therefore, the PHE model is properly a promising model for single cell mechanics studies. Moreover, it has been found that the hydraulic permeability of living chondrocytes reduced with decreasing of strain-rates. It might be due to the intracellular fluid volume fraction and the fluid pore pressure gradients of chondrocytes were higher when higher strain-rates applied.