Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

A high-resolution multiproxy study of sediment core GeoB7186-3 from SE Pacific off southern South America

A high-resolution multiproxy study performed on a marine record from SE Pacific off southern South America was used to reconstruct past regional environmental changes and their relation to global climate, particularly to El Niño/Southern Oscillation (ENSO) phenomenon during the last 2200 years. Our results suggest a sustained northward shift in the position of the zonal systems, i.e. the Southern Westerly Wind belt and the Antarctic Circumpolar Current, which occurred between 1300 and 750 yr BP. The synchrony of the latitudinal shift with cooling in Antarctica and reduced ENSO activity observed in several marine and terrestrial archives across South America suggests a causal link between ENSO and the proposed displacement of the zonal systems. This shift might have acted as a positive feedback to more La Niña-like conditions between 1300 and 750 yr BP by steepening the hemispheric and tropical Pacific zonal sea surface temperature gradient. This scenario further suggests different boundary conditions for ENSO before 1300 and after 750 yr BP.

Fly-through movie for the reconstructed folds in the onset region of the Petermann Glacier, North Greenland

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier.

Stable carbon and oxygen isotope ratios of foraminifera from Holocene sediments of the North Atlantic

The Denmark Strait Overflow (DSO) today compensates for the northward flowing Norwegian and Irminger branches of the North Atlantic Current that drive the Nordic heat pump. During the Last Glacial Maximum (LGM), ice sheets constricted the Denmark Strait aperture in addition to ice eustatic/isostatic effects which reduced its depth (today ~630 m) by ~130 m. These factors, combined with a reduced north-south density gradient of the water-masses, are expected to have restricted or even reversed the LGM DSO intensity. To better constrain these boundary conditions, we present a first reconstruction of the glacial DSO, using four new and four published epibenthic and planktic stable-isotope records from sites to the north and south of the Denmark Strait. The spatial and temporal distribution of epibenthic delta18O and delta13C maxima reveals a north-south density gradient at intermediate water depths from sigma0 ~28.7 to 28.4/28.1 and suggests that dense and highly ventilated water was convected in the Nordic Seas during the LGM. However, extremely high epibenthic delta13C values on top of the Mid-Atlantic Ridge document a further convection cell of Glacial North Atlantic Intermediate Water to the south of Iceland, which, however, was marked by much lower density (sigma0 ~28.1). The north-south gradient of water density possibly implied that the glacial DSO was directed to the south like today and fed Glacial North Atlantic Deep Water that has underthrusted the Glacial North Atlantic Intermediate Water in the Irminger Basin.

Geochemistry of sediments cores in the Mediterranean Sea

We investigated five time-equivalent core sections (180-110 kyr BP) from the Balearic Sea (Menorca Rise), the easternmost Levantine Basin and southwest, south, and southeast of Crete to reconstruct spatial patterns of productivity during deposition of sapropels S5 and S6 in the Mediterranean Sea. Our indicators are Ba, total organic carbon and carbonate contents. We found no indications of Ba remobilization within the investigated core intervals, and used the accumulation rate of biogenic Ba to compute paleoproductivity. Maximum surface water productivity (up to 350 g C/m2/yr) was found during deposition of S5 (isotope stage 5e) but pronounced spatial variability is evident. Coeval sediment intervals in the Balearic Sea show very little productivity change, suggesting that chemical and biological environments in the eastern and western Mediterranean basins were decoupled in this interval. We interpret the spatial variability as the result of two different modes of nutrient delivery to the photic zone: riverderived nutrient input and shoaling of the pycnocline/nutricline to the photic zone. The productivity increase during the formation of S6 was moderate compared to S5 and had a less marked spatial variability within the study area of the eastern Mediterranean Sea. Given that S6 formed during a glacial interval, glacial boundary conditions such as high wind stress and/or cooler surface water temperatures apparently favored lateral and vertical mixing and prevented the development of the spatial gradients within the Eastern Mediterranean Sea (EMS) observed for S5. A non-sapropel sediment interval with elevated Ba content and depleted 18O/16O ratios in planktonic foraminifer calcite was detected between S6 and S5 that corresponds to the weak northern hemisphere insolation maximum at 150 kyr. At this time, productivity apparently increased up to five times over surrounding intervals, but abundant benthic fauna show that the deep water remained oxic. Following our interpretation, the interval denotes a failed sapropel, when a weaker monsoon did not force the EMS into permanent stratification. The comparison of interglacial and glacial sapropels illustrates the relevance of climatic boundary conditions in the northern catchment in determining the facies and spatial variability of sapropels within the EMS.

