Settlement history of a mountain oasis in Northern Oman

To unravel the settlement history of oases in northern Oman, data on topography, the agricultural setting, water and soil parameters and archaeological findings were collected in the Wadi Bani Awf with its head oasis Balad Seet. Data collection lasted from April 2000 to April 2003 and was based on the establishment of a 3D-georeferenced map of the oasis comprising all its major infrastructural and agronomic features. At today's Balad Seet, a total of 8.8 ha are planted to 2,800 date palms and 4.6 ha are divided into 385 small fields dedicated to wheat, barley, sorghum, oats, alfalfa, garlic, onion, lime and banana. Radiocarbon dating of charcoal in the lower part of the main terrace system determined its age to 911 ± 43 years. Monthly flow measurements of four major aflaj systems showed a total maximum flow of 32 m^3 h^-1 with the largest falaj contributing 78% of the total flow. During drought periods, average water flow decreased by 3% per month, however, with significant differences between the spring systems. The analysis of the tritium/^3helium ratio in the water led to an estimated water age of up to 10 years. In combination with the flow data, this provided insights into the elasticity of the spring flow over time. The use of the natural resources of the Wadi Bani Awf by a pastoral population started probably in the early 3rd millennium BC. The first permanent settlement might have been established at Balad Seet during the first part of the 1st millennium BC. Presumably it was initiated by settlers from al-Hamra, a village at the southern foot of the Hajar mountains. Given an abundant und stable flow of springs, even in periods of drought, the construction of Balad Seet's first irrigation systems may have occurred at this early time. The combination of topographic, agricultural, hydro-pedological and archaeological data allowed assessment of the carrying capacity of this oasis over the three millennia of its likely existence. The changing scarcity of land and water and the eventual optimisation of their use by different aflaj constructions have been major driving forces for the development and apparent relativeley stable existence of this oasis.

A High Order Finite Volume Scheme for the 2D Shallow Water Equations Including Topography

Inhalt dieser Arbeit ist ein Verfahren zur numerischen Lösung der zweidimensionalen Flachwassergleichung, welche das Fließverhalten von Gewässern, deren Oberflächenausdehnung wesentlich größer als deren Tiefe ist, modelliert. Diese Gleichung beschreibt die gravitationsbedingte zeitliche Änderung eines gegebenen Anfangszustandes bei Gewässern mit freier Oberfläche. Diese Klasse beinhaltet Probleme wie das Verhalten von Wellen an flachen Stränden oder die Bewegung einer Flutwelle in einem Fluss. Diese Beispiele zeigen deutlich die Notwendigkeit, den Einfluss von Topographie sowie die Behandlung von Nass/Trockenübergängen im Verfahren zu berücksichtigen. In der vorliegenden Dissertation wird ein, in Gebieten mit hinreichender Wasserhöhe, hochgenaues Finite-Volumen-Verfahren zur numerischen Bestimmung des zeitlichen Verlaufs der Lösung der zweidimensionalen Flachwassergleichung aus gegebenen Anfangs- und Randbedingungen auf einem unstrukturierten Gitter vorgestellt, welches in der Lage ist, den Einfluss topographischer Quellterme auf die Strömung zu berücksichtigen, sowie in sogenannten \glqq lake at rest\grqq-stationären Zuständen diesen Einfluss mit den numerischen Flüssen exakt auszubalancieren. Basis des Verfahrens ist ein Finite-Volumen-Ansatz erster Ordnung, welcher durch eine WENO Rekonstruktion unter Verwendung der Methode der kleinsten Quadrate und eine sogenannte Space Time Expansion erweitert wird mit dem Ziel, ein Verfahren beliebig hoher Ordnung zu erhalten. Die im Verfahren auftretenden Riemannprobleme werden mit dem Riemannlöser von Chinnayya, LeRoux und Seguin von 1999 gelöst, welcher die Einflüsse der Topographie auf den Strömungsverlauf mit berücksichtigt. Es wird in der Arbeit bewiesen, dass die Koeffizienten der durch das WENO-Verfahren berechneten Rekonstruktionspolynome die räumlichen Ableitungen der zu rekonstruierenden Funktion mit einem zur Verfahrensordnung passenden Genauigkeitsgrad approximieren. Ebenso wird bewiesen, dass die Koeffizienten des aus der Space Time Expansion resultierenden Polynoms die räumlichen und zeitlichen Ableitungen der Lösung des Anfangswertproblems approximieren. Darüber hinaus wird die wohlbalanciertheit des Verfahrens für beliebig hohe numerische Ordnung bewiesen. Für die Behandlung von Nass/Trockenübergangen wird eine Methode zur Ordnungsreduktion abhängig von Wasserhöhe und Zellgröße vorgeschlagen. Dies ist notwendig, um in der Rechnung negative Werte für die Wasserhöhe, welche als Folge von Oszillationen des Raum-Zeit-Polynoms auftreten können, zu vermeiden. Numerische Ergebnisse die die theoretische Verfahrensordnung bestätigen werden ebenso präsentiert wie Beispiele, welche die hervorragenden Eigenschaften des Gesamtverfahrens in der Berechnung herausfordernder Probleme demonstrieren.

