The interaction of two co-rotating quasi-geostrophic vortices in the vicinity of a surface buoyancy filament

In this paper, we investigate the interaction between two like-signed quasi-geostrophic uniform potential vorticity internal vortices in the vicinity of a surface buoyancy anomaly filament in a three dimensional, stably stratified and rapidly rotating fluid. The surface buoyancy distribution locally modifies the pressure fields and generates a shear flow. We start the study by first considering the effects of a uniform linear horizontal shear on the binary vortex interaction. We confirm that a cooperative shear facilitates the merger of a pair of vortices while an adverse shear has the opposite effect. We next investigate the binary vortex interaction in the vicinity of the surface buoyancy filament explicitly. Here, not only the filament generates a shear flow, but it also responds dynamically to the forcing by the vortex pair. The filament destabilises and forms buoyancy billows at the surface. These billows interact with the internal vortices. In particular, a surface billow may pair with one of the internal vortices. In such cases, the like-signed internal vortex pair may separate if they are initially moderately distant from each other.     

Interaction between a surface quasi-geostrophic buoyancy filament and an internal vortex

This paper focuses on the nonlinear interaction between a surface quasi-geostrophic buoyancy filament and an internal vortex. We first revisit the stability of an isolated buoyancy filament. The buoyancy profile considered is continuous and leads to a continuous velocity field, albeit one with infinite shear just outside its edge. The stability properties of an isolated filament help to interpret the unsteady interaction with a sub-surface (internal) vortex studied next. We find that, in all cases, the filament breaks into billows, analogous in form to those occurring in Kelvin–Helmholtz shear instability. For intense buoyancy filaments, the vortex itself may undergo strong deformations, including being split into several pieces. Generally, the nonlinear interaction causes both the filament and the vortex to lose their respective “self”-energies to the energy of interaction. The flow evolution depends sensitively on whether the vertical vorticity of the filament and the vortex have the same or opposite signs – termed “cooperative” and “adverse” shear respectively. In cooperative shear, the filament rolls up into a coherent surface eddy above a vortex initially placed below it, whereas in adverse shear, buoyancy is expelled above the vortex. Although sufficiently great shear induced by the buoyancy filament may split the vortex in both cases, adverse shear is significantly more destructive.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

Impacts of model initialization on an integrated surface water–groundwater model

Integrated hydrologic models characterize catchment responses by coupling the subsurface flow with land surface processes. One of the major areas of uncertainty in such models is the specification of the initial condition and its influence on subsequent simulations. A key challenge in model initialization is that it requires spatially distributed information on model states, groundwater levels and soil moisture, even when such data are not routinely available. Here, the impact of uncertainty in initial condition was explored across a 208 km² catchment in Denmark using the ParFlow.CLM model. The initialization impact was assessed under two meteorological conditions (wet vs dry) using five depth to water table and soil moisture distributions obtained from various equilibrium states (thermal, root zone, discharge, saturated and unsaturated zone equilibrium) during the model spin‐up. Each of these equilibrium states correspond to varying computation times to achieve stability in a particular aspect of the system state. Results identified particular sensitivity in modelled recharge and stream flow to the different initializations, but reduced sensitivity in modelled energy fluxes. Analysis also suggests that to simulate a year that is wetter than the spin‐up period, an initialization based on discharge equilibrium is adequate to capture the direction and magnitude of surface water–groundwater exchanges. For a drier or hydrologically similar year to the spin‐up period, an initialization based on groundwater equilibrium is required. Variability of monthly subsurface storage changes and discharge bias at the scale of a hydrological event show that the initialization impacts do not diminish as the simulations progress, highlighting the importance of robust and accurate initialization in capturing surface water–groundwater dynamics. Copyright © 2015 John Wiley & Sons, Ltd.     

Intercomparison across scales between remotely-sensed land surface temperature and representative equilibrium temperature from a distributed energy water balance model

Abstract

