Murray Basin Research Articles

Groundwater exhibits spatially opposing trends during the Australian Millennium Drought

The adverse impacts of Australia’s Millennium Drought on both surface and groundwater hydrological systems are extensively documented. During the Millennium Drought, the Murray Basin experienced a severe rainfall deficit. Our study revisited groundwater table trends in 451 wells within the Murray Basin during the drought from 1997 to 2009. These trends varied, 70% showed significant downward shifts, 19% were insignificant, and 11% even displayed upward trends. The results from K-means clustering analysis indicate a markedly slow recuperation of groundwater levels post-drought. We used multiple regression models to link interannual groundwater dynamics with climate variables, revealing climate as the primary driver of declining groundwater levels. This connection is influenced by land cover and thickness of the vadose zone, resulting in hysteresis effects and spatial variations. In cases with a thick vadose zone and minimal evapotranspiration, the influence of the Millennium Drought on the groundwater system is reduced. The increasing trends may also be related to lateral recharge from mountainous areas, human activities in adjacent irrigation districts, and east-west geostress. Our findings reveal the complex interactions between climate, land characteristics, and groundwater behavior during and after the Millennium Drought, holding significant implications for understanding hydrological processes under extreme drought conditions and for the sustainable management of water resources.

Sediment thickness across Australia from passive seismic methods

SUMMARY Around the world the Earth’s crust is blanketed to various extents by sediment. For continental regions, knowledge of the distribution and thickness of sediments is crucial for a wide range of applications including seismic hazard, resource potential and our ability to constrain the deeper crustal geology. Excellent constraints on the sediment thickness can be obtained from borehole drilling or active seismic surveys. However, these approaches are expensive and impractical in remote continental interiors such as central Australia. Recently, a method for estimating the sediment thickness using passive seismic data, the collection of which is relatively simple and low cost, was developed and applied to seismic stations in South Australia. This method uses receiver functions, specifically the time delay of the P-to-S converted phase generated at the sediment-basement interface, relative to the direct-P arrival, to generate a first-order estimate of the thickness of sediments. In this work, we expand the analysis to the vast array of over 1500 seismic stations across Australia, covering an entire continent and numerous sedimentary basins that span the entire range from Precambrian to present day. We compare with an established yet separate method to estimate the sediment thickness, which utilizes the autocorrelation of the radial receiver functions to ascertain the two-way traveltime of shear waves reverberating in a sedimentary layer. Across the Australian continent the new results match the broad pattern of expected sedimentary features based on the various geological provinces. We are able to delineate the boundaries of many sedimentary basins, such as the Eucla and Murray Basins, which are Cenozoic, and the boundary between the Karumba Basin and the mineral rich Mount Isa Province. Contrasts in seismic delay time across these boundaries are upwards of 0.4 $\, \mathrm{s}$. The delay signal is found to diminish to <0.1 $\, \mathrm{s}$ for older Proterozoic basins, likely due to compaction and metamorphism of the sediments over time. As an application of the method, a comparison with measurements of sediment thickness from local boreholes allows for a straightforward predictive relationship between the delay time and the cover thickness to be defined. This offers future widespread potential, providing a simple and cheap way to characterize the sediment thickness in underexplored areas from passive seismic data.

Provenance of the Copi North heavy mineral sands deposit, Murray Basin, Australia, based on the geochronology and geochemistry of detrital zircon, rutile and monazite

