Early Last Interglacial Greenland Ice Sheet melting and a sustained period of meridional overturning weakening: a model analysis of the uncertainties

Proxy-data suggest that the Last Interglacial (LIG; ~130–116 ka BP) climate was characterized by higher temperatures, a partially melted Greenland Ice Sheet (GIS) and a changed Atlantic meridional overturning circulation (AMOC). Notwithstanding the uncertainties in LIG palaeoclimatic reconstructions, this setting potentially provides an opportunity to evaluate the relation between GIS melt and the AMOC as simulated by climate models. However, first we need to assess the extent to which a causal relation between early LIG GIS melt and the weakened AMOC is plausible. With a series of transient LIG climate simulations with the LOVECLIM earth system model, we quantify the importance of the major known uncertainties involved in early LIG GIS melt scenarios. Based on this we construct a specific scenario that is within the parameter space of uncertainties and show that it is physically consistent that early LIG GIS melting kept the AMOC weakened. Notwithstanding, this scenario is at the extreme end of the parameter space. Assuming that proxy-based reconstructions of early LIG AMOC weakening offer a realistic representation of its past state, this indicates that either (1) the AMOC weakening was caused by other forcings than early LIG GIS melt or (2) the early LIG AMOC was less stable than indicated by our simulations and a small amount of GIS melt was sufficient to keep the AMOC in the weak state of a bi-stable regime. We argue that more intensive research is required because of the high potential of the early LIG to evaluate model performance in relation to the AMOC response to GIS melt.     

A double slip non-coaxial flow rule for viscous-plastic Cosserat materials

We propose a double slip non-coaxial plastic model within the framework of a Cosserat continuum theory. In a Cosserat continuum, a material point possesses the degrees of freedom of an infinitesimal rigid body: two translations and one rotation in 2D. We formulate the plastic model into viscous-plastic constitutive relationships and illustrate the viscous-plastic behaviour of the model by means of numerical solution of a simple shear problem.     

A non-hydrostatic algorithm for free-surface ocean modelling

