Palaeolimnological records of shallow lake biodiversity change: exploring the merits of single versus multi-proxy approaches

Shallow lakes are among the most threatened ecosystems in the world and many contemporary studies have demonstrated declines in biodiversity due to anthropogenic forcing. Mostly, however, these studies have not covered the full period of human-induced diversity change in lakes which is typically over decades-centuries. Here we provide two examples of palaeoecological studies focussed on reconstructing biodiversity changes in contrasting shallow lake environments that demonstrate the efficacy of the approach—a shallow UK lake and a suite of floodplain lakes (called billabongs) in the Murray-Darling basin, Australia. In the Murray-Darling billabongs, complex sedimentary processes operate, sediment chronologies are less certain and replication of sites is needed to confirm patterns. The combination of sediment records from 10 billabongs showed that diatom diversity changes pre- and post-European (>1850) disturbance were inconsistent; however, reductions in diversity were more common and appear to reflect reductions in macrophyte abundance. At Felbrigg Lake, a multi-proxy study with strong chronological control demonstrated divergent responses of macrophyte, diatom, cladoceran and chironomid richness and diversity to a century of eutrophication. Eutrophication of the site was qualitatively inferred from changes in the macrophyte community in turn reconstructed from plant macrofossils. Benthic cladocerans showed a consistent decline in richness through the record, reflecting the gradual reduction in their macrophyte associated habitat over the past century. Diatom richness and diversity responses were complex, with increases in diversity and richness linked to both increases and decreases in macrophyte species richness and abundance. Chironomid richness and diversity patterns were less consistently linked to eutrophication. The loss of the dominant zooplanktivore (perch) in the 1970s was reflected in the richness and diversity profiles for all groups. Our study reveals clear potential for using sediment cores to infer biodiversity change in shallow lakes and shallow lake regions. However, given the contrasting patterns of diversity change for the different biological groups both in Felbrigg Lake and between Felbrigg and the billabongs, caution is required when interpreting whole-ecosystem biodiversity changes (or the absence of change) based on single as opposed to multi-proxy studies.     

Internal structure of Europa and Callisto

Models of the internal structure of completely differentiated Europa and partially differentiated Callisto have been constructed on the basis of Galileo gravity measurements, geochemical constraints on composition of ordinary and carbonaceous chondrites, and thermodynamic data on the equations of state of water, high-pressure ices, and meteoritic material. We assume thermal and mechanical equilibrium for the interiors of the satellites. A geophysically and geochemically permissible thickness of Europa's outer water-ice shell lies between 105 and 160 km (6.2-9.2% of total mass). Our results show that the bulk composition of the rock-iron core of Europa may be described by material approaching the L/LL-type chondrites in composition, but cannot be correlated either with the material of CI chondrites or H chondrites. For Europa's L/LL-chondritic models, core radii are estimated to be 470-640 km (5.3-12.5% of total mass). The allowed thickness of Europa's H₂O layer ranges from 115±10 km for a differentiated L/LL-type chondritic mantle with a crust to 135±10 km for an undifferentiated mantle. We show that Callisto must only be partially differentiated into an outer ice-I layer, a water ocean, a rock-ice mantle, and a rock-iron core (mixture of anhydrous silicates and/or hydrous silicates + FeFeS alloy). We accept that the composition of the rock-iron material of Callisto is similar to the bulk composition of L/LL-type chondritic material containing up to 10-15% of iron and iron sulfide. Assuming conductive heat transfer through the ice-I crust [Ruiz, 2001. The stability against freezing of an internal liquid-water ocean on Gallisto. Nature, 412, 409-411], heat flows were estimated and the possibility of the existence of a water ocean in Callisto was evaluated. The liquid phase is stable (not freezing) beneath the ice crust, if the heat flow is between 3.3 and 3.7 mW m⁻², which corresponds to the heat flow from radiogenic sources. The thickness of the ice-I crust is 135-150 km, and that of the underlying water layer, 120-180 km. The results of modeling support the hypothesis that Callisto may have an internal liquid-water ocean. The allowed total (maximum) thickness of the outer water-ice shell is up to 270-315 km. Rock-iron core radii, depending on the presence or absence of hydrous silicates, do not exceed 500-700 km, the thickness of an intermediate ice-rock mantle is not less than 1400 km, and its density is in the range of 1960-2500 kg m⁻³. The surface temperature of Callisto is expected to be 100-112 K. The total amount of H₂O in Callisto is found to be 49-55 wt%. The correspondence between the density and moment of inertia values for bulk ice-free Io, rock-iron core of ice-poor Europa, and rock-iron cores of Ganymede and Callisto shows that their bulk compositions may be, in general, similar and may be described by the composition close to a material of the L/LL-type chondrites with the (Fe_tot/Si) weight ratios ranging from 0.9 to 1.3. Planetesimals composed of these types of ordinary chondrites could be considered as analogues of building material for the rock-iron cores of the Galilean satellites. Similarity of bulk composition of the rock-iron cores of the inner and outer satellites implies the absence of iron-silicon fractionation in the protojovian nebula.     

Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model