The Murray Basin hosts one of the world’s largest strandplain systems and has enormous potential for heavy mineral sand resources of titanium, zircon and rare earth elements. Further discovery of heavy mineral sand deposits can be aided via improved understanding of the source and transport of detrital heavy minerals of the basin. In this study, the provenance of two heavy mineral sand units of the Copi North deposit (western NSW) was determined using U–Pb geochronology and geochemistry of detrital zircon, rutile and monazite. Results from the two units are similar, indicating they derive from the same sources. Approximately 42% of the zircon population is aged <470 Ma and is interpreted to have been derived from the Tasmanides. Older zircon populations, up to ca 2.5 Ga, were largely derived from the Adelaide Superbasin and Delamerian Orogen, or from metasedimentary rocks of the Lachlan Orogen. Over 90% of rutile U–Pb ages are between 570 and 450 Ma, correlating either to Delamerian metamorphism or to detrital rutile from the Darling Basin; both sources are supported by rutile geochemistry that is indicative of a metasedimentary source. Most monazites have U–Pb ages of ca 600 Ma, which is consistent with an ultimate source from the Mt Arrowsmith Volcanics and/or Antarctica. A minor proportion of the monazite was sourced from igneous rocks of the New England Orogen. Temporal changes in sediment source regions can be inferred when considering the detrital mineral record for Copi North with other mineral sand deposits of the basin. Older (>ca 6.7 Ma) deposits such as Copi North feature a significant proportion of post-350 Ma detritus, reflecting the important role of the ancestral Darling River in supplying sediment from the north at the earliest stages of Murray Basin development. Post-6.7 Ma deposits lack a major component of detritus for far northern sources, and instead consist of sediment derived dominantly from western (Adelaide Superbasin) and eastern (Lachlan Orogen) sources, the latter delivered via the Murray River.

Multi-source precipitation products assessment on drought monitoring across global major river basins

Increased availability of precipitation observational datasets allows for accurate drought detection, monitoring and forecast. Most studies have shown accuracy of various precipitation datasets on regional drought detection and/or monitoring, how they perform across different global river basins, however, is less understood. Here, we evaluate five popular precipitation products including GPM-IMERG, ERA5, GLADS, CPC and PREC/L for meteorological drought detection in 46 large river basins globally, by employing the Standardized Precipitation Index (SPI) at 3-month time scale (SPI-3). The results show that there is a large performance difference among these products. On the whole, GPM-IMERG and PREC/L are ideal precipitation datasets for drought monitoring in North America, South America and Asia, while GLDAS is better in Europe. In terms of specific mega drought events over major basins, all precipitation products perform well in capturing mega drought events overall, especially for the Yangtze, Danube, and Colorado river basins, with PREC/L having the highest accuracy. Specifically, PREC/L well captures the centroid and extent of the typical mega drought events over the Orange and Murray basins. The outcomes of this study can provide a deeper insight into applicability of various precipitation datasets in drought monitoring for different global river basins, which guides scientists and engineers to select suitable precipitation products and further better manage drought disasters.

Past and future rainfall change in sub-regions of Victoria, Australia

We examine rainfall variability and change in three sub-regions of the state of Victoria in Australia: the Murray Basin Victoria (MBVic), southeast Victoria (SEVic), and southwest Victoria (SWVic). These sub-regions represent three different hydrological super-catchments over Victoria and received average cool season rainfall for the 1997–2018 period, about 15%, 11%, and 8% less, respectively, than the 1900–1959 average. All three observed declines are shown to be very unusual in terms of historical variability. On analysing CMIP5 models under different forcing conditions (preindustrial, historical-all, historical-GHGs-only, historical-natural-only, RCP2.6, RCP4.5, and RCP8.5), we estimate that external forcing caused 30% of the observed drying in SWVic, 18% in MBVic, and 17% in SEVic. The external forcing contributions to the observed trend for the 1900–2018 period are estimated to be 56%, 17%, and 24% for SWVic, MBVic, and SEVic, respectively. Taken at face value, these figures suggest that only the 1900–2018 trend in SWVic was dominated by external forcing. Nearly all models underestimate the magnitude of the observed drying. This arises because models underestimate the magnitude of decadal variability, and because models might also underestimate externally forced drying, and/or the contribution of internal variability in the real world to the observed event was unusually large. By 2037, approximately 90% of the models simulate drying in SWVic, even under a low emissions scenario. Under a high emission scenario, the anthropogenically forced drying towards the late twenty-first century is so large in all three sub-regions that internal variability appears too small to offset it.