An original implementation of a non-hydrostatic, free-surface algorithm based on a pressure correction method is proposed for ocean modelling. The free surface is implemented through an explicit scheme combined with a mode-spitting method but the depth-averaged velocity and the position of the free surface are updated at each non-hydrostatic iteration. The vertical momentum equation is also integrated up to the surface enabling a natural and accurate treatment of the surface layer. The consistent specification of the numerical schemes provides balanced transfers of potential and kinetic energy. This algorithm is well-suited for implementation as a non-hydrostatic kernel on originally hydrostatic free-surface ocean models such as Symphonie (http://poc.obs-mip.fr/pages/research_topics/modelling/symphonie/symphonie.htm) for which it has originally been developed.Energy balances associated with the propagation of short surface waves and solitary waves are presented for two dedicated well-documented configurations over closed domains. The buoyancy flux, the work rate of the pressure force together with the power of the advective terms are evaluated and discussed for the generation and the propagation of these two types of waves. The dissipation rate is in particular shown to be several orders of magnitude smaller than the work rates of the hydrostatic and non-hydrostatic pressure forces confirming the necessity for the exchanges of energy to be numerically balanced. The algorithm is subsequently applied to the complex generation of non-linear solitary internal waves by surface tides over Georges Bank, in the Gulf of Maine. The generation and the propagation of the observed non-linear and non-hydrostatic features in this region are correctly reproduced.     

In situ observation of crystal growth in a basalt melt and the development of crystal size distribution in igneous rocks

The speciation of CO₂ in dacite, phonolite, basaltic andesite, and alkali silicate melt was studied by synchrotron infrared spectroscopy in diamond anvil cells to 1,000 °C and more than 200 kbar. Upon compression to 110 kbar at room temperature, a conversion of molecular CO₂ into a metastable carbonate species was observed for dacite and phonolite glass. Upon heating under high pressure, molecular CO₂ re-appeared. Infrared extinction coefficients of both carbonate and molecular CO₂ decrease with temperature. This effect can be quantitatively modeled as the result of a reduced occupancy of the vibrational ground state. In alkali silicate (NBO/t = 0.98) and basaltic andesite (NBO/t = 0.42) melt, only carbonate was detected up to the highest temperatures studied. For dacite (NBO/t = 0.09) and phonolite melts (NBO/t = 0.14), the equilibrium CO₂ + O^2? = CO₃ ^2? in the melt shifts toward CO₂ with increasing temperature, with ln K = ?4.57 (±1.68) + 5.05 (±1.44) 10³ T ^?1 for dacite melt (ΔH = ?42 kJ mol^?1) and ln K = ?6.13 (±2.41) + 7.82 (±2.41) 10³ T ^?1 for phonolite melt (ΔH = ?65 kJ mol^?1), where K is the molar ratio of carbonate over molecular CO₂ and T is temperature in Kelvin. Together with published data from annealing experiments, these results suggest that ΔS and ΔH are linear functions of NBO/t. Based on this relationship, a general model for CO₂ speciation in silicate melts is developed, with ln K = a + b/T, where T is temperature in Kelvin and a = ?2.69 ? 21.38 (NBO/t), b = 1,480 + 38,810 (NBO/t). The model shows that at temperatures around 1,500 °C, even depolymerized melts such as basalt contain appreciable amounts of molecular CO₂, and therefore, the diffusion coefficient of CO₂ is only slightly dependent on composition at such high temperatures. However, at temperatures close to 1,000 °C, the model predicts a much stronger dependence of CO₂ solubility and speciation on melt composition, in accordance with available solubility data.     

New type of tsunami barrier

Tsunami impact on a coast can be reduced by applying a submerged vertical barrier to reflect tsunami before the catastrophic waves are built up near the coast. However, construction of such long walls by conventional submarine technology is difficult. In this paper, the construction of extended submarine walls at a depth of between 50 and 500 m below sea level by a relatively simple and efficient technology is described. The submarine walls can consist either of high-strength steel fences with anchors or two parallel steel fences with distance holders, lowered into the sea and fixed with rocks inserted from top. Alternatively, the barriers could be built from gabions (prefabricated steel net baskets filled with rocks) and lowered into the sea. The space between these tsunami barriers and the coast can be filled with solid material, thus allowing reclamation of new land, or this gap can be used for fish farming. These barriers can contribute to preservation of beaches and natural ecosystems at the coast.     

Hydrologic Variability and Its Control of Phytoplankton Community Structure and Function in Two Shallow,Coastal, Lagoonal Ecosystems: The Neuse and New River Estuaries,North Carolina,USA

Hydrologic conditions, especially changes in freshwater input, play an important, and at times dominant, role in determining the structure and function of phytoplankton communities and resultant water quality of estuaries. This is particularly true for microtidal, shallow water, lagoonal estuaries, where water flushing and residence times show large variations in response to changes in freshwater inputs. In coastal North Carolina, there has been an increase in frequency and intensity of extreme climatic (hydrologic) events over the past 15 years, including eight hurricanes, six tropical storms, and several record droughts; these events are forecast to continue in the foreseeable future. Each of the past storms exhibited unique hydrologic and nutrient loading scenarios for two representative and proximate coastal plain lagoonal estuaries, the Neuse and New River estuaries. In this synthesis, we used a 13-year (1998–2011) data set from the Neuse River Estuary, and more recent 4-year (2007–2011) data set from the nearby New River Estuary to examine the effects of these hydrologic events on phytoplankton community biomass and composition. We focused on the ability of specific taxonomic groups to optimize growth under hydrologically variable conditions, including seasonal wet/dry periods, episodic storms, and droughts. Changes in phytoplankton community composition and biomass were strongly modulated by the amounts, duration, and seasonality of freshwater discharge. In both estuaries, phytoplankton total and specific taxonomic group biomass exhibited a distinctive unimodal response to varying flushing rates resulting from both event-scale (i.e., major storms, hurricanes) and more chronic seasonal changes in freshwater input. However, unlike the net negative growth seen at long flushing times for nano-/microphytoplankton, the pigments specific to picophytoplankton (zeaxanthin) still showed positive net growth due to their competitive advantage under nutrient-limited conditions. Along with considerations of seasonality (temperature regimes), these relationships can be used to predict relative changes in phytoplankton community composition in response to hydrologic events and changes therein. Freshwater inputs and droughts, while not manageable in the short term, must be incorporated in water quality management strategies for these and other estuarine and coastal ecosystems faced with increasing frequencies and intensities of tropical cyclones, flooding, and droughts.     

Anticlockwise pressure–temperature paths record Variscan upper‐plate exhumation: Example from micaschists of the Porto Vecchio region,Corsica

To better understand the evolution of deep‐seated crust of the Variscan orogen in the Sardinia‐Corsica region, we studied garnet‐bearing micaschists which were sampled 3 km east and 15 km northeast of Porto Vecchio, south‐eastern Corsica. After a careful investigation of the textural relations and compositions of minerals, especially of zoned garnet, a P–T path was reconstructed using contoured P–T pseudosections. U–Th–Pb dating of monazite in the micaschists was undertaken with the electron microprobe. The micaschists from both localities were formed along similar anticlockwise P–T paths. The prograde branch of these paths starts at 3 kbar close to 600°C in the P–T field of sillimanite and reaches peak conditions at 7 kbar and 600 (15 km NE of Porto Vecchio) to 630°C (3 km E of Porto Vecchio). The metamorphism at peak P–T conditions happened c. 340 Ma based on low‐Y (<0.65 wt% Y₂O₃) monazite. Ages of monazite with high‐Y contents (>2 wt% Y₂O₃), which probably have formed before garnet, scatter around 362 Ma. The retrograde branch of the P–T paths passes through 4 kbar at ~550°C. We conclude that the micaschists belong to a common metasedimentary sequence, which extends over the Porto Vecchio region and is separated from other metamorphic rock sequences in the north and the south by major tectonic boundaries. This sequence had experienced peak pressures which are lower than those determined for metamorphic rocks, such as micaschist and gneiss, from north‐eastern Sardinia. At present, we favour a continent–continent collisional scenario with the studied metasedimentary sequence buried during the collisional event as part of the upper plate. The contemporaneous high‐P metamorphic rocks from NE Sardinia were part of the upper portion of the lower plate. The addressed rocks from both plates were exhumed in an exhumation channel.