Hydraulic modeling of runoff over a rough surface under partial inundation

This study presents a numerical method to derive the Darcy- Weisbach friction coefficient for overland flow under partial inundation of surface roughness. To better account for the variable influence of roughness with varying levels of emergence, we model the flow over a network which evolves as the free surface rises. This network is constructed using a height numerical map, based on surface roughness data, and a discrete geometry skeletonization algorithm. By applying a hydraulic model to the flows through this network, local heads, velocities, and Froude and Reynolds numbers over the surface can be estimated. These quantities enable us to analyze the flow and ultimately to derive a bulk friction factor for flow over the entire surface which takes into account local variations in flow quantities. Results demonstrate that although the flow is laminar, head losses are chiefly inertial because of local flow disturbances. The results also emphasize that for conditions of partial inundation, flow resistance varies nonmonotonically but does generally increase with progressive roughness inundation.

The Effects of Stratification on Flow Separation

Separation of stratified flow over a two-dimensional hill is inhibited or facilitated by acceleration or deceleration of the flow just outside the attached boundary layer. In this note, an expression is derived for this acceleration or deceleration in terms of streamline curvature and stratification. The expression is valid for linear as well as nonlinear deformation of the flow. For hills of vanishing aspect ratio a linear theory can be derived and a full regime diagram for separation can be constructed. For hills of finite aspect ratio scaling relationships can be derived that indicate the presence of a critical aspect ratio, proportional to the stratification, above which separation will occur as well as a second critical aspect ratio above which separation will always occur irrespective of stratification.

The influence of hilly terrain on canopy-atmosphere carbon dioxide exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO2 (N-ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO2 sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO2 exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO2. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N-ee = F-c, the vertical turbulent flux of CO2. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F-c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F-c-at canopy top differs by more than 20% compared to the flat-terrain case.

Factors governing cellular convection in orographic precipitation

The development of shallow cellular convection in warm orographic clouds is investigated through idealized numerical simulations of moist flow over topography using a cloud-resolving numerical model. Buoyant instability, a necessary element for moist convection, is found to be diagnosed most accurately through analysis of the moist Brunt–Väisälä frequency (N_m) rather than the vertical profile of θ_e. In statically unstable orographic clouds (N_m^2) < 0), additional environmental and terrain-related factors are shown to have major effects on the amount of cellularity that occurs in 2D simulations. One of these factors, the basic-state wind shear, may suppress convection in 2D yet allow for longitudinal convective roll circulations in 3D. The presence of convective structures within an orographic cloud substantially enhanced the maximum rainfall rates, precipitation efficiencies, and precipitation accumulations in all simulations.

Further investigation of the impact of idealised continents and SST distributions on the Northern Hemisphere storm tracks

Building on studies by Brayshaw et al. (2009, 2011) of the basic ingredients of the North Atlantic storm track (land-sea contrast, orography and SST), this article investigates the impact of Eurasian topography and Pacific SST anomalies on the North Pacific and Atlantic storm tracks through a hierarchy of atmospheric GCM simulations using idealised boundary conditions in the HadGAM1 model. The Himalaya-Tibet mountain complex is found to play a crucial role in shaping the North Pacific storm track. The northward deflection of the westerly flow around northern Tibet generates an extensive pool of very cold air in the north-eastern tip of the Asian continent, which strengthens the meridional temperature gradient and favours baroclinic growth in the western Pacific. The Kuroshio SST front is also instrumental in strengthening the Pacific storm track through its impact on near-surface baroclinicity, while the warm waters around Indonesia tend to weaken it through the impact on baroclinicity of stationary Rossby waves propagating poleward from the convective heating regions. Three mechanisms by which the Atlantic storm track may be affected by changes in the boundary conditions upstream of the Rockies are discussed. In the model configuration used here, stationary Rossby waves emanating from Tibet appear to weaken the North Atlantic storm track substantially, whereas those generated over the cold waters off Peru appear to strengthen it. Changes in eddy-driven surface winds over the Pacific generally appear to modify the flow over the Rocky Mountains, leading to consistent modifications in the Atlantic storm track. The evidence for each of these mechanisms is, however, ultimately equivocal in these simulations.