A Comparison of Granite Genesis in the Adelaide Fold Belt and Glenelg River Complex Using U–Pb, Hf and O Isotopes in Zircon

Abstract We present new U–Pb ages and Hf and O isotope data for zircon from I-, S- and A-type granites from both the western and eastern edges of the Delamerian Orogen in southeastern Australia. The I-type Tanunda Creek Gneiss contains zircon populations of 507 ± 4 and 492 ± 6 Ma inferred to reflect igneous and metamorphic ages, respectively. The I-type Palmer Granite yielded an age of 509 ± 3 Ma, and the Port Elliot S-type Granite has a magmatic age of 508 ± 7 Ma. Inherited zircon in these granites range from 1092 to 3343 Ma, probably derived from assimilation of Adelaide Group sediments. The Murray Bridge A-type Granite is 490 ± 2 Ma in age and lacks inherited zircon. In the Glenelg River Complex, an S-type migmatite from near Harrow contains a complex zircon population. It is most likely ~500 Ma in age and has inherited zircon of 550–700, 1000–1100 and 2437 Ma, hence matching those from the Kanmantoo Group. From this and detrital zircon ages, we infer that only the Kanmantoo Group extends across the Murray Basin into the Glenelg River Complex. The Wando Tonalite and Loftus Creek I-type granites yielded ages of 501 ± 2 and 486 ± 3 Ma, respectively. Zircon from the Dergholm Granite has suffered Pb loss, and the best age estimate for this granite is 488 ± 5 Ma. Combining all the granite data together, εHft and δ18O in the magmatic zircon range from 5.6 to −10.3 and from 5.8 to 8.1, respectively, and are well correlated. The zircon indicates the same temporal and compositional evolution of granitic petrogenesis across ~300 km of strike, reaffirming the notion that these terranes form part of the same orogen. Westward-directed subduction caused orogenic thickening, heating and increasing amounts of crustal contribution. This was followed by convective thinning of the thickened mantle lithosphere and a return to more primitive magmas lacking significant crustal contributions. It contrasts significantly with inferred granite petrogenesis and tectonic style in the younger Lachlan and New England Fold Belts further east that were not built upon extended cratonic lithosphere.

A new volute, Ericusa ngayawang sp. nov. (Gastropoda: Volutidae), from the Miocene of South Australia.

Ericusa ngayawang sp. nov. is described from shells preserved in the Middle Miocene Cadell Formation in the western Murray Basin of South Australia. At the time the Murray Basin was part of the Southeastern Australian Marine Biogeographic Province. Ericusa ngayawang is a small heavily costate species of Ericusa with clear affinities to the Early Miocene E. atkinsoni of Victoria and Tasmania but can be distinguished from it by its smaller size, more slender proportions and its heavily costate body whorl. Ericusa atkinsoni and its relative, E. macroptera, inhabited the basins to the east of the Murray Basin during the Late Oligocene and Early Miocene but were extinct there before the end of the Burdigalian Stage of the Early Miocene. The persistence of E. ngayawang into the Langhian Stage of the Middle Miocene is another piece of evidence for partial biogeographic isolation of the western Murray Basin from the rest of the Southeastern Australian Province during the Miocene.

Latest Miocene ostracods from the Bookpurnong Formation in the Murray Basin of southeastern Australia: shallow marine migrants into an epicontinental sea

