Large-scale pesticide monitoring across Great Barrier Reef catchments--Paddock to Reef Integrated Monitoring, Modelling and Reporting Program

The transport and potential toxicity of pesticides in Queensland (QLD) catchments from agricultural areas is a key concern for the Great Barrier Reef (GBR). In 2009, a pesticide monitoring program was established as part of the Australian and QLD Governments' Reef Plan (2009). Samples were collected at eight End of System sites (above the tidal zone) and three sub-catchment sites. At least two pesticides were detected at every site including insecticides, fungicides, herbicides, and the Reef Plan's (2009) five priority photosystem II (PSII) herbicides (diuron, atrazine, hexazinone, tebuthiuron and ametryn). Diuron, atrazine and metolachlor exceeded Australian and New Zealand water quality guideline trigger values (TVs) at eight sites. Accounting for PSII herbicide mixtures increased the estimated toxicity and led to larger exceedances of the TVs at more sites. This study demonstrates the widespread contamination of pesticides, particularly PSII herbicides, across the GBR catchment area which discharges to the GBR.     

Biogeomorphology in the marine landscape: Modelling the feedbacks between patches of the polychaete worm Lanice conchilega and tidal sand waves

Tidal sand waves are dynamic bedforms found in coastal shelf seas. Moreover, these areas are inhabited by numerous benthic species, of which the spatial distribution is linked to the morphological structure of sand waves. In particular, the tube-building worm Lanice conchilegais of interest as this organism forms small mounds on the seabed, which provide shelter to other organisms. We investigate how the interactions between small-scale mounds (height ∼dm) and large-scale sand waves (height ∼m) shape the bed of the marine environment. To this end, we present a two-way coupled process-based model of sand waves and tube-building worm patches in Delft3D. The population density evolves according to a general law of logistic growth, with the bed shear stress controlling the carrying capacity. Worm patches are randomly seeded and the tubes are mimicked by small cylinders that influence flow and turbulence, thereby altering sediment dynamics. Model results relate the patches with the highest worm densities to the sand wave troughs, which qualitatively agrees with field observations. Furthermore, the L. conchilegatubes trigger the formation of sandy mounds on the seabed. Because of the population density distribution, the mounds in the troughs can be several centimetres higher than on the crests. Regarding sand wave morphology, the combination of patches and mounds are found to shorten the time-to-equilibrium. Also, if the initial bed comprised small sinusoidal sand waves, the equilibrium wave height decreased with a few decimetres compared to the situation without worm patches. As the timescale of mound formation (years) is shorter than that of sand wave evolution (decades), the mounds induce (and accelerate) sand wave growth on a similar spatial scale to the mounds. Initially, this leads to shorter sand waves than they would be in an abiotic environment. However, near equilibrium the wavelengths tend towards their abiotic counterparts again. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Modelling sand wave migration in shallow shelf seas

Sand waves form a prominent regular pattern in the offshore seabed of sandy shallow seas. The positions of sand-wave crests and troughs slowly change in time. Sand waves are usually assumed to migrate in the direction of the residual current. This paper considers the physical mechanisms that may cause sand waves to migrate and methods to quantify the associated migration rates. We carried out a theoretical study based on the assumption that sand waves evolve as free instabilities of the system. A linear stability analysis was then performed on a 2DV morphological model describing the interaction between the vertically varying water motion and an erodible bed in a shallow sea. Here, we disrupted the basic tidal symmetry by choosing a combination of a steady current (M₀) and a sinusoidal tidal motion (M₂) as the basic flow. We allowed for two different physical mechanisms to generate the steady current: a sea surface wind stress and a pressure gradient. The results show that similar sand waves develop for both flow conditions and that these sand waves migrate slowly in the direction of the residual flow. The rates of migration and wavelengths found in this work agree with theoretical and empirical values reported in the literature.     

Formation of offshore tidal sand banks triggered by a gasmined bed subsidence