Current theories of F-layer storms are discussed using numerical simulations with the Upper Atmosphere Model, a global self-consistent, time dependent numerical model of the thermosphere-ionosphere-plasmasphere-magnetosphere system including electrodynamical coupling effects. A case study of a moderate geomagnetic storm at low solar activity during the northern winter solstice exemplifies the complex storm phenomena. The study focuses on positive ionospheric storm effects in relation to thermospheric disturbances in general and thermospheric composition changes in particular. It investigates the dynamical effects of both neutral meridional winds and electric fields caused by the disturbance dynamo effect. The penetration of short-time electric fields of magnetospheric origin during storm intensification phases is shown for the first time in this model study. Comparisons of the calculated thermospheric composition changes with satellite observations of AE-C and ESRO-4 during storm time show a good agreement. The empirical MSISE90 model, however, is less consistent with the simulations. It does not show the equatorward propagation of the disturbances and predicts that they have a gentler latitudinal gradient. Both theoretical and experimental data reveal that although the ratio of [O]/[N₂] at high latitudes decreases significantly during the magnetic storm compared with the quiet time level, at mid to low latitudes it does not increase (at fixed altitudes) above the quiet reference level. Meanwhile, the ionospheric storm is positive there. We conclude that the positive phase of the ionospheric storm is mainly due to uplifting of ionospheric F2-region plasma at mid latitudes and its equatorward movement at low latitudes along geomagnetic field lines caused by large-scale neutral wind circulation and the passage of travelling atmospheric disturbances (TADs). The calculated zonal electric field disturbances also help to create the positive ionospheric disturbances both at middle and low latitudes. Minor contributions arise from the general density enhancement of all constituents during geomagnetic storms, which favours ion production processes above ion losses at fixed height under day-light conditions.     

Comparisons between the reproductive biology of black bream Acanthopagrus butcheri (Teleostei: Sparidae) in four estuaries with widely differing characteristics

The reproductive biology of Acanthopagrus butcheri has been studied in the permanently open Swan River and intermittently open Moore River estuaries on the lower west coast of Australia (31--32 °S) and in the permanently open Nornalup Walpole and normally closed Wellstead estuaries on the southern coast of Western Australia (34--35 °S). Trends exhibited by gonadosomatic indices, gonadal maturity stages and the sizes and developmental stages of the oocytes demonstrate that A. butcher typically spawns in spring and early summer. However, spawning occurred in salinities ranging from as low as 3.5-8 g L^–1 in the Moore River Estuary to as high as 41-45 g L^–1 in the Wellstead Estuary. Furthermore, water temperatures during spawning were greater in the two northern estuaries (19.7--28.5 °C) than in the two southern and cooler estuaries (17.5--23.4 °C). Histological studies strongly indicate that A. butcheri spawn more than once in a breeding season and demonstrate that the development of its oocytes exhibits group synchrony sensu de Vlaming (1983). The ages and total lengths at which, on average, female and male A. butcheri both first attain maturity in the Swan River Estuary were ca 2 years and ca 215 mm. However, the age at which individual fish in that system reach maturity was influenced by body size. This suggests that the attainment of first maturity at an older age but smaller length in the Moore River and Nornalup Walpole estuaries than is the case in the Swan River Estuary is a consequence of the slower growth rates of A. butcheri in those estuaries. The combination of the young age (ca 2 years) but small length (ca 145 mm) at which maturity is first attained in the Wellstead Estuary could have resulted from selection pressures brought about by high mortality rates and/or heavy fishing pressure in this estuary. The mean fecundity of A. butcheri, based on the combined number of yolk vesicle and yolk granule oocytes found in ovaries just prior to the onset of spawning, was 1580 × 10³. The significance of the sizes at first maturity, minimum legal length for capture, mesh selectivity data and closure of certain regions of estuaries to fishing for the management of the recreational and commercial fishery for A. butcheri is discussed.     