The uppermost Miocene Bookpurnong Formation within the Murray Basin of southeastern Australia overlies a regionally extensive subaerial unconformity formed by relatively low late Miocene eustatic sea levels, and the initial phase of the Kosciuszko Uplift tectonic event. A diverse marine fossil ostracod fauna has been recovered from the Bookpurnong Formation, and is associated with a marine transgression that flooded inland regions of southeastern Australia to form a shallow epicontinental sea. Many of the Bookpurnong Formation ostracods represent immigrant taxa, with species such as Puriana lubbockiana, evidencing a subtropical range expansion of thermophilic warm water forms into southern mid-latitudes. We attribute this to warm plumes from the East Australian Current, which would have impacted southeastern Australia at that time. In general, the Bookpurnong Formation ostracod assemblages indicate low to moderate energy shallow offshore palaeoenvironments subject to warm-temperate and subtropical conditions. One new genus and five new species are described: Fortistriginglymus gen. nov., Bradyleberis praecristatella sp. nov., Callistocythere bookpurnongensis sp. nov., Callistocythere mchenryi sp. nov., Callistocythere zigzaga sp. nov., and Parakeijia notoreticularis sp. nov.

A strontium isoscape of inland southeastern Australia

Abstract. The values and distribution patterns of the strontium (Sr) isotope ratio 87Sr/86Sr in Earth surface materials are of use in the geological, environmental, and social sciences. Ultimately, the 87Sr/86Sr ratios of soils and everything that lives in and on them are inherited from the rocks that are the parent materials of the soil's components. In Australia, there are few large-scale surveys of 87Sr/86Sr available, and here we report on a new, low-density dataset using 112 catchment outlet (floodplain) sediment samples covering 529 000 km2 of inland southeastern Australia (South Australia, New South Wales, Victoria). The coarse (&lt;2 mm) fraction of bottom sediment samples (depth ∼ 0.6–0.8 m) from the National Geochemical Survey of Australia were milled and fully digested before Sr separation by chromatography and 87Sr/86Sr determination by multicollector-inductively coupled plasma mass spectrometry. The results show a wide range of 87Sr/86Sr values from a minimum of 0.7089 to a maximum of 0.7511 (range 0.0422). The median 87Sr/86Sr (± median absolute deviation) is 0.7199 (± 0.0071), and the mean (± standard deviation) is 0.7220 (± 0.0106). The spatial patterns of the Sr isoscape observed are described and attributed to various geological sources and processes. Of note are the elevated (radiogenic) values (≥∼ 0.7270; top quartile) contributed by (1) the Palaeozoic sedimentary country rock and (mostly felsic) igneous intrusions of the Lachlan geological region to the east of the study area; (2) the Palaeoproterozoic metamorphic rocks of the central Broken Hill region; both these sources contribute radiogenic material mainly by fluvial processes; and (3) the Proterozoic to Palaeozoic rocks of the Kanmantoo, Adelaide, Gawler, and Painter geological regions to the west of the area; these sources contribute radiogenic material mainly by aeolian processes. Regions of low 87Sr/86Sr (≤∼ 0.7130; bottom quartile) belong mainly to (1) a few central Murray Basin catchments; (2) some Darling Basin catchments in the northeast; and (3) a few Eromanga geological region-influenced catchments in the northwest of the study area; these sources contribute unradiogenic material mainly by fluvial processes. The new spatial Sr isotope dataset for the DCD (Darling–Curnamona–Delamerian) region is publicly available (de Caritat et al., 2022; https://dx.doi.org/10.26186/146397)​​​​​​​.

The impact of effective tortuosity on carbonate acidizing and the validation of Damköhler and Péclet dimensionless phase space

Acid stimulation is one of the main methodologies to enhance well productivity and increase permeability lost by near-wellbore formation damage through matrix acidizing. The desired outcome of acidizing carbonate rocks is to form deep conductive flow channels called “wormholes”. The productivity enhancement depends on how deep these wormholes penetrate into the reservoir. The porous medium is characterized by three fundamental parameters which are anisotropic or isotropic porosity, permeability, and percolation. Here we test the hypothesis that effective tortuosity - a key parameter for microstructure characterization in porous media - has a significant effect on wormhole formation. Significant progress has been made by empirical work on carbonate acidizing to increase the stimulation efficiency, but further improvements are possible by developing a deeper understanding of the impact of effective tortuosity on wormhole initiation and development. We use Indiana limestone core-flooding experiments to inject diluted hydrochloric acid (HCl) at different acid concentrations and injection rates. Effective tortuosity was calculated for the tight Indiana limestone and was compared with an earlier study on the highly porous Mount Gambier limestone to evaluate the effect of tortuosity on wormhole formation. Our results suggest that the inverse Damköhler and Péclet phase space, traditionally used for evaluating criteria for wormhole formation, needs to be extended by tortuosity as a third dimension to capture the effect of microstructure for the initiation of wormhole phenomenon. Incorporating this extension will improve the prediction of wormhole initiation conditions for all carbonate rock types by using a universal dimensionless phase space.