Hydroclimate records and radiocarbon dating of sediment cores SO189/2 _39KL, SO189/2 _119KL, and SO189/2 _144KL from the eastern tropical Indian Ocean

The response of the tropical climate in the Indian Ocean realm to abrupt climate change events in the North Atlantic Ocean is contentious. Repositioning of the intertropical convergence zone is thought to have been responsible for changes in tropical hydroclimate during North Atlantic cold spells1, 2, 3, 4, 5, but the dearth of high-resolution records outside the monsoon realm in the Indian Ocean precludes a full understanding of this remote relationship and its underlying mechanisms. Here we show that slowdowns of the Atlantic meridional overturning circulation during Heinrich stadials and the Younger Dryas stadial affected the tropical Indian Ocean hydroclimate through changes to the Hadley circulation including a southward shift in the rising branch (the intertropical convergence zone) and an overall weakening over the southern Indian Ocean. Our results are based on new, high-resolution sea surface temperature and seawater oxygen isotope records of well-dated sedimentary archives from the tropical eastern Indian Ocean for the past 45,000 years, combined with climate model simulations of Atlantic circulation slowdown under Marine Isotope Stages 2 and 3 boundary conditions. Similar conditions in the east and west of the basin rule out a zonal dipole structure as the dominant forcing of the tropical Indian Ocean hydroclimate of millennial-scale events. Results from our simulations and proxy data suggest dry conditions in the northern Indian Ocean realm and wet and warm conditions in the southern realm during North Atlantic cold spells.

Stable carbon and oxygen isotope record of IODP Site 306-U1313, MIS 23-19

Since the seminal work by Hays et al. (1976), a plethora of studies has demonstrated a correlation between orbital variations and climatic change. However, information on how changes in orbital boundary conditions affected the frequency and amplitude of millennial-scale climate variability is still fragmentary. The Marine Isotope Stage (MIS) 19, an interglacial centred at around 785 ka, provides an opportunity to pursue this question and test the hypothesis that the long-term processes set up the boundary conditions within which the short-term processes operate. Similarly to the current interglacial, MIS 19 is characterised by a minimum of the 400-kyr eccentricity cycle, subdued amplitude of precessional changes, and small amplitude variations in insolation. Here we examine the record of climatic conditions during MIS 19 using high-resolution stable isotope records from benthic and planktonic foraminifera from a sedimentary sequence in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313) in order to assess the stability and duration of this interglacial, and evaluate the climate system's response in the millennial band to known orbitally induced insolation changes. Benthic and planktonic foraminiferal d18O values indicate relatively stable conditions during the peak warmth of MIS 19, but sea-surface and deep-water reconstructions start diverging during the transition towards the glacial MIS 18, when large, cold excursions disrupt the surface waters whereas low amplitude millennial scale fluctuations persist in the deep waters as recorded by the oxygen isotope signal. The glacial inception occurred at ~779 ka, in agreement with an increased abundance of tetra-unsaturated alkenones, reflecting the influence of icebergs and associated meltwater pulses and high-latitude waters at the study site. After having combined the new results with previous data from the same site, and using a variety of time series analysis techniques, we evaluate the evolution of millennial climate variability in response to changing orbital boundary conditions during the Early-Middle Pleistocene. Suborbital variability in both surface- and deep-water records is mainly concentrated at a period of ~11 kyr and, additionally, at ~5.8 and ~3.9 kyr in the deep ocean; these periods are equal to harmonics of precession band oscillations. The fact that the response at the 11 kyr period increased over the same interval during which the amplitude of the response to the precessional cycle increased supports the notion that most of the variance in the 11 kyr band in the sedimentary record is nonlinearly transferred from precession band oscillations. Considering that these periodicities are important features in the equatorial and intertropical insolation, these observations are in line with the view that the low-latitude regions play an important role in the response of the climate system to the astronomical forcing. We conclude that the effect of the orbitally induced insolation is of fundamental importance in regulating the timing and amplitude of millennial scale climate variability.

