Flow over Hills: A Large-Eddy Simulation of the Bolund Case

Simulation of local atmospheric flows around complex topography is important for several applications in wind energy (short-term wind forecasting and turbine siting and control), local weather prediction in mountainous regions and avalanche risk assessment. However, atmospheric simulation around steep mountain topography remains challenging, and a number of different approaches are used to represent such topography in numerical models. The immersed boundary method (IBM) is particularly well-suited for efficient and numerically stable simulation of flow around steep terrain. It uses a homogenous grid and permits a fast meshing of the topography. Here, we use the IBM in conjunction with a large-eddy simulation (LES) and test it against two unique datasets. In the first comparison, the LES is used to reproduce experimental results from a wind-tunnel study of a smooth three-dimensional hill. In the second comparison, we simulate the wind field around the Bolund Hill, Denmark, and make direct comparisons with field measurements. Both cases show good agreement between the simulation results and the experimental data, with the largest disagreement observed near the surface. The source of error is investigated by performing additional simulations with a variety of spatial resolutions and surface roughness properties.     

Identification and assessment of groundwater flow and storage components of the relict Schöneben Rock Glacier,Niedere Tauern Range,Eastern Alps (Austria)

More than 2,600 relict rock glaciers are known in the Austrian Alps but the knowledge of their hydraulic properties is severely limited. The relict Schöneben Rock Glacier (Niedere Tauern Range, Austria), with an extension of 0.17 km², was investigated based on spring data (2006–2014) and seismic refraction survey. Spring-discharge hydrographs and natural and artificial tracer data suggest a heterogeneous aquifer with a layered internal structure for the relict rock glacier. The discharge behavior exhibits a fast and a delayed flow component. The spring discharge responds to recharge events within a few hours but a mean residence time of several months can also be observed. The internal structure of the rock glacier (up to several tens of meters thick) consists of: an upper blocky layer with a few meters of thickness, which lacks fine-grained sediments; a main middle layer with coarse and finer-grained sediments, allowing for fast flow; and an approximately 10-m-thick basal till layer as the main aquifer body responsible for the base flow. The base-flow component is controlled by (fine) sandy to silty sediments with low hydraulic conductivity and high storage capacity, exhibiting a difference in hydraulic conductivity to the upper layer of about three orders of magnitude. The high storage capacity of relict rock glaciers has an impact on water resources management in alpine catchments and potentially regulates the risk of natural hazards such as floods and related debris flows. Thus, the results highlight the importance of such aquifer systems in alpine catchments.     

A method to reconstruct flood scenarios using field interviews and hydrodynamic modelling: application to the 2017 Suleja and Tafa,Nigeria flood

Natural Hazards - The scarcity of model input and calibration data has limited efforts in reconstructing scenarios of past floods in many regions globally. Recently, the number of studies that use...     

Low-temperature thermochronology of the flanks of the southern Upper Rhine Graben

The Upper Rhine Graben (URG) is the most perceptible part of the European Cenozoic Rift System. Uplifted Variscan basement of the Black Forest and the Vosges forms the flanks of the southern part of the graben. Apatite and zircon fission-track (FT) analyses indicate a complex low-temperature thermal history of the basement that was deciphered by inverse modelling of FT parameters. The models were tested against the observed data and independent geological constraints. The zircon FT ages of 28 outcrop samples taken along an E–W trending transect across the Black Forest and the Vosges range from 136 to 312 Ma, the apatite FT ages from 20 to 83 Ma. The frequency distributions of confined track lengths are broad and often bimodal in shape indicating a complex thermal history. Cooling below 120°C in the Early Cretaceous to Palaeogene was followed by a discrete heating episode during the late Eocene and subsequent cooling to surface temperature. The modelled time–temperature (t–T) paths point to a total denudation of the flanks of URG in the range of 1.0–1.7 km for a paleogeothermal gradient of 60°C/km, and 1.3–2.2 km for a paleogeothermal gradient of 45°C/km since the late Eocene.     

Trace elements and cathodoluminescence of quartz in stockwork veins of Mongolian porphyry-style deposits