Comparison of Petrophysical Properties of Porous Rocks Using NMR, Micro-CT, and Fluid Flow Simulations

Micro-computed tomography (micro-CT) is increasingly utilized to image the pore network and to derive petrophysical properties in combination with modelling software. The effect of micro-CT image resolution and size on the accuracy of the derived petrophysical properties is addressed in this study using a relatively homogenous sandstone and a heterogenous, highly porous bioclastic limestone. Standard laboratory procedures including NMR (nuclear magnetic resonance) analysis, micro-CT analysis at different image resolutions and sizes and pore-scale flow simulations were used to determine and compare petrophysical properties. NMR-derived pore-size distribution (PSD) was comparable to the micro-CT-derived PSD at a resolution of 7 µm for both the rock types. Porosity was higher using the water saturation method as compared to the NMR method in both rocks. The resolution did not show a significant effect on the porosity of the homogeneous sandstone, but porosity in the heterogeneous limestone varies depending on the location of the sub-sample. The transport regime in the sandstone was derived by simulations and changed with the resolution of the micro-CT image. The transport regime in the sandstone was advection-dominated at higher image resolution and diffusion-dominated when using a lower image resolution. In contrast, advection was the dominant transport regime for the limestone based on simulations using higher and lower image resolutions. Simulation-derived permeability for a 400 Voxel3 image at 7 µm resolution in the Berea sandstone matched laboratory results, although local heterogeneity within the rock plays an integral role in the permeability estimation within the sub-sampled images. The simulation-derived permeability was highly variable in the Mount Gambier limestone depending on the image size and resolution with the closest value to a laboratory result simulated with an image resolution of 2.5 µm and a size of 300 Voxel3. Overall, the study demonstrates the need to decide on micro-CT parameters depending on the type of petrophysical property of interest and the degree of heterogeneity within the rock types.

Revisiting the abandoned shorelines of Lake George, Australia: a refined optical dating framework

ABSTRACTAbandoned shorelines are an important archive used to constrain past fluctuations in the hydrological balance of lakes around the globe. Within Australia, the shorelines preserved at Lake George, NSW, form one of the few shoreline archives in the south‐east of the continent that record palaeoenvironmental conditions throughout the late Quaternary. Here, we examined and tested the lake‐level record for Lake George constructed in the 1970s by dating a well‐preserved shoreline sequence at Luckdale, on the lake's eastern shore, using single‐grain optically stimulated luminescence (OSL) dating. Ten stratigraphic units were identified, and these suggest a late Quaternary highstand for Lake George in MIS 3, with fluctuations superimposed upon an overall drying trend throughout MIS 2 and into the present. At Luckdale, the highest four shoreline‐associated units were deposited ~13 to 19 m above lake base and date to between 39 ± 2 and 29 ± 1 ka ago. Our study pushes back the timing of maximum lake depth at Lake George to at least MIS 3, rather than MIS 2. The overall drying trend is supported by similar reductions in both Riverine Plain fluvial activity and other associated lake‐level records from within the Murray basin.

3‐D Sedimentary Structures Beneath Southeastern Australia Constrained by Passive Seismic Array Data

