Assessment of interactions between alluvial, Great Artesian Basin (GAB) and volcanic aquifers using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia

A detailed 3D lithological model framework was developed using GOCAD software to understand interactions between alluvial, volcanic and GAB aquifers and the spatial and temporal distribution of groundwater recharge to the alluvium of the Lockyer Valley. Groundwater chemistry, isotope data (H20-δ2H and δ18O , 87Sr/86Sr, 3H and 14C) and groundwater level time-series data from approximately 550 observation wells were integrated into the catchment-wide 3D model to assess the recharge processes involved. This approach enabled the identification of zones where recharge to the alluvium primarily occurs from stream water during episodic flood events. Importantly, the study also demonstrates that in some sections of the alluvium recharge is also from storm rainfall and seepage discharge from the underlying GAB aquifers. These other sources of recharge are indicated by (a) the absence of a response of groundwater levels to flooding in some areas, (b) old radiocarbon ages, and (c) distinct bedrock water chemistry and δ2H and δ18O signatures in alluvial groundwater at these locations. Integration of isotopes, water chemistry and time-series displays of groundwater levels before and after the 2010/2011 flood into the 3D model suggest that the spatial variations in the alluvial groundwater response are mostly controlled by valley morphology and lithological (i.e. permeability) variations within the alluvium. Examination of the groundwater level variations in the 3D model also enabled quantification of the volumetric change of groundwater stored in the unconfined alluvial aquifer prior to and post-flood events.

Visualisation modelling in 3D of an alluvial aquifer system at a valley-wide scale to understand groundwater status : Lockyer, QLD

Visualisation provides a method to efficiently convey and understand the complex nature and processes of groundwater systems. This technique has been applied to the Lockyer Valley to aid in comprehending the current condition of the system. The Lockyer Valley in southeast Queensland hosts intensive irrigated agriculture sourcing groundwater from alluvial aquifers. The valley is around 3000 km2 in area and the alluvial deposits are typically 1-3 km wide and to 20-35 m deep in the main channels, reducing in size in subcatchments. The configuration of the alluvium is of a series of elongate “fingers”. In this roughly circular valley recharge to the alluvial aquifers is largely from seasonal storm events, on the surrounding ranges. The ranges are overlain by basaltic aquifers of Tertiary age, which overall are quite transmissive. Both runoff from these ranges and infiltration into the basalts provided ephemeral flow to the streams of the valley. Throughout the valley there are over 5,000 bores extracting alluvial groundwater, plus lesser numbers extracting from underlying sandstone bedrock. Although there are approximately 2500 monitoring bores, the only regularly monitored area is the formally declared management zone in the lower one third. This zone has a calibrated Modflow model (Durick and Bleakly, 2000); a broader valley Modflow model was developed in 2002 (KBR), but did not have extensive extraction data for detailed calibration. Another Modflow model focused on a central area river confluence (Wilson, 2005) with some local production data and pumping test results. A recent subcatchment simulation model incorporates a network of bores with short-period automated hydrographic measurements (Dvoracek and Cox, 2008). The above simulation models were all based on conceptual hydrogeological models of differing scale and detail.

Use of 3D visualisation to enhance groundwater stable isotope data and its interpretation, Lockyer Valley, southeast Queensland

The use of stable isotope ratios δ18O and δ2H are well established in assessment of groundwater systems and their hydrology. The conventional approach is based on x/y plots and relation to various MWL’s, and plots of either ratio against parameters such as Clor EC. An extension of interpretation is the use of 2D maps and contour plots, and 2D hydrogeological vertical sections. An enhancement of presentation and interpretation is the production of “isoscapes”, usually as 2.5D surface projections. We have applied groundwater isotopic data to a 3D visualisation, using the alluvial aquifer system of the Lockyer Valley. The 3D framework is produced in GVS (Groundwater Visualisation System). This format enables enhanced presentation by displaying the spatial relationships and allowing interpolation between “data points” i.e. borehole screened zones where groundwater enters. The relative variations in the δ18O and δ2H values are similar in these ambient temperature systems. However, δ2H better reflects hydrological processes, whereas δ18O also reflects aquifer/groundwater exchange reactions. The 3D model has the advantage that it displays borehole relations to spatial features, enabling isotopic ratios and their values to be associated with, for example, bedrock groundwater mixing, interaction between aquifers, relation to stream recharge, and to near-surface and return irrigation water evaporation. Some specific features are also shown, such as zones of leakage of deeper groundwater (in this case with a GAB signature). Variations in source of recharging water at a catchment scale can be displayed. Interpolation between bores is not always possible depending on numbers and spacing, and by elongate configuration of the alluvium. In these cases, the visualisation uses discs around the screens that can be manually expanded to test extent or intersections. Separate displays are used for each of δ18O and δ2H and colour coding for isotope values.