Estimation of the friction velocity in stably stratified boundary-layer flows over hills

A method is suggested for the calculation of the friction velocity for stable turbulent boundary-layer flow over hills. The method is tested using a continuous upstream mean velocity profile compatible with the propagation of gravity waves, and is incorporated into the linear model of Hunt, Leibovich and Richards with the modification proposed by Hunt, Richards and Brighton to include the effects of stability, and the reformulated solution of Weng for the near-surface region. Those theoretical results are compared with results from simulations using a non-hydrostatic microscale-mesoscale two-dimensional numerical model, and with field observations for different values of stability. These comparisons show a considerable improvement in the behaviour of the theoretical model when the friction velocity is calculated using the method proposed here, leading to a consistent variation of the boundary-layer structure with stability, and better agreement with observational and numerical data.

Pulsatory andesite lava flow at Bagana Volcano

Using a time series of TerraSAR-X spaceborne radar images we have measured the pulsatory motion of an andesite lava flow over a 14-month period at Bagana volcano, Papua New Guinea. Between October 2010 and December 2011, lava flowed continuously down the western flank of the volcano forming a 3 km-long blocky lava flow with a channel, levees, overflows and branches. We captured four successive pulses of lava advancing down the channel system, the first such behaviour of an andesite flow to be recorded using radar. Each pulse had a volume of the order of 107 m3 emplaced over many weeks. The average extrusion rate estimated from the radar data was 0.92 ± 0.35 m3 s-1 , and varied between 0.3 and 1.8 m3 s-1, with higher rates occurring earlier in each pulse. This, together with observations of sulphur dioxide emissions, explosions and incandescence suggest a variable supply rate of magma through Bagana’s conduit as the most likely source of the pulsatory behaviour.

The effects of rotation and ice shelf topography on frazil-laden ice shelf water plumes

A model of the dynamics and thermodynamics of a plume of meltwater at the base of an ice shelf is presented. Such ice shelf water plumes may become supercooled and deposit marine ice if they rise (because of the pressure decrease in the in situ freezing temperature), so the model incorporates both melting and freezing at the ice shelf base and a multiple-size-class model of frazil ice dynamics and deposition. The plume is considered in two horizontal dimensions, so the influence of Coriolis forces is incorporated for the first time. It is found that rotation is extremely influential, with simulated plumes flowing in near-geostrophy because of the low friction at a smooth ice shelf base. As a result, an ice shelf water plume will only rise and become supercooled (and thus deposit marine ice) if it is constrained to flow upslope by topography. This result agrees with the observed distribution of marine ice under Filchner–Ronne Ice Shelf, Antarctica. In addition, it is found that the model only produces reasonable marine ice formation rates when an accurate ice shelf draft is used, implying that the characteristics of real ice shelf water plumes can only be captured using models with both rotation and a realistic topography.

Curl free pressure gradients over orography in a solution of the fully compressible Euler equations with implicit treatment of acoustic and gravity waves

Steep orography can cause noisy solutions and instability in models of the atmosphere. A new technique for modelling flow over orography is introduced which guarantees curl free gradients on arbitrary grids, implying that the pressure gradient term is not a spurious source of vorticity. This mimetic property leads to better hydrostatic balance and better energy conservation on test cases using terrain following grids. Curl-free gradients are achieved by using the co-variant components of velocity over orography rather than the usual horizontal and vertical components. In addition, gravity and acoustic waves are treated implicitly without the need for mean and perturbation variables or a hydrostatic reference profile. This enables a straightforward description of the implicit treatment of gravity waves. Results are presented of a resting atmosphere over orography and the curl-free pressure gradient formulation is advantageous. Results of gravity waves over orography are insensitive to the placement of terrain-following layers. The model with implicit gravity waves is stable in strongly stratified conditions, with N∆t up to at least 10 (where N is the Brunt-V ̈ais ̈al ̈a frequency). A warm bubble rising over orography is simulated and the curl free pressure gradient formulation gives much more accurate results for this test case than a model without this mimetic property.

Synoptic versus orographic control on stationary convective banding

Quasi-stationary convective bands can cause large localised rainfall accumulations and are often anchored by topographic features. Here, the predictability of and mechanisms causing one such band are determined using ensembles of the Met Office Unified Model at convection-permitting resolution (1.5 km grid length). The band was stationary over the UK for 3 h and produced rainfall accumulations of up to 34 mm. The amount and location of the predicted rainfall was highly variable despite only small differences between the large-scale conditions of the ensemble members. Only three of 21 members of the control ensemble produced a stationary rain band; these three had the weakest upstream winds and hence lowest Froude number. Band formation was due to the superposition of two processes: lee-side convergence resulting from flow around an upstream obstacle and thermally forced convergence resulting from elevated heating over the upstream terrain. Both mechanisms were enhanced when the Froude number was lower. By increasing the terrain height (thus reducing the Froude number), the band became more predictable. An ensemble approach is required to successfully predict the possible occurrence of such quasi-stationary convective events because the rainfall variability is largely modulated by small variations of the large-scale flow. However, high-resolution models are required to accurately resolve the small-scale interactions of the flow with the topography upon which the band formation depends. Thus, although topography provides some predictability, the quasi-stationary convective bands anchored by it are likely to remain a forecasting challenge for many years to come.