AbstractA high‐resolution sedimentary map of southeastern Australia is of great importance to explore natural resources as well as to evaluate seismic hazards. In this study, we build a 3‐D Vs model of the upper crust beneath southeastern Australia by jointly inverting Rayleigh wave phase velocity, ellipticity, and teleseismic P‐coda waves. At depths less than 3 km, sedimentary basins are revealed as prominent low‐velocity anomalies. The average sediment depth in the Murray Basin is ∼713 m with the maximum depth reaching 1,500 m, while, in the Sydney‐Gunnedah‐Bowen Basin (SGBB), the average sediment depth is 600 m with a maximum depth of ∼3,000 m. Vs of the sedimentary layer in the Murray Basin is lower than that in the SGBB, suggesting that the sediments in the Murray Basin are less consolidated compared with those in the SGBB. The unconsolidated sediments make the Murray Basin more vulnerable to seismic hazards. In the upper and middle crystalline crust, we image two low‐velocity zones: one beneath the Adelaide Fold Belt and the other beneath the Begargo Hill volcanic zone. The low‐velocity zone beneath the Adelaide Fold Belt is likely caused by high temperature and/or a weak granitic composition, which lead to the concentration of intraplate seismic activities in this area. The low‐velocity zone beneath the Begargo Hill Volcanic zone indicates that there exists magma in the upper to middle crust or the temperature there is still very high after the latest eruption.

Identifying Periods of Forecast Model Confidence for Improved Subseasonal Prediction of Precipitation

AbstractSubseasonal forecast skill is not homogeneous in time, and prior assessment of the likely forecast skill would be valuable for end-users. We propose a method for identifying periods of high forecast confidence using atmospheric circulation patterns, with an application to southern Australia precipitation. In particular, we use archetypal analysis to derive six patterns, called archetypes, of daily 500-hPa geopotential height (Z500) fields over Australia. We assign Z500 reanalysis fields to the closest-matching archetype and subsequently link the archetypes to precipitation for three key regions in the Australian agriculture and energy sectors: the Murray Basin, southwest Western Australia, and western Tasmania. Using a 20-yr hindcast dataset from the European Centre for Medium-Range Weather Forecasts subseasonal-to-seasonal prediction system, we identify periods of high confidence as when hindcast Z500 fields closely match an archetype according to a distance criterion. We compare the precipitation hindcast accuracy during these confident periods compared to normal. Considering all archetypes, we show that there is greater skill during confident periods for lead times of less than 10 days in the Murray Basin and western Tasmania, and for greater than 6 days in southwest Western Australia, although these conclusions are subject to substantial uncertainty. By breaking down the skill results for each archetype individually, we highlight how skill tends to be greater than normal for those archetypes associated with drier-than-average conditions.

Divergent lineages in a semi-arid mallee species, Eucalyptus behriana, correspond to a major geographic break in southeastern Australia.

AimTo infer relationships between populations of the semi‐arid, mallee eucalypt, Eucalyptus behriana, to build hypotheses regarding evolution of major disjunctions in the species' distribution and to expand understanding of the biogeographical history of southeastern Australia.LocationSoutheastern Australia.Taxon Eucalyptus behriana (Myrtaceae, Angiospermae).MethodsWe developed a large dataset of anonymous genomic loci for 97 samples from 11 populations of E. behriana using double digest restriction site‐associated DNA sequencing (ddRAD‐seq), to determine genetic relationships between the populations. These relationships, along with species distribution models, were used to construct hypotheses regarding environmental processes that have driven fragmentation of the species’ distribution.ResultsGreatest genetic divergence was between populations on either side of the Lower Murray Basin. Populations west of the Basin showed greater genetic divergence between one another than the eastern populations. The most genetically distinct population in the east (Long Forest) was separated from others by the Great Dividing Range. A close relationship was found between the outlying northernmost population (near West Wyalong) and those in the Victorian Goldfields despite a large disjunction between them.ConclusionsPatterns of genetic variation are consistent with a history of vicariant differentiation of disjunct populations. We infer that an early disjunction to develop in the species distribution was that across the Lower Murray Basin, an important biogeographical barrier separating many dry sclerophyll plant taxa in southeastern Australia. Additionally, our results suggest that the western populations fragmented earlier than the eastern ones. Fragmentation, both west and east of the Murray Basin, is likely tied to climatic changes associated with glacial‐interglacial cycles although it remains possible that major geological events including uplift of the Mount Lofty Ranges and basalt flows in the Newer Volcanics Province also played a role.

