ALBERT

Description: Transformation of refractory cratonic mantle into more fertile lithologies is the key to the fate of cratonic lithosphere. This process has been extensively studied in the eastern North China Craton (NCC) while that of its western part is still poorly constrained. A comprehensive study of newly-found pyroxenite xenoliths from the Langshan area, in the northwestern part of this craton is integrated with a regional synthesis of pyroxenite and peridotite xenoliths to constrain the petrogenesis of the pyroxenites and provide an overview of the processes involved in the modification of the deep lithosphere. The Langshan pyroxenites are of two types, high-Mg# [Mg2+/(Mg2++Fe2+)*100 = ∼ 90, atomic ratios] olivine-bearing websterites with high equilibration temperatures (880 ∼ 970 oC), and low-Mg# (70 ∼ 80) plagioclase-bearing websterites with low equilibration temperatures (550 ∼ 835 oC). The high-Mg# pyroxenites show trade-off abundances of olivine and orthopyroxene, highly depleted bulk Sr-Nd (ƐNd = +11.41, 87Sr/86Sr = ∼0.7034) and low clinopyroxene Sr isotopic ratios (mean 87Sr/86Sr = ∼0.703). They are considered to reflect the reaction of mantle peridotites with silica-rich silicate melts derived from the convective mantle. Their depletion in fusible components (e.g., FeO, TiO2 and Na2O) and progressive exhaustion of incompatible elements suggest melt extraction after their formation. The low-Mg# pyroxenites display layered structures, convex-upward rare earth element patterns, moderately enriched bulk Sr-Nd isotopic ratios (ƐNd = -14.20 ∼ -16.74, 87Sr/86Sr = 0.7070 ∼ 0.7078) and variable clinopyroxene Sr-isotope ratios (87Sr/86Sr = 0.706-0.711). They are interpreted to be crustal cumulates from hypersthene-normative melts generated by interaction between the asthenosphere and heterogeneous lithospheric mantle. Combined with studies on regional peridotite xenoliths, it is shown that the thinning and refertilization of the lithospheric mantle was accompanied by crustal rejuvenation and that such processes occurred ubiquitously in the northwestern part of the NCC. A geodynamic model is proposed for the evolution of the deep lithosphere, which includes long-term mass transfer through a mantle wedge into the deep crust from the Paleozoic to the Cenozoic, triggered by subduction of the Paleo-Asian ocean and the Late Mesozoic lithospheric extension of eastern Asia.

Description: The orogenic development after the continental collision between Laurussia and Gondwana, led to two contrasting associations of mantle-derived magmatic rocks on the territory of the Bohemian Massif: (i) a 340–310 Ma lamprophyre-lamproite orogenic association and (ii) a 300–275 Ma lamprophyre association of anorogenic affinity. Major types of potassic mantle-derived magmatic rocks recognised in the orogenic and anorogenic associations include: (i) calc-alkaline to alkaline lamprophyres, (ii) alkaline “orthopyroxene minettes” (and geochemically related rocks), and (iii) peralkaline lamproites. These three types significantly differ with respect to mineral, whole-rock and Sr–Nd–Pb–Li isotope composition, and spatial distribution. The calc-alkaline lamprophyres occur throughout the entire Saxo-Thuringian and Moldanubian zones, whereas the different types of malte-derived potassic rocks are spatially restricted to particular zones. Rocks of the Carboniferous lamprophyre-lamproite orogenic association are characterised by variable negative εNd(i) and variably radiogenic Sr(i), whereas the rocks of the Permian lamprophyre association of anorogenic affinity are characterised by positive εNd(i) and relatively young depleted-mantle Nd-model ages reflecting increasing input from upwelling asthenospheric mantle. The small variation in the Pb isotopic composition of post-collisional potassic mantle-derived magmatic rocks (of both the orogenic and anorogenic series) implies that the Pb budget of the mantle beneath the Bohemian Massif is dominated by the same crust-derived material, which itself may include material derived from several sources. The source rocks of “orthopyroxene minettes” are characterised by isotopically light (“eclogitic”) Li and strongly radiogenic (crustal) Sr and may have been metasomatised by high-pressure fluids along the edge of a subduction zone. In contrast, the strongly Al2O3 and CaO depleted mantle source of the lamproites is characterised by isotopically heavy Li and high SiO2 and extreme K2O contents. This mantle source may have been metasomatised predominantly by melts. The mantle source of the lamprophyres may have undergone metasomatism by both fluids and melts.

Description: The Lu–Hf isotope system and Sr–Nd–Hf–Os isotope systematics of mantle rocks are capable of unravelling the early processes in collision belts, especially in a hot subduction context where the Sm–Nd and U–Pb systems in crustal rocks are prone to resetting owing to high temperatures and interaction with melts during exhumation. To improve models of the Devonian–Carboniferous evolution of the Bohemian Massif, we investigated in detail mafic and ultramafic rocks (eclogite, pyroxenite, and peridotite) from the ultrahigh-pressure and ultrahigh-temperature Kutná Hora Crystalline Complex (KHCC: Úhrov, Bečváry, Doubrava, and Spačice localities). Petrography, multiphase solid inclusions, major and trace element compositions of rocks and minerals, and radiogenic isotopic data document contrasting sources and protoliths as well as effects of subduction-related processes for these rocks. The Úhrov peridotite has a depleted composition corresponding to the suboceanic asthenospheric mantle, whereas Bečváry and Doubrava peridotites represent lithospheric mantle that underwent melt refertilization by basaltic and SiO2-undersaturated melts, respectively. Multiphase solid inclusions enclosed in garnet from Úhrov and Bečváry peridotites represent trapped H2O ± CO2-bearing metasomatizing agents and Fe–Ti-rich melts. The KHCC eclogites either formed by high-pressure crystal accumulation from mantle-derived basaltic melts (Úhrov) or represent a fragment of mid-ocean ridge basalt-like gabbroic cumulate (Spačice) and crustal-derived material (Doubrava) both metamorphosed at high P–T conditions. The Lu–Hf age of 395 ± 23 Ma obtained for the Úhrov peridotite reflects garnet growth related to burial of the asthenospheric mantle during subduction of the oceanic slab. By contrast, Spačice and Doubrava eclogites yield younger Lu–Hf ages of ∼350 and 330 Ma, respectively, representing mixed ages as demonstrated by the strong granulite-facies overprint and trace element zoning in garnet grains. We propose a refined model for the Early Variscan evolution of the Bohemian Massif starting with the subduction of the oceanic crust (Saxothuringian ocean) and associated oceanic asthenospheric mantle (Úhrov) beneath the Teplá–Barrandian at ≥380 Ma, which was responsible for melt refertilization of the associated mantle wedge (Bečváry, Doubrava). This was followed by continental subduction (∼370–360 Ma?) accompanied by the oceanic slab break-off and incorporation of the upwelling asthenospheric mantle into the Moldanubian lithospheric mantle and subsequent coeval exhumation of mantle and crustal rocks at ∼350–330 Ma.

Description: The Changning–Menglian orogenic belt (CMOB) in the southeastern Tibetan Plateau is an important link between the Longmu Co–Shuanghu suture (LCSS) in the northern Tibetan Plateau and the Chiang Mai–Inthanon and Bentong–Raub sutures in Thailand and Peninsular Malaysia. These belts and sutures are generally regarded as containing the remnants of the oceanic crust of the Palaeo-Tethys that formed by seafloor spreading as a result of the separation of Gondwana- and Eurasia-derived blocks during the Middle Cambrian. In this paper we report the first discovery of abundant unaltered and retrograde eclogites that occur as irregular lenses and blocks in metasedimentary rocks of the CMOB, and these eclogites form an elongate and almost north–south-trending high-pressure (HP)–ultrahigh-pressure (UHP) metamorphic belt that is ∼200 km long and ∼50 km wide. The newly discovered phengite/talc/epidote–glaucophane eclogites, lawsonite–talc–phengite eclogites, dolomite/magnesite–kyanite eclogites and phengite–kyanite-bearing retrograde eclogites have enriched (E-) and normal mid-ocean ridge basalt (N-MORB)-like affinities and mainly positive as well as some negative whole-rock εNd values (–4·34 to +7·89), which suggest an enriched and depleted oceanic lithosphere source for their protoliths. Magmatic zircons separated from the epidote–glaucophane, magnesite–kyanite and (phengite–kyanite-bearing) retrograde eclogites gave protolith ages of 317–250 Ma, which fit well within the time frame of the opening of the Palaeo-Tethys during the Middle Cambrian and its closure during the Triassic. Abundant metamorphic zircons in the eclogites indicate a Triassic metamorphic event related to the subduction of the Palaeo-Tethys oceanic crust from 235 to 227 Ma. Taking into account previous isotopic age data, we now establish the periods of Early–Middle Triassic (246–227 Ma) and Late Triassic (222–209 Ma) as the ages of subduction and exhumation of the Palaeo-Tethyan oceanic crust, respectively. Thermodynamic modelling revealed that the eclogites record distinct HP–UHP peak metamorphic conditions of 23·0–25·5 kbar and 582–610 °C for the phengite–glaucophane eclogites, 24·0–25·5 kbar and 570–586 °C for the talc–glaucophane eclogites, 29·0–31·0 kbar and 675–712 °C for the dolomite–kyanite eclogites, and 30·0–32·0 kbar and 717–754 °C for the magnesite–kyanite eclogites. These P–T estimates and geochronological data indicate that the Palaeo-Tethys oceanic slab was subducted to different mantle depths from 75 km down to 95 km, forming distinct types of eclogite with a variety of peak eclogite-facies mineral assemblages. The eclogites consistently record clockwise metamorphic P–T–t paths characterized by a heating–compression prograde loop under a low geothermal gradient of 5–10 °C km–1, indicating the rapid subduction of cold oceanic crust at a rate of 4·5–6·0 km Ma–1, followed by isothermal or cooling–decompressive retrogression and exhumation at an average rate of 3·2–4·2 km Ma–1. The newly discovered eclogites of the CMOB with their signatures of ocean-crust subduction are petrologically, geochemically and geochronologically comparable with those of the LCSS, providing powerful support for the idea that a nearly 2000 km long HP–UHP eclogite belt extends from the northern Tibetan Plateau to the southeastern Tibetan Plateau, and that it represents the main boundary suture of the Palaeo-Tethyan domain. These results have far-reaching implications for the tectonic framework and complex metamorphic evolution of the Palaeo-Tethyan domain.

Description: Replacive symplectites (vermicular intergrowths of two or more minerals) are an important feature of layered igneous intrusions, recording evidence of late-stage reactions between interstitial liquid and crystals. They are common throughout the Layered Series of the 564 Ma Sept Iles layered intrusion in Quebec, Canada, and fall into three types: oxy-symplectites, ‘Type I’ symplectites, and ‘Type II’ symplectites. Oxy-symplectites are comprised of magnetite and orthopyroxene, nucleate on olivine primocrysts, and form via the reaction Olivine + O2 → Orthopyroxene + Magnetite; Type I symplectites (of which there are 3 distinct categories) are comprised of anorthitic plagioclase with pyroxene, amphibole, or olivine vermicules, grow from primocryst oxide grains, and replace primary plagioclase; and Type II symplectites (of which there are 2 distinct categories) are comprised of anorthitic plagioclase with orthopyroxene ± amphibole vermicules, grow from primocryst olivine grains, and replace primocryst plagioclase. Rare symplectites composed of biotite and plagioclase are also present. Symplectite growth occurred at 700-1030 °C with pressure constraints of 1-2 kbar. We propose that Type I symplectites, and some Type II symplectites, formed from interaction of primocrysts with residual Fe-rich liquid as a consequence of differential loss of an immiscible Si-rich liquid conjugate from the crystal mush. However, redistribution and concentration of hydrous fluids in incompletely solidified rock, or an increase in water activity of the interstitial melt, may be more plausible processes responsible for the formation of replacive symplectites comprising abundant hydrous mineral assemblages.

Description: Although copepods have been considered tolerant against the direct influence of the ocean acidification (OA) projected for the end of the century, some recent studies have challenged this view. Here, we have examined the direct impact of short-term exposure to a pCO 2 / pH level relevant for the year 2100 ( pH NBS , control: 8.18, low pH : 7.78), on the physiological performance of two representative marine copepods: the calanoid Acartia grani and the cyclopoid Oithona davisae . Adults of both species, from laboratory cultures, were preconditioned for four consecutive days in algal suspensions ( Akashiwo sanguinea ) prepared with filtered sea water pre-adjusted to the targeted pH values via CO 2 bubbling. We measured the feeding and respiratory activity and reproductive output of those pre-conditioned females. The largely unaffected fatty acid composition of the prey offered between OA treatments and controls supports the absence in the study of indirect OA effects (i.e. changes of food nutritional quality). Our results show no direct effect of acidification on the vital rates examined in either copepod species. Our findings are compared with results from previous short- and long-term manipulative experiments on other copepod species.

Description: We analyzed the genetic structure of the radiolarian morphospecies Larcopyle buetschlii from the surface to deep waters (up to 2000 m) in the Japan Sea using the internal transcribed spacer (ITS) regions of ribosomal RNA genes. Each individual had several ITS variants, but these polymorphisms show no vertical phylogeographic structure, suggesting a single biological species. Its rapid clonal reproduction suggested by high ITS variation likely plays a pivotal role in maintaining its wide vertical distribution.

Description: Phronima sedentaria is a hyperiid amphipod that diel migrates into a pronounced oxygen minimum zone (OMZ) in the Eastern Tropical North Pacific. In this study, oxygen consumption and lactate production were measured in P. sedentaria to estimate the aerobic and anaerobic contributions to total metabolism under conditions that mimic its day- (1% oxygen, 10°C) and night-time (20% oxygen, 20°C) habitat. When exposed to hypoxia and low temperature, the total metabolism of P. sedentaria was depressed by 78% compared with normoxic conditions. The metabolic enzymes citrate synthase (CS) and lactate dehydrogenase (LDH) were also measured as indicators of aerobic and anaerobic metabolism, and compared with specimens collected from the California Current and the North Atlantic to assess potential adaptations to low oxygen. LDH activity was not significantly different between regions. Significant differences in CS activity may be due to variation in food availability. Climate change is predicted to increase surface temperatures and cause the expansion of OMZs. This will result in vertical compression of the night-time range for P. sedentaria and is likely to have the same impact on other diel migrators. Habitat compression will reduce zooplankton contribution to carbon cycling and alter oceanic ecology, including predator–prey interactions.

Description: Many phytoplankton exploit phosphorus (P) from organic sources when dissolved inorganic P (DIP) is depleted. This process is, however, rarely considered in ecological and biogeochemical models. We present a mechanistic model describing explicitly the ability of phytoplankton to use dissolved organic P (DOP) when DIP is limiting, by synthesizing alkaline phosphatase (AP) that releases DIP from DOP. This model, applicable to any phytoplankton species expressing AP, is here specifically developed for the colony-forming Phaeocystis globosa. It describes the main processes related to P metabolism, including DIP transport, intracellular accumulation and assimilation. Model behaviour is explored in DIP-limiting batch-type conditions for different DOP ranging between 0 and 1.5 mmol P m –3 . Simulations show that the DOP-derived DIP increases the maximum biomass reached and extends the period of net growth. The magnitude of the enhanced biomass production is controlled by the DOP initially present as well as the released DOP, the latter being recycled by lysis of P. globosa cells. We also present a simplified model version derived from the mechanistic model, which involves fewer state variables and parameters. The latter is directly usable in both variable (quota-type) and fixed stoichiometry descriptions of phytoplankton growth.