Impact of non-hydrostatic effects and trapped lee waves on mountain wave drag in directionally sheared flow

The orographic gravity wave drag produced in flow over an axisymmetric mountain when both vertical wind shear and non-hydrostatic effects are important was calculated using a semi-analytical two-layer linear model, including unidirectional or directional constant wind shear in a layer near the surface, above which the wind is constant. The drag behaviour is determined by partial wave reflection at the shear discontinuity, wave absorption at critical levels (both of which exist in hydrostatic flow), and total wave reflection at levels where the waves become evanescent (an intrinsically non-hydrostatic effect), which produces resonant trapped lee wave modes. As a result of constructive or destructive wave interference, the drag oscillates with the thickness of the constant-shear layer and the Richardson number within it (Ri), generally decreasing at low Ri and when the flow is strongly non-hydrostatic. Critical level absorption, which increases with the angle spanned by the wind velocity in the constant-shear layer, shields the surface from reflected waves, keeping the drag closer to its hydrostatic limit. While, for the parameter range considered here, the drag seldom exceeds this limit, a substantial drag fraction may be produced by trapped lee waves, particularly when the flow is strongly non-hydrostatic, the lower layer is thick and Ri is relatively high. In directionally sheared flows with Ri = O(1), the drag may be misaligned with the surface wind in a direction opposite to the shear, a behaviour which is totally due to non-trapped waves. The trapped lee wave drag, whose reaction force on the atmosphere is felt at low levels, may therefore have a distinctly different direction from the drag associated with vertically propagating waves, which acts on the atmosphere at higher levels.

Tide-topography interaction along the eastern Brazilian shelf

In the Tropics, continental shelves governed by western boundary currents are considered to be among the least productive ocean margins in the world, unless eddy-induced shelf-edge upwelling becomes significant. The eastern Brazilian shelf in the Southwest Atlantic is one of these, and since the slight nutrient input from continental sources is extremely oligotrophic. It is characterized by complex bathymetry with the presence of shallow banks and seamounts. In this work, a full three-dimensional nonlinear primitive equation ocean model is used to demonstrate that the interaction of tidal currents and the bottom topography of the east Brazil continental shelf is capable of producing local upwelling of South Atlantic Central Water, bringing nutrients up from deep waters to the surface layer. Such upper layer enrichment is found to be of significance in increasing local primary productivity. (c) 2005 Elsevier Ltd. All rights reserved.

Amazon Rainforest Exchange of Carbon and Subcanopy Air Flow: Manaus LBA Site-A Complex Terrain Condition

On the moderately complex terrain covered by dense tropical Amazon Rainforest (Reserva Biologica do Cuieiras-ZF2-02 degrees 36'17.1 '' S, 60 degrees 12'24.4 '' W), subcanopy horizontal and vertical gradients of the air temperature, CO2 concentration and wind field were measured for the dry and wet periods in 2006. We tested the hypothesis that horizontal drainage flow over this study area is significant and can affect the interpretation of the high carbon uptake rates reported by previous works at this site. A similar experimental design as the one by Tota et al. (2008) was used with a network of wind, air temperature, and CO2 sensors above and below the forest canopy. A persistent and systematic subcanopy nighttime upslope (positive buoyancy) and daytime downslope (negative buoyancy) flow pattern on a moderately inclined slope (12%) was observed. The microcirculations observed above the canopy (38 m) over the sloping area during nighttime presents a downward motion indicating vertical convergence and correspondent horizontal divergence toward the valley area. During the daytime an inverse pattern was observed. The microcirculations above the canopy were driven mainly by buoyancy balancing the pressure gradient forces. In the subcanopy space the microcirculations were also driven by the same physical mechanisms but probably with the stress forcing contribution. The results also indicated that the horizontal and vertical scalar gradients (e. g., CO2) were modulated by these micro-circulations above and below the canopy, suggesting that estimates of advection using previous experimental approaches are not appropriate due to the tridimensional nature of the vertical and horizontal transport locally. This work also indicates that carbon budget from tower-based measurement is not enough to close the system, and one needs to include horizontal and vertical advection transport of CO2 into those estimates.