Analysis of productive performance of crop and animal production systems : an integrated analytical framework

This article presents a two-stage analytical framework that integrates ecological crop (animal) growth and economic frontier production models to analyse the productive efficiency of crop (animal) production systems. The ecological crop (animal) growth model estimates "potential" output levels given the genetic characteristics of crops (animals) and the physical conditions of locations where the crops (animals) are grown (reared). The economic frontier production model estimates "best practice" production levels, taking into account economic, institutional and social factors that cause farm and spatial heterogeneity. In the first stage, both ecological crop growth and economic frontier production models are estimated to calculate three measures of productive efficiency: (1) technical efficiency, as the ratio of actual to "best practice" output levels; (2) agronomic efficiency, as the ratio of actual to "potential" output levels; and (3) agro-economic efficiency, as the ratio of "best practice" to "potential" output levels. Also in the first stage, the economic frontier production model identifies factors that determine technical efficiency. In the second stage, agro-economic efficiency is analysed econometrically in relation to economic, institutional and social factors that cause farm and spatial heterogeneity. The proposed framework has several important advantages in comparison with existing proposals. Firstly, it allows the systematic incorporation of all physical, economic, institutional and social factors that cause farm and spatial heterogeneity in analysing the productive performance of crop and animal production systems. Secondly, the location-specific physical factors are not modelled symmetrically as other economic inputs of production. Thirdly, climate change and technological advancements in crop and animal sciences can be modelled in a "forward-looking" manner. Fourthly, knowledge in agronomy and data from experimental studies can be utilised for socio-economic policy analysis. The proposed framework can be easily applied in empirical studies due to the current availability of ecological crop (animal) growth models, farm or secondary data, and econometric software packages. The article highlights several directions of empirical studies that researchers may pursue in the future.

Assessment of flood-related recharge to alluvial aquifers of an irrigated catchment following extended drought conditions using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia.

The Lockyer Valley in southeast Queensland, Australia, hosts an economically significant alluvial aquifer system which has been impacted by prolonged drought conditions (~1997 to ~ 2009). Throughout this time, the system was under continued groundwater extraction, resulting in severe aquifer depletion. By 2008, much of the aquifer was at <30% of storage but some relief occurred with rains in early 2009. However, between December 2010 and January 2011, most of southeast Queensland experienced unprecedented flooding, which generated significant aquifer recharge. In order to understand the spatial and temporal controls of groundwater recharge in the alluvium, a detailed 3D lithological property model of gravels, sands and clays was developed using GOCAD software. The spatial distribution of recharge throughout the catchment was assessed using hydrograph data from about 400 groundwater observation wells screened at the base of the alluvium. Water levels from these bores were integrated into a catchment-wide 3D geological model using the 3D geological modelling software GOCAD; the model highlights the complexity of recharge mechanisms. To support this analysis, groundwater tracers (e.g. major and minor ions, stable isotopes, 3H and 14C) were used as independent verification. The use of these complementary methods has allowed the identification of zones where alluvial recharge primarily occurs from stream water during episodic flood events. However, the study also demonstrates that in some sections of the alluvium, rainfall recharge and discharge from the underlying basement into the alluvium are the primary recharge mechanisms of the alluvium. This is indicated by the absence of any response to the flood, as well as the observed old radiocarbon ages and distinct basement water chemistry signatures at these locations. Within the 3D geological model, integration of water chemistry and time-series displays of water level surfaces before and after the flood suggests that the spatial variations of the flood response in the alluvium are primarily controlled by the valley morphology and lithological variations within the alluvium. The integration of time-series of groundwater level surfaces in the 3D geological model also enables the quantification of the volumetric change of groundwater stored in the unconfined sections of this alluvial aquifer during drought and following flood events. The 3D representation and analysis of hydraulic and recharge information has considerable advantages over the traditional 2D approach. For example, while many studies focus on singular aspects of catchment dynamics and groundwater-surface water interactions, the 3D approach is capable of integrating multiple types of information (topography, geological, hydraulic, water chemistry and spatial) into a single representation which provides valuable insights into the major factors controlling aquifer processes.

Poorly characterized phases in C2M carbonaceous chondrites : proposed structures and significance