Description: Zooplankton diel vertical migration (DVM) is often explained as a balance between predator avoidance and resource acquisition. However, recent studies suggest that ultraviolet radiation (UV) may also be important in driving zooplankton DVM in some systems. Williamson et al. ( Williamson et al ., 2011 ) proposed the "transparency-regulator hypothesis," which integrates UV into our current understanding of the drivers of DVM and predicts that the relative roles of UV and visual predation pressure will vary systematically across a gradient of lake transparency. To assess this hypothesis, we conducted in situ mesocosm experiments in five different lakes: two lakes without fish and three lakes with fish that spanned a range of UV and visible light transparency. We used an open-bottomed mesocosm design that allowed for the direct manipulation of UV that did not constrain visual predators or the amplitude or timing of natural migrations. Consistent with the transparency-regulator hypothesis, we found that UV is an important driver of Daphnia DVM in highly UV transparent lakes with and without fish but not in low transparency systems. Our results also suggest that UV and visual predation pressure may interact in systems of intermediate transparency.

Description: Competition for resources can lead to species exclusion. However, this exclusion may be avoided if species show differential adaptation to physical environment. Empirical studies on competition are difficult when species are phylogenetically close and have complex life cycles. This is the case of B. plicatilis and B. manjavacas , two cryptic rotifer species differing in their salinity niches and in life-history traits related to sex and diapause. These differences have been suggested to promote the stable co-occurrence observed in natural populations of these species. However, in a previous empirical study, the outcome of competition between both species was always exclusion. Here, we theoretically explored the effect of complex life-history traits and salinity fluctuations on the long-term competitive outcome of B. plicatilis and B. manjavacas . We developed a model and simulated ecological scenarios combining different growing period lengths, levels of crossed induction of sex between species and salinity regimes. Results show that a fluctuating salinity regime, an intermediate length of growing season and a low level of crossed induction of sex are essential conditions to take into account to explain coexistence.

Description: Colonization of new habitats through dispersal of phytoplankton cysts might be limited, if resident populations outcompete invaders during germination. We reciprocally transferred Gonyostomum semen (Raphidophyceae) cysts from three lakes into native and foreign waters originating from the respective habitats. Germination rate and germling growth were impacted by water origin, but there was no preference for native water. Gonyostomum semen 's ability to germinate in different conditions might explain its expansion in northern Europe.

Description: Measuring zooplankton biomass and physiological rates is of paramount importance in biological oceanography in order to assess the role of this community in, e.g. carbon fluxes. Classical methods (incubations) are very time-consuming and cannot match the frequency of physical and chemical measurements. Attempting to solve this, a variety of methods (e.g. egg production, RNA/DNA ratio or enzyme activities) have been developed over the last decades. These methods also show uncertainties and hitherto only incubation methods have been widely accepted. Predictive equations relating physiological processes and body weight (bw) and temperature are a rough alternative, normally used to ascertain the role of these organisms in the oceanic ecosystem. However, using imaging systems and empirical relationships to determine bw allows the application of physiological models to each individual, obtaining reliable estimates for taxonomic groups and size classes. In this study, we developed predictive equations suitable for growth and respiration estimations in subtropical regions. In addition, biomass and physiological rates assessed from empirical equations in combination with an image-based system (ZooImage) were compared with standard and enzymatic methods, respectively. We observed a consistent agreement between methodologies, the former resulting in an inexpensive and faster procedure for the appraisal of biomass and community carbon fluxes at large spatial and temporal scales.

Description: Carbon-specific prey clearance and ingestion rates of 1.5-mm tentaculate larvae of the ctenophore Mnemiopsis leidyi increased linearly between 6 and 25°C but declined between 25 and 30°C. Both absolute (length) and carbon-specific growth rate increased linearly with increasing temperature. The latter was 0.87 d –1 at 25°C. Extremely low or negative growth rates observed at 6 and 30°C help define the thermal limits to population growth of this successful biological invader.

Description: We have supplemented available, concurrent measurements of fresh weight ( W , g) and body carbon (C, g) (46 individuals, 14 species) and nitrogen (N, g) (11 individuals, 9 species) of marine gelatinous animals with data obtained during the global ocean MALASPINA 2010 Expedition (totalling 267 individuals and 33 species for the W versus C data; totalling 232 individuals and 31 species for the N versus C data). We then used those data to test the allometric properties of the W versus C and N versus C relationships. Overall, gelatinous organisms contain 1.13 ± 1.57% of C (by weight, mean ± SD) in their bodies and show a C:N of 4.56 ± 2.46, respectively, although estimations can be improved by using separate conversion coefficients for the carnivores and the filter feeders. Reduced major axis regression indicates that W increases isometrically with C in the carnivores (cnidarians and ctenophores), implying that their water content can be described by a single conversion coefficient of 173.78 gW(g C) –1 , or a C content of 1.17 ± 1.90% by weight, although there is much variability due to the existence of carbon-dense species. In contrast, W increases more rapidly than C in the filter feeders (salps and doliolids), according to a power relationship W = 446.68C 1.54 . This exponent is not significantly different from 1.2, which is consistent with the idea that the watery bodies of gelatinous animals represent an evolutionary response towards increasing food capture surfaces, i.e. a bottom-up rather than a top-down mechanism. Thus, the available evidence negates a bottom-up mechanism in the carnivores, but supports it in the filter feeders. Last, N increases isometrically with C in both carnivores and filter feeders with C:N ratios of 3.89 ± 1.34 and 4.38 ± 1.21, respectively. These values are similar to those of compact, non-gelatinous organisms and reflect a predominantly herbivorous diet in the filter feeders, which is confirmed by a difference of one trophic level between filter feeders and carnivores, according to stable N isotope enrichment data.

Description: I explored mortality estimation for stage-structured populations, building on previous work that applied vertical life-table methods to populations of copepods. A new Bayesian approach for estimating mortality rates accounts for uncertainties in stage duration and number counted by stage, which have not been fully incorporated into previous analyses. This method assumes that mortality is similar among similar life stages. Results using simulated data show that realistic values of the standard deviation of stage duration and number of individuals counted result in reliable mortality estimates, though with wide confidence intervals. This uncertainty obscures variation in estimated mortality between successive stages and can also obscure bias due to violation of underlying assumptions such as that of a stable stage distribution. More importantly, the uncertainty calls into question many previous mortality estimates across pairs of life stages that do not account for these sources of uncertainty. The method was applied to an introduced population of the brackish-water cyclopoid copepod Limnoithona tetraspina in the San Francisco Estuary. Despite the uncertainties, results were interpretable: mortality was highest in nauplii and lowest in adults, probably because of high vulnerability of nauplii to invertebrate predators and low vulnerability of adults to fish predation.

Description: The strong La Niña of 2010–2011 provided an opportunity to investigate the ecological impacts of El Niño-Southern Oscillation on coastal plankton communities using the nine national reference stations around Australia. Based on remote sensing and across the entire Australian region 2011 (La Niña) was only modestly different from 2010 (El Niño) with the average temperature declining 0.2%, surface chlorophyll a up 3% and modelled primary production down 14%. Other changes included a poleward shift in Prochlorococcus and Synechococcus . Along the east coast, there was a reduction in salinity, increase in nutrients, Chlorophytes and Prasinophytes (taxa with chlorophyll b , neoxanthin and prasinoxanthin). The southwest region had a rise in the proportion of 19-hexoyloxyfucoxanthin; possibly coccolithophorids in eddies of the Leeuwin Current and along the sub-tropical front. Pennate diatoms increased, Ceratium spp. decreased and Scrippsiella spp. increased in 2011. Zooplankton biomass declined significantly in 2011. There was a reduction in the abundance of Calocalanus pavo and Temora turbinata and increases in Clausocalanus farrani , Oncaea scottodicarloi and Macrosetella gracilis in 2011. The changes in the plankton community during the strong La Niña of 2011 suggest that this climatic oscillation exacerbates the tropicalization of Australia.

Description: Jellyfish are effective predators on mesozooplankton and release large amounts of dissolved organic matter. Nevertheless, jellyfish initiated trophic cascades and bottom-up influences impacting lower trophic levels have received limited attention. We conducted a mesocosm experiment to quantify simultaneous top-down and bottom-up effects of a common jellyfish, Cyanea capillata , in a natural plankton community during autumn. Treatments were 0, 2 or 5 jellyfish per 2.5 m 3 mesocosm, four replicates each, with initial additions of inorganic nutrients. Primary and bacterial production, species abundance and composition of several trophic levels and nutrient and carbon dynamics were followed during the 8-day experiment. Multivariate statistics and generalized additive mixed modelling were applied to test whether jellyfish carbon concentration (0–1.26 mg jellyC L –1 ) in the mesocosms affected the variables monitored. Unexpected negligible predatory impact of jellyfish on mesozooplankton was observed, potentially related to jellyfish senescence. Community compositions of bacteria, phytoplankton and mesozooplankton changed with time, but did not differ between treatments. However, nutrient regeneration by jellyfish was evident, and jellyfish had a positive impact on total and specific bacterial production, total primary production and the 〉10 µm chlorophyll a fraction. Bottom-up influences from abundant jellyfish could thus stimulate productivity in nutrient depleted autumnal surface waters.

Description: The feeding ecology of Blackfordia virginica was evaluated concurrently with their ecophysiological condition in a temperate estuary. The diet of B. virginica is composed not only of metazooplankton, as commonly observed for other jellyfish species, but also of phytoplankton, ciliates and detritus. This feeding behavior might explain their good nutritional condition and sustainable growth during bloom peaks, when zooplankton abundance has already decreased significantly.

Description: The pelagic dynamics of the cosmopolitan scyphozoan Aurelia sp . was investigated in three French Mediterranean lagoons, Thau, Berre and Bages-Sigean, which harbour resident populations. The annual cycles showed a common univoltine pattern in all lagoons where the presence of pelagic stages in the water column lasted ~8 months. Field observations showed a release of ephyrae in winter time followed by pronounced growth between April and July, when individuals reached the largest sizes, before disappearing from the water column. Maximum abundance of ephyrae and medusae were registered in Thau. Medusae abundance attained a maximum of 331 ind 100 m –3 in Thau, 18 ind 100 m –3 in Berre and 7 ind 100 m –3 in Bages-Sigean lagoons. Temperature and zooplankton abundance appeared as leading factors of growth, where Bages-Sigean showed the population with higher growth rates (2.66 mm day –1 ) and maximum size (32 cm), followed by Thau (0.57–2.56 mm day –1 ; 22.4 cm) and Berre (1.57–2.22 mm day –1 ; 17 cm). The quantification of environmental windows used by the species showed wider ranges than previously reported in the Mediterranean Sea, which suggests a wide ecological plasticity of Aurelia spp. populations in north-western Mediterranean lagoons.

Description: Two lobate ctenophores, Bolinopsis infundibulum and Mnemiopsis leidyi , occur in the North Sea. Stomach contents of field-collected B. infundibulum were recorded and clearance rates for cladocerans and copepods calculated. In starvation experiments, daily body carbon losses of 2.2 and 1.2% and total carbon content losses of 76 and 63% were observed for B. infundibulum (after 68 days) and M. leidyi (after 67 days), respectively.

Description: The ctenophore Mnemiopsis leidyi is characterized by high growth rates and a large reproductive capacity. However, reproductive dynamics are not yet well understood. Here, we present laboratory data on food-dependent egg production in M. leidyi and egg hatching time and success. Further, we report on the reproduction of laboratory-reared and field-caught animals during starvation. Our results show that the half-saturation zooplankton prey concentration for egg production is reached at food levels of 12–23 µgC L –1 , which is below the average summer food concentration encountered in invaded areas of northern Europe. Furthermore, starved animals continue to produce eggs for up to 12 days after cessation of feeding with high overall hatching success of 65–90%. These life history traits allow M. leidyi to thrive and reproduce in environments with varying food conditions and give it a competitive advantage under unfavourable conditions. This may explain why recurrent population blooms are observed and sustained in localized areas in invaded northern Europe, where water exchange is limited and zooplankton food resources are quickly depleted by M. leidyi . We suggest that these reproductive life history traits are key to its invasion success.

Description: Surface-dwelling colonies of Velella velella occur throughout tropical to cold-temperate oceans of the world and sometimes are stranded in masses along hundreds of kilometers of beaches. Large-scale blooms in the Western Mediterranean Sea in 2013 and 2014 allowed the study of diet, prey digestion times and predation rates. Gastrozooid content analyses showed that 59% of the 769 identified prey were euphausiid larvae (calytopsis and furcilia) captured at night. Copepods (41%), fish eggs (2.2%) and larvae (0.5%) were captured both at day and night. Digestion times at ambient temperature (~17°C) of calytopsis, furcilia and copepods were estimated to be 〉6.5, 4.4 and 3.9 h, respectively. Estimated prey consumption was substantially lower in 2014 than in 2013 (41 vs. 75 prey day –1 colony –1 ). Velella velella and other gelatinous species bloomed in the Mediterranean Sea and the northeastern Atlantic and Pacific oceans in 2013 and 2014. Because of the wide distribution of V. velella colonies, their mass occurrences, potential importance as predators and competitors of fish, additional production from symbiotic zooxanthellae and stranding on beaches, they could be important in open-ocean carbon cycling and in transport of pelagic production to landmasses.

Description: The Matachewan Large Igneous Province (LIP) is interpreted to have formed during the early stages of mantle plume-induced continental break-up in the early Proterozoic. When the Matachewan LIP is reconstructed to its original configuration with units from the Superior Craton and other formerly adjacent blocks (Karelia, Kola, Wyoming and Hearne), the dyke swarms, layered intrusions and flood basalts, emplaced over the lifetime of the province, form one of the most extensive magmatic provinces recognized in the geological record. New geochemical data allow, for the first time, the Matachewan LIP to be considered as a single, coherent entity and show that Matachewan LIP rocks share a common tholeiitic composition and trace element geochemistry, characterized by enrichment in the most incompatible elements and depletion in the less incompatible elements. This signature, ubiquitous in early Proterozoic continental magmatic rocks, may indicate that the Matachewan LIP formed through contamination of the primary magmas with lithospheric material or that the early Proterozoic mantle had a fundamentally different composition from the modern mantle. In addition to the radiating geometry of the dyke swarms, a plume origin for the Matachewan LIP is consistent with the geochemistry of some of the suites; these suites are used to constrain a source mantle potential temperature of c. 1500–1550°C. Comparison of these mantle potential temperatures with estimated temperatures for the early Proterozoic upper mantle indicates that they are consistent with a hot mantle plume source for the magmatism. Geochemical data from coeval intrusions suggest that the plume head was compositionally heterogeneous and sampled material from both depleted and enriched mantle. As has been documented with less ancient but similarly vast LIPs, the emplacement of the Matachewan LIP probably had a significant impact on the early Proterozoic global environment. Compilation of the best age estimates for various suites shows that the emplacement of the Matachewan LIP occurred synchronously with the Great Oxidation Event. We explore the potential for the eruption of this LIP and the emission of its associated volcanic gases to have been a driver of the irreversible oxygenation of the Earth.