Offshore gasmining is an example of a human intervention with a morphological impact. On land, it is usually attended with a dish-like bed depression. We show that, if located at sea, such a bed depression can become morphodynamically active by triggering mechanisms related to tidal sand bank formation. To that end, a simple morphological model is considered which describes an erodible bed subject to a tidal wave in a shallow sea. The continuous subsidence is modelled by a sink term in the sediment balance. Then, a linear approximation is carried out to describe the bed evolution after the onset of subsidence. The results, presented in physical space, show that the subsidence triggers the formation of a sand bank pattern that gradually spreads around the centre of subsidence, at a rate that may go up to 160 m year⁻¹, depending on the tidal transport rate and the tidal eccentricity. The dimension of the depression does not affect the spreading rate nor the orientation of the sand banks, but it does influence their spacing. The main conclusion is that the horizontal extent of the area influenced by the bed depression by far exceeds that of the direct subsidence, thus showing that bed depressions on land and at sea indeed behave in fundamentally different ways. The results suggest that nonlinear effects are worthwhile to be investigated in order to describe finite amplitude development of sand banks as well as the interaction between subsidence and bed forms.     

Some genetic implications of silica and aluminosilicates in peat and coal

The inorganic content of Okefenokee peat, in the absence of a significant detrital contribution, is a function of sponge spicule, phytolith, and diatom abundance in addition to the degree of organic matter decomposition and the frequency of fires. The detrital component (quartz) dominates ash composition where it is present.No clear correlation has been observed between the inorganic content of coals examined in this study and their maceral composition. Vitrinite-rich coals have a wide range of ash yield. These observations suggest that the minerals have been introduced detritally rather than accumulated as a result of peat degradation. Quartz in the coals appears to be detrital, and the presence of 2m₁ mica indicates that there is also a detrital contribution to the clay minerals present.     

The timing of regional Lateglacial events and post-glacial sedimentation rates from Lake Superior

We analyze both new and previously published paleomagnetic records of secular variation (PSV) from Lake Superior sediment cores and compare these records to correlated rhythmite (varve) thickness records to determine post-glacial sedimentation rates and to reassess the termination of glaciolacustrine varves in the basin. The results suggest that offshore sedimentation rates have exhibited considerable spatial variation over the past 8000 years, particularly during the mid-Holocene. We attribute offshore, mid-Holocene sedimentation changes to alterations in whole basin circulation, perhaps precipitated by a greater dominance of the Gulf of Mexico air mass during the summer season. Nearshore bays are characterized by high sedimentation rates for at least 1000 years after varve cessation and during a period between around 4500 and 2000 cal. BP. After 2000 cal. BP, sedimentation rates subsided to earlier rates. The increases between 4500 and 2000 cal. BP are probably due to lake level fall after the Nipissing II highstand.The older glaciolacustrine varve thickness records suggest that the influx of glacially derived sediment ended abruptly everywhere in the lake, except near the Lake Nipigon inlets. Multiple sediment cores reveal 36 anomalously thick varves, previously ascribed to the formation of the Nakina moraine, which were deposited just prior to varve cessation in the open lake. The PSV records support the observation that the cessation of these thick varves is a temporally correlative event, occurring at 9035±170 cal. BP (calibrated years before 1950, ca 7950–8250 ¹⁴C BP). This date would correlate to the eastern diversion of Lake Agassiz and glacial meltwater into Lake Ojibway.     

The Siberian lithosphere traverse: mantle terranes and the assembly of the Siberian Craton

