Estimation Of Sensible Heat Flux From The Surface Temperature Wave And One-Time-Of-Day Air Temperature Observation, Part Ii – Heterogeneous Surfaces

By means of the algorithm presented in Part I of this study, the temporal course H(t) and the daily mean H of the sensible heat flux H can be determined from measurements of the thermodynamic surface temperature (as a function of time) and from a one-time-of-day air temperature observation. Inaddition to these temperatures, one needs estimates of daily mean wind speed,of the roughness lengths of momentum and heat transfer, and of the displacementheight. In Part I, the algorithm was derived for areas with homogeneous surfaces,i.e., with uniform surface temperature, and the method was verified with measurements taken during several field campaigns. The root mean square error for the temperature difference between surface and air, in the comparison between measurement and model, amounted to one or two kelvin, and the error of H was 10 to 25 per cent. The method can be used to determine the sensible heat flux from measurements of surface temperatures by satellites, but can also be applied to ground based measurements.In Part II, the procedure is generalized for areas that consist of various surface types (sub-regions) with different surface temperatures, and can be usedwhen only a few (at least one) air temperature measurements per day are available over only one of the different sub-regions. This generalization should allow improvements to the estimates for H(t) by means of temperature measurements from, e.g., NOAA/AVHRR or LANDSAT/TM, taking into account the heterogeneity of the area contained in one METEOSAT pixel. Criteria are given as to whether effective (areal mean) surface temperatures and roughness lengths may be used for the computation of H or if the above mentioned generalized procedure has to be applied. The new algorithm is verified by measurements sampled during the field campaigns EFEDA 91 (Spain) and HIBE 89 (Hildesheimer Börde in Germany), and by using synthetic data (due to the lack of measured data) for one further combined surfacetype [soil and water (lakes)].     

On the Parameterisation of the Effective Roughness Length for Momentum Transfer over Heterogeneous Terrain

We develop a parameterisation for the effective roughness length of terrain that consists of a repeating sequence of patches, in which each patch is composed of strips of two roughness types. A numerical model with second-order closure in the turbulent stress is developed and used to show that: (i) the normalised Reynolds stress develops as a self-similar profile; (ii) the mixing-length parameterisation is a good first-order approximation to the Reynolds stress. These findings are used to characterise the blending layer, where the stress adjusts smoothly from its local surface value to its effective value aloft. Previous studies have assumed that this adjustment occurs abruptly at a single level, often called the blending height. The blending layer is shown to be characterised by height scales that arise naturally in linear models of surface layer flow over roughness changes, and calculations with the numerical model show that these height scales remain appropriate in the nonlinear regime. This concept of the blending layer allows the development of a new parameterisation of the effective roughness length, which gives values for the effective roughness length that are shown to compare well with both atmospheric measurements and values determined from the second-order model.     

Signatures in flowing fluid electric conductivity logs

Flowing fluid electric conductivity logging provides a means to determine hydrologic properties of fractures, fracture zones, or other permeable layers intersecting a borehole in saturated rock. The method involves analyzing the time-evolution of fluid electric conductivity (FEC) logs obtained while the well is being pumped and yields information on the location, hydraulic transmissivity, and salinity of permeable layers. The original analysis method was restricted to the case in which flows from the permeable layers or fractures were directed into the borehole (inflow). Recently, the method was adapted to permit treatment of both inflow and outflow, including analysis of natural regional flow in the permeable layer. A numerical model simulates flow and transport in the wellbore during flowing FEC logging, and fracture properties are determined by optimizing the match between simulation results and observed FEC logs. This can be a laborious trial-and-error procedure, especially when both inflow and outflow points are present. Improved analyses methods are needed. One possible tactic would be to develop an automated inverse method, but this paper takes a more elementary approach and focuses on identifying the signatures that various inflow and outflow features create in flowing FEC logs. The physical insight obtained provides a basis for more efficient analysis of these logs, both for the present trial and error approach and for a potential future automated inverse approach. Inflow points produce distinctive signatures in the FEC logs themselves, enabling the determination of location, inflow rate, and ion concentration. Identifying outflow locations and flow rates typically requires a more complicated integral method, which is also presented in this paper.     

Dating of lower terrace sediments from the Middle Rhine area, Germany

Modern and known-age Pleistocene fluvial sediments were investigated by optical dating of quartz to test the suitability of the approach for dating deposits from the deeply incised Middle Rhine Valley. Samples from modern flood sediments revealed skewed distributions indicating different residual levels of equivalent dose (D_e) within the different aliquots. Nevertheless, a substantial number of aliquots from the modern deposits reflect D_e values close to zero. For the Pleistocene samples, optical ages are in general consistent with age control given by the presence of the Laacher See Tephra and radiocarbon dating. However, some samples overestimate the known age by a few thousand years when using the arithmetic mean. This is apparently explained by including aliquots in the determination of mean D_e where the optical signal was incompletely bleached at deposition. The most difficult issue in this context is identifying a suitable approach that can distinguish between the variability of D_e due to partial bleaching and microdosimetry. However, even when considering these limitations it appears that optical dating will by a quite suitable method to date Pleistocene sediments from such a complex fluvial environment, especially when focusing on a precision scale beyond a few thousand years.     