Effective stress, porosity, p-wave velocity and mineral composition of ODP Hole 174A-1073A sediments

Porosity, permeability, and compressional (P-wave) velocity were measured as a function of stress on sediments from Ocean Drilling Program Site 1073, U.S. Mid-Atlantic continental slope. Thin sections, scanning electron microscopy, and X-ray diffraction analyses provided mineralogical characteristics of the samples. Uniaxial strain boundary conditions were imposed on the samples during consolidation tests with the maximum effective axial stress reaching 13 MPa. The maximum effective radial stress necessary to maintain uniaxial strain was 7.6 MPa. Over an effective axial stress interval of 0 to 5.2 MPa, Sample 174A-1073A-26X-2, 82-89 cm (226.65 meters below seafloor [mbsf]), exhibited the largest decrease in porosity (51% to 41%), whereas Sample 71X-1, 2-8 cm (644.70 mbsf), exhibited the smallest decrease in porosity (48% to 45%). All samples showed negligible porosity increases during unloading. The permeability (on the order of 1 x 10-17 m**2) of Sample 174A-1073A-71X-1, 2-8 cm, was twice that measured on Sample 8H-1, 23-26 cm (63.75 mbsf), even though the former was considerably deeper and older. The differences in porosity-stress behavior and permeability between shallow and deep samples is related to lithologic, mineralogic, and diagenetic differences between the sediments above and below the Pliocene-Pleistocene to Miocene unconformity. P-wave velocity for Samples 174A-1073A-41X-5, 97-103 cm (372.35 mbsf), and 71X-1, 2-8 cm, increased with decreasing porosity, but did not change significantly during unloading.

The influence of seawater pH on U/ Ca ratios in the scleractinian cold-water coral Lophelia pertusa

The increasing pCO2 in seawater is a serious threat for marine calcifiers and alters the biogeochemistry of the ocean. Therefore, the reconstruction of past-seawater properties and their impact on marine ecosystems is an important way to investigate the underlying mechanisms and to better constrain the effects of possible changes in the future ocean. Cold-water coral (CWC) ecosystems are biodiversity hotspots. Living close to aragonite undersaturation, these corals serve as living laboratories as well as archives to reconstruct the boundary conditions of their calcification under the carbonate system of the ocean. We investigated the reef-building CWC Lophelia pertusa as a recorder of intermediate ocean seawater pH. This species-specific field calibration is based on a unique sample set of live in situ collected L. pertusa and corresponding seawater samples. These data demonstrate that uranium speciation and skeletal incorporation for azooxanthellate scleractinian CWCs is pH dependent and can be reconstructed with an uncertainty of ±0.15. Our Lophelia U / Ca-pH calibration appears to be controlled by the high pH values and thus highlighting the need for future coral and seawater sampling to refine this relationship. However, this study recommends L. pertusa as a new archive for the reconstruction of intermediate water mass pH and hence may help to constrain tipping points for ecosystem dynamics and evolutionary characteristics in a changing ocean.

Reduced order models for industrial thermo-fluid problems