Poorly characterized phases (PCP's) may constitute up to 30 volume percent of some C2M carbonaceous chondrite matrices [1] and are an important key to an understanding of matrix evolution. PCPs are usually fine-grained (

Proposed structures for poorly characterized phases in C2M carbonaceous chondrite meteorites

Structure and chemistry of poorly characterized phases (PCP). We suggest here that approximately 10 angstrom PCP, a dominant matrix variety, has a structure equivalent to iron-rich tochilinite [6Fe (sub 0.9) S 5(Fe, Mg) (OH) (sub 2) ] which consists of coherently interstratified mackinawite and brucite sheets. approximately 17 angstrom PCP, previously described as an SBB-type mixed-layer structure, is a commensurate intergrowth of serpentine and tochilinite layers. A wide range of cation substitutions is possible within both tochilinite and serpentine-tochilinite structural types. Various forms of PCP observed in carbonaceous chondrites are intergrowths of tochilinite, serpentine, serpentine-tochilinite and/or valleriite-type minerals.--Modified journal abstract.

Poorly graphitized carbon as a new cosmothermometer for primitive extraterrestrial materials

The presence of carbon in primitive extraterrestrial materials has long been considered a useful indicator of prevailing geochemical conditions early in the formation of the Solar System. A recent addition to the suite of primitive materials available for study by cosmochemists includes particles collected from the stratosphere called chondritic porous (CP) aggregates1. Carbon-rich CP aggregates are less abundant in stratospheric collections and contain many low-temperature phases (such as layer silicates) as minor components2,3. We describe here the nature of the most abundant carbon phase in a carbon-rich CP aggregate (sample no. W7029* A) collected from the stratosphere as part of the Johnson Space Center (JSC) Cosmic Dust Program4. By comparison with experimental and terrestrial studies of poorly graphitized carbon (PGC), we show that the graphitization temperature, or the degree of ordering in the PGC, may provide a useful cosmothermometer for primitive extraterrestrial materials.

Planning for transit system reliability using productive performance and risk assessment

Urban transit system performance may be quantified and assessed using transit capacity and productive capacity for planning, design and operational management. Bunker (4) defines important productive performance measures of an individual transit service and transit line. Transit work (p-km) captures transit task performed over distance. Transit productiveness (p-km/h) captures transit work performed over time. This paper applies productive performance with risk assessment to quantify transit system reliability. Theory is developed to monetize transit segment reliability risk on the basis of demonstration Annual Reliability Event rates by transit facility type, segment productiveness, and unit-event severity. A comparative example of peak hour performance of a transit sub-system containing bus-on-street, busway, and rail components in Brisbane, Australia demonstrates through practical application the importance of valuing reliability. Comparison reveals the highest risk segments to be long, highly productive on street bus segments followed by busway (BRT) segments and then rail segments. A transit reliability risk reduction treatment example demonstrates that benefits can be significant and should be incorporated into project evaluation in addition to those of regular travel time savings, reduced emissions and safety improvements. Reliability can be used to identify high risk components of the transit system and draw comparisons between modes both in planning and operations settings, and value improvement scenarios in a project evaluation setting. The methodology can also be applied to inform daily transit system operational management.

Planning for transit system reliability using productive performance and risk assessment

Urban transit system performance may be quantified and assessed using transit capacity and productive capacity for planning, design and operational management. Bunker (4) defines important productive performance measures of an individual transit service and transit line. Transit work (p-km) captures transit task performed over distance. Transit productiveness (p-km/h) captures transit work performed over time. This paper applies productive performance with risk assessment to quantify transit system reliability. Theory is developed to monetize transit segment reliability risk on the basis of demonstration Annual Reliability Event rates by transit facility type, segment productiveness, and unit-event severity. A comparative example of peak hour performance of a transit sub-system containing bus-on-street, busway, and rail components in Brisbane, Australia demonstrates through practical application the importance of valuing reliability. Comparison reveals the highest risk segments to be long, highly productive on street bus segments followed by busway (BRT) segments and then rail segments. A transit reliability risk reduction treatment example demonstrates that benefits can be significant and should be incorporated into project evaluation in addition to those of regular travel time savings, reduced emissions and safety improvements. Reliability can be used to identify high risk components of the transit system and draw comparisons between modes both in planning and operations settings, and value improvement scenarios in a project evaluation setting. The methodology can also be applied to inform daily transit system operational management.

Analysis of productive performance of crop production systems: an integrated analytical framework

Recent literature has argued that environmental efficiency (EE), which is built on the materials balance (MB) principle, is more suitable than other EE measures in situations where the law of mass conversation regulates production processes. In addition, the MB-based EE method is particularly useful in analysing possible trade-offs between cost and environmental performance. Identifying determinants of MB-based EE can provide useful information to decision makers but there are very few empirical investigations into this issue. This article proposes the use of data envelopment analysis and stochastic frontier analysis techniques to analyse variation in MB-based EE. Specifically, the article develops a stochastic nutrient frontier and nutrient inefficiency model to analyse determinants of MB-based EE. The empirical study applies both techniques to investigate MB-based EE of 96 rice farms in South Korea. The size of land, fertiliser consumption intensity, cost allocative efficiency, and the share of owned land out of total land are found to be correlated with MB-based EE. The results confirm the presence of a trade-off between MB-based EE and cost allocative efficiency and this finding, favouring policy interventions to help farms simultaneously achieve cost efficiency and MP-based EE.

Five commercial citrus rate poorly as hosts of the polyphagous fruit flyBactrocera tryoni(Froggatt) (Diptera: Tephritidae) in laboratory studies

Queensland fruit fly, Bactrocera tryoni (Froggatt), is a polyphagous pest, and many citrus types are included among its hosts. While quantification of citrus host use by B. tryoni is lacking, citrus is generally considered a ‘low pressure’ crop. This paper investigates B. tryoni female oviposition preference and offspring performance in five citrus types; Murcott mandarin (Citrus reticulata), Navel orange and Valencia orange (Citrus sinensis), Eureka lemon (Citrus limon) and yellow grapefruit (Citrus paradisi). Oviposition preference was investigated in laboratory-based choice and no-choice experiments, while immature survival and offspring performance were investigated by infesting fruits in the laboratory and evaluating pupal recovery, pupal emergence and F1 fecundity. Fruit size, Brix level and peel toughness were also measured for correlation with host use. Bactrocera tryoni demonstrated an oviposition preference hierarchy among the citrus fruits tested; Murcott and grapefruit were most preferred for oviposition and lemon the least, while preference for Navel and Valencia was intermediate. Peel toughness was negatively correlated with B. tryoni oviposition preference, while no significant correlations were detected between oviposition and Brix level or fruit size. Immature survival in the tested fruit was very low. Murcott was the best host (21% pupal recovery), while all other citrus types that showed pupal recovery of 6% or lower and no pupae were recovered from Valencia orange. In pupae recovered from Navel orange and lemon, adult eclosion was greatly reduced, while in grapefruit and lemon, no eggs were recovered from F1 adults. Based on these laboratory results, many commercial citrus varieties appear to be poor hosts for B. tryoni and may pose a low post-harvest and quarantine risk. These findings need to be confirmed in the field, as they impact on both pre-harvest and post-harvest countermeasures.

Modelling sea water intrusion in coastal aquifers using heterogeneous computing

The objective of this PhD research program is to investigate numerical methods for simulating variably-saturated flow and sea water intrusion in coastal aquifers in a high-performance computing environment. The work is divided into three overlapping tasks: to develop an accurate and stable finite volume discretisation and numerical solution strategy for the variably-saturated flow and salt transport equations; to implement the chosen approach in a high performance computing environment that may have multiple GPUs or CPU cores; and to verify and test the implementation. The geological description of aquifers is often complex, with porous materials possessing highly variable properties, that are best described using unstructured meshes. The finite volume method is a popular method for the solution of the conservation laws that describe sea water intrusion, and is well-suited to unstructured meshes. In this work we apply a control volume-finite element (CV-FE) method to an extension of a recently proposed formulation (Kees and Miller, 2002) for variably saturated groundwater flow. The CV-FE method evaluates fluxes at points where material properties and gradients in pressure and concentration are consistently defined, making it both suitable for heterogeneous media and mass conservative. Using the method of lines, the CV-FE discretisation gives a set of differential algebraic equations (DAEs) amenable to solution using higher-order implicit solvers. Heterogeneous computer systems that use a combination of computational hardware such as CPUs and GPUs, are attractive for scientific computing due to the potential advantages offered by GPUs for accelerating data-parallel operations. We present a C++ library that implements data-parallel methods on both CPU and GPUs. The finite volume discretisation is expressed in terms of these data-parallel operations, which gives an efficient implementation of the nonlinear residual function. This makes the implicit solution of the DAE system possible on the GPU, because the inexact Newton-Krylov method used by the implicit time stepping scheme can approximate the action of a matrix on a vector using residual evaluations. We also propose preconditioning strategies that are amenable to GPU implementation, so that all computationally-intensive aspects of the implicit time stepping scheme are implemented on the GPU. Results are presented that demonstrate the efficiency and accuracy of the proposed numeric methods and formulation. The formulation offers excellent conservation of mass, and higher-order temporal integration increases both numeric efficiency and accuracy of the solutions. Flux limiting produces accurate, oscillation-free solutions on coarse meshes, where much finer meshes are required to obtain solutions with equivalent accuracy using upstream weighting. The computational efficiency of the software is investigated using CPUs and GPUs on a high-performance workstation. The GPU version offers considerable speedup over the CPU version, with one GPU giving speedup factor of 3 over the eight-core CPU implementation.