The kimberlite fields scattered across the NE part of the Siberian Craton have been used to map the subcontinental lithospheric mantle (SCLM), as it existed during Devonian to Late Jurassic time, along a 1000-km traverse NE–SW across the Archean Magan and Anabar provinces and into the Proterozoic Olenek Province. 4100 garnets and 260 chromites from 65 kimberlites have been analysed by electron probe (major elements) and proton microprobe (trace elements). These data, and radiometric ages on the kimberlites, have been used to estimate the position of the local (paleo)geotherm and the thickness of the lithosphere, and to map the detailed distribution of specific rock types and mantle processes in space and time. A low geotherm, corresponding approximately to the 35 mW/m² conductive model of Pollack and Chapman [Tectonophysics 38, 279–296, 1977], characterised the Devonian lithosphere beneath the Magan and Anabar crustal provinces. The Devonian geotherm beneath the northern part of the area was higher, rising to near a 40 mW/m² conductive model. Areas intruded by Mesozoic kimberlites are generally characterised by this higher, but still ‘cratonic' geotherm. Lithosphere thickness at the time of kimberlite intrusion varied from ca. 190 to ca. 240 km beneath the Archean Magan and Anabar provinces, but was less (150–180 km) beneath the Proterozoic Olenek Province already in Devonian time. Thinner Devonian lithosphere (140 km) in parts of this area may be related to Riphean rifting. Near the northern end of the traverse, differences in geotherm, lithosphere thickness and composition between the Devonian Toluopka area and the nearby Mesozoic kimberlite fields suggest thinning of the lithosphere by ca. 50–60 km, related to Devonian rifting and Triassic magmatism. A major conclusion of this study is that the crustal terrane boundaries defined by geological mapping and geophysical data (extended from outcrops in the Anabar Shield) represent major lithospheric sutures, which continue through the upper mantle and juxtapose lithospheric domains that differ significantly in composition and rock-type distribution between 100 and 250 km depth. The presence of significant proportions of harzburgitic and depleted lherzolitic garnets beneath the Magan and Anabar provinces is concordant with their Archean surface geology. The lack of harzburgitic garnets, and the chemistry of the lherzolitic garnets, beneath most of the other fields are consistent with the Proterozoic surface rocks. Mantle sections for different terranes within the Archean portion of the craton show pronounced differences in bulk composition, rock-type distribution, metasomatic overprint and lithospheric thickness. These observations suggest that individual crustal terranes, of both Archean and Proterozoic age, had developed their own lithospheric roots, and that these differences were preserved during the Proterozoic assembly of the craton. Data from kimberlite fields near the main Archean–Proterozoic suture (the Billyakh Shear Zone) suggest that reworking and mixing of Archean and Proterozoic mantle was limited to a zone less than 100 km wide.     

Geochemical characteristics of lava-field basalts from eastern Australia and inferred sources: Connections with the subcontinental lithospheric mantle?

 A large new database of major, trace elements and Sr-Nd isotopic ratios from 11 lava-field provinces in New South Wales and Queensland, eastern Australia allows detailed interpretation of the origin of these basaltic magmas. Isotopic signatures and trace element patterns identify an OIB-type (oceanic island basalt) source as a dominant component for most of these and some provinces appear to have additional significant components derived from the subcontinental lithospheric mantle (SCLM). The SCLM components have geochemical characteristics that overlap those observed in spinel lherzolite xenoliths (samples of shallow lithospheric mantle) from eastern Australia. These SCLM components show geochemical provinciality that indicates the occurrence of distinct geochemical lithospheric domains reflecting the timing and style of tectonic evolution of different regions. One component reflects modification by subduction-related processes during the late Paleozoic and Mesozoic, one records enrichment by fluids during old metasomatic events and another suggests a metasomatic event involving a distinctive amphibole and apatite-style enrichment. The composition and age distribution of volcanic lava-field provinces older than 10 Ma are consistent with a model involving a regional upwelling (elongated N–S along eastern Australia) of deep hot mantle related to marginal rifting and with OIB-type source geochemical characteristics. Thermal inhomogeneities within this plume swath resulted in small diapirs which may have undergone melt segregation at about 100 km and incorporated varying amounts of SCLM components there or from higher levels of the SCLM during ascent. Subsequent hot-spot generated central volcanoes overprinted this lava-field volcanism, tapped a similar OIB-type source component and truncated the thermal events. Accepted: 15 March 1995