Boussinesq-type modelling using an unstructured finite element technique

A model for solving the two-dimensional enhanced Boussinesq equations is presented. The model equations are discretised in space using an unstructured finite element technique. The standard Galerkin method with mixed interpolation is applied. The time discretisation is performed using an explicit three-step Taylor–Galerkin method. The model is extended to the surf and swash zone by inclusion of wave breaking and a moving boundary at the shoreline. Breaking is treated by an existing surface roller model, but a new procedure for the detection of the roller thickness is devised. The model is verified using four test cases and the results are compared with experimental data and results from an existing finite difference Boussinesq model.     

Successful application of multiscale methods in a real reservoir simulator environment

For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs.     

Estimating environmentally relevant fixed nitrogen demand in the 21st century

Human activities affect the impact of the nitrogen cycle on both the environment and climate. The rate of anthropogenic nitrogen fixation from atmospheric N₂ may serve as an indicator to the magnitude of this impact, acknowledging that relationship to be effect-dependent and non-linear. Building on the set of Representative Concentration Pathway (RCP) scenarios developed for climate change research, we estimate anthropogenic industrial nitrogen fixation throughout the 21st century. Assigning characteristic key drivers to the four underlying scenarios we arrive at nitrogen fixation rates for agricultural use of 80 to 172 Tg N/yr by 2100, which is slightly less to almost twice as much compared with the fixation rate for the year 2000. We use the following key drivers of change, varying between scenarios: population growth, consumption of animal protein, agricultural efficiency improvement and additional biofuel production. Further anthropogenic nitrogen fixation for production of materials such as explosives or plastics and from combustion are projected to remain considerably smaller than that related to agriculture. While variation among the four scenarios is considerable, our interpretation of scenarios constrains the option space: several of the factors enhancing the anthropogenic impact on the nitrogen cycle may occur concurrently, but never all of them. A scenario that is specifically targeted towards limiting greenhouse gas emissions ends up as the potentially largest contributor to nitrogen fixation, as a result of large amounts of biofuels required and the fertilizer used to produce it. Other published data on nitrogen fixation towards 2100 indicate that our high estimates based on the RCP approach are rather conservative. Even the most optimistic scenario estimates that nitrogen fixation rate will remain substantially in excess of an estimate of sustainable boundaries by 2100.     

Downscaling a twentieth century global climate simulation to the North Sea

The regional ocean model system (ROMS) is used to downscale a 26-year period of the twentieth century 20C3M experiment from the global coupled Bergen climate model (BCM) for the North Sea. Compared to an observational-based climatology, BCM have good results on the mean temperature, except for too low winter temperature. This is connected to a too weak inflow of Atlantic water. The downscaling gives added value to the BCM results by providing regional details, doubling the Atlantic inflow, and improving the mean winter temperature. For mean salinity, BCM has values very close to the climatology, whereas the downscaling becomes too fresh. The downscaling, however, improves the sea surface salinity, the vertical structure, and the Norwegian Coastal Current. It is concluded that the downscaling procedure as presented here is a suitable tool for assessing the future Atlantic inflow and sea temperature in the North Sea based on a global climate projection.     

Biogeochemistry of dissolved carbon,major, and trace elements during spring flood periods on the Ob River