Description: Olivine major and trace element compositions from 12 basalts from the southern Payenia volcanic province in Argentina have been analyzed by electron microprobe and laser ablation inductively coupled plasma mass spectrometry. The olivines have high Fe/Mn and low Ca/Fe and many fall at the end of the global olivine array, indicating that they were formed from a pyroxene-rich source distinct from typical mantle peridotite. The olivines with the highest Fe/Mn have higher Zn/Fe, Zn and Co and lower Co/Fe than the olivines with lower Fe/Mn, also suggesting contributions from a pyroxene-rich source. Together with whole-rock radiogenic isotopes and elemental concentrations, the samples indicate mixing between two mantle sources: (1) a pyroxene-rich source with EM-1 ocean island basalt type trace element and isotope characteristics; (2) a peridotitic source with more radiogenic Pb that was metasomatized by subduction-zone fluids and/or melts. The increasing contributions from the pyroxene-rich source in the southern Payenia basalts are correlated with an increasing Fe-enrichment, which caused the olivines to have lower forsterite contents at a given Ni content. Al-in-olivine crystallization temperatures measured on olivine–spinel pairs are between 1155 and 1243°C and indicate that the magmas formed at normal upper mantle (asthenospheric) temperatures of ~1350°C. The pyroxene-rich material is interpreted to have been brought up from the deeper parts of the upper mantle by vigorous asthenospheric upwelling caused by break-off of the Nazca slab south of Payenia during the Pliocene and roll-back of the subducting slab beneath Payenia. The pyroxene-rich mantle mixed with peridotitic metasomatized South Atlantic mantle in the mantle wedge beneath Payenia.

Description: To reconstruct the magmatic–hydrothermal processes leading to porphyry Mo ore formation at the Climax Mo mine, Colorado, four magma units that were emplaced before, during and shortly after the mineralization events were investigated: (1) a pre-mineralization white dike of the Alma district; (2) the syn-mineralization Chalk Mountain Rhyolite; (3) a late- to post-mineralization rhyolite porphyry dyke; (4) a mafic enclave within the productive Bartlett stock. Melt inclusions, mineral inclusions and fluid inclusions in quartz phenocrysts were investigated by means of laser ablation inductively coupled plasma mass spectrometry, electron microprobe and microthermometry. Based on melt inclusion data both the Chalk Mountain Rhyolite and the rhyolite porphyry were ~10 times more fractionated than average granite and show geochemical characteristics of topaz rhyolites. They were saturated in magnetite, Mn-rich ilmenite, fluorite, aeschynite, monazite, pyrrhotite and thorite, and crystallized predominantly at 710–730°C, 1·2–2·6 kbar and log f O 2 FMQ + 2·2 (where FMQ is fayalite–magnetite–quartz). The silicate melt of the Chalk Mountain Rhyolite contained 3·5 ± 0·4 wt % F, 0·09 ± 0·03 wt % Cl, ≥ 3·0 wt % H 2 O, 15–90 µg g –1 Cs, 500–1500 µg g –1 Rb and 5–7 µg g –1 Mo, whereas that of the rhyolite porphyry contained 1·1 ± 0·3 wt % F and 4·9 ± 1·2 wt % H 2 O, but otherwise had a virtually identical major and trace element composition. The fluid exsolving from the latter melt had a bulk salinity of 10 ± 2 wt % NaCl equiv and contained of the order of 100 µg g –1 Mo. After emplacement of the Chalk Mountain Rhyolite magma at subvolcanic levels, extremely fractionated silicate melts coexisting with hypersaline brines (salt melts) and low-density vapor percolated at near-solidus conditions through the rock. These silicate melts contained 6·6 ± 0·4 wt % F, ≥ 7·5 ± 0·6 wt % H 2 O, 0·51 ± 0·05 wt % Cl, and up to 0·5 wt % Cs and 100 µg g –1 Mo, whereas the hypersaline brines contained 1–2 wt % Cs and 0·3–0·6 wt % Mo. However, owing to their negligible masses these liquids are unlikely to have played a major role in the mineralization process. The majority of Mo in the Climax deposit appears to have been derived from melts containing 5–7 µg g –1 Mo and bulk fluids containing ~100 µg g –1 Mo. These concentrations are similar to those found in similarly fractionated melts and fluids in barren and sub-economically mineralized intrusions. However, whereas in the latter intrusions fractionated melts occurred in a rather dispersed state, they seem to have been present as large, coherent masses in the apical parts of Climax-type porphyry Mo-forming magma systems. Efficient segregation of fractionated melts and fluids into the top of mineralizing magma chambers appears to have been promoted by high fluorine concentrations in the silicate melt, which was partly a primary feature, and partly an indirect consequence of other characteristics of within-plate magmatism.

Description: The Skaergaard intrusion, Greenland, is the type locality for Skaergaard-type mineralizations. Mineralization levels are perfectly concordant with igneous layering, up to 5 m thick, internally fractionated, and contain crystallized sulphide droplets and precious metal alloys, sulphides, arsenides and telluride. Immiscible Cu-rich sulphide droplets, formed in a mush zone below the roof, scavenged precious metals. They were subsequently dissolved and transported to the floor in late-formed, immiscible, Fe-rich mush melts. Mineralized stratigraphic intervals of floor gabbro formed in ‘proto-macrolayers‘, owing to local sulphide saturation in melt concentrated between floating plagioclase and sinking clinopyroxene. The floor mineralization is divided into four stratigraphic sections. Formation of the Lower Platinum Group Element Mineralization (LPGEM) involved: (1) crystallization of the bulk liquid liquidus paragenesis and in situ fractionation; (2) sulphide saturation and formation of sulphide droplets in melt in the upper part of ‘proto-macrolayers‘. After further in situ fractionation, the following steps occurred: (3) the onset of silicate–silicate immiscibility and the consequent loss of buoyant and immiscible Si-rich melt; (4) dissolution of unprotected droplets of sulphide melt present in the Fe-rich mush melt; (5) compaction-driven upwards loss of residual mush melt enriched in, for example, Au. The LPGEM preserves upward increasing bulk Pd/Pt (~6–13) owing to a continued supply of PGE and Au, with high Pd/Pt. The further development of the LPGEM ceased as the supply of precious metals to the floor waned. The Upper PGE Mineralization (UPGEM) subsequently formed from precious metals recycled in the floor. The UPGEM is characterized by increasing Au substitution in PGE phases, and a decrease in total PGE and Pd/Pt owing to upward fractionation in migrating mush melts and exhaustion of Pd and Pt. An upper Au-rich mineralization level (UAuM) was caused by late remobilization of Au and deposition on grain boundaries in fully crystallized gabbro. Cu concentrations (~150 ppm) are not correlated with PGE and Au. Repeated Cu mineralization levels (CuM), attaining 〉1000 ppm, occur above the Au levels, caused by local mush layer sulphide saturation. PGE, Au and Cu distributions in the floor mineralization reflect sub-liquidus, but supra-solidus, processes and reactions in mushes at the roof, wall and floor. Constraints provided by a new model for the mineralization provide the basis for re-evaluation of the solidification processes in the Skaergaard intrusion. We have identified the importance of extensive in situ fractionation and intrusion-wide elemental redistributions in immiscible Fe- and Si-rich silicate melts. Our model characterizes the floor cumulates as bulk liquid orthocumulates containing an upwards-increasing proportion crystallized from Fe-rich, immiscible mush melt. The roof-rocks are complementary to the floor, with downwards increasing proportions crystallized from the conjugate Si-rich melt. Petrographic observations and the relative timing of crystallization support the hypothesis that crystallization was restricted to marginal mush zones. Bulk melt remaining in the magma chamber evolved not, as generally assumed, as a result of loss of crystals grown from the bulk melt, but as the consequence of mixing with recycled and evolved melt expelled from the mush by compaction. Redistribution of Fe in immiscible melts may be common to mafic intrusions and puts into question the validity of petrogenetic modelling of bulk liquids in mafic intrusions based only on consideration of floor cumulates.

Description: We have determined experimentally the hydrous phase relations and trace element partitioning behaviour of ocean floor basalt protoliths at pressures and temperatures (3 GPa, 750–1000°C) relevant to melting in subduction zones. To avoid potential complexities associated with trace element doping of starting materials we have used natural, pristine mid-ocean ridge basalt (MORB from Kolbeinsey Ridge) and altered oceanic crust (AOC from Deep Sea Drilling Project leg 46, ~20°N Atlantic). Approximately 15 wt % water was added to starting materials to simulate fluid fluxing from dehydrating serpentinite underlying the oceanic crust. The vapour-saturated solidus is sensitive to basalt K 2 O content, decreasing from 825 ± 25°C in MORB (~0·04 wt % K 2 O) to 750°C in AOC (~0·25 wt % K 2 O). Textural evidence indicates that near-solidus fluids are sub-critical in nature. The residual solid assemblage in both MORB and AOC experiments is dominated by garnet and clinopyroxene, with accessory kyanite, epidote, Fe–Ti oxide and rutile (plus quartz–coesite, phengite and apatite below the solidus). Trace element analyses of quenched silica-rich melts show a strong temperature dependence of key trace elements. In contrast to the trace element-doped starting materials of previous studies, we do not observe residual allanite. Instead, abundant residual epidote provides the host for thorium and light rare earth elements (LREE), preventing LREE from being released (LREE 〈3 ppm at 750–900°C). Elevated Ba/Th ratios, characteristic of many arc basalts, are found to be generated within a narrow temperature field above the breakdown temperature of phengite, but below exhaustion of epidote. Melts with Ba/Th 〉1500 and La/Sm PUM (where PUM indicates primitive upper mantle) ~1, most closely matching the geochemical signal of arc lavas worldwide, were generated from AOC at 800–850°C.

Description: We constrain the physical nature of the magma reservoir and the mechanisms of rhyolite generation at Yellowstone caldera via detailed characterization of zircon and sanidine crystals hosted in three rhyolites erupted during the ( c . 170–70 ka) Central Plateau Member eruptive episode—the most recent post-caldera magmatism at Yellowstone. We present 238 U– 230 Th crystallization ages and trace-element compositions of the interiors and surfaces (i.e. unpolished rims) of single zircon crystals from each rhyolite. We compare these zircon data with 238 U– 230 Th crystallization ages of bulk sanidine separates coupled with chemical and isotopic data from single sanidine crystals. Zircon age and trace-element data demonstrate that the magma reservoir that sourced the Central Plateau Member rhyolites was long-lived (150–250 kyr) and genetically related to the preceding episode of magmatism, which occurred c . 256 ka. The interiors of most zircons in each rhyolite were inherited from unerupted material related to older stages of Central Plateau Member magmatism or the preceding late Upper Basin Member magmatism (i.e. are antecrysts). Conversely, most zircon surfaces crystallized near the time of eruption from their host liquids (i.e. are autocrystic). The repeated recycling of zircon interiors from older stages of magmatism demonstrates that sequentially erupted Central Plateau Member rhyolites are genetically related. Sanidine separates from each rhyolite yield 238 U– 230 Th crystallization ages at or near the eruption age of their host magmas, coeval with the coexisting zircon surfaces, but are younger than the coexisting zircon interiors. Chemical and isotopic data from single sanidine crystals demonstrate that the sanidines in each rhyolite are in equilibrium with their host melts, which considered along with their near-eruption crystallization ages suggests that nearly all Central Plateau Member sanidines are autocrystic. The paucity of antecrystic sanidine crystals relative to antecrystic zircons requires a model in which eruptible rhyolites are generated by extracting melt and zircons from a long-lived mush of immobile crystal-rich magma. In this process the larger sanidine crystals remain trapped in the locked crystal network. The extracted melts (plus antecrystic zircon) amalgamate into a liquid-dominated (i.e. eruptible) magma body that is maintained as a physically distinct entity relative to the bulk of the long-lived crystal mush. Zircon surfaces and sanidines in each rhyolite crystallize after melt extraction and amalgamation, and their ages constrain the residence time of eruptible magmas at Yellowstone. Residence times of the large-volume rhyolites (~40–70 km 3 ) are ≤1 kyr (conservatively 〈6 kyr), which suggests that large volumes of rhyolite can be generated rapidly by extracting melt from a crystal mush. Because the lifespan of the crystal mush that sourced the Central Plateau Member rhyolites is two orders of magnitude longer than the residence time of eruptible magma bodies within the reservoir, it is apparent that the Yellowstone magma reservoir spends most of its time in a largely crystalline (i.e. uneruptible) state, similar to the present-day magma reservoir, and that eruptible magma bodies are ephemeral features.

Description: We report and interpret new geochemical and Pb–Sr–Nd isotopic data from 325 samples of shield, late-shield, postshield, and rejuvenated stage lavas from Kauai and Niihau, the two most northwesterly islands in the Hawaiian island chain. Kauai is unique in the Hawaiian chain in that it exhibits a near-continuous geochemical transition from shield to postshield to rejuvenated stage volcanism between 4·4 and 3·6 Ma and has been continuously active over ~6 Myr. From c . 5·7 to 4·3 Ma, the shield stage of both islands produced tholeiitic basalts typical of other Hawaiian shield volcanoes. The Niihau basalts are more evolved and have high Gd/Yb compared with Kauai, indicating a higher residual garnet content in the source. Both Kauai and Niihau shield basalts have Kea-like trace element ratios, but isotopic ratios are transitional between Kea- and Loa-like compositions. The geochemical similarity of the two shields indicates that mantle sources in different regions of the plume source were similar, and that the 〈2 Ma Loa and Kea trends of the southeastern Hawaiian volcanoes are not observed. More Loa-like compositions are evident in shield lavas from eastern Kauai, where the enhanced Loa composition may reflect melting of low-melting temperature plume components as the island migrates off the hotter, more Kea-like, center of the Hawaiian plume. Postshield lavas and intrusive rocks on both islands are rare: Kauai includes alkalic basalts, hawaiites and mugearites that are isotopically homogeneous and include a significant depleted mantle component compared with the shield basalts, whereas the Niihau late-shield and postshield rocks consist of highly contrasting transitional tholeiites or basanites with variable but shield-like isotopic compositions. The Niihau postshield rocks represent variable, but lower degrees of melting of the shield mantle source. Large volumes of rejuvenated stage lavas cover both islands and also form submarine cone fields, but lava compositions are different. On Kauai, rejuvenated lavas range from melilitite to alkalic basalt with trace element, Nd isotope, and Pb isotope ratios that vary as a function of Th and SiO 2 content. Low-degree (high-Th) melts are dominated by a mixed Kea–Loihi component and high-degree (low-Th) melts include more of a depleted rejuvenated component (DRC) typified by rejuvenated stage lavas and xenoliths from nearby Kaula Island. With the exception of a single basanite, the Niihau rejuvenated stage lavas are uniformly alkalic basalt, with Sr and Ba excesses combined with depleted Th and Nb abundances relative to the light rare earth elements. Rejuvenated stage alkalic basalts from both islands are dominated by contributions from the DRC, which have high Sr/Ce and 87 Sr/ 86 Sr but low 206 Pb/ 204 Pb. The Sr-rich, possibly carbonate-bearing, DRC component may be distributed patchily in the rejuvenated stage mantle source such that, where present, the degree of partial melting was enhanced compared with the degree of partial melting of the Sr-poor, mixed Kea–Loihi component. Given the lack of a hiatus between postshield and rejuvenated stages on Kauai, the rejuvenated mantle source is already able to melt at the tail end of shield stage activity and no secondary melting mechanism is required to explain the rejuvenated stage.

Description: Several hypotheses exist that describe phytoplankton spring blooms in temperate and subpolar oceans: the critical depth, shoaling mixed layer (ML), critical turbulence, onset of stratification and disturbance-recovery hypotheses. These theories appear to be mutually exclusive and none of them describe the annual cycle of phytoplankton biomass. Here, we present a model of the annual cycle in phytoplankton that recognizes that phytoplankton are not always mixed throughout the so-called ML, and that it is important to distinguish between the surface biomass and depth-integrated phytoplankton. Once these important distinctions are made, the annual cycles and blooms in surface and depth-integrated phytoplankton can be described straightforwardly in terms of the physical drivers and biotic responses.