Quantifying the reliability of distributed hydrological models is an important task in hydrology to understand their ability to estimate energy and water fluxes at the agricultural district scale as well the basin scale for water resources management in drought monitoring and flood forecasting. In this context, the paper presents an intercomparison of simulated representative equilibrium temperature (RET) derived from a distributed energy water balance model and remotely-sensed land surface temperature (LST) at spatial scales from the agricultural field to the river basin. The main objective of the study is to evaluate the use of LST retrieved from operational remote sensing data at different spatial and temporal resolutions for the internal validation of a distributed hydrological model to control its mass balance accuracy as a complementary method to traditional calibration with discharge measurements at control river cross-sections. Modelled and observed LST from different radiometric sensors located on the ground surface, on an aeroplane and a satellite are compared for a maize field in Landriano (Italy), the agricultural district of Barrax (Spain) and the Upper Po River basin (Italy). A good ability of the model in reproducing the observed LST values in terms of mean bias error, root mean square error, relative error and Nash-Sutcliffe index is shown.

Editor Z.W. Kundzewicz; Associate editor D. Gerten     

Relationship between thermally forced surface wind and sea surface temperature gradient

An important part of the influence of the oceans on the atmosphere is through direct radiation, sensible heat flux and release of latent heat of evaporation, whereby all of these processes are directly related to the surface temperature of the oceans. A main effect of the atmosphere on the oceans is through momentum exchange at the air-ocean interface, and this process is directly related to the surface wind stress. The sea surface temperature (SST) and the surface wind stress are the two important components in the air-ocean system. If SST is given, a thermally forced boundary layer atmospheric circulation can be simulated. On the other hand, if the surface wind stress is given, the wind-driven ocean waves and ocean currents can be computed.The relationship between SST and surface wind is a coupling of the atmosphere and the oceans. It changes a one-way effect (ocean mechanically driven by atmosphere, or atmosphere thermally forced by oceans) into two-way air-sea interactions. Through this coupling the SST distribution, being an output from an ocean model, leads to the thermally forced surface winds, which feeds back into the ocean model as an additional forcing.Based on Kuo's planetary boundary layer model a linear algebraic equation is established to link the SST gradient with the thermally forced surface wind. The surface wind blows across the isotherms from cold to warm region with some deflection angle to the right (left) in the Northern (Southern) Hemisphere. Results from this study show that the atmospheric stratification reduces both the speed and the deflection angle of the thermally forced wind, however, the Coriolis' effect increases the wind speed in stable atmosphere (Ri>10^–4) and increases the deflection angle.     

Pilot-point based multi-objective calibration in a surface–subsurface distributed hydrological model

Abstract

In this study, a fully-coupled surface–subsurface, distributed, physics-based hydrological model was calibrated using the pilot-point method. A minimum variance field rule was included in the objective function to regularize the extensive calibration exercise that included 74 parameters (72 associated with pilot points and two spatially-invariant channel parameters). Because the overland and vadose zone systems are not in permanent hydrological connection, the information contained in the observation points may not be accessible by the pilot points at all times, rendering them insensitive to the observations and hindering the calibration process. An analysis of the spatial and temporal variability of parameter sensitivities was done to explore how the information contained in local observations spreads from the observation points to the pilot points, where parameter values are identified. The results show that the channel flow time series is valuable to identify the parameters at all pilot-point locations, indicating that the information in channel flow propagates to the entire basin. However, information in soil moisture measurements is of local extent and thus only valuable to identify the parameters at locations close to the observation point.

Editor D. Koutsoyiannis; Associate editor I. Nalbantis

Citation Maneta, M.P. and Wallender, W.W., 2013. Pilot-point based multi-objective calibration in a surface–subsurface distributed hydrological model. Hydrological Sciences Journal, 58 (2), 390–407.     

Coseismic and postseismic subsurface displacements and strains for a vertical strike-slip fault in a three-layer elastic medium