Size-dependent Behaviour of Weak Intact Rocks

Characterising the size-dependent behaviour of rock has been a significant challenge in rock engineering particularly during the design of structures on or within a rock mass. Generally, the mechanical characterisation of rock starts at the laboratory scale where the intact rock is tested and then the resulting data is extrapolated to the field conditions using a suitable size effect model. Despite extensive research on the size effect of intact rock, very few have included the weak rocks with low uniaxial compressive strength (UCS). Thus, in this study, the size-dependent behaviour of two different weak rocks namely, Gambier limestone and artificial rock with uniaxial compressive strength of less than 10 MPa were investigated experimentally and analytically. The diameter of the cylindrical Gambier limestone samples varied from 26 to 285 mm while the diameter of artificial rock samples ranged between 26 and 139 mm. For Gambier limestone, the uniaxial compressive, Brazilian and point load experiments were carried out while for artificial rock, only the uniaxial compressive tests were conducted. In both rock types, the ascending and then descending size effect trend was a pronounced behaviour for UCS and Young’s modulus data while the size effect behaviour of Poisson’s ratio was inconclusive. Also, the tensile strength and point load index data obtained from Gambier limestone revealed only descending size effect trends. The unified size effect law and its improved version were fitted to the UCS and Young’s modulus data leading to a very good agreement between the data and the model predictions. It was confirmed that an improved version of unified size effect law can predict a more realistic strength value for a sample with an infinite size. Finally, the applicability of two descending size effect models to the tensile strength and point load data was assessed and concluded that the multifractal scaling law is a suitable model for the point load data while the size effect law can better predict the tensile strength data.

Modelling Groundwater Returns to Streams From Irrigation Areas with Perched Water Tables

Quantifying the magnitude and timing of groundwater returns to streams from irrigation is important for the management of natural resources in irrigation districts where the quantity or quality of surface water can be affected. Deep vadose zones and perched water tables can complicate the modelling of these fluxes, and model outputs may be biased if these factors are misrepresented or ignored. This study was undertaken in the Murray Basin in southern Australia to develop and test an integrated modelling method that links irrigation activity to surface water impacts by accounting for all key hydrological processes, including perching and vadose zone transmission. The method incorporates an agronomic water balance to simulate root zone processes, semi-analytical transfer functions to simulate the deeper vadose zone, and an existing numerical groundwater model to simulate irrigation returns to the Murray River and inform the management of river salinity. The integrated modelling can be calibrated by various means, depending on context, and has been shown to be beneficial for management purposes without introducing an unnecessary level of complexity to traditional modelling workflows. Its applicability to other irrigation settings is discussed.

Modelling Recharge from Irrigation Developments with a Perched Water Table and Deep Unsaturated Zone

Modelling of recharge under irrigation zones for input to groundwater modelling is important for assessment and management of environmental risks. Deep vadose zones, when coupled with perched water tables, affect the timing and magnitude of recharge. Despite the temporal and spatial complexities of irrigation areas; recharge in response to new developments can be modelled semi-analytically, with most outputs comparing well with numerical models. For parameter ranges relevant to the western Murray Basin in southern Australia, perching can reduce the magnitude of recharge relative to irrigation accessions and will cause significant time lags for changes to move through vadose zone. Recharge in the vicinity of existing developments was found to be similar to that far from existing developments. This allows superposition to be implemented spatially for new developments, thus simplifying estimation of recharge. Simplification is further aided by the use of exponential approximants for recharge responses from individual developments.

Bluecap: A geospatial model to assess regional economic-viability for mineral resource development