Description: Impacts of climate change on marine ecosystems have become increasingly apparent during the past decades. In consequence, it is necessary to study how these alterations can affect the habitat and population dynamics of key organisms. Here we used a video plankton recorder (VPR) to investigate the effect of climate-induced habitat changes on the copepod Pseudocalanus acuspes , a key species in the Baltic Sea. The VPR allowed the observation of reproducing copepod females, identified by attached egg sacs, usually lost during traditional net sampling. We compared the small-scale distribution of our target species during non-inflow and inflow periods. Our study showed a large increase in the availability of suitable habitat after the inflow event due to improved oxygen and salinity conditions. Furthermore, increased copepod abundance and a deeper and wider vertical distribution was apparent. Applying a new approach to estimate in situ egg production rates from VPR-derived images revealed no changes. However, we observed increased offspring survival with improved hydrographic conditions pointing toward the importance of salinity and oxygen for the population dynamics of Baltic P. acuspes . Our observations illustrate the strong impact that climate change can have on the habitat of key marine ecosystem species, important for overall ecosystem dynamics.

Description: Research on nutrient controls of planktonic productivity tends to focus on a few standard fractions of inorganic or total nitrogen (N) and phosphorus (P). However, there is a wide range in the degree to which land-derived dissolved organic nutrients can be assimilated by biota. Thus, in systems where such fractions form a majority of the macronutrient resource pool, including many boreal inland waters and estuaries, our understanding of bacterio- and phytoplankton production dynamics remains limited. To adequately predict aquatic productivity in a changing environment, improved standard methods are needed for determining the sizes of active (bioavailable) pools of N, P and organic carbon (C). A synthesis of current knowledge suggests that variation in the C:N:P stoichiometry of bioavailable resources is associated with diverse processes that differentially influence the individual elements across space and time. Due to a generally increasing organic nutrient bioavailability from C to N to P, we hypothesize that the C:N and N:P of bulk resources often vastly overestimates the corresponding ratios of bioavailable resources. It is further proposed that basal planktonic production is regulated by variation in the source, magnitude and timing of terrestrial runoff, through processes that have so far been poorly described.

Description: Understanding the factors that regulate the abundance, size structure and community structure of cladocerans is an important goal of aquatic ecologists. While both top-down and bottom-up factors help to structure cladoceran communities, there may be interactions between these factors. We conducted a mesocosm study to determine how alien large-bodied Daphnia , zebra mussels and fish affected cladoceran community and size structure. We found that large-bodied Daphnia reduced algal resources and the fecundity of smaller bodied cladocerans. Fish removed the large-bodied Daphnia magna from the mesocosms and shifted the cladoceran community to a smaller body size. Fish also appeared to promote increases in cladoceran diversity through the coexistence of several smaller bodied taxa. In contrast, zebra mussels increased cyanobacteria and helped to promote the success of the alien Daphnia , but reduced the biomass of small-bodied cladocerans. Zebra mussels reduced the carbon (C):phosphorus (P) ratio of the phytoplankton in the mesocosms which may have favored the relatively P-limited Daphnia . Combined, our results highlight the complex interactions of multiple factors that help to regulate cladoceran community and size structure.

Description: In freshwater systems, Daphnia has been demonstrated to show adaptive responses following the light–dark cycle. The adjustment of these responses to the change of day and night is probably transmitted via the hormone melatonin. The rate-limiting enzyme in melatonin synthesis is the arylalkylamine N-transferase (AANAT). We identified three genes coding for insect-like AANATs in Daphnia , of which we measured the gene expression in an ecologically relevant light–dark cycle. We demonstrated that Daphnia 's insect-like AANAT gene expression oscillated in a daily manner, and that the highest peak of expression after the onset of darkness was followed by a peak of melatonin production at midnight. Moreover, we could show an oscillation of endogenous melatonin synthesis in Daphnia . In most organisms, melatonin synthesis is due to rhythmic expression of genes of the circadian clock, since transcription of aanat s is directly linked to a circadian transcription factor. We could demonstrate that putative clock genes and insect-like AANAT genes of Daphnia were equally expressed. Therefore, we propose that melatonin synthesis is coupled to the expression of Daphnia clock genes, and that insect-like AANATs of crustaceans have a similar function as AANATs of vertebrates: The initiation of melatonin synthesis. In future studies with Daphnia , it will be necessary to take the time of day into account since melatonin concentrations might influence stress responses.

Description: It has been hypothesized that terrestrial particulate organic matter (t-POM) makes important contributions to Daphnia production in some lakes. We conducted a series of feeding experiments to explore the fatty acid responses in Daphnia to diets comprised of different terrestrial resources (i.e. Alnus rubra , Phragmites australis , Betula nana and Betula pendula ) and mixed diets with terrestrial and phytoplankton ( Scenedesmus or Cryptomonas ) resources. When fed 100% phytoplankton, Daphnia had very similar ( r 2 〉 0.80) fatty acid profiles to their diets, whereas Daphnia that consumed t-POM diets had weak correlations ( r 2 = 0.002–0.56) with the corresponding diet sources. Unusual 16 carbon chain polyunsaturated fatty acids (16:26, 16:33 and 16:43), linoleic acid (18:26) and α-linolenic acid (18:33) were diagnostic fatty acids for Scenedesmus and Daphnia that consumed this alga. Stearidonic acid (18:43) and eicosapentaenoic acid (20:53) were diagnostic for Cryptomonas and Daphnia that consumed this diet. All of the t-POM resources were characterized by a high content of saturated fatty acids (SAFA; 79 ± 12%), especially the diagnostic long-chain SAFA (20:0, 22:0, 24:0, 26:0, 28:0). Daphnia that consumed t-POM assimilated very little of these terrestrial biomarkers, but the shorter chain SAFA 16:0 and 18:0 were very prevalent in juvenile and adult Daphnia that consumed terrestrial plant matter. The -3:-6 ratios were distinctive between terrestrial (0.3–1.6) and phytoplankton resources (3–15), and this ratio in Daphnia was strongly associated with their diets ( r 2 = 0.88). These results suggest that Daphnia , and perhaps zooplankton in general, preferentially retain algae-derived 3 fatty acids, and low -3:-6 ratios in Daphnia indicate a mainly terrestrial diet or poor nutritional condition.

Description: In order to trace community dynamics and reticulate evolution in hybrid species complexes, long-term comparative studies of natural populations are necessary. Such studies require the development of tools for fine-scale genetic analyses. In the present study, we developed species-diagnostic SNP-based markers for hybridizing freshwater crustaceans: the multispecies Daphnia longispina complex. Specifically, we took advantage of transcriptome data from a key species of this hybrid complex, the annotated genome of a related Daphnia species and well-defined reference genotypes from three parental species. Altogether eleven nuclear loci with several species-specific SNP sites were identified in sequence alignments of these reference genotypes from three parental species and their interspecific hybrids. A PCR-RFLP assay was developed for cost-efficient large population screening by SNP-based genotyping. Taxon assignment by RFLP patterns was nearly perfectly concordant with microsatellite genotyping across several screened populations from Europe. Finally, we were able to amplify two short regions of these loci in formaldehyde-preserved samples dating back to the year 1960. The species-specific SNP-based markers developed here provide valuable tools to study hybridization over time, including the long-term impact of various environmental factors on hybridization and biodiversity changes. SNP-based genotyping will finally allow eco-evolutionary dynamics to be revealed at different time scales.

Description: A new thermodynamic model is presented for calculating phase relations in peridotite, from 0·001 to 60 kbar and from 800°C to liquidus temperatures, in the system NCFMASOCr. This model system is large enough to simulate phase relations and melting of natural peridotite and basaltic liquids. Calculations in the program thermocalc illustrate mantle phase relationships and melting conditions, specifically for the peridotite composition KLB-1. The garnet–spinel transition zone intersects the solidus at 21·4–21·7 kbar, where both Fe 3+ and Cr increase spinel stability, expanding the width of the transition. Orthopyroxene is lost at the solidus at 42 kbar in KLB-1, although this pressure is very sensitive to bulk composition. Calculated oxidation states are in excellent agreement with measured log f O 2 for xenolith suites with mantle Fe 2 O 3 contents in the range 0·1–0·3 wt %. It appears that mantle oxidation state is not just a simple function of P and T , but depends on phase assemblage, and may vary in a complex way within a single assemblage. The liquid model performs well, such that calculated solidus, melt productivity and liquid compositions compare favourably with those of experimental studies, permitting its use in interpolating between, and extrapolating from, experimental P–T conditions. Experimentally challenging but geologically useful regimes can be explored, such as subsolidus samples and very low melt fractions, with application to both mantle xenoliths and the origin of basalt.

Description: The anorthositic members of the Mealy Mountains Intrusive Suite (MMIS; Labrador, Canada) are host to 0·5–5 m diameter pegmatitic, pod-like segregations, originally described as graphic granite pods. U–Pb zircon geochronology confirms that the pods are coeval with the 1650–1630 Ma emplacement age range for the MMIS, yielding ages of 1654 ± 8 to 1628 ± 3·5 Ma. Petrographic and geochemical analysis of five pods from anorthositic rocks of the MMIS reveals that the pods have a diverse compositional range from monzodiorite to granite, varying from Fe-rich and Si-poor, to Fe-poor and Si-rich compositions. Fe-rich, Si-poor pods in the MMIS and other massifs (e.g. Laramie Anorthosite Complex) tend to be hosted by olivine-bearing anorthosites, whereas Si-rich, Fe-poor pods are hosted by pyroxene-bearing anorthosites. Each pod shows a range of graphic, myrmekitic and symplectitic textures, along with distinctive mineral assemblages (e.g. apatite and zircon) and highly enriched trace-element compositions. Evolved mineral assemblages, high concentrations of Fe, Ti and P (and in some cases SiO 2 ), and 10–1000 x chondrite enrichment in light rare earth elements, U, Th and Rb indicate that many of the pods are highly fractionated. The array of textural intergrowths provides clues about the final stages of crystallization in the pods and, by extension, the anorthosites. Macroscopic quartz–K-feldspar graphic intergrowths indicate high-viscosity, fluid-bearing and significantly undercooled magmatic conditions, whereas microscopic myrmekitic (plagioclase–quartz) and symplectitic (plagioclase–orthopyroxene) intergrowths on primary grain boundaries indicate replacement of phases in the presence of reactive fluids. In assessing the nature of these pegmatitic pods based on field, petrographic and geochemical evidence, we conclude that they represent the fluid-bearing, late-stage crystallization products of a residual liquid in the massif anorthosite system. The Fe and Si compositional variations observed in these late-stage pods can be linked to a fundamental olivine–pyroxene dichotomy observed in most Proterozoic anorthosite massifs, suggesting that pulses of magma experience variable contamination (in amount and/or composition) leading to varying differentiation paths. A range of lithologies (monzonites, monzonorites, ferrodiorites and jotunites) observed in similar pod-like structures, as well as dykes and plutons, has been observed in other Proterozoic anorthosite massifs and all have, at one time or another, been interpreted as the residual liquids of anorthosite crystallization. Our observation of in situ pods with similar compositions to all of the aforementioned lithologies, and displaying textures indicative of late-stage crystallization, supports the notion that all of these associated lithologies can be interpreted as comagmatic with, but variably contaminated and isolated residual liquids of, anorthosite crystallization. However, using isotopic evidence we cannot support the notion that the far larger granitic plutons associated with Proterozoic anorthosites are also residual liquids of anorthositic magma fractionation.

Description: We report the results of experiments on two natural marine sediments with different carbonate contents (calcareous clay: CO 2 = 6·1 wt %; marl: CO 2 = 16·2 wt %) at subduction-zone conditions (3 GPa, 750–1200°C). Water (7–15 wt %) was added to the starting materials to simulate the effects of external water addition from within the subducting slab. The onset of melting is at 760°C in water-rich experiments; melt becomes abundant by 800°C. In contrast, the onset of melting in published, water-poor experiments occurs at variable temperatures with the production of significant melt fractions being restricted to more than 900°C (phengite-out). The different solidus temperatures ( T solidus ) can be ascribed to variable fluid X H2O [H 2 O/(CO 2 + H 2 O)], which, in turn, depends on bulk K 2 O, H 2 O and CO 2 . Partial melts in equilibrium with residual garnet, carbonate, quartz/coesite, epidote, rutile, kyanite, phengite, and clinopyroxene are granitic in composition, with substantial dissolved volatiles. Supersolidus runs always contain both silicate melt and solute-rich fluid, indicating that experimental conditions lie below the second critical endpoint in the granite–H 2 O–CO 2 system. Carbonatite melt coexists with silicate melt and solute-rich fluid above 1100°C in the marl. The persistence of carbonate to high temperature, in equilibrium with CO 2 -rich hydrous melts, provides a mechanism to both supply CO 2 to arc magmas and recycle carbon into the deep Earth. The trace element compositions of the experimental glasses constrain the potential contribution of calcareous sediment to arc magmas. The presence of residual epidote and carbonate confers different trace element characteristics when compared with the trace element signal of Ca-poor marine sediments (e.g. pelagic clays). Notably, epidote retains Th and light rare earth elements, such that some melts derived from calcareous sediments have elevated Ba/Th and U/Th, and low La/Sm PUM , thereby resembling fluids conventionally ascribed to altered oceanic crust. Our results emphasize the importance of residual mineralogy, rather than source lithology, in controlling the trace element characteristics of slab-derived fluids.

Description: Three crystal-poor obsidian samples (one dacite, 67 wt % SiO 2 ; two rhyolites, 73 and 75 wt % SiO 2 ), which erupted effusively from monogenetic vents, contain sparse (〈2%) plagioclase phenocrysts that span a remarkably wide and continuous range in composition (≤30 mol % An). Many, but not all, of the plagioclase crystals display diffusion-limited growth textures (e.g. swallow-tails, skeletal, vermiform). Hypotheses to explain the paradox of a wide compositional range despite a low abundance of plagioclase include (1) incorporation of xenocrysts and/or magma mingling, (2) slow crystallization of plagioclase driven by slow cooling in a magma chamber, (3) slow crystallization of plagioclase followed by a resorption (e.g. heating) event, and (4) crystallization driven by rapid degassing (i.e. loss of melt H 2 O) ± rapid cooling during ascent. To test these hypotheses, a series of phase equilibrium experiments were conducted under pure-H 2 O fluid-saturated conditions in a cold-seal pressure vessel between 30 and 300 MPa and 750 and 950°C. The results show that the plagioclase population in each obsidian sample could have grown from their respective melts, with the exception of a single calcic core (An 60–63 ) in one sample. The results additionally rule out slow cooling in a magma chamber, because this would lead to equilibrium abundances of plagioclase (≤20%), which are far higher than what is observed in the samples (〈2%). Nor can resorption (i.e. heating) explain the low abundance of plagioclase, because this would eliminate the more sodic plagioclase crystals and hence the wide compositional range of plagioclase that is observed. The most viable hypothesis is that the sparse plagioclase phenocrysts grew relatively rapidly during magma ascent to the surface; this is consistent with the results of isothermal (850°C) continuous decompression experiments (2·9, 1·0, 0·8, and 0·1 MPa h –1 ), under pure-H 2 O fluid-saturated conditions, which were performed on one of the rhyolites (MLV-36; 73 wt % SiO 2 ) and quenched at P H2O = 89, 58 and 40 MPa. The four decompression rates correspond to degassing rates of 1·6, 0·56, 0·45 and 0·06 wt % H 2 O per day. Decompressions ≥1·0 MPa( P H2O ) h –1 , initiated above the liquidus, quenched to 100% glass at all final P H2O . Decompressions at 0·8 MPa( P H2O ) h –1 , also initiated above the liquidus, led to plagioclase crystals nearly five times larger than those grown in runs decompressed at the same rate, but initiated just below the plagioclase-in curve. It is the kinetic hindrance to nucleation that permits crystal growth to be concentrated on relatively few crystals, leading to larger crystals. Plagioclase growth rates from these experiments show that the largest phenocrysts (~1 mm) in the MLV-36 obsidian could have grown in 〈42 h. A cooling rate of ~1·2°C h –1 closely matches both the increase in melt viscosity with time and the effective undercooling with time that occurs during the 0·8 MPa( P H2O ) h –1 decompression over the first 50 h. The combined results point to crystallization of sparse plagioclase driven by relatively rapid rates of degassing ± cooling during ascent to the surface of melts that were initially above their liquidus. The obsidian samples must have been efficiently segregated as nearly 100% liquids from their respective source regions at H 2 O-fluid undersaturated conditions to attain a degree of superheating upon ascent before reaching fluid saturation.