Modeling patterns of coral bleaching at a remote Central Pacific atoll

A mild bleaching event (9.2% prevalence) at Palmyra Atoll occurred in response to the 2009 ENSO, when mean water temperature reached 29.8-30.1 °C. Prevalence among both abundant and sparse taxa varied with no clear pattern in susceptibility relating to coral morphology. Seven taxon-specific models showed that turbidity exacerbated while prior exposure to higher background temperatures alleviated bleaching, with these predictors explaining an average 16.3% and 11.5% variation in prevalence patterns, respectively. Positive associations occurred between bleaching prevalence and both immediate temperature during the bleaching event (average 8.4% variation explained) and increased sand cover (average 3.7%). Despite these associations, mean unexplained variation in prevalence equalled 59%. Lower bleaching prevalence in areas experiencing higher background temperatures suggests acclimation to temperature stress among several coral genera, while WWII modifications may still be impacting the reefs via shoreline sediment re-distribution and increased turbidity, exacerbating coral bleaching susceptibility during periods of high temperature stress.     

Statistical characteristics of flow as indicators of channeling in heterogeneous porous and fractured media

We introduce two new channeling indicators D_ic and D_cc based on the Lagrangian distribution of flow rates. On the basis of the participation ratio, these indicators characterize the extremes of both the flow-tube width distribution and the flow rate variation along flow lines. The participation ratio is an indicator biased toward the larger values of a distribution and is equal to the normalized ratio of the square of the first-order moment to the second-order moment. Compared with other existing indicators, they advantageously provide additional information on the flow channel geometry, are consistently applicable to both porous and fractured media, and are generally less variable for media generated using the same parameters than other indicators. Based on their computation for a broad range of porous and fracture permeability fields, we show that they consistently characterize two different geometric properties of channels. D_ic gives a characteristic scale of low-flow zones in porous media and a characteristic distance between effectively flowing structures in fractured cases. D_cc gives a characteristic scale of the extension of high-flow zones in porous media and a characteristic channel length in fractured media. D_ic is mostly determined by channel density and permeability variability. D_cc is, however, more affected by the nature of the correlation structure like the presence of permeability channels or fractures in porous media and the length distribution in fracture networks.     

Density-dependent dispersion in heterogeneous porous media Part I: A numerical study

In this paper, we describe carefully conducted numerical experiments, in which a dense salt solution vertically displaces fresh water in a stable manner. The two-dimensional porous media are weakly heterogeneous at a small scale. The purpose of these simulations, conducted for a range of density differences, is to obtain accurate concentration profiles that can be used to validate nonlinear models for high-concentration-gradient dispersion. In this part we focus on convergence of the computations, in numerical and statistical sense, to ensure that the uncertainty in the results is small enough.Concentration variances are computed, which give estimates of the uncertainty in local concentration values. These local variations decrease with increasing density contrast. For tracer transport, obtained longitudinal dispersivities are in accordance with analytical findings. In the case of high-density contrasts, stabilizing gravity forces counteract the growth of dispersive fingers, decreasing the effective width of the transition zone. For small log-permeability variances, the decrease of the apparent dispersivity that is found is in agreement with laboratory results for homogeneous columns.     

Testing for sufficient signal resetting during sediment transport using a polymineral multiple-signal luminescence approach

We present a multiple luminescence signal measurement procedure that simultaneously measures six different luminescence signals from a single polymineral aliquot (i.e. multiple-signal, short MS-SAR approach). The six signals show different bleaching rates in bleaching experiments, ranging from rapid bleaching for the quartz dominated blue stimulated luminescence signal (measured at 125 °C, BSL-125), to the slow-bleaching polymineral thermoluminescence signal. The bleaching rate of the infrared stimulated luminescence (IRSL) measured at room temperature (IR-25) and elevated temperature post-IR IRSL (pIRIR-90, pIRIR-155, pIRIR-225) signals decrease with increasing measurement temperature. Owing to these different bleaching rates, the MS-SAR approach allows inference of the degree of bleaching, and thereby information on the transport history of sediments. We test this approach by applying the MS-SAR to four coastal samples from a well-monitored sand-nourishment site at the Dutch coast. Our results show that the proposed MS-SAR approach can be utilised to construct bleaching plateaus which provide an independent and time-effective measure of the degree of poor bleaching in a sediment sample based on the measurement of only a few large aliquots. We propose that the MS-SAR protocol can be used to profile the age, luminescence properties and degree of bleaching of minimal prepared polymineral. This pre-profiling will allow the selection of suitable samples for full luminescence dating analysis in a target-orientated and time-effective manner.     

Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.