Inertial effects during irreversible meniscus reconfiguration in angular pores

In porous media, the dynamics of the invading front between two immiscible fluids is often characterized by abrupt reconfigurations caused by local instabilities of the interface. As a prototype of these phenomena we consider the dynamics of a meniscus in a corner as it can be encountered in angular pores. We investigate this process in detail by means of direct numerical simulations that solve the Navier–Stokes equations in the pore space and employ the Volume of Fluid method (VOF) to track the evolution of the interface. We show that for a quasi-static displacement, the numerically calculated surface energy agrees well with the analytical solutions that we have derived for pores with circular and square cross sections. However, the spontaneous reconfigurations are irreversible and cannot be controlled by the injection rate: they are characterized by the amount of surface energy that is spontaneously released and transformed into kinetic energy. The resulting local velocities can be orders of magnitude larger than the injection velocity and they induce damped oscillations of the interface that possess their own time scales and depend only on fluid properties and pore geometry. In complex media (we consider a network of cubic pores) reconfigurations are so frequent and oscillations last long enough that increasing inertial effects leads to a different fluid distribution by influencing the selection of the next pore to be invaded. This calls into question simple pore-filling rules based only on capillary forces. Also, we demonstrate that inertial effects during irreversible reconfigurations can influence the work done by the external forces that is related to the pressure drop in Darcy’s law. This suggests that these phenomena have to be considered when upscaling multiphase flow because local oscillations of the menisci affect macroscopic quantities and modify the constitutive relationships to be used in macro-scale models. These results can be extrapolated to other interface instabilities that are at the origin of fast pore-scale events, such as Haines jumps, snap-off and coalescence.     

Embedding complex hydrology in the regional climate system – Dynamic coupling across different modelling domains

To improve our understanding of the impacts of feedback between the atmosphere and the terrestrial water cycle including groundwater and to improve the integration of water resource management modelling for climate adaption we have developed a dynamically coupled climate–hydrological modelling system. The OpenMI modelling interface is used to couple a comprehensive hydrological modelling system, MIKE SHE running on personal computers, and a regional climate modelling system, HIRHAM running on a high performance computing platform. The coupled model enables two-way interaction between the atmosphere and the groundwater via the land surface and can represent the lateral movement of water in both the surface and subsurface and their interactions, not normally accounted for in climate models. Meso-scale processes are important for climate in general and rainfall in particular. Hydrological impacts are assessed at the catchment scale, the most important scale for water management. Feedback between groundwater, the land surface and the atmosphere occurs across a range of scales. Recognising this, the coupling was developed to allow dynamic exchange of water and energy at the catchment scale embedded within a larger meso-scale modelling domain. We present the coupling methodology used and describe the challenges in representing the exchanges between models and across scales. The coupled model is applied to one-way and two-way coupled simulations for a managed groundwater-dominated catchment, the Skjern River, Denmark. These coupled model simulations are evaluated against field observations and then compared with uncoupled climate and hydrological model simulations. Exploratory simulations show significant differences, particularly in the summer for precipitation and evapotranspiration the coupled model including groundwater and the RCM where groundwater is neglected. However, the resulting differences in the net precipitation and the catchment runoff in this groundwater dominated catchment were small. The need for further decadal scale simulations to understand the differences and insensitivity is highlighted.     

Improved numerical solvers for implicit coupling of subsurface and overland flow

Due to complex dynamics inherent in the physical models, numerical formulation of subsurface and overland flow coupling can be challenging to solve. ParFlow is a subsurface flow code that utilizes a structured grid discretization in order to benefit from fast and efficient structured solvers. Implicit coupling between subsurface and overland flow modes in ParFlow is obtained by prescribing an overland boundary condition at the top surface of the computational domain. This form of implicit coupling leads to the activation and deactivation of the overland boundary condition during simulations where ponding or drying events occur. This results in a discontinuity in the discrete system that can be challenging to resolve. Furthermore, the coupling relies on unstructured connectivities between the subsurface and surface components of the discrete system, which makes it challenging to use structured solvers to effectively capture the dynamics of the coupled flow. We present a formulation of the discretized algebraic system that enables the use of an analytic form of the Jacobian for the Newton–Krylov solver, while preserving the structured properties of the discretization. An effective multigrid preconditioner is extracted from the analytic Jacobian and used to precondition the Jacobian linear system solver. We compare the performance of the new solver against one that uses a finite difference approximation to the Jacobian within the Newton–Krylov approach, previously used in the literature. Numerical results explores the effectiveness of using the analytic Jacobian for the Newton–Krylov solver, and highlights the performance of the new preconditioner and its cost. The results indicate that the new solver is robust and generally outperforms the solver that is based on the finite difference approximation to the Jacobian, for problems where the overland boundary condition is activated and deactivated during the simulation. A parallel weak scaling study highlights the efficiency of the new solver.     

Distribution of the low velocity bulk in the middle-to-lower crust of southern Tibet: implications for formation of the north–south trending rift zones

We conducted the ambient noise tomography to image the shallow crustal structure of southern Tibet. The 2D maps of phase velocity anomalies at the periods of 10–16 s show that the low velocities are mainly confined along or near some of the rift zones. While the maps at the periods of 18–25 s show that the coherent patterns that the low velocities expand outside of the rift zones. It means that the low velocities are prevailing in the middle crust of southern Tibet. According to the previous study of surface wave tomography with teleseismic data, we find that the low velocities from the lower crust to the lithospheric mantle are also restricted to the same rift zones. Thus, the integrated knowledge of the distribution of the low velocities in southern Tibet provides some new insight on the formation of the north–south trending rift zones. Compiling the multidiscipline evidences, we conclude that the rifting was an integrated process of the entire lithosphere in the early stage (~26–10 Ma), but mainly occurred within the upper crust due to the weakening a decoupling in the low velocity middle crust in the late stage (later than ~8 Ma).     