The combination of scanning electron microscope–cathodoluminescence (CL), fluid inclusion analysis and high-resolution electron probe microanalysis of Al, Ti, K and Fe in vein quartz has yielded results permitting a greater understanding of the complex mineralisation of the Central Oyu Tolgoi and Zesen Uul porphyry-style deposits, southern Mongolia. These data demonstrate the relationship between quartz precipitation, dissolution and ore deposition as the mineralising fluid chemistry changed through time. Four major quartz generations are identified in the A-type veins from the stockworks of both the Central Oyu Tolgoi (OTi to OTiv) and Zesen Uul deposits (ZUi to ZUiv). Despite differences in the associated alteration and mineralisation style, the observed CL textures and trace element signatures of the quartz generations are comparable between deposits. The OTi and ZUi stage formed both the primary network of A-type veins and pervasive silicification of the host rock. Using the Ti-in-quartz geothermometer, crystallisation temperatures for OTi and ZUi of between 598°C and 880°C are indicated. The main stage of sulphide mineralisation was accompanied by the dissolution of pre-existing quartz (OTi and ZUi) and precipitation of a weakly luminescent generation of quartz (OTii and ZUii) with a low Ti content, reflected in a calculated temperature drop from approximately 700°C to 340°C in Central Oyu Tolgoi and 445°C in Zesen Uul. OTii and ZUii stage quartzes show high and variable Al concentrations. The next stage of quartz in both deposits (OTiii and ZUiii) forms a fine network of veins in cracks formed in pre-existing quartz. OTiii and ZUiii quartz contain measurable fluid inclusions of moderate salinity (3–17.1 wt.% NaCl eq.), entrapped in the temperature range 256°C to 385°C. OTiii and ZUiii are not related to any sulphide mineralisation. The final OTiv and ZUiv stages are characterised by quartz–calcite micro-breccias that penetrate the A-type veins. Based on the calculated entrapment temperatures, the OTiv/ZUiv stage crystallised between 212°C and 335°C, and the quartz is characterised by elevated but variable Al and Fe contents. The CL and trace element signatures of the OTi to OTiii and ZUi to ZUiii stages of the two Mongolian porphyries show similar features to those observed in porphyry-style deposits from other regions. This suggests that a common sequence of quartz crystallisation occurs during the formation of early veins in many porphyry copper systems.     

Intercomparison of Two Hydroxyl Radical Measurement Techniques at the Atmosphere Simulation Chamber SAPHIR

At the atmosphere simulation chamber SAPHIR in Jülich both Laser-Induced Fluorescence Spectroscopy (LIF) and Long-Path Differential Optical Laser Absorption Spectroscopy (DOAS) are operational for the detection of OH radicals at tropospheric levels. The two different spectroscopic techniques were compared within the controlled environment of SAPHIR based on all simultaneous measurements acquired in 2003 (13 days). Hydroxyl radicals were scavenged by added CO during four of these days in order to experimentally check the calculated precisions at the detection limit. LIF measurements have a higher precision (σ= 0.88×10⁶ cm^–3) and better time resolution (Δt = 60 s), but the DOAS method (σ= 1.24×10⁶ cm^–3, Δt = 135 s) is regarded as primary standard for comparisons because of its good accuracy. A high correlation coefficient of r = 0.95 was found for the whole data set highlighting the advantage of using a simulation chamber. The data set consists of two groups. The first one includes 3 days, where the LIF measurements yield (1 – 2) ×10⁶ cm^–3 higher OH concentrations than observed by the DOAS instrument. The experimental conditions during these days are characterized by increased NO_x concentration and a small dynamic range in OH. Excellent agreement is found within the other group of 6 days. The regression to the combined data of this large group yields unity slope without a significant offset.     

Multi-variable error correction of regional climate models

Climate change impact research needs regional climate scenarios of multiple meteorological variables. Those variables are available from regional climate models (RCMs), but affected by considerable biases. We evaluate the application of an empirical-statistical error correction method, quantile mapping (QM), for a small ensemble of RCMs and six meteorological variables. Annual and monthly biases are reduced to close to zero by QM for all variables in most cases. Exceptions are found, if non-stationarity of the model’s error characteristics occur. Even in the worst cases of non-stationarity, QM clearly improves the biases of raw RCMs. In addition, QM successfully adjusts the distributions of the analysed variables. To approach the question whether time series and inter-variable relationships are still plausible after correction, we evaluate the root-mean-square error (RMSE), autocorrelation and inter-variable correlation. We found improvement or no clear effect in RMSE and autocorrelation, and no clear effect on the correlation between meteorological variables. These results demonstrate that QM retains the quality of the temporal structure in time series and the inter-variable dependencies of RCMs. It has to be emphasised that this cannot be interpreted as an improvement and that deficiencies of the RCMs in those features are retained as well. Our results give some indication for the performance of QM applied to future scenarios, since our evaluation relies on independent calibration and evaluation periods, which are affected by climate variability and change. The effect of non-stationarity, however, can be expected to be larger in far future. We demonstrate the retainment of the RCM’s temporal structure and inter-variable dependencies, and large improvements in biases. This qualifies QM as a valuable, though not perfect, method in the interface between climate models and climate change impact research. Nonetheless, in case of no correlation between re-analysis driven RCM and observation, one should consider that QM does not correct this correlation.     