Tracking Solar Active Region Outflow Plasma from Its Source to the Near-Earth Environment

Seeking to establish whether active-region upflow material contributes to the slow solar wind, we examine in detail the plasma upflows from Active Region (AR) 10978, which crossed the Sun’s disc in the interval 8 to 16 December 2007 during Carrington rotation (CR) 2064. In previous work, using data from the Hinode/EUV Imaging Spectrometer, upflow velocity evolution was extensively studied as the region crossed the disc, while a linear force-free-field magnetic extrapolation was used to confirm aspects of the velocity evolution and to establish the presence of quasi-separatrix layers at the upflow source areas. The plasma properties, temperature, density, and first ionisation potential bias [FIP-bias] were measured with the spectrometer during the disc passage of the active region. Global potential-field source-surface (PFSS) models showed that AR 10978 was completely covered by the closed field of a helmet streamer that is part of the streamer belt. Therefore it is not clear how any of the upflowing AR-associated plasma could reach the source surface at 2.5 R_⊙ and contribute to the slow solar wind. However, a detailed examination of solar-wind in-situ data obtained by the Advanced Composition Explorer (ACE) spacecraft at the L₁ point shows that increases in O⁷⁺/O⁶⁺, C⁶⁺/C⁵⁺, and Fe/O – a FIP-bias proxy – are present before the heliospheric current-sheet crossing. These increases, along with an accompanying reduction in proton velocity and an increase in density are characteristic of both AR and slow-solar-wind plasma. Finally, we describe a two-step reconnection process by which some of the upflowing plasma from the AR might reach the heliosphere.     

Astrometry from space: GAIA and planet detection

The proposed baseline GAIA mission will be able to detect the astrometric signature of Jupiter-size planets around of the order of a million stars, using either global or narrow-angle astrometry. If the mission can realize the higher astrometric accuracy that photon statistics allows for bright stars, lower-mass planets (from Earth size to ten times larger) can be found around ten to a few hundred stars.     

Aerial electromagnetic sounding of the lithosphere of Venus

Electromagnetic (EM) investigation depths are larger on Venus than Earth due to the dearth of water in rocks, in spite of higher temperatures. Whistlers detected by Venus Express proved that lightning is present, so the Schumann resonances ～10–40 Hz may provide a global source of electromagnetic energy that penetrates ～10–100 km. Electrical conductivity will be sensitive at these depths to temperature structure and hence thermal lithospheric thickness. Using 1D analytic and 2D numerical models, we demonstrate that the Schumann resonances—transverse EM waves in the ground-ionosphere waveguide—remain sensitive at all altitudes to the properties of the boundaries. This is in marked contrast to other EM methods in which sensitivity to the ground falls off sharply with altitude. We develop a 1D analytical model for aerial EM sounding that treats the electrical properties of the subsurface (thermal gradient, water content, and presence of conductive crust) and ionosphere, and the effects of both random errors and biases that can influence the measurements. We initially consider specified 1D lithospheric thicknesses 100–500 km, but we turn to 2D convection models with Newtonian temperature-dependent viscosity to provide representative vertical and lateral temperature variations. We invert for the conductivity-depth structure and then temperature gradient. For a dry Venus, we find that the error on temperature gradient obtained from any single local measurement is ～100%—perhaps enough to distinguish “thick” vs. “thin” lithospheres. When averaging over thousands of kilometers, however, the standard deviation of the recovered thermal gradient is within the natural variability of the convection models, <25%. A “wet” interior (hundreds of ppm H₂O) limits EM sounding depths using the Schumann resonances to <20 km, and errors are too large to estimate lithospheric properties. A 30-km conductive crust has little influence on the dry-interior models because the Schumann penetration depths are significantly larger. We conclude that EM sounding of the interior of Venus is feasible from a 55-km high balloon. Lithospheric thickness can be measured if the upper-mantle water content is low. If H₂O at hundreds of ppm is present, the deeper, temperature-sensitive structure is screened, but the “wet” nature of the upper mantle, as well as structure of the upper crust, is revealed.     