Description: The Jurassic Vestfjella dyke swarm at the volcanic rifted margin of western Dronning Maud Land represents magmatism related to the incipient Africa–Antarctica rift zone; that is, rift-assemblage magmatism of the Karoo continental flood basalt (CFB) province. Geochemical and Nd–Sr isotopic data for basaltic and picritic dyke samples indicate diverse low-Ti and high-Ti tholeiitic compositions with Nd (180 Ma) ranging from +8 to –17. Combined with previously reported data on a subcategory of ferropicritic dykes, our new data facilitate grouping of the Vestfjella dyke swarm into seven geochemically distinct types. The majority of the dykes exhibit geochemical affinity to continental lithosphere and can be correlated with two previously identified chemical types (CT) of the wall-rock CFB lavas and are accordingly referred to as the CT1 and CT3 dykes. The less abundant Low-Nb and High-Nb dykes, a relatively enriched subtype of CT3 (CT3-E) dykes, and dykes belonging to the depleted and enriched ferropicrite suites represent magma types found only as intrusions. The chemically mid-ocean ridge basalt (MORB)-like Low-Nb and the depleted ferropicrite suite dykes represent, respectively, relatively high- and low-degree partial melting of the same overall depleted mantle (DM)-affinity source in the sublithospheric mantle. In contrast, we ascribe the chemically ocean island basalt (OIB)-like High-Nb dykes and the enriched ferropicrite suite dykes to melting of enriched components in the sublithospheric mantle. Geochemical modelling suggests that the low-Ti affinity CT1 and CT3, and high-Ti affinity CT3-E magma types of Vestfjella dyke may predominantly result from mixing of DM-sourced Low-Nb type magmas with 〈10 wt % of crust- and lithospheric mantle-derived melts. U/Pb zircon dating confirms synchronous emplacement of CT1 dykes and Karoo main-stage CFBs at 182·2 ± 0·9 and 182·2 ± 0·8 Ma, whereas two 40 Ar/ 39 Ar plagioclase plateau ages of 189·2 ± 2·3 Ma (CT1) and 185·5 ± 1·8 Ma (depleted ferropicrite suite), and a mini-plateau age of 186·9 ± 2·8 Ma (CT3-E) for the Vestfjella dykes raise the question of whether the onset of rift-zone magmatism could predate the province-wide c. 179–183 Ma main stage of Karoo magmatism. Notably variable Ca/K spectra suggest that younger 40 Ar/ 39 Ar plagioclase plateau ages of 173, 170, 164, and 154 Ma are related to crystallization of secondary minerals during the late-stage tectono-magmatic development of the Antarctic rifted margin. The occurrence of rare MORB- and OIB-like magma types in Vestfjella and along the African and Antarctic rifted margins suggests melting of geochemically variable depleted and enriched sublithospheric mantle beneath the Africa–Antarctica rift zone. Our models for the Vestfjella dyke swarm indicate that the voluminous lithosphere-affinity low-Ti and high-Ti rift-assemblage tholeiites could have been derived from MORB-like parental magmas by contamination, which implies sublithospheric depleted mantle as the principal source of the CFB magmas of the Africa–Antarctica rift zone.

Description: Monazite laser ablation–split-stream inductively coupled plasma–mass spectrometry (LASS) was used to date monazite in situ in Barrovian-type micaschists of the Moravian zone in the Thaya window, Bohemian Massif. Petrography and garnet zoning combined with pseudosection modelling show that rocks from staurolite–chlorite, staurolite, kyanite and kyanite–sillimanite zones record burial in the S 1 fabric under a moderate geothermal gradient from 4–4·5 kbar and ~530–540°C to 5 kbar and 570°C, 6–7 kbar and 600–640°C, 7·5–8 kbar and 630–650°C, and 8 kbar and 650°C, respectively. In the kyanite and kyanite–sillimanite zones, garnet rim chemistry and local syntectonic replacement of garnet by sillimanite–biotite aggregates point to re-equilibration at 5·5–6 kbar and 630–650°C in the S 2 fabric. Heterogeneously developed retrograde shear zones (S 3 ) are marked by widespread chloritization, but minor chlorite is present in the studied samples. Monazite abundance and size increase with metamorphic grade from 5 µm in the staurolite–chlorite zone to 〉100 µm in the kyanite and kyanite–sillimanite zones. Irrespective of the monazite-forming reaction, this is interpreted as the onset of limited prograde monazite growth at staurolite grade, and continued prograde monazite growth after the kyanite-in reaction, compatible with conditions of about 5·5 kbar and 570°C and 7·5 kbar and 630°C from pseudosection modelling. Monazite is zoned, showing embayments and sharp boundaries between zones, with low Y in the staurolite zone, high-Y cores and low-Y rims in the kyanite zone, and high-Y cores, a low-Y mantle and a high-Y rim in the sillimanite zone. The 207 Pb-corrected 238 U/ 206 Pb ages from three samples range from 344 ± 7 to 330 ± 7 Ma, irrespective of metamorphic grade. The dates from monazite inclusions are interpreted as the ages of the staurolite- and kyanite-in reactions along the prograde path at 340 and 337 ± 7 Ma, respectively. The monazite in the matrix (and some inclusions) is interpreted as dating the prograde crystallization at (340–337) ± 7 Ma within the S 1 fabric, and then being affected by recrystallization at or down to 332 ± 7 Ma in the S 2 and S 3 fabrics. The two groups of data, for 340–337 and 332 Ma, are significantly different when only their in-run uncertainties (±1–3 Myr) are compared and indicate a 9 ± 3 Myr period of monazite (re)crystallization. A systematic increase in heavy rare earth element (HREE) content with decreasing monazite age from 344 to 335 Ma is correlated with growth on the prograde P–T path; the drop in HREE of monazite at 335–328 Ma is assigned to recrystallization. The presence of chlorite even in the least retrogressed samples witnesses limited external fluid availability on the retrograde P–T path. Migration of this fluid was probably responsible for heterogeneous fluid-assisted recrystallization and resetting of original prograde monazite, even where included in garnet, staurolite or kyanite. It is suggested that the rocks passed the chlorite-in reaction on the retrograde path at 332 ± 7 Ma. The timing of burial in the Thaya window, a deformed part of the underthrust Brunia microcontinent, was coeval with exhumation of granulites and migmatites of the Moldanubian orogenic root at c. 340 Ma.

Description: The Fanshan intrusion in the North China Craton (NCC) is concentrically zoned with syenite in the core (Unit 1), surrounded by ultramafic rocks (clinopyroxenite and biotite clinopyroxenite; Unit 2), and an outer rim of garnet-rich clinopyroxenite and orthoclase clinopyroxenite and syenite (Unit 3). The intrusive rocks are composed of variable amounts of Ca-rich augite, biotite, orthoclase, melanite, garnet, magnetite and apatite, with minor primary calcite. Monomineralic apatite rocks, nelsonite and glimmerite exclusively occur in Unit 2. Geochemically, the Fanshan rocks are highly enriched in light rare earth elements (LREE) and large ion lithophile elements (LILE), moderately depleted in high field strength elements (HFSE), and have a limited range of Sr–Nd–O isotopic compositions. The similar mineralogy, mineral compositions, and trace element characteristics of the three units suggest that all the rocks are co-magmatic. The parental magma is ultrapotassic and is akin to kamafugite. Very low-degree partial melting of metasomatized lithospheric mantle best explains the geochemistry and petrogenesis of the parental magmas of the Fanshan intrusion. We propose that the mantle source may have been metasomatized by a hydrous carbonate-bearing melt, which has imprinted the enriched Sr–Nd isotopic signature and incompatible element enrichment with conspicuous negative Nb–Ta–Zr–Hf–Ti anomalies and LREE enrichments. The mantle source enrichment may be correlated with oceanic sediment recycling during southward subduction of the Paleo-Asian oceanic plate during the Carboniferous and Permian. We propose that crystal settling and mechanical sorting combined with repeated primitive magma replenishment and mixing with previously fractionated magma is the predominant process responsible for the formation of the apatite ores.

Description: Pointwise error estimates for the first-order div least-squares (LS) finite element method for second-order elliptic partial differential equations are presented. Direct flux approximation is considered as an important advantage of the LS method. However, there are no known pointwise error estimates for the direct flux approximation. In this paper, we provide optimal pointwise estimates which show local dependence of the error at a point and weak dependence of the global norm. As an elementary consequence of these estimates, we provide an asymptotic error expansion inequality. The inequality has applications to superconvergence and a posteriori estimates.

Description: The first-order and higher-order derivatives of a function can be viewed as the solutions of Volterra integral equations of the first kind. In this paper we propose a fast multiscale solver for the numerical solution of the Tikhonov regularization of the Volterra equations. In association with the special form of the kernels, the matrices resulting from the discretization by multiscale bases are sparse. Moreover, they can be truncated using proper strategies with only a minor loss of accuracy. In the best case, the number of nonzero entries of the truncated matrices is linear with respect to the dimensions of the matrices. The accuracy of the solution from the solver is analysed theoretically and verified by numerical experiments.

Description: The finite element method with $\mathscr {Q}_p$ elements is applied to a singularly perturbed convection–diffusion problem on an L-shaped domain. As an effect of corner singularities the exact solution is not $H^2$ -regular. Therefore, we combine a layer-adapted Shishkin mesh with a special grading adapted to the corner singularity. On such meshes we prove error estimates and estimates for the closeness error which explicitly show the influence of the grading parameter $\mu$ . Hence, $\mu$ can be chosen such that optimal error bounds are obtained. Thereby, it turns out that in the problem studied the influence of the corner singularity becomes small if the perturbation parameter $\varepsilon$ decreases. Moreover, we conduct numerical experiments that verify the theoretical results.

Description: Pták's method of nondiscrete induction is based on the idea that in the analysis of iterative processes one should aim at rates of convergence as functions rather than just numbers, because functions may give convergence estimates that are tight throughout the iteration rather than just asymptotically. In this paper we motivate and prove a theorem on nondiscrete induction, originally due to Potra and Pták, and we apply it to the Newton iterations for computing the matrix polar decomposition and the matrix square root. Our goal is to illustrate the application of the method of nondiscrete induction in the finite-dimensional numerical linear algebra context. We show the sharpness of the resulting convergence estimate analytically for the polar decomposition iteration and on some examples for the square root iteration. We also discuss some of the method's limitations and possible extensions.

Description: We consider the numerical solution, by a Petrov–Galerkin finite-element method, of a singularly perturbed reaction–diffusion differential equation posed on the unit square. In Lin & Stynes (2012, A balanced finite element method for singularly perturbed reaction-diffusion problems. SIAM J. Numer. Anal. , 50 , 2729–2743), it is argued that the natural energy norm, associated with a standard Galerkin approach, is not an appropriate setting for analysing such problems, and there they propose a method for which the natural norm is ‘balanced’. In the style of a first-order system least squares method, we extend the approach of Lin & Stynes (2012, A balanced finite element method for singularly perturbed reaction-diffusion problems. SIAM J. Numer. Anal. , 50 , 2729–2743) by introducing a constraint which simplifies the associated finite-element space and the method's analysis. We prove robust convergence in a balanced norm on a piecewise-uniform (Shishkin) mesh, and present supporting numerical results. Finally, we demonstrate how the resulting linear systems are solved optimally using multigrid methods.

Description: The construction of tensor-product surface patches with a family of Pythagorean-hodograph (PH) isoparametric curves is investigated. The simplest nontrivial instances, interpolating four prescribed patch boundary curves, involve degree $(5,4)$ tensor-product surface patches $\bf{x}(u,v)$ whose $v=\hbox {constant}$ isoparametric curves are all spatial PH quintics. It is shown that the construction can be reduced to solving a novel type of quadratic quaternion equation, in which the quaternion unknown and its conjugate exhibit left and right coefficients, while the quadratic term has a coefficient interposed between them. A closed-form solution for this type of equation is derived, and conditions for the existence of solutions are identified. The surfaces incorporate three residual scalar freedoms which can be exploited to improve the interior shape of the patch. The implementation of the method is illustrated through a selection of computed examples.

Description: Interior eigenvalues of bounded scattering objects can be rigorously characterized from multi-static and multi-frequency far field data, that is, from the behaviour of scattered waves far away from the object. This characterization, the so-called inside–outside duality, holds for various types of penetrable and impenetrable scatterers and is based on the behaviour of a particular eigenvalue of the far field operator. It naturally leads to a numerical algorithm for computing interior eigenvalues of a scatterer that does not require shape or physical properties of the scatterer as input. Since the nonlinear inverse problem to compute such interior eigenvalues from far field data is ill-posed, we propose a regularizing algorithm that is shown to converge as the noise level of the far field data tends to zero. We illustrate feasibility and accuracy of our algorithm by numerical experiments where we compute interior transmission eigenvalues and Robin eigenvalues of the Laplacian in three-dimensional domains from scattering data of these domains due to plane incident waves.

Description: Experimental studies of mantle petrology find that small concentrations of water and carbon dioxide have a large effect on the solidus temperature and distribution of melting in the upper mantle. However, it has remained unclear what effect small fractions of deep, volatile-rich melts have on melt transport and reactive melting in the shallow asthenosphere. Here we present theory and computations indicating that low-degree, reactive, volatile-rich melts cause channelization of magmatic flow at depths approximately corresponding to the anhydrous solidus temperature. These results are obtained with a novel method to simulate the thermochemical evolution of the upper mantle in the presence of volatiles. The method uses a thermodynamically consistent framework for reactive, disequilibrium, multi-component melting. It is coupled with a system of equations representing conservation of mass, momentum, and energy for a partially molten grain aggregate. Application of this method in two-phase, three-component upwelling-column models demonstrates that it reproduces leading-order features of hydrated and carbonated peridotite melting; in particular, it captures the production of low-degree, volatile-rich melt at depths far below the volatile-free solidus. The models predict that segregation of volatile-rich, deep melts promotes a reactive channelling instability that creates fast and chemically isolated pathways of melt extraction. Reactive channelling occurs where volatile-rich melts flux the base of the silicate melting region, enhancing dissolution of fusible components from the ambient mantle. We find this effect to be similarly expressed for models of both hydrated and carbonated mantle melting. These findings indicate that despite their small concentrations, water and carbon dioxide have an important control on the extent and style of magma genesis, as well as on the dynamics of melt transport.

Description: We derive optimal-order a posteriori error estimates for fully discrete approximations of initial and boundary value problems for linear parabolic equations. For the discretization in time we apply the fractional-step $\vartheta $ -scheme, and for the discretization in space the finite element method with finite element spaces that are allowed to change with time. The first optimal-order a posteriori error estimates for the norms of $L^\infty (0,T;L^2(\varOmega ))$ and $L^2(0,T;H^1(\varOmega ))$ are derived by applying the reconstruction technique.