A new method is presented to generate reduced order models (ROMs) in Fluid Dynamics problems of industrial interest. The method is based on the expansion of the flow variables in a Proper Orthogonal Decomposition (POD) basis, calculated from a limited number of snapshots, which are obtained via Computational Fluid Dynamics (CFD). Then, the POD-mode amplitudes are calculated as minimizers of a properly defined overall residual of the equations and boundary conditions. The method includes various ingredients that are new in this field. The residual can be calculated using only a limited number of points in the flow field, which can be scattered either all over the whole computational domain or over a smaller projection window. The resulting ROM is both computationally efficient(reconstructed flow fields require, in cases that do not present shock waves, less than 1 % of the time needed to compute a full CFD solution) and flexible(the projection window can avoid regions of large localized CFD errors).Also, for problems related with aerodynamics, POD modes are obtained from a set of snapshots calculated by a CFD method based on the compressible Navier Stokes equations and a turbulence model (which further more includes some unphysical stabilizing terms that are included for purely numerical reasons), but projection onto the POD manifold is made using the inviscid Euler equations, which makes the method independent of the CFD scheme. In addition, shock waves are treated specifically in the POD description, to avoid the need of using a too large number of snapshots. Various definitions of the residual are also discussed, along with the number and distribution of snapshots, the number of retained modes, and the effect of CFD errors. The method is checked and discussed on several test problems that describe (i) heat transfer in the recirculation region downstream of a backwards facing step, (ii) the flow past a two-dimensional airfoil in both the subsonic and transonic regimes, and (iii) the flow past a three-dimensional horizontal tail plane. The method is both efficient and numerically robust in the sense that the computational effort is quite small compared to CFD and results are both reasonably accurate and largely insensitive to the definition of the residual, to CFD errors, and to the CFD method itself, which may contain artificial stabilizing terms. Thus, the method is amenable for practical engineering applications. Resumen Se presenta un nuevo método para generar modelos de orden reducido (ROMs) aplicado a problemas fluidodinámicos de interés industrial. El nuevo método se basa en la expansión de las variables fluidas en una base POD, calculada a partir de un cierto número de snapshots, los cuales se han obtenido gracias a simulaciones numéricas (CFD). A continuación, las amplitudes de los modos POD se calculan minimizando un residual global adecuadamente definido que combina las ecuaciones y las condiciones de contorno. El método incluye varios ingredientes que son nuevos en este campo de estudio. El residual puede calcularse utilizando únicamente un número limitado de puntos del campo fluido. Estos puntos puede encontrarse dispersos a lo largo del dominio computacional completo o sobre una ventana de proyección. El modelo ROM obtenido es tanto computacionalmente eficiente (en aquellos casos que no presentan ondas de choque reconstruir los campos fluidos requiere menos del 1% del tiempo necesario para calcular una solución CFD) como flexible (la ventana de proyección puede escogerse de forma que evite contener regiones con errores en la solución CFD localizados y grandes). Además, en problemas aerodinámicos, los modos POD se obtienen de un conjunto de snapshots calculados utilizando un código CFD basado en la versión compresible de las ecuaciones de Navier Stokes y un modelo de turbulencia (el cual puede incluir algunos términos estabilizadores sin sentido físico que se añaden por razones puramente numéricas), aunque la proyección en la variedad POD se hace utilizando las ecuaciones de Euler, lo que hace al método independiente del esquema utilizado en el código CFD. Además, las ondas de choque se tratan específicamente en la descripción POD para evitar la necesidad de utilizar un número demasiado grande de snapshots. Varias definiciones del residual se discuten, así como el número y distribución de los snapshots,el número de modos retenidos y el efecto de los errores debidos al CFD. El método se comprueba y discute para varios problemas de evaluación que describen (i) la transferencia de calor en la región de recirculación aguas abajo de un escalón, (ii) el flujo alrededor de un perfil bidimensional en regímenes subsónico y transónico y (iii) el flujo alrededor de un estabilizador horizontal tridimensional. El método es tanto eficiente como numéricamente robusto en el sentido de que el esfuerzo computacional es muy pequeño comparado con el requerido por el CFD y los resultados son razonablemente precisos y muy insensibles a la definición del residual, los errores debidos al CFD y al método CFD en sí mismo, el cual puede contener términos estabilizadores artificiales. Por lo tanto, el método puede utilizarse en aplicaciones prácticas de ingeniería.

Theoretical analysis and numerical verification of the consistency of viscous smoothed-particle-hydrodynamics formulations in simulating free-surface flows

The theoretical formulation of the smoothed particle hydrodynamics (SPH) method deserves great care because of some inconsistencies occurring when considering free-surface inviscid flows. Actually, in SPH formulations one usually assumes that (i) surface integral terms on the boundary of the interpolation kernel support are neglected, (ii) free-surface conditions are implicitly verified. These assumptions are studied in detail in the present work for free-surface Newtonian viscous flow. The consistency of classical viscous weakly compressible SPH formulations is investigated. In particular, the principle of virtual work is used to study the verification of the free-surface boundary conditions in a weak sense. The latter can be related to the global energy dissipation induced by the viscous term formulations and their consistency. Numerical verification of this theoretical analysis is provided on three free-surface test cases including a standing wave, with the three viscous term formulations investigated.