Detailed knowledge of the flood period of Arctic rivers remains one of the few factors impeding rigorous prediction of the effect of climate change on carbon and related element fluxes from the land to the Arctic Ocean. In order to test the temporal and spatial variability of element concentration in the Ob River (western Siberia) water during flood period and to quantify the contribution of spring flood period to the annual element export, we sampled the main channel year round in 2014–2017 for dissolved C, major, and trace element concentrations. We revealed high stability (approximately ≤10% relative variation) of dissolved C, major, and trace element concentrations in the Ob River during spring flood period over a 1‐km section of the river channel and over 3 days continuous monitoring (3‐hr frequency). We identified two groups of elements with contrasting relationship to discharge: (a) DIC and soluble elements (Cl, SO₄, Li, B, Na, Mg, Ca, P, V, Cr, Mn, As, Rb, Sr, Mo, Ba, W, and U) negatively correlated (p < 0.05) with discharge and exhibited minimal concentrations during spring flood and autumn high flow and (b) DOC and particle‐reactive elements (Al, Fe, Ti, Y, Zr, Nb, Cs, REEs, Hf, Tl, Pb, and Th), some nutrients (K), and metalloids (Ge, Sb, and Te), positively correlated (p < 0.05) with discharge and showed the highest concentrations during spring flood. We attribute the decreased concentration of soluble elements with discharge to dilution by groundwater feeding and increased concentration of DOC and particle‐reactive metals with discharge to leaching from surface soil, plant litter, and suspended particles. Overall, the present study provides first‐order assessment of fluxes of major and trace elements in the middle course of the Ob River, reveals their high temporal and spatial stability, and characterizes the mechanism of river water chemical composition acquisition.     

How Well Does Chlorophyll Explain the Seasonal Variation in Phytoplankton Activity?

The seasonal variation in phytoplankton activity is determined by analysing 1385 primary production (PP) profiles, chlorophyll a (Chl) concentration profiles and phytoplankton carbon biomass concentrations (C) from the period 1998–2012. The data was collected at six different stations in the Baltic Sea transition zone (BSTZ) which is a location with strong seasonal production patterns with light as the key parameter controlling this productivity. We show that the use of Chl as a proxy for phytoplankton activity strongly overestimates the contribution from the spring production to annual pelagic carbon flow. Spring (February and March) Chl comprised 16–30% of the total annual Chl produced, whereas spring C was much lower (8–23%) compared to the annual C. Spring PP accounted for 10–18% of the total annual PP, while the July–August production contributed 26–33%, i.e. within the time frame when zooplankton biomass and grazing pressure are highest. That is, Chl failed in this study to reflect the importance of the high summer PP. A better proxy for biomass may be C, which correlated well with the seasonal pattern of PP (Pearson correlation, p < 0.05). Thus, this study suggests to account for the strong seasonal pattern in C/Chl ratios when considering carbon flow in coastal systems. Seasonal data for PP were fitted to a simple sinusoidal wave model describing the seasonal distribution of PP in the BSTZ and were proposed to present a better parameterizaton of PP in shallow stratified temperate regions than more commonly applied proxies.     

Sulphur-containing compounds in sulphur-rich crude oils from hypersaline lake sediments and their geochemical implications

Three sulphur-rich commercial crude oils have been studied, which contain sulphur as high as up to 4–12%. These samples were collected from Tertiary hypersaline lake sediments of the Jianghan Basin, Hubei Province at different depths, but above the oil generation threshold (2200m). FPD-GC and GC-MS data show that aromatic fractions of the crude oils are composed of different homologues of sulphur-containing compounds, including long-chain normal alkyl-thiophenes and-thiolanes, long-chain isoprenoid-thiophenes and -thiolanes, and benzothiophenes. It is worth noting that the distribution patterns of long-chain alkyl-thiophenes and -thiolanes from two shallow-seated crude oils are quite similar to those of normal alkanes showing marked even-odd predominance. It seems that the even-odd predominance of sulphur-containing compounds decreases with increasing burial depth of the crude oils. The major component of aliphatic fraction is phytane, and similarly the major peaks of aromatic fractions also represent C₂₀ isoprenoid thiophenes. Some preliminary conclusions have been drawn from the above discussion: (1) Abundant sulphur-containing compounds may be used as an indicator of low mature or immature crude oils produced from hypersaline lake sediments; (2) Sulphur-containing compounds are considered to be early diagenetic products of reactions between elemental sulphur or sulfides and alkanes or their precursors (phytols, fatty acids, alcohols, etc.), or of bacterial activities, but not direct inputs of organisms.