Description: In this paper, we propose a fast and accurate numerical method based on Fourier transforms to solve Kolmogorov forward equations of symmetric scalar Lévy processes. The method is based on the accurate numerical formulas for Fourier transforms proposed by Ooura. These formulas are combined with nonuniform fast Fourier transforms (FFT) and fractional FFT to speed up the numerical computations. Moreover, we propose a formula for numerical indefinite integration on equispaced grids as a component of the method. The proposed integration formula is based on the sinc-Gauss sampling formula, which is a function approximation formula. This integration formula is also combined with the FFT. Therefore, all steps of the proposed method are executed using the FFT and its variants. The proposed method allows us to be free from some special treatments for a nonsmooth initial condition and numerical time integration. The numerical solutions obtained by the proposed method appear to be exponentially convergent on the interval if the corresponding exact solutions do not have sharp cusps. Furthermore, the real computational times are approximately consistent with the theoretical estimates.

Description: In this paper, we consider the heat equation coupled with Darcy's law with a nonlinear source term describing heat production due to an exothermic chemical reaction. The existence and uniqueness of a solution are established. Next, a spectral discretization of the problem is presented and thoroughly analysed. Finally, we present some numerical experiments which confirm the interest of the discretization.

Description: The convergence of the so-called quadratic method for computing eigenvalue enclosures of general self-adjoint operators is examined. Explicit asymptotic bounds for convergence to isolated eigenvalues are found. These bounds turn out to improve significantly upon those determined in previous investigations. The theory is illustrated by means of several numerical experiments performed on particularly simple benchmark models of one-dimensional Schrödinger operators.

Description: A study of whole-rock major and trace element and Sr–Nd–Pb–Hf isotope compositions, combined with zircon U–Pb ages and Hf–O isotopes, for postcollisional intermediate volcanic rocks from the Dabie orogen, China provides constraints on the origin of andesitic magmas and insights into slab–mantle interaction in continental subduction channels. The volcanic rocks exhibit variable contents of SiO 2 (50·28–63·86 wt %), MgO (1·18–4·65 wt %), (Fe 2 O 3 ) T (3·60–8·53 wt %), Al 2 O 3 (12·92–18·95 wt %), Na 2 O (2·08–6·30 wt %) and K 2 O (0·73–5·25 wt %). They are mainly trachyandesites, with lesser amounts of basaltic trachyandesite, andesite, dacite and trachyte, characterized by subduction-like trace element distribution patterns showing enrichment of large ion lithophile elements and light rare earth elements but depletion of high field strength elements. The volcanic rocks have relatively enriched Sr–Nd–Pb–Hf isotope compositions, with high initial 87 Sr/ 86 Sr ratios of 0·7075–0·7110, highly negative Nd ( t ) values of –23·1 to –15·0, Hf ( t ) values of –29·8 to –18·3 and elevated 207 Pb/ 204 Pb and 208 Pb/ 204 Pb ratios at given 206 Pb/ 204 Pb ratios. Zircon Hf–O isotope analyses yield negative Hf ( t ) values of –31·0 to –17·8 and 18 O values of 4·4–6·8 for syn-magmatic domains. Zircon U–Pb dating yields consistent Early Cretaceous ages of 124 ± 3 to 130 ± 2 Ma for magma emplacement. Residual zircon cores yield Triassic, Neoproterozoic and Paleoproterozoic U–Pb ages, consistent with the ages of tectonothermal events involving ultrahigh-pressure metamorphism and protolith formation in the Dabie orogen. The zircon cores also yield a few low to negative 18 O values. An integrated interpretation of these geochemical characteristics is that the andesitic magmas were derived by partial melting of metasomatized zones in the orogenic lithospheric mantle. The metasomatites were generated by reaction of subcontinental lithospheric mantle wedge peridotite with felsic melts that originated from deeply subducted continental crust during continental collision in the Triassic. Melt–peridotite reaction in a subduction channel is therefore a key to the origin of the mantle sources of andesitic magmas in collisional orogens.

Description: The late Proterozoic Ntaka Ultramafic Complex is a body of dominantly pyroxenitic cumulate rocks containing cyclic alternations of olivine–orthopyroxene cumulates. Chemical zoning in the pyroxenes has been imaged at 25–40 µm resolution using desktop microbeam X-ray fluorescence mapping followed up with laser ablation–inductively coupled plasma mass spectrometry analysis for minor and trace elements on selected samples. Poikilitic and granular harzburgites are finely intermingled, in some cases on a centimetre scale in the same thin section. Poikilitic varieties display spectacular textures, ranging from isolated equant orthopyroxene oikocrysts within olivine-rich heteradcumulate harzburgites to rocks composed entirely of interlocking centimetre-sized anhedral orthopyroxene oikocrysts containing sharply bounded idiomorphic Cr-enriched cores. The poikilitic harzburgites are interlayered with cumulate pyroxenites in which orthopyroxene grains show a variety of zoning patterns: Cr-rich cores similar to those in the oikocrysts; sharply bounded oscillatory zoned cores; and reverse zoning with Cr-poor cores and Cr-enriched rims. A further variation is the presence of a mingled harzburgite lithology in which dunite or poikilitic harzburgite is invaded on a centimetre scale by diffuse vein networks or patches of coarse orthopyroxenite. This range of textures and lithologies attests to a more complex set of processes than implied by the standard cumulus theory model in which oikocrysts are considered to have crystallized from intercumulus liquid within a permeable crystal mush. A range of hypotheses is proposed, including infiltration metasomatism of original olivine cumulates by migrating orthopyroxene-saturated pore fluid; however, the textural relationships, whole-rock chemistry and Cr zoning within the grains can best be explained by a model in which the orthopyroxene oikocrysts form in part or whole as mechanically accumulated cumulus grains. The complexity of zoning patterns is attributed to stirring of entrained olivine and orthopyroxene crystals within a heterogeneous flowing crystal mush, where the transporting magma has a wide range of silica contents owing to poorly stirred incorporation of siliceous country-rock material. The Cr-rich orthopyroxenite component grew from Si-enriched chromite-saturated magma. Mingled lithologies developed after accumulation as a result of percolation and infiltration metasomatism by Si-enriched liquid derived by melting of xenoliths within the crystal pile. The model may be more generally applicable: dunite–harzburgite cycles, common in many layered intrusions, may reflect variable degrees of contamination rather than cycles of fractional crystallization and replenishment.

Description: New thermodynamic data for skiagite garnet (Fe 3 Fe 2 3+ Si 3 O 12 ) are derived from experimental phase-equilibrium data that extend to 10 GPa and are applied to oxybarometry of mantle peridotites using a revised six-component garnet mixing model. Skiagite is more stable by 12 kJ mol –1 than in a previous calibration of the equilibrium 2 skiagite = 4 fayalite + ferrosilite + O 2 , and this leads to calculated oxygen fugacities that are higher (more oxidized) by around 1–1·5 logfO2 units. A new calculation method and computer program incorporates four independent oxybarometers (including 2 pyrope + 2 andradite + 2 ferrosilite = 2 grossular + 4 fayalite + 3 enstatite + O 2 ) for use on natural peridotite samples to yield optimum logfO2 estimates by the method of least squares. These estimates should be more robust than those based on any single barometer and allow assessment of possible disequilibrium in assemblages. A new set of independent oxybarometers for spinel-bearing peridotites is also presented here, including a new reaction 2 magnetite + 3 enstatite = 3 fayalite + 3 forsterite + O 2 . These recalibrations combined with internally consistent PT determinations for published analyses of mantle peridotites with analysed Fe 2 O 3 data for garnets, from both cratonic (Kaapvaal, Siberia and Slave) and circumcratonic (Baikal Rift) regions, provide revised estimates of oxidation state in the lithospheric mantle. Estimates of logfO2 for spinel assemblages are more reduced than those based on earlier calibrations, whereas garnet-bearing assemblages are more oxidized. Importantly, this lessens considerably the difference between garnet and spinel oxybarometry that was observed with previous published calibrations.

Description: The Kidnappers [~1200 km 3 dense rock equivalent (DRE)] and Rocky Hill (~200 km 3 DRE) caldera-forming events in the Taupo Volcanic Zone, New Zealand, were erupted in close succession from the Mangakino volcanic centre. They have identical radiometric ages at ~1 Ma, yet erosion along the contact between the two deposits suggests that some years to decades separated the two eruptions. Field constraints and the similarities of crystal textures and compositions and glass chemistries of both eruption deposits demonstrate that they came from one overall magmatic system with a common crystal mush source. However, second-order variations in these parameters confirm that the Kidnappers and Rocky Hill deposits represent distinct events and are not the products of a single zoned magma chamber. The systematically zoned Kidnappers fall deposits provide evidence for the tapping of three discrete magma bodies, whereas the succeeding Kidnappers ignimbrite is compositionally more diverse. The transition from fall to flow deposition marks a change in the style of caldera collapse and the simultaneous evacuation of discrete but compositionally diverse melts, each of which underwent a distinct evolution and was held at slightly different P–T conditions prior to eruption. Contrasting plagioclase and orthopyroxene zonation patterns are present in pumices originating from three discrete magma bodies. Less evolved mafic melts interacted with the system, which mobilized portions of the final erupted melt through heating and volatile or chemical exchange in the mush. The two largest Kidnappers melt-dominant bodies were re-tapped in modified form, or re-established from their common mush source, prior to the Rocky Hill event. Rocky Hill pumices contain common, fluid-affected antecrystic crystal clots derived from chamber wall material. Amphibole compositions from each eruption reflect melt evolution processes and, in particular, the contemporaneous crystallization of biotite and breakdown of orthopyroxene. Plagioclase and orthopyroxene from Rocky Hill pumices share common zonation patterns with those from the two largest magma bodies in the Kidnappers. The rapid production of new melt-dominant bodies and the triggering of the Rocky Hill eruption reflect the ability of the magmatic system to rejuvenate on a geologically short timescale. The Mangakino centre did not follow a typical cycle of decreased activity after the supervolcanic Kidnappers event, instead producing a second caldera-forming eruption, within years to decades from the same system.

Description: Ambrym, a basaltic volcano in the Vanuatu Arc, has displayed variable eruptive behaviour throughout the past century, with major eruptions occurring both on the volcano flanks and at multiple vents within its caldera. These have been interspersed with periods of relative quiescence marked by extensive passive degassing at active, intra-caldera lava lakes, which experience occasional Strombolian explosions. Volcanic rocks from all vents and eruptive styles display similar isotope and incompatible trace element compositions, suggesting that all are derived from the same primary melt by fractional crystallization. Major eruptions are commonly responsible for effusion of the least evolved lavas examined (SiO 2 ~ 50 wt %; MgO ~ 5 wt %). Although all are geochemically similar, petrological differences discriminate between lavas erupted during flank and intra-caldera eruptions. Phyric basalts with homogeneous mineral compositions are common to flank eruptions, whereas crystal-rich basalts with variable mineral compositions, many not in equilibrium with their host liquid, are a feature of intra-caldera lavas. Lava lake samples are slightly more evolved than those from effusive eruptions (SiO 2 ~ 51–52 wt %; MgO ~ 4 wt %), as a result of additional crystallization during periods of relative quiescence. The diverse petrology of the intra-caldera lavas can be explained by mixing of replenishment magmas similar to those erupted from the volcano flanks with residual magma from lava lake activity. Flank eruptions exploit dykes that bypass the shallow reservoir involved with lava lake activity, limiting their interaction with this component.

Description: The generation and evolution of basaltic magmas at Usu volcano, located at the junction between the NE Japan arc and the Kuril arc, have been investigated. The mafic products, which form the somma edifice of the volcano, consist of basalt (49·6–51·3 wt % SiO 2 ) and basaltic andesite (52·0–54·9 wt % SiO 2 ) lavas. The basaltic lavas show relatively tight compositional trends, and 87 Sr/ 86 Sr ratios tend to decrease with increasing whole-rock SiO 2 content. The water content of the basaltic magmas was determined to be ~4·8 wt % based on plagioclase–melt thermodynamic equilibrium. Using this information and an olivine maximum fractionation model, the water content of the primary Usu magma was estimated to be 3·9 wt %. Multi-component thermodynamic calculations suggest that the primary magma was generated by ~23% melting of the source mantle with ~0·94 wt % H 2 O at ~1300°C and ~1·4 GPa. The 0·94 wt % water content of the source mantle is significantly higher than that beneath volcanoes in the main NE Japan arc (generally 〈0·7 wt % H 2 O); this implies that the wedge mantle at the arc–arc junction is intensively hydrated. The temperature of the wedge mantle of ~1300°C at ~1·4 GPa is also significantly higher than that of the mantle in the main NE Japan arc. Unlike the basaltic lavas, the whole-rock compositions of the basaltic andesite lavas are scattered in Harker variation diagrams. This observation suggests that the compositional diversity was produced by at least two independent processes. To elucidate the processes responsible for this compositional diversity, principal component analysis was applied to the major element compositions of the samples. This suggests that 47% of the diversity of the whole-rock compositions can be explained by mixing with partial melts of lower crustal materials, 25% is explained by redistribution of plagioclase phenocrysts, and 16% is explained by fractionation of accessory minerals.

Description: Numerous models have been developed to simulate the reaction of magmas to changes of thermodynamic variables, such as pressure, temperature, oxygen fugacity, and water activity. However, the extensive experimental database still lacks information on the distinct effect of small amounts of H 2 O on olivine + plagioclase + clinopyroxene cotectic crystallization in tholeiitic basalt. We present an experimental study addressing the effects of pressure (at 100, 200, 400, and 700 MPa) and small amounts of H 2 O on phase relations and liquid lines of descent in three tholeiitic basalts representing different evolutionary stages of the Shatsky Rise oceanic plateau magmatic system (compositions AH6, AH3, and AH5 with 8·6, 8·0, and 6·4 wt % MgO, respectively). Two experimental approaches (dry and low H 2 O) are designed to maintain contrasting H 2 O activities during crystallization using (1) graphite–platinum double capsules to perform nearly anhydrous experiments (〈0·15 wt % H 2 O in the melt) and (2) Fe pre-saturated Au 20 Pd 80 capsules to obtain low melt H 2 O contents ranging from 0·4 to 1·1 wt % H 2 O. Under dry conditions, at lower pressures (≤400 MPa), the crystallization in the MgO-rich AH6 and intermediate AH3 basalts follows the typical sequence of tholeiitic differentiation with olivine crystallization at the liquidus followed by olivine + plagioclase and olivine + plagioclase + clinopyroxene. Both basalts are close to multiple saturation at pressures between 400 and 700 MPa. At high pressure (700 MPa) the crystallization sequence is reversed, starting with clinopyroxene at the liquidus. Under low-H 2 O conditions, AH6 and AH3 are very close to multiple saturation, even at the low pressures of 100 and 200 MPa, and the reversed crystallization sequence (clinopyroxene, plagioclase + clinopyroxene, olivine + plagioclase + clinopyroxene) is observed already at 400 MPa. In contrast to the two more MgO-rich basalts, in the most evolved AH5 basalt, clinopyroxene is the liquidus phase at all investigated pressures and under both dry and low-H 2 O conditions, followed by crystallization of plagioclase + clinopyroxene and olivine + plagioclase + clinopyroxene. The most striking observation in our experiments is that the stability of clinopyroxene increases not only with pressure increase but also in the presence of small amounts of H 2 O (when compared with dry counterparts at similar pressures). Small amounts of H 2 O increase the proportion of clinopyroxene in the olivine + plagioclase + clinopyroxene phase assemblage. Our experiments clearly show that the effect of adding 0·4 wt % H 2 O to cotectic melt compositions (e.g. CaO/Al 2 O 3 ratio at a given MgO) is similar to that caused by an increase of pressure from 100 to ~ 300 MPa. This implies that small amounts of H 2 O can lead to significant overestimation of cotectic crystallization pressures (by up to 300 MPa) and that H 2 O contents need to be taken into account in geobarometric models. Our new experiments emphasize the role of low melt H 2 O contents in stabilizing clinopyroxene and provide some new insights into the problem of the ‘pyroxene paradox’. The apparent mantle pressures obtained for some mid-ocean ridge basalts using ‘dry’ geobarometric approaches can actually represent depths within the lower crust, if small amounts of H 2 O are present. The application of our experimental data to natural Shatsky Rise basalts implies that the magmas record partial crystallization processes occurring mainly at low pressure (100 MPa), corresponding to depths of ~3 km beneath the former spreading center, although the more primitive lavas show evidence of differentiation in a deeper reservoir at ~14 km depth (400 MPa).