Basic analysis of climate and urban bioclimate of Dar es Salaam, Tanzania

Better understanding of urban microclimate and bioclimate of any city is imperative today when the world is constrained by both urbanisation and global climate change. Urbanisation generally triggers changes in land cover and hence influencing the urban local climate. Dar es Salaam city in Tanzania is one of the fast growing cities. Assessment of its urban climate and the human biometeorological conditions was done using the easily available synoptic meteorological data covering the period 2001–2011. In particular, the physiologically equivalent temperature (PET) was calculated using the RayMan software and results reveal that the afternoon period from December to February (DJF season) is relatively the most thermal stressful period to human beings in Dar es Salaam where PET values of above 35 °C were found. Additionally, the diurnal cycle of the individual meteorological elements that influence the PET index were analysed and found that air temperature of 30–35 °C dominate the afternoon period from 12:00 to 15:00 hours local standard time at about 60 % of occurrence. The current results, though considered as preliminary to the ongoing urban climate study in the city, provide an insight on how urban climate research is of significant importance in providing useful climatic information for ensuring quality of life and wellbeing of city dwellers.     

Topography-adjusted solar radiation indices and their importance in hydrology

Abstract

Solar radiation, direct and diffuse, is affected by surface characteristics, such as slope, aspect, altitude and shading. The paper examines the effects of topography on radiation, at multiple spatiotemporal scales, using suitable geometric methods for the direct and diffuse components. Two indices are introduced for comparing the direct radiation received by areas at the same and different latitudes. To investigate the profile of direct radiation across the whole of Greece, these are evaluated from an hourly to annual basis, via GIS techniques. Moreover, different approaches are examined for estimating the actual global radiation at operational spatial scales (sub-basin and terrain), according to the available meteorological data. The study indicates that the errors of typical hydrometeorological modelling formulas, which ignore the topographic effects and the seasonal allocation of direct and diffuse radiation, depend on the spatial scale and are non-uniformly distributed in time. In all cases, the estimations are improved by applying the proposed adjustment approaches. In particular, the adjustment of the measured global radiation ensures up to 10% increase of efficiency, while the modified Angström formula achieves slight (i.e. 2–4%) increase of efficiency and notable reduction of bias.

Editor Z.W. Kundzewicz

Citation Mamassis, N., Efstratiadis, A. and Apostolidou, I.-G., 2012. Topography-adjusted solar radiation indices and their importance in hydrology. Hydrological Sciences Journal, 57 (4), 756–775.     

Robust flood statistics: comparison of peak over threshold approaches based on monthly maxima and TL-moments

ABSTRACT

Flood quantile estimation based on partial duration series (peak over threshold, POT) represents a noteworthy alternative to the classical annual maximum approach since it enlarges the available information spectrum. Here the POT approach is discussed with reference to its benefits in increasing the robustness of flood quantile estimations. The classical POT approach is based on a Poisson distribution for the annual number of exceedences, although this can be questionable in some cases. Therefore, the Poisson distribution is compared with two other distributions (binomial and Gumbel-Schelling). The results show that only rarely is there a difference from the Poisson distribution. In the second part we investigate the robustness of flood quantiles derived from different approaches in the sense of their temporal stability against the occurrence of extreme events. Besides the classical approach using annual maxima series (AMS) with the generalized extreme value distribution and different parameter estimation methods, two different applications of POT are tested. Both are based on monthly maxima above a threshold, but one also uses trimmed L-moments (TL-moments). It is shown how quantile estimations based on this “robust” POT approach (rPOT) become more robust than AMS-based methods, even in the case of occasional extraordinary extreme events.

Editor M.C. Acreman Associate editor A. Viglione