A sensitivity analysis of lake water level response to changes in climate and river regimes

Lake water level regimes are influenced by climate, hydrology and land use. Intensive land use has led to a decline in lake levels in many regions, with direct impacts on lake hydrology, ecology and ecosystem services. This study examined the role of climate and river flow regime in controlling lake regimes using three different lakes with different hydraulic characteristics (volume-inflow ratio, CIR). The regime changes in the lakes were determined for five different river inflows and five different climate patterns (hot-arid, tropical, moderate, cold-arid, cold-wet), giving 75 different combinations of governing factors in lake hydrology. The input data were scaled to unify them for lake comparisons. By considering the historical lake volume fluctuations, the duration (number of months) of lake volume in different ‘wetness’ regimes from ‘dry’ to ‘wet’ was used to develop a new index for lake regime characterisation, ‘Degree of Lake Wetness’ (DLW). DLW is presented as two indices: DLW₁, providing a measure of lake filling percentage based on observed values and lake geometry, and DLW₂, providing an index for lake regimes based on historical fluctuation patterns. These indices were used to classify lake types based on their historical time series for variable climate and river inflow. The lake response time to changes in hydrology or climate was evaluated. Both DLW₁ and DLW₂ were sensitive to climate and hydrological changes. The results showed that lake level in high CIR systems depends on climate, whereas in systems with low CIR it depends more on river regime.     

Phytoplankton and periphyton production and its relation to temperature in a humic lagoon

To analyse the existence of interactive competition between phytoplankton and periphyton, we studied their photosynthesis–irradiance (P–E) response during one year in a humic lagoon. Lake production was dominated by phytoplankton, which followed seasonal changes in temperature. Periphyton primary production and algae biomass increased in winter, when phytoplankton biomass and production were lower. In this study we show that even in conditions of phytoplankton dominance, the habitat coupling between phytoplankton and periphyton can still be noticed.     

Transitional and freshwater bioassessments: One site,two perspectives?

The freshwater–saltwater-transition-zone was analysed using two different sampling protocols and assessment methodologies, developed for freshwater and estuaries, to compare their agreement level in terms of community composition and quality assessments. The use of different protocols resulted in significant differences in macroinvertebrate communities, in index scores and initially in quality classes. After modifications in the sensitivity scores of the IBMWP and AMBI indices (average scores or the use of a score of the other index when both were present), the differences were largely reduced and quality classes became coincident for the assessments provided by IPtIs and BAT tools. Such harmonisation of quality assessments for adjacent water categories (e.g., large rivers vs. transitional waters), exemplified here as an harmonisation in one of the metrics comprised in the assessment tools, is essential as it has direct implications on the expansion and accomplishment of River Basin Management Plans committed by the Water Framework Directive.     

The isolated elliptical galaxy NGC 5812–MOND or dark matter?

There exist isolated elliptical galaxies, whose dynamics can be modeled without resorting to dark matter or MOND, for example, NGC 7507. Such objects lack understanding within the current framework of galaxy formation. The isolated elliptical NGC 5812 is another object to investigate a possible role of isolation. We use globular clusters (GCs) and the galaxy light itself as dynamical tracers to constrain its mass profile. We employ Gemini/GMOS mask spectroscopy, apply the GMOS reduction procedures provided within IRAF, measure GC velocities by cross correlation methods and extract the line-of-sight kinematics of galaxy spectra using the tool pPXF. We identify 28 GCs with an outermost galactocentric distance of 20 kpc, for which velocities could be obtained. Furthermore, 16 spectra of the integrated galaxy light out to 6 kpc have been used to model the central kinematics. These spectra provide evidence for a disturbed velocity field, which is plausible given the disturbed morphology of the galaxy. We construct spherical Jeans models for the galaxy light and apply tracer mass estimators for the globular clusters. With the assumptions inherent to the mass estimators, MOND is compatible with the mass out to 20 kpc. However, a dark matter free galaxy is not excluded, given the uncertainties related to a possible nonsphericity and a possible nonequilibrium state. We find one globular cluster with an estimated mass of $1.6\times 10^7{M}_{\odot }$, the first Ultra Compact Dwarf in an isolated elliptical. We put NGC 5812 into the general context of dark matter or alternative ideas in elliptical galaxies. The case for a MONDian phenomenology also among early-type galaxies has become so strong that deviating cases appear astrophysically more interesting than agreements. The baryonic Tully Fisher relation (BTFR) as predicted by MOND is observed in some samples of early-type galaxies, in others not. However, in cases of galaxies that deviate from the MONDian prediction, data quality and data completeness are often problematic.     

The influence of quark anomalous magnetic moment in the Nambu–Jona-Lasinio model with different regularizations

The Nambu–Jona-Lasinio model is known for its simplicity and capacity to reproduce some of the basic characteristics of the quantum chromodynamics phase diagram. However, since it is a nonrenormalizable model, there are regularization issues that should be treated conveniently. This is the case when considering the quark anomalous magnetic moment (AMM) when external constant magnetic fields are present. Regularization procedures based on entangled functions between the magnetic field and the cutoff of the model can predict first-order phase transitions for chiral symmetry restoration at finite values of magnetic fields and inverse magnetic catalysis. The strengths of magnetic fields explored in NJL model and lattice QCD do not show first-order phase transition. In the present work, we show that some of the previous results are regularization-dependent effects and how to handle the divergences using the vacuum magnetic regularization scheme.