A proposed parameterization of interface discontinuity factors depending on neighborhood for pin-by-pin diffusion computations for LWR

There exists an interest in performing full core pin-by-pin computations for present nuclear reactors. In such type of problems the use of a transport approximation like the diffusion equation requires the introduction of correction parameters. Interface discontinuity factors can improve the diffusion solution to nearly reproduce a transport solution. Nevertheless, calculating accurate pin-by-pin IDF requires the knowledge of the heterogeneous neutron flux distribution, which depends on the boundary conditions of the pin-cell as well as the local variables along the nuclear reactor operation. As a consequence, it is impractical to compute them for each possible configuration. An alternative to generate accurate pin-by-pin interface discontinuity factors is to calculate reference values using zero-net-current boundary conditions and to synthesize afterwards their dependencies on the main neighborhood variables. In such way the factors can be accurately computed during fine-mesh diffusion calculations by correcting the reference values as a function of the actual environment of the pin-cell in the core. In this paper we propose a parameterization of the pin-by-pin interface discontinuity factors allowing the implementation of a cross sections library able to treat the neighborhood effect. First results are presented for typical PWR configurations.

The estimation of truncation error by tau-estimation revisited

The aim of this paper was to accurately estimate the local truncation error of partial differential equations, that are numerically solved using a finite difference or finite volume approach on structured and unstructured meshes. In this work, we approximated the local truncation error using the @t-estimation procedure, which aims to compare the residuals on a sequence of grids with different spacing. First, we focused the analysis on one-dimensional scalar linear and non-linear test cases to examine the accuracy of the estimation of the truncation error for both finite difference and finite volume approaches on different grid topologies. Then, we extended the analysis to two-dimensional problems: first on linear and non-linear scalar equations and finally on the Euler equations. We demonstrated that this approach yields a highly accurate estimation of the truncation error if some conditions are fulfilled. These conditions are related to the accuracy of the restriction operators, the choice of the boundary conditions, the distortion of the grids and the magnitude of the iteration error.

A three-dimensional dynamical model of Amundsenisen icefield, Svalbard.

Amundsenisen is an ice field, 80 km2 in area, located in Southern Spitsbergen, Svalbard. Radio-echo sounding measurements at 20 MHz show high intensity returns from a nearly flat basal reflector at four zones, all of them with ice thickness larger than 500m. These reflections suggest possible subglacial lakes. To determine whether basal liquid water is compatible with current pressure and temperature conditions, we aim at applying a thermo mechanical model with a free boundary at the bed defined as solution of a Stefan problem for the interface ice-subglaciallake. The complexity of the problem suggests the use of a bi-dimensional model, but this requires that well-defined flowlines across the zones with suspected subglacial lakes are available. We define these flow lines from the solution of a three-dimensional dynamical model, and this is the main goal of the present contribution. We apply a three-dimensional full-Stokes model of glacier dynamics to Amundsenisen icefield. We are mostly interested in the plateau zone of the icefield, so we introduce artificial vertical boundaries at the heads of the main outlet glaciers draining Amundsenisen. At these boundaries we set velocity boundary conditions. Velocities near the centres of the heads of the outlets are known from experimental measurements. The velocities at depth are calculated according to a SIA velocity-depth profile, and those at the rest of the transverse section are computed following Nye’s (1952) model. We select as southeastern boundary of the model domain an ice divide, where we set boundary conditions of zero horizontal velocities and zero vertical shear stresses. The upper boundary is a traction-free boundary. For the basal boundary conditions, on the zones of suspected subglacial lakes we set free-slip boundary conditions, while for the rest of the basal boundary we use a friction law linking the sliding velocity to the basal shear stress,in such a way that, contrary to the shallow ice approximation, the basal shear stress is not equal to the basal driving stress but rather part of the solution.