Description: Mantle-derived xenoliths hosted by melilitite lavas from In Teria (Ahaggar, SE Algeria) include garnet and spinel peridotites, pyroxenite and phlogopite megacrysts. The spinel and garnet peridotites record an early deformation event, which formed porphyroclastic microstructures and olivine crystal preferred orientations, followed by static infiltration of hydrous alkaline melts. This metasomatic stage (stage 1) is characterized by the crystallization of phlogopite in the garnet and spinel peridotites, amphibole in the spinel peridotites and clinopyroxene in the garnet peridotite, which record chemical equilibration with an alkaline silicate melt. These early events were largely overprinted by carbonatitic metasomatism (stage 2), which is observed only in the spinel peridotites. Spinel peridotite major and trace element compositions, as well as the compositions of newly formed minerals, are characteristic of interaction with carbonate melt, associated with strong enrichment in incompatible trace elements in clinopyroxene. This second stage was followed by crystallization of pyroxenites (stage 3) in vein conduits, probably segregated from alkaline melts. We propose a scenario in which the different metasomatic imprints record successive stages of interaction between lithospheric mantle and sublithospheric melts throughout the Cenozoic. In Sr–Nd isotope space, the host melilitites and several xenoliths are clustered and plot close to the HIMU mantle end-member. However, some peridotite xenoliths are shifted towards more radiogenic 87 Sr/ 86 Sr values. In 207 Pb/ 204 Pb– 206 Pb/ 204 Pb and 208 Pb/ 204 Pb– 06 Pb/ 204 Pb space the In Teria samples define a relatively large domain characterized by high 206 Pb/ 204 Pb and 208 Pb/ 204 Pb, consistent with a contribution of an HIMU component, considered to represent a sublithospheric signature. The highest 87 Sr/ 86 Sr values are comparable with those ascribed to the EM1 mantle end-member, representing the signature of the lower continental lithosphere, and are probably inherited from the pre-metasomatic lithospheric mantle beneath In Teria. Numerical modelling of porous percolation of melt of sublithospheric origin through an EM1-like lithospheric mantle protolith reproduces the In Teria peridotite compositions, using moderately sub-chondritic Sr/Nd values for the peridotite (e.g. In Teria garnet peridotite) and moderately super-chondritic Sr/Nd values in the melt (approximately ocean island basalt values). A few spinel peridotites require a component characterized by a 143 Nd/ 144 Nd signature higher than both the EM1 end-member and the local Ahaggar basalts; the 208 Pb/ 204 Pb compositions of several samples point to a component with a depleted mid-ocean ridge basalt (MORB) mantle (DMM) signature. Thus the lithospheric mantle beneath In Teria probably did not have a uniform EM1 signature before the onset of metasomatism; it included a DMM peridotite component as well as some peridotites with elevated 143 Nd/ 144 Nd values recording long-term LREE depletion.

Description: Modelled primary magma compositions of Palaeogene basalts from the North Atlantic Igneous Province (NAIP) require melting at mantle potential temperatures ( T P ) in the range 1480–1550°C. Modern lavas from Icelandic rift zones require T P ~ 1500°C and those from the rift flanks T P ~ 1450°C. Secular cooling of the NAIP thermal anomaly was therefore of the order of ~50°C over the past 61 Myr. There were systematic variations in T P of 50–100°C from the centre of the thermal anomaly to its margins at any one time, although limits on the stratigraphical distribution of T P determinations do not rule out thermal pulsing on a timescale of millions of years. Variation in extent of melting at similar T P was controlled by local variability in lithospheric thickness. In the west of the NAIP, lithosphere thickness varied from ~90 km at Disko Island to ~65 km at Baffin Island, with similar thickness variations being evident for magmatism in the Faroe Islands, Faroe–Shetland Basin and the British Palaeogene Igneous Province (BPIP). Mean pressure of melting was greater than or equal to the final pressure of melting; the two values converge for melting columns with a melting interval of 〈1·5 GPa, regardless of T P . The majority of the BPIP magmas were generated in the garnet–spinel transition in the upper mantle. Calculated and observed rare earth element distributions in NAIP lavas are entirely consistent with the melting regimes derived from major element melting models. This allows a calibration of rare earth element fractionation and melting conditions that can be applied to other flood basalt provinces.

Description: We consider an elliptic optimal control problem where the objective functional contains evaluations of the state at a finite number of points. In particular, we use a fidelity term that encourages the state to take certain values at these points, which means that our problem is related to ones with state constraints at points. The analysis and numerical analysis differs from when the fidelity is in the $L^2$ -norm because we need the state space to embed into the space of continuous functions. In this paper, we discretize the problem using two different piecewise linear finite element methods. For each discretization we use two different approaches to prove a priori $L^2$ -error estimates for the control. We discuss the differences between these methods and approaches, and present numerical results that agree with our analytical results.

Description: A numerical scheme for the approximation of large vibration deformations of inextensible elastic curves is devised. Its unconditional stability and convergence under a regularity assumption on the velocities are demonstrated and numerical experiments are provided.

Description: We present a new scaleable algorithm for approximating the $H_{\infty }$ norm, an important robust stability measure for linear dynamical systems with input and output. Our spectral-value-set-based method uses a novel hybrid expansion–contraction scheme that, under reasonable assumptions, is guaranteed to converge to a stationary point of the optimization problem defining the $H_{\infty }$ norm, and, in practice, typically returns local or global maximizers. We prove that the hybrid expansion–contraction method has a quadratic rate of convergence that is also confirmed in practice. In comprehensive numerical experiments, we show that our new method is not only robust but exceptionally fast, successfully completing a large-scale test set 25 times faster than an earlier method by Guglielmi, Gürbüzbalaban & Overton (2013, SIAM J. Matrix Anal. Appl. , 34 , 709–737), which occasionally breaks down far from a stationary point of the underlying optimization problem.

Description: The first part of this paper enfolds a medius analysis for mixed finite element methods (FEMs) and proves a best-approximation result in $L^2$ for the stress variable independent of the error of the Lagrange multiplier under stability, compatibility and efficiency conditions. The second part applies the general result to the FEM of Arnold and Winther for linear elasticity: the stress error in $L^2$ is controlled by the $L^2$ best-approximation error of the true stress by any discrete function plus data oscillations. The analysis is valid without any extra regularity assumptions on the exact solution and also covers coarse meshes and Neumann boundary conditions. Further applications include Raviart–Thomas finite elements for the Poisson and the Stokes problems. The result has consequences for nonlinear approximation classes related to adaptive mixed FEMs.

Description: We propose a numerical method to approximate the solution of a nonlocal diffusion problem on a general setting of metric measure spaces. These spaces include, but are not limited to, fractals, manifolds and Euclidean domains. We obtain error estimates in $L^\infty (L^p)$ for $p=1,\infty $ under the sole assumption of the initial datum being in $L^p$ . An improved bound for the error in $L^\infty (L^1)$ is obtained when the initial datum is in $L^2$ . We also derive some qualitative properties of the solutions like stability, comparison principles and study the asymptotic behaviour as $t\to \infty $ . We finally present two examples on fractals: the Sierpinski gasket and the Sierpinski carpet, which illustrate on the effect of nonlocal diffusion for piece-wise constant initial datum.

Description: This paper presents a Riemannian trust region algorithm for unconstrained optimization problems with locally Lipschitz objective functions defined on complete Riemannian manifolds. To this end we define a function $\Phi :TM\rightarrow \mathbb {R}$ on the tangent bundle $TM$ , and at the $k$ th iteration, using the restricted function $\Phi |_{T_{x_k}M}$ , where $T_{x_k}M$ is the tangent space at $x_k$ , a local model function $Q_k$ that carries both first- and second-order information for the locally Lipschitz objective function $f:M\rightarrow \mathbb {R}$ on a Riemannian manifold $M$ , is defined and minimized over a trust region. We establish the global convergence of the proposed algorithm. Moreover, using the Riemannian $\varepsilon $ -subdifferential, a suitable model function is defined. Numerical experiments illustrate our results.

Description: Variational inequalities play an important role in many applications and are an active research area. Optimal a priori error estimates in the natural energy norm do exist, but only very few results are known for different norms. Here, we consider as prototype a simple Signorini problem, and provide new optimal order a priori error estimates for the trace and the flux on the Signorini boundary. The a priori analysis is based on a continuous and a discrete Steklov–Poincaré operator, as well as on Aubin–Nitsche-type duality arguments. Numerical results illustrate the convergence rates of the finite-element approach.

Description: We show for adaptive triangulations in two dimensions, which are generated by the newest vertex bisection, an optimal grading estimate. Roughly speaking, we construct from the piecewise constant mesh-size function a regularized one with the following two properties. First, the two functions are equivalent, and second, the regularized mesh-size function differs at most by a factor of 2 on neighbouring elements. In combination with Bank & Yserentant (2014, Numer. Math. 126 , 361–381), this optimal grading estimate enables us to show that the $L_2$ -orthogonal projections onto the space of continuous Lagrange finite elements up to order 12 is $H^1$ -stable. We extend these results to a modified red–green refinement.

Description: There exist at least three first-order finite volume element methods with totally different dual meshes and conforming or nonconforming $P_1$ finite element trial functions and the question arises of whether they are comparable. The fact that the underlying norms are very different does not prevent the proof that the errors are equivalent on any mesh in some norm. This equivalence is independent of the regularity of the exact solution and holds for any coarse or fine mesh with or without local mesh refining, but up to equivalence constants and additional explicit data-oscillation terms of higher order. The equivalence constants depend on the minimal angle of the shape-regular triangulation and the penalization parameter of the discontinuous Galerkin scheme. This also implies quasi-optimality in the sense that the error is bounded by the best approximation of the flux by piecewise constants. An a posteriori error analysis for the discontinuous Galerkin finite volume element scheme is also discussed. The analysis is exemplified for a boundary value model problem for some second-order elliptic partial differential equation in two dimensions.

Description: The dihedral angle formed at junctions between two plagioclase grains and a grain of augite is only very rarely in textural equilibrium in gabbros from kilometre-scale crustal layered intrusions. The median of a population of these disequilibrium angles, cpp , varies systematically within a single layered intrusion, remaining constant over large stretches of stratigraphy with significant increases and decreases associated with the addition or reduction respectively of the number of phases on the liquidus of the bulk magma. The stepwise changes in cpp are present in the Upper Zone of the Bushveld Complex, the Megacyclic Unit I of the Sept Iles Intrusion, and the Layered Series of the Skaergaard intrusion. The plagioclase-bearing cumulates of Rum have a bimodal distribution of cpp , dependent on whether the cumulus assemblage includes clinopyroxene. The presence of the stepwise changes is independent of the order of arrival of cumulus phases and of the composition of either the cumulus phases or the inferred composition of the interstitial liquid. The only parameter that behaves in an exactly analogous manner to cpp is the rate of change in enthalpy with temperature ( H / T ) during crystallization. Both H / T and cpp increase with the addition of a liquidus phase, and decrease with the removal of a liquidus phase. The replacement of one phase by another has little effect on H / T and no discernible effect on cpp . An increase of H / T results in an increase in the fraction of the total enthalpy budget that is the latent heat of crystallization (the fractional latent heat). It also increases the mass crystallized in each incremental temperature drop (the crystal productivity). These increases of both fractional latent heat and crystal productivity are likely to cause an increase in the time taken to form three-grain junctions in the mush via thermal buffering of a thickened mushy layer. We suggest these are the underlying causes of stepwise increases in cpp . Stepwise changes in the geometry of three-grain junctions in fully solidified gabbros thus provide a clear microstructural marker for the progress of fractionation down the liquid line of descent in layered intrusions.

Description: High-resolution sampling in monogenetic fields has the potential to reveal fine-scale heterogeneity of the mantle, a feature that may be overwhelmed by larger fluxes of magma, or missed by under-sampling. The Quaternary Auckland Volcanic Field (AVF) in northern New Zealand is a basaltic field of 51 small-volume volcanic centres, and is one of the best-sampled examples of a monogenetic volcanic field. We present data for 12 centres in the volcanic field. These show the large compositional variations between volcanoes as well as through single eruptive sequences. Whole-rock compositions range from subalkaline basalt in the larger centres, through alkali basalt to nephelinite in the smallest centres. Fractional crystallization has had a limited effect in many of the centres, but high-pressure clinopyroxene crystallization may have occurred in others. Three end-members are observed in Pb isotope space, indicating that distinct mantle source components are involved in the petrogenesis of the magmas. Whole-rock multi-element patterns show that the larger centres have prominent positive Sr anomalies and lack K anomalies, whereas the smaller centres have prominent negative K anomalies and lack Sr anomalies. The melting parameters and compositions of the sources involved are modelled using trace element ratios and multi-element patterns, and three components are characterized: (1) fertile peridotite with a Pb-isotope composition similar to Pacific mid-ocean ridge basalt; (2) eclogite domains with a HIMU-like isotope composition dispersed within the fertile peridotite; (3) slightly depleted subduction-metasomatized peridotitic lithospheric mantle (containing c . 3% subduction fluids). Modelling shows that melting in the AVF begins in garnet-bearing fertile asthenosphere (with preferential melting of eclogite domains) and that melts are variably diluted by melts of the lithospheric source. The U–Th isotope compositions of the end-members in the AVF show 230 Th excess [( 230 Th/ 232 Th) ratios of 1·11–1·38], with the samples of lower ( 230 Th/ 232 Th) exhibiting higher ( 238 U/ 232 Th), which we attribute to the dilution effect of the melts from the lithospheric mantle source. Modelling reveals a correlation between melting in the asthenosphere, the degree of melting and incorporation of the metasomatized lithospheric mantle source, and the resultant size of the volcanic centre. This suggests that the scale of the eruption may essentially be controlled by asthenospheric mantle dynamics.