Multipoint observations of low latitude ULF Pc3 waves in south-east Australia

Geomagnetic pulsations recorded on the ground are the signatures of the integrated signals from the magnetosphere. Pc3 geomagnetic pulsations are quasi-sinusoidal variations in the earth’s magnetic field in the period range 10–45 seconds. The magnitude of these pulsations ranges from fraction of a nT (nano Tesla) to several nT. These pulsations can be observed in a number of ways. However, the application of ground-based magnetometer arrays has proven to be one of the most successful methods of studying the spatial structure of hydromagnetic waves in the earth’s magnetosphere. The solar wind provides the energy for the earth’s magnetospheric processes. Pc3–5 geomagnetic pulsations can be generated either externally or internally with respect to the magnetosphere. The Pc3 studies undertaken in the past have been confined to middle and high latitudes. The spatial and temporal variations observed in Pc3 occurrence are of vital importance because they provide evidence which can be directly related to wave generation mechanisms both inside and external to the magnetosphere. At low latitudes (L < 3) wave energy predominates in the Pc3 band and the spatial characteristics of these pulsations have received little attention in the past. An array of four low latitude induction coil magnetometers were established in south-east Australia over a longitudinal range of 17 degrees at L = 1.8 to 2.7 for carrying out the study of the effect of the solar wind velocity on these pulsations. Digital dynamic spectra showing Pc3 pulsation activity over a period of about six months have been used to evaluate Pc3 pulsation occurrence. Pc3 occurrence probability at low latitudes has been found to be dominant for the solar wind velocity in the range 400–700 km/s. The results suggest that solar wind controls Pc3 occurrence through a mechanism in which Pc3 wave energy is convected through the magnetosheath and coupled to the standing oscillations of magnetospheric field lines.     

Crystal field stabilization energies of almandine-pyrope and almandine-spessartine garnets determined by FTIR near infrared measurements

The band positions of three partially overlapping Fe²⁺ spin-allowed transitions located between 4000 and 9000 cm^–1 in almandine-pyrope and almandinespessartine garnets solid solutions were measured using near-infrared (NIR) spectroscopy. The crystal field stabilization energies (CFSE) along both binaries were calculated assuming a splitting of 1100 cm^–1 for the lower orbitals. The CFSE show a slight increase along the almandine-spessartine binary from 3730 to 3810 cm^–1 and a larger increase from 3730 to 3970 cm^–1 for the almandine-pyrope binary. Dodecahedral Fe²⁺-site distortion increases with an increase in spessartine component and decreases with increasing pyrope component, in agreement with average dodecahedral site distortions determined from diffraction experiments. The excess CFSE's along both joins are negative. For the almandinespessartine binary they are small, but are about 3.5 times larger in magnitude along the join almandine-pyrope, where an interaction parameter of W= -2.9 KJ/mole has been derived from a symmetric mixing model. The excess CFSE are relatively small compared to the magnitudes of the excess enthalpies of mixing that have been assigned to garnet solid solutions. Moreover, they give no indication which could support the positive and asymmetric excess enthalpies of mixing that have been proposed for almandine-pyrope solid solutions.