Frontier mineral exploration is often exclusively focused on assessing geological potential without consideration for the economic viability of resource development. This strategy may overlook potentially prosperous zones for more geologically-favoured but financially-disadvantageous regions, or conversely, may introduce implicit biases against potential developments without due regard to underlying economies of scale or proximity to infrastructure. Accordingly, in this paper, we introduce a numerical model aimed at identifying economic fairways, i.e. areas permissive to mineral development from an economic perspective. The model, Bluecap, combines large-scale infrastructure and geological datasets to conduct geospatial analysis of the economic-viability of mining operations across Australia.We provide a detailed description of the inputs and assumptions that underlie the cost models employed in Bluecap, outlining the methods used to evaluate mining, processing, administrative and infrastructure expenses. We also describe the databases used by the model to evaluate available infrastructure, transportation distances and depth of cover. Finally, we present examples that demonstrate the use of the Bluecap model on regions around Mount Isa and the Murray Basin to verify its ability to evaluate commercially feasible mineral prospects. While the immediate utility of this model stands to benefit mineral explorers, its ability to map mineral economic fairways also provides an objective, evidence base to underpin government decision making with respect to position of new infrastructure and consideration of competing land use claims.

Mineral distribution and provenance of heavy mineral sands (zircon, ilmenite, rutile) deposits from the NW Murray Basin, far western NSW, Australia

There is significant economic interest in the Murray Basin of southeastern Australia as it is proving to be a major heavy mineral sands (HMS) province that will be one of Australia’s major source for production of rutile, zircon and ilmenite. The distribution and provenance of HMS resources in the Murray Basin is poorly understood because of its huge size, limited exploration and the complex depositional, structural and weathering mechanisms in their development. In this paper, we focus on the Copi North and Magic deposits some 130–180 km south of Broken Hill, NSW. The heavy mineral assemblages of the Copi North and Magic deposits are very similar, with the main economic minerals being ilmenite, leucoxene, pseudorutile, rutile and zircon. Intensive fracturing and brecciation are identified in many samples and are inferred to have been initially caused by multi-stage deformational processes associated with metamorphism and then further developed through alluvial and eolian transportation. Both deposits are classified as ‘medium sands, symmetrical, mesokurtic and moderately well-sorted’. The majority of economic minerals are of low to medium sphericity and subrounded, along with abundant polished eolian quartz grains. The Copi North deposit has coarser and more poorly sorted sediments with higher HMS grades and magnetics content than the Magic deposit, reflecting a higher energy depositional environment. The main source for the HMS for the Copi North and Magic deposits is largely ascribed to the Broken Hill Block. Previous studies have shown that the Broken Hill orebody underwent substantial sub-aerial weathering over hundreds of millions of years. In addition, the complex metamorphic events experienced by the Broken Hill Block were capable of forming the broad series of minerals identified within the Copi North and Magic deposits. The HMS were believed to have been transported through paleovalleys near the Mulculca Fault in a southeast direction representing a feeder system into the NW Murray Basin. Both deposits feature a relatively linear geometry (roughly parallel to the strike of the paleoshoreline), with high HMS grades, modest tonnages, and coarser sediments when compared to WIM-style offshore deposits. Compared to other strandline HMS deposits of the Murray Basin, they are smaller in size although have similar high grades of 3.7–6.9% THM and similar proportions of the HMS assemblage of ilmenite, leucoxene, rutile and zircon. Deposits of similar size and grade are likely to occur throughout the northern part of the Murray Basin.KEY POINTSBoth the Copi North and Magic deposits contain similar mineral assemblages with the provenance of the heavy minerals ascribed to the Broken Hill Block.A relatively high energy inshore environment is inferred for the Copi North deposit while a lower energy environment associated with either a foreshore or backshore environment is inferred for the Magic deposit.Deposits of similar mineralogy, grades and size are likely to occur elsewhere throughout the northern Murray Basin.