Description: The origin of mafic and ultramafic sills exhibiting different whole-rock compositional profiles (e.g. I-, C-, D-, M- and S-shaped profiles) remains controversial. We have addressed this issue by revisiting three ~100 m thick Siberian dolerite sills (Vavukansky, Kuz’movsky and Vilyuysky) that display remarkable internal differentiation. The Vavukansky sill has an M-shaped profile with prominent basal and top reversals showing inward increases in whole-rock MgO, Mg-number [100Mg/(Mg + Fe)] and normative An content [100An/(An + Ab)], followed by the Layered and Upper Border Series with inward decreases in these indices. The Kuz’movsky and Vilyuysky sills both show S-shaped profiles similar to the Vavukansky sill, but lack a top reversal. These whole-rock M- and S-shaped profiles are accompanied by similar profiles in mineral compositions. Plagioclase and, to a lesser extent, olivine show systematic inward increases in An content and Mg-number, respectively, across basal and top reversals. These compositional trends are followed by inward decreases in these ratios in the interiors of the Vavukansky and Kuz’movsky sills. Currently accepted models attribute whole-rock M- and S-shaped compositional profiles to crystal settling, compositional convection or compaction operating in closed systems. Our observations challenge these traditional interpretations because variations in mineral compositions observed in marginal reversals cannot result from closed-system fractionation. We suggest instead that initially the sills evolved as open systems that were slowly inflated by magmas that became gradually more primitive with time. The inflation was accompanied by in situ crystallization that preserved the preceding fractionation history of the injected magmas by forming basal and top reversals with minerals becoming more primitive inwards. This process culminated with rapid inflation of the sills to their current size owing to a major influx of primitive magma. Subsequently, magma flow through the sills ceased and they evolved as closed systems by fractional crystallization. This resulted in the Layered and Upper Border Series with minerals becoming more evolved inwards. This model can be extended to explain other compositional profiles and petrological features in mafic and ultramafic sills.

Description: We show the existence and uniqueness of a solution for the nonlocal vector-valued Allen–Cahn variational inequality in a formulation involving Lagrange multipliers for local and nonlocal constraints. Furthermore, we propose and analyse a primal–dual active set (PDAS) method for local and nonlocal vector-valued Allen–Cahn variational inequalities. The local convergence behaviour of the PDAS algorithm is studied by interpreting the approach as a semismooth Newton method and numerical simulations are presented demonstrating its efficiency.

Description: In this paper, we will present a generalized convolution quadrature for solving linear parabolic and hyperbolic evolution equations. The original convolution quadrature method by Lubich works very nicely for equidistant time steps while the generalization of the method and its analysis to nonuniform time stepping is by no means obvious. We will introduce the generalized convolution quadrature allowing for variable time steps and develop a theory for its error analysis. This method opens the door for further development towards adaptive time stepping for evolution equations. As the main application of our new theory, we will consider the wave equation in exterior domains which is formulated as a retarded boundary integral equation.

Description: We address the error control of Galerkin discretization (in space) of linear second-order hyperbolic problems. More specifically, we derive a posteriori error bounds in the L ( L 2 ) norm for finite element methods for the linear wave equation, under minimal regularity assumptions. The theory is developed for both the space-discrete case and for an implicit fully discrete scheme. The derivation of these bounds relies crucially on carefully constructed space and time reconstructions of the discrete numerical solutions, in conjunction with a technique introduced by Baker (1976, Error estimates for finite element methods for second-order hyperbolic equations. SIAM J. Numer. Anal. , 13 , 564–576) in the context of a priori error analysis of Galerkin discretization of the wave problem in weaker-than-energy spatial norms.

Description: The surface finite element method can be used to approximate curvatures on embedded hypersurfaces and to discretize geometric partial differential equations. In this paper, we present a definition of discrete Ricci curvature on polyhedral hypersurfaces of arbitrary dimension based on the discretization of a weak formulation with isoparametric finite elements. We prove that for a piecewise quadratic approximation of a two- or three-dimensional hypersurface R n +1 , this definition approximates the Ricci curvature of with a linear order of convergence in the L 2 ( ) norm. By using a smoothing scheme in the case of a piecewise linear approximation of , we still get a convergence of order 2/3 in the L 2 ( ) norm and of order 1/3 in the W 1, 2 ( ) norm.

Description: We give general conditions which guarantee that the sequence generated by a descent algorithm converges to an equilibrium point. The convergence result is based on the Lojasiewicz gradient inequality; optimal convergence rates are also derived, as well as a stability result. We show how our results apply to a large variety of standard time discretizations of gradient-like flows. Schemes with variable time step are considered and optimal conditions on the maximal step size are derived. Applications to time and space discretizations of the Allen–Cahn equation, the sine–Gordon equation and a damped wave equation are given.

Description: This paper presents quadratic finite-volume methods for elliptic and parabolic problems on quadrilateral meshes that use Barlow points (optimal stress points) for dual partitions. Introducing Barlow points into the finite-volume formulations results in better approximation properties at the cost of loss of symmetry. The novel ‘symmetrization’ technique adopted in this paper allows us to derive optimal-order error estimates in the H 1 - and L 2 -norms for elliptic problems and in the L ( H 1 )- and L ( L 2 )-norms for parabolic problems. Superconvergence of the difference between the gradients of the finite-volume solution and the interpolant can also be derived. Numerical results confirm the proved error estimates.

Description: A linear parabolic differential equation on a moving surface is discretized in space by evolving-surface finite elements and in time by backward difference formulas (BDFs). Using results from Dahlquist's G-stability theory and Nevanlinna & Odeh's multiplier technique together with properties of the spatial semidiscretization, stability of the full discretization is proved for BDF methods up to order 5 and optimal-order convergence is shown. Numerical experiments illustrate the behaviour of the fully discrete method.

Description: In this article, we develop the a priori and a posteriori error analysis of hp -version interior penalty discontinuous Galerkin finite element methods for strongly monotone quasi-Newtonian fluid flows in a bounded Lipschitz domain R d , d = 2, 3. In the latter case, computable upper and lower bounds on the error are derived in terms of a natural energy norm, which are explicit in the local mesh size and local polynomial degree of the approximating finite element method. A series of numerical experiments illustrate the performance of the proposed a posteriori error indicators within an automatic hp -adaptive refinement algorithm.

Description: In recent years, it has become increasingly clear that the critical issue in gradient methods is the choice of the step length, whereas using gradient as the search direction may lead to very effective algorithms, whose surprising behaviour has only been partially explained, mostly in terms of the spectrum of the Hessian matrix. On the other hand, the convergence of the classical Cauchy steepest descent (SD) method has been analysed extensively and related to the spectral properties of the Hessian matrix, but the connection with the spectrum of the Hessian has not been exploited much to modify the method in order to improve its behaviour. In this work, we show how, for convex quadratic problems, moving from some theoretical properties of the SD method, second-order information provided by the step length can be exploited to dramatically improve the usually poor practical behaviour of this method. This allows us to achieve computational results comparable with those of the Barzilai and Borwein algorithm, with the further advantage of monotonic behaviour.

Description: Magma dynamics and time scales during the VEI 5, 2000 bp eruption of El Misti volcano, southern Peru (EM2000BP) are investigated to address cyclic explosive activity at this hazardous volcano. The 1·4 km 3 of pumice falls and flows have abundant mingled pumice of high-K, calc-alkaline rhyolite and andesite composition. Phenocryst zoning and compositions reveal mutual exchange of plagioclase between the two magmas; amphibole in the rhyolite was derived from the andesite. Amphiboles in the andesite are predominantly unrimmed crystals whereas those in the rhyolite mostly exhibit reaction rims. Phase equilibria indicate that the andesite formed at ~900–950°C and 2–3 kbar pressure and was water-saturated with 5·1–6·0 wt % H 2 O, broadly similar to El Misti magmas overall. Amphibole, plagioclase, Ti-magnetite, and two pyroxenes were the crystallizing phases. A separate rhyolite magma existed higher in the crust at a temperature of 816 ± 30°C and ~5% H 2 O in which only plagioclase and Fe–Ti oxides were stable. The lack of cognate amphibole in the rhyolite despite H 2 O saturation requires that it staged above the stability limit of amphibole (〈100 MPa). Exchange reactions in amphibole (dominantly pargasitic) and trace element partitioning in plagioclase indicate that both andesite and rhyolite magmas were broadly constant in temperature and H 2 O content. These constraints suggest that the initially separate rhyolite and deeper andesite magmas interacted by an initial andesite recharge event that resulted in mingling and crystal exchange. A period of 50–60 days is required for amphibole introduced into the rhyolite to develop reaction rims owing to decompression. These rims are dominated by plagioclase, a consequence of the Al-rich nature of the amphibole. The lack of reaction rims on amphibole in the andesite implicates a second, more-forceful and voluminous eruption-triggering recharge event during which andesite rose rapidly from source to surface in ≤5 days at ascent rates of at least 0·023 m s –1 . Further decompression-driven crystallization is recorded in plagioclase rims and microlite growth that may have contributed to a rapid increase in viscosity leading to explosive eruption. This VEI 5 plinian eruption shares characteristics with other explosive events at El Misti on a time scale of 2000–4000 years, suggesting periodic recharge-driven explosive activity.

Description: We report structural evidence of ductile strain localization in mantle pyroxenite from the spinel to plagioclase websterite transition in the Ronda Peridotite (southern Spain). Mapping shows that, in this domain, small-scale shear zones occurring at the base of the lithospheric section are systematically located within thin pyroxenite layers, suggesting that the pyroxenite was locally weaker than the host peridotite. Strain localization is associated with a sudden decrease of grain size and increasing volume fractions of plagioclase and amphibole as a result of a spinel to plagioclase phase transformation reaction during decompression. This reaction also fostered hydrogen extraction (‘dehydroxylation’) from clinopyroxene producing effective fluid saturation that catalyzed the synkinematic net-transfer reaction. This reaction produced fine-grained olivine and plagioclase, allowing the onset of grain-size sensitive creep and further strain localization in these pyroxenite bands. The strain localization in the pyroxenites is thus explained by their more fertile composition, which allowed earlier onset of the phase transition reactions. Geothermobarometry undertaken on compositionally zoned constituent minerals suggests that this positive feedback between reactions and deformation is associated with cooling from at least 1000°C to 700°C and decompression from 1·0 to 0·5 GPa.

Description: In recent years, there has been an enormous interest in developing methods for the approximation of manifold-valued functions. In this paper, we focus on the manifold of symmetric positive-definite (SPD) matrices. We investigate the use of SPD-matrix means to adapt linear positive approximation methods to SPD-matrix-valued functions. Specifically, we adapt corner-cutting subdivision schemes and Bernstein operators. We present the concept of admissible matrix means and study the adapted approximation schemes based on them. Two important cases of admissible matrix means are treated in detail: the exp–log and the geometric matrix means. We derive special properties of the approximation schemes based on these means. The geometric mean is found to be superior in the sense of preserving more properties of the data, such as monotonicity and convexity. Furthermore, we give error bounds for the approximation of univariate SPD-matrix-valued functions by the adapted operators.

Description: We present a mass-preserving scheme for the stochastic nonlinear Schrödinger equation with multiplicative noise of Stratonovich type. It is a splitting scheme and we present an explicit formula for solving the sub-step related to the nonlinear part. The scheme is unconditionally stable in the L 2 norm. For the linear stochastic Schrödinger equation, we prove that the scheme has a strong convergence rate in time equal to 1, which is not common for stochastic partial differential equations with noise depending on space and time.

Description: A numerical scheme for the approximation of the elastic flow of inextensible curves is devised and convergence of approximations to exact solutions of the nonlinear time-dependent partial differential equation is proved. The nonlinear, pointwise constraint of local length preservation is linearized about a previous solution in each time step which leads to a sequence of linear saddle-point problems. The spatial discretization is based on piecewise Bézier curves and the resulting semiimplicit scheme is unconditionally stable and convergent.

Description: We study the coercivity properties and the norm dependence on the wave-number k of certain regularized combined field boundary integral operators that we recently introduced for the solution of two- and three-dimensional acoustic scattering problems with Neumann boundary conditions. We show that in the case of circular and spherical boundaries, our regularized combined field boundary integral operators are L 2 coercive for large enough values of the coupling parameter, and that the norms of these operators are bounded by constant multiples of the coupling parameter. We establish that the norms of the regularized combined field boundary integral operators grow modestly with the wave-number k for smooth boundaries and we provide numerical evidence that these operators are L 2 coercive for two-dimensional starlike boundaries. We present and analyse a fully discrete collocation (Nyström) method for the solution of two-dimensional acoustic scattering problems with Neumann boundary conditions based on regularized combined field integral equations. In particular, for analytic boundaries and boundary data, we establish pointwise superalgebraic convergence rates of the discrete solutions.

Description: The classical theory of Gaussian quadrature assumes a positive weight function. We will show that in some cases Gaussian rules can be constructed with respect to an oscillatory weight, yielding methods with complex quadrature nodes and positive weights. These rules are well suited to highly oscillatory integrals because they attain optimal asymptotic order. We show that, for the Fourier oscillator, this approach yields the numerical method of steepest descent, a method with optimal asymptotic order that has previously been proposed for this class of integrals. However, the approach readily extends to more general kernels, such as Bessel functions that appear as the kernel of the Hankel transform.

Description: We consider anisotropic Allen–Cahn equations with interfacial energy induced by an anisotropic surface energy density . Assuming that is positive, positively homogeneous of degree 1, strictly convex in tangential directions to the unit sphere and sufficiently smooth, we show the stability of various time discretizations. In particular, we consider a fully implicit and a linearized time discretization of the interfacial energy combined with implicit and semiimplicit time discretizations of the double-well potential. In the semiimplicit variant, concave terms are taken explicitly. The arising discrete spatial problems are solved by globally convergent truncated nonsmooth Newton multigrid methods. Numerical experiments show the accuracy of the different discretizations. We also illustrate that pinch-off under anisotropic mean curvature flow is no longer invariant under rotation of the initial configuration for a fixed orientation of the anisotropy.

Description: In this paper, we define a new finite element method for numerically approximating the solution of a partial differential equation in a bulk region coupled with a surface partial differential equation posed on the boundary of the bulk domain. The key idea is to take a polyhedral approximation of the bulk region consisting of a union of simplices, and to use piecewise polynomial boundary faces as an approximation of the surface. Two finite element spaces are defined, one in the bulk region and one on the surface, by taking the set of all continuous functions which are also piecewise polynomial on each bulk simplex or boundary face. We study this method in the context of a model elliptic problem; in particular, we look at well-posedness of the system using a variational formulation, derive perturbation estimates arising from domain approximation and apply these to find the optimal-order error estimates. A numerical experiment is described which demonstrates the order of convergence.

Description: As a model of more general contour integration problems we consider the numerical calculation of high-order derivatives of holomorphic functions using Cauchy's integral formula. Bornemann (2011, Accuracy and stability of computing high-order derivatives of analytic functions by Cauchy integrals. Found. Comput. Math. , 11 , 1–63) showed that the condition number of the Cauchy integral strongly depends on the chosen contour and solved the problem of minimizing the condition number for circular contours. In this paper, we minimize the condition number within the class of grid paths of step size h using Provan's algorithm for finding a shortest enclosing walk in weighted graphs embedded in the plane. Numerical examples show that optimal grid paths yield small condition numbers even in those cases where circular contours are known to be of limited use, such as for functions with branch-cut singularities.

Description: Anisotropic meshes are important for efficiently resolving incompressible flow problems that include boundary layer or corner singularity phenomena. Unfortunately, the stability of standard inf–sup stable mixed approximation methods is prone to degeneracy whenever the mesh aspect ratio becomes large. As an alternative, a stabilized mixed approximation method is considered here. Specifically, a robust a priori error estimate for the local jump stabilized Q 1 – P 0 approximation introduced by Kechkar & Silvester (1992, Analysis of locally stabilized mixed finite element methods for the Stokes problem. Math. Comp. , 58 , 1–10) is established for anisotropic meshes. Our numerical results demonstrate that the stabilized Q 1 – P 0 method is competitive with the nonconforming, nonparametric, rotated approximation method introduced by Rannacher & Turek (1992, Simple nonconforming quadrilateral Stokes element. Numer. Meth. Partial Differential Equations , 8 , 97–111).