A three-layer elastic-gravitational fault displacement model using dislocation theory has been developed and used to examine the effect of layering of earth elastic moduli on surface and subsurface displacement fields for a vertical strike-slip fault. The model has been used to examine the effect of depth variation of elastic properties at coseismic and postseismic time scales. For pure strike-slip motion the effect of gravity on coseismic and postseismic horizontal deformation is negligible. For coseismic deformation the model predicts that (for constant Poisson's ratio) an increase in elastic moduli with depth attenuates the displacements within the upper layers with respect to displacement distribution for a uniform half-space, while an inclusion of a soft layer between the top layer and lower half-space amplifies upper layer displacements. The effect of variation in Poisson's ratio on surface and subsurface displacements has also been examined.The effect of postseismic stress relaxation on surface and subsurface displacements for a three-layer model has been calculated and compared with that of a uniformly relaxed half-space model. Layer 1 is assumed to correspond to the upper crust, layer 2 the lower crust and layer 3 the upper mantle. The effect of postseismic stress relaxation within a uniform half-space and within just the lower crust and upper mantle has been examined. Stress relaxation within the whole half-space decreases the amplitude and shortens the wavelength of displacements, while stress relaxation within the lower two layers increases the amplitude and broadens the wavelength of displacements. The difference between uniform and layered postseismic relaxation is particularly pronounced at the base of the crust.Coseismic and postseismic normal and volumetric strains for a vertical strike-slip fault have also been examined. For a uniformly relaxed half-space model, an increase in normal strains is shown with respect to the coseismic elastic solution, whereas the postseismic volumetric strain is effectively zero. For a three-layer model with stress relaxation in the lower layers only, the normal and volumetric strains within the top elastic layer resemble coseismic strains, while in the lower layers which suffer a rigidity decrease, the postseismic volumetric strain is effectively zero.     

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Relationship between sea surface salinity and ocean circulation and climate change

Based on Argo sea surface salinity(SSS) and the related precipitation(P), evaporation(E), and sea surface height data sets, the climatological annual mean and low-frequency variability in SSS in the global ocean and their relationship with ocean circulation and climate change were analyzed. Meanwhile, together with previous studies, a brief retrospect and prospect of seawater salinity were given in this work. Freshwater flux(E-P) dominated the mean pattern of SSS, while the dynamics of ocean circulation modulated the spatial structure and low-frequency variability in SSS in most regions. Under global warming, the trend in SSS indicated the intensification of the global hydrological cycle, and featured a decreasing trend at low and high latitudes and an increasing trend in subtropical regions. In the most recent two decades, global warming has slowed down, which is called the"global warming hiatus". The trend in SSS during this phase, which was different to that under global warming, mainly indicated the response of the ocean surface to the decadal and multi-decadal variability in the climate system, referring to the intensification of the Walker Circulation. The significant contrast of SSS trends between the western Pacific and the southeastern Indian Ocean suggested the importance of oceanic dynamics in the cross-basin interaction in recent decades. Ocean Rossby waves and the Indonesian Throughflow contributed to the freshening trend in SSS in the southeastern Indian Ocean, while the increasing trend in the southeastern Pacific and the decreasing trend in the northern Atlantic implied a long-term linear trend under global warming. In the future, higher resolution SSS data observed by satellites, together with Argo observations, will help to extend our knowledge on the dynamics of mesoscale eddies, regional oceanography, and climate change.     

Interaction between two sliders in a system with rate- and state-dependent friction

This study examines slip recurrence patterns in a two-block spring-slider model with rate- and state-dependent friction. Both weak and strong heterogeneities are considered with different settings of coupling stiffness. The results show that the recurrence pattern of slips strongly depends on the degree of coupling between the two blocks. With strong coupling between the two blocks (e.g., kc/ki max >~1), the slip pattern of the system is simple and characterized by periodical stick-slips, with the two blocks slipping together. With strong heterogeneity in friction strength, period-2 motion is found for moderate coupling stiffness (kc/ki max=0.4) between the two blocks. More complicated patterns are found with weak coupling stiffness (kc/ki max=0.2) and strong heterogeneity. With strong heterogeneity, very weak coupling leads to chaotic slip patterns. With coupling stiffness kc=5 ki max and strong heterogeneity, chaotic slip patterns are not found, in contrast with the results by Huang and Turcotte who employed the classical static/kinetic friction law.     

A coupled model for simulating surface water and groundwater interactions in coastal wetlands

Coastal wetlands are characterized by strong, dynamic interactions between surface water and groundwater. This paper presents a coupled model that simulates interacting surface water and groundwater flow and solute transport processes in these wetlands. The coupled model is based on two existing (sub) models for surface water and groundwater, respectively: ELCIRC (a three‐dimensional (3‐D) finite‐volume/finite‐difference model for simulating shallow water flow and solute transport in rivers, estuaries and coastal seas) and SUTRA (a 3‐D finite‐element/finite‐difference model for simulating variably saturated, variable‐density fluid flow and solute transport in porous media). Both submodels, using compatible unstructured meshes, are coupled spatially at the common interface between the surface water and groundwater bodies. The surface water level and solute concentrations computed by the ELCIRC model are used to determine the boundary conditions of the SUTRA‐based groundwater model at the interface. In turn, the groundwater model provides water and solute fluxes as inputs for the continuity equations of surface water flow and solute transport to account for the mass exchange across the interface. Additionally, flux from the seepage face was routed instantaneously to the nearest surface water cell according to the local sediment surface slope. With an external coupling approach, these two submodels run in parallel using time steps of different sizes. The time step (Δt_g) for the groundwater model is set to be larger than that (Δt_s) used by the surface water model for computational efficiency: Δt_g = M × Δt_s where M is an integer greater than 1. Data exchange takes place between the two submodels through a common database at synchronized times (e.g. end of each Δt_g). The coupled model was validated against two previously reported experiments on surface water and groundwater interactions in coastal lagoons. The results suggest that the model represents well the interacting surface water and groundwater flow and solute transport processes in the lagoons. Copyright © 2011 John Wiley & Sons, Ltd.     

Thermal and isotopic evidence for surface and subsurface water contributions to baseflow in a high Arctic river

Seven longitudinal water temperature tow surveys were conducted to attempt to identify the location of surface and subsurface river water exchanges along the length of the West River at the Cape Bounty Arctic Watershed Observatory, Melville Island, Nunavut, Canada (74°55′ N, 109°35′ W). Water temperature data were collected using a calibrated thermistor with an accuracy of ±0.002 °C (resolution <0.00005 °C) along the river during July 2014 in conjunction with stable water isotope sampling to support the thermal data and to determine the extent of surface water mixing from different sources such as precipitation, snowmelt, and surface/subsurface water contributions to the river. Atmospheric conditions were found to be the main contributor to seasonal temperature variance in the river, whereas tributary inflows and residual channel snow also had important thermal effects to river temperatures. Residual channel snow was a sustained source of cold water during much of the 2014 summer season (June–August) and substantially offset downstream warming. The longitudinal temperature profiles indicate notable changes to the thermal state of the river, which are interpreted to be indicative of subsurface and surface water exchange through inputs of relatively cold or warm water. Broadly, surface inflows were found to provide warmer water relative to the West River, and contributed to downstream warming of the river, along with downstream enrichment of δD and δ¹⁸O. Subsurface inflows provided cooler water relative to the river, and contributed to downstream depletion of δD and δ¹⁸O and downstream cooling of river temperatures. These results demonstrate that localized changes in river temperature, in conjunction with isotopic tracers, can be used to track channel–slope water interactions in Arctic hydrological systems, work previously limited to alpine and temperate settings.     

Quantifying watershed surface depression storage: determination and application in a hydrologic model

Hydrologic models often require correct estimates of surface macro‐depressional storage to accurately simulate rainfall–runoff processes. Traditionally, depression storage is determined through model calibration or lumped with soil storage components or on an ad hoc basis. This paper investigates a holistic approach for estimating surface depressional storage capacity (DSC) in watersheds using digital elevation models (DEMs). The methodology includes implementing a lumped DSC model to extract geometric properties of storage elements from DEMs of varying grid resolutions and employing a consistency zone criterion to quantify the representative DSC of an isolated watershed. DSC obtained using the consistency zone approach is compared to DSC estimated by “brute force” (BF) optimization method. The BF procedure estimates optimal DSC by calibrating DRAINMOD, a quasi‐process based hydrologic model, with observed streamflow under different climatic conditions. Both methods are applied to determine the DSC for relatively low‐gradient coastal plain watersheds on forested landscape with slopes less than 3%. Results show robustness of the consistency zone approach for estimating depression storage. To test the adequacy of the calculated DSC values obtained, both methods are applied in DRAINMOD to predict the daily watershed flow rates. Comparison between observed and simulated streamflow reveals a marginal difference in performance between BF optimization and consistency zone estimated DSCs during wet periods, but the latter performed relatively better in dry periods. DSC is found to be dependent on seasonal antecedent moisture conditions on surface topography. The new methodology is beneficial in situations where data on depressional storage is unavailable for calibrating models requiring this input parameter. Copyright © 2012 John Wiley & Sons, Ltd.     

A study of interaction between surface water and groundwater using environmental isotope in Huaisha River basin

The surface water and groundwater are important components of water cycle, and the interaction between surface water and groundwater is the important part in water cycle research. As the effective tracers in water cycle research, environmental isotope and hydrochemistry can reveal the interrelationships between surface water and groundwater effectively. The study area is the Huaisha River basin, which is located in Huairou district, Beijing. The field surveying and sampling for spring, river and well water were finished in 2002 and 2003. The hydrogen and oxygen isotopes and water quality were measured at the laboratory. The spatial characteristics in isotope and evolution of water quality along river lines at the different area were analyzed. The altitude effect of oxygen isotope in springs was revealed, and then using this equation, theory foundation for deducing recharge source of spring was estimated. By applying the mass balance method, the annual mean groundwater recharge rate at the catchment was estimated. Based on the groundwater recharge analysis, combining the hydrogeological condition analysis, and comparing the rainfall-runoff coefficients from the 1960s to 1990s in the Huaisha River basin and those in the Chaobai River basin, part of the runoff in the Huaisha River basin is recharged outside of this basin, in other words, this basin is an un-enclosed basin. On the basis of synthetically analyses, combining the compositions of hydrogen and oxygen isotopes and hydrochemistry, geomorphology, geology, and watershed systems characteristics, the relative contributions between surface water and groundwater flow at the different areas at the catchments were evaluated, and the interaction between surface water and groundwater was re- vealed lastly.     

三种下垫面温度及结冰预报模型研究

应用地表热量平衡方程,考虑太阳短波辐射、大气和地面的长波辐射、潜热、感热传输等能量之间的平衡,并考虑水汽、气溶胶、浮尘以及云等对太阳短波辐射的吸收和散射,建立了一种较实用的下垫面温度预报模型.应用湖北省恩施和金沙2009年冬季2月对土壤、水泥、沥青三种不同下垫面温度和自动气象站的常规气象要素观测进行模拟分析,并与该时段...     

Statistical equilibrium distributions of baroclinic vortices in a rotating two-layer model at low Froude numbers

The point-vortex equilibrium statistical model of two-layer baroclinic quasigeostrophic vortices in an unbounded f-plane is examined. A key conserved quantity, angular momentum, serves to confine the vortices to a compact domain, thereby justifying the statistical mechanics model, and also eliminating the need for boundary conditions in a practical method for its resolution. The Metropolis method provides a fast and efficient algorithm for solving the mean field non-linear elliptic PDEs of the equilibrium statistical theory. A verification of the method is done by comparison with the exact Gaussian solution at the no interaction limit of zero inverse temperature. The numerical results include a geophysically and computationally relevant power law for the radii at which the most probable vortex distribution is non-vanishing: For fixed total circulation, and fixed average angular momentum, the radii of both layers are proportional to the square root of the inverse temperature β. By changing the chemical potentials μ of the runs, one is able to model the most probable vorticity distributions for a wide range of total circulation and energy. The most probable vorticity distribution obtained at low positive temperatures are consistently close to a radially symmetric flat-top profiles. At high temperatures, the radially symmetric vorticity profiles are close to the Gaussian distribution.     

Sea surface height and transport stream function of the South China Sea from a variable-grid global ocean circulation model

The study on the South China Sea (SCS) circulation has a history of more than 40 years. Nevertheless, the SCS circulation is not fully understood compared with the Bohai Sea, Yellow Sea and East China Sea (ECS). Many numerical studies on the SCS circulati…     

Evaluating the interactions between surface water and groundwater in the arid mid-eastern Yanqi Basin,northwestern China

The environment of Bosten Lake in the Mid-Eastern Yanqi Basin (MEYB), an arid inland area in northwest China, has deteriorated greatly due to increasing groundwater exploitation and changes in the interactions between groundwater and surface water. This study intended to simulate the spatio-temporal variability of groundwater and surface water across the entire MEYB over the period 2000–2013. The applicable groundwater flow model and mass balance calculation method for river water were constructed to evaluate the change in groundwater recharged by and discharged to different segments of the Kaidu River. Simulation results show that the entire river seepage in the MEYB increased from 1.05 to 6.17 × 10⁸ m³/year between 2000 and 2013. The increasing river seepage, induced by increasing groundwater exploitation, plays the most important role in the water level decline in the downstream reaches of the Kaidu River and in Bosten Lake. This implies that the current utilization of groundwater resources in the MEYB is unsustainable.