Calculated wind climatology of the South-Saxonian/North-Czech mountain topography including improved resolution of mountains

A mesoscale model has been applied to calculate climatological means of the surface wind. A reliable average requires more than 40 model runs, which are differentiated by the direction and speed of the geostrophic wind under the assumption of neutral stratification. The frequency distributions of the geostrophic wind have been taken from observations of the 850-hPa winds at the radiosonde station in Prague for a 10-year period. The simulation results have been averaged over all sectors and speed classes of the geostrophic wind according to their frequencies. A comparison of the calculated mean wind speeds with observed ones shows deviations of about 0.4 ms⁻¹ outside the mountains. The representation of steep topography and isolated mountains on the basis of a 3-km horizontal resolution of the simulations needs special treatment in order to reduce the gap of up to 4 ms⁻¹ between observed and simulated mean wind speeds over mountains. Therefore, an empiric speed-up formula has been applied to the isolated mountains that otherwise would fall through the 3-km meshes. The corresponding deviations have been reduced to 1.5 ms⁻¹.     

Dynamic identification of moisture sources in the Orinoco basin in equatorial South America

Abstract

The main areas of net moisture uptake are examined in air masses over the Orinoco River basin, located in equatorial South America, north of the Amazon basin. Although the Orinoco River has the third largest average annual discharge in the world (with 5.4 × 10¹¹ m³ year^?1 draining into the Atlantic Ocean), the sources of moisture that feed it have not previously been studied in any detail. The results are presented from analyses of back-tracking of all the air masses over the Orinoco basin over a period of five years (2000–2004) using the diagnostic Lagrangian tool FLEXPART. The input data for the model were obtained from the European Centre for Medium-range Weather Forecasts (ECMWF). Air transported into the Orinoco basin experiences a large uptake of water over the tropical North Atlantic within the three days prior to its arrival over the basin. The Tropical South Atlantic and the eastern coast of the Pacific become significant moisture sources for about 5–10 days before arriving over the Orinoco basin. Contrary to what might be expected, large areas of the Amazon basin, along with the Gulf of Mexico, do not provide significant moisture to the study area. Interestingly, over these zones the air experiences net moisture loss. Preliminary analysis of the processes that occur leads to the conclusion that most of the water observed over the Orinoco basin derives from advective fluxes into the area, while recycling of moisture is negligible.     

Erosion modelling for land management in the Tahoe basin,USA: scaling from plots to forest catchments

Abstract

Recent developments in hydrological modelling of river basins are focused on prediction in ungauged basins, which implies the need to improve relationships between model parameters and easily-obtainable information, such as satellite images, and to test the transferability of model parameters. A large-scale distributed hydrological model is described, which has been used in several large river basins in Brazil. The model parameters are related to classes of physical characteristics, such as soil type, land use, geology and vegetation. The model uses two basin space units: square grids for flow direction along the basin and GRU—group response units—which are hydrological classes of the basin physical characteristics for water balance. Expected ranges of parameter values are associated with each of these classes during the model calibration. Results are presented of the model fitting in the Taquari-Antas River basin in Brazil (26 000 km² and 11 flow gauges). Based on this fitting, the model was then applied to the Upper Uruguay River basin (52 000 km²), having similar physical conditions, without any further calibration, in order to test the transferability of the model. The results in the Uruguay basin were compared with recorded flow data and showed relatively small errors, although a tendency to underestimate mean flows was found.     

The equatorial dynamics of a shallow,homogeneous ocean

Abstract

It is shown that the linear equatorial dynamics of a shallow ocean is characterized by two boundary layers of width γ^? L and γL (γ is the Ekman number of the flow, assumed small, and L is a horizontal dimension of the basin). In the γ^? layer stress in the bottom Ekman layer is comparable to that in the surface Ekman layer. In the γ layer vertical friction is important throughout the depth of the ocean. Should the Rossby number ? be so large as to invalidate a linear theory (? > γ^5/3), then inertial effects become important at a distance ?^2/5 L from the equator. The role played in the circulation of the basin by the non-linear equatorial current first studied by Charney (1960) is shown to be similar to that of the γ layer of the linear theory. Though lateral friction is unimportant in a linear model of the flow, shear layers at the equator are found to be a necessary feature of non-linear flow.     

Simulating strong ground motion from complex sources by reciprocal Green functions

We have developed a method to calculate site and path effects for complex heterogeneous media using synthetic Green’s functions. The Green’s functions are calculated numerically by imposing body forces at the site of interest and then storing the reciprocal Green’s functions along arbitrary finite-fault surfaces. By using reciprocal Green’s functions, we can then simulate many source scenarios for those faults because the primary numerical calculations need be done only once. The advantage of the proposed method is shown by evaluation of the site and path effects for three sites in the vicinity of the Los Angeles basin using the Southern California Velocity Model (version 2.2, Magistrale et al., 2000). In this example, we have simulated 300 source scenarios for 5 major southern California faults and compared their responses for period longer then 3 seconds at the selected sites. However, a more detailed comparison with strong motion records will be necessary before a particular hazard assessment can be made. For the tested source scenarios the results show that the variations in the peak velocity amplitudes and durations due to a source scenarios are as large as variations due to a heterogeneous velocity model.     

Motions at core surface in the geostrophic approximation

Motions at the top of the core which generate the observed Secular Variation (S.V.) field are computed. To reduce the well known ambiguity of the solution, two constraints are added: the flow is a large scale one and is geostrophic. The computed flow then has a very simple geometry; its poloïdal part is roughly axisymmetrical with respect to an equatorial diameter. This geometry is almost unchanged from 1970 to 1980 while the intensity of the velocity is nearly doubled.     

Evolution of a bifurcation in a meandering river with adjustable channel widths,Rhine delta apex,The Netherlands

Rivers may dramatically change course on a fluvial plain. Such an avulsion temporarily leads to two active channels connected at a bifurcation. Here we study the effect of dynamic meandering at the bifurcation and the effect of channel width adjustment to changing discharge in both downstream branches on the evolution of a bifurcation and coexisting channels. As an example, we reconstructed the last major avulsion at the Rhine delta apex. We combined historical and geological data to reconstruct a slowly developing avulsion process spanning 2000 years and involving channel width adjustment and meandering at the bifurcation. Based on earlier idealised models, we developed a one‐dimensional model for long‐term morphodynamic prediction of upstream channel and bifurcates connected at the bifurcation node. The model predicts flow and sediment partitioning at the node, including the effect of migrating meanders at the bifurcation and channel width adjustment. Bifurcate channel width adaptation to changing discharge partitioning dramatically slows the pacing of bifurcation evolution because the sediment balance for width adjustment and bed evolution are coupled. The model further shows that meandering at the bifurcation modulates channel abandonment or enlargement periodically. This explains hitherto unrecognised reactivation signals in the sedimentary record of the studied bifurcation meander belts, newly identified in our geological reconstruction. Historical maps show that bifurcation migration due to meander bend dynamics increases the bifurcation angle, which increases the rate of closure of one bifurcate. The combination of model and reconstruction identifies the relevant timescales for bifurcation evolution and avulsion duration. These are the time required to fill one downstream channel over one backwater length, the time to translate one meander wavelength downstream and, for strong river banks, the adaptation timescale to adjust channel width. The findings have relevance for all avulsions where channel width can adjust to changing discharge and where meandering occurs. Copyright © 2011 John Wiley & Sons, Ltd.     

Inertial oscillations and frontogenesis driven by a quadratic horizontal density variation

We consider inviscid rotating flow driven by a horizontally quadratic density variation in a horizontally unbounded slab. This configuration permits a similarity solution, removing the dependence on the horizontal coordinate from the vorticity and temperature equations, which are then solved by numerical integration along characteristics. At large values of Rossby number, the flow proceeds to a singularity in a similar manner to the non-rotating flow with the same initial conditions. At small values of Rossby number there are inertial oscillations of growing amplitude, which have been analysed using the method of multiple scales. The oscillations become desynchronised between the upper and lower parts of the domain, and static instability appears for a small fraction of each oscillation period. Eventually the oscillations give way to the rapid formation of a singularity, in contrast to geostrophic adjustment theory which predicts that a singularity will form only if the Rossby number is sufficiently large.     

旋转流场中的格子波耳兹曼模型

在大尺度的旋转流场中,由于存在哥氏力,使流体的流动出现了一系列复杂的动力学现象.在原格子Boltzmann模型研究的基础上,引入了哥氏力效应,发展了一个旋转流场中的格子Boltzmann模型.从该模型出发可导出地球流体力学方程,用这一模型对理想边界条件下的北半球大气环流进行了数值计算.数值结果很好地再现了大尺度地转流的流动特征.从理论和数值实验上验证了该模型的适用性.     

DYNAMIC AND KINEMATIC WAVE MODELS FOR 1-D RILL FLOW: A COMPARISON

1INTRODUCTIONInprocess-basedoverlandflowandsoilerosionmodels,surfacerunoffonahillslopeisoftenrepresentedaseitherbroadsheetfloworflowinrillswithassumedrectangularchannelcrosssections(e.g.,Bairdelal1992,NSERL1995).Inmostcasesthehydraulicsofoverlandflowiscalculatedbyusingthekinematicwavemodel,whichisasirnplificationofthedynamicwavemodel(theequationsystemofSaint-Venantequationandequationofcontinuity).Forabroadsheetoverlandflowonhillslopeduetorainfallexcess,theequationofcontinuityiswherehis…     

Assessment of flooding in urbanized ungauged basins: a case study in the Upper Tiber area,Italy

The reliability of a procedure for investigation of flooding into an ungauged river reach close to an urban area is investigated. The approach is based on the application of a semi‐distributed rainfall–runoff model for a gauged basin, including the flood‐prone area, and that furnishes the inlet flow conditions for a two‐dimensional hydraulic model, whose computational domain is the urban area. The flood event, which occurred in October 1998 in the Upper Tiber river basin and caused significant damage in the town of Pieve S. Stefano, was used to test the approach. The built‐up area, often inundated, is included in the gauged basin of the Montedoglio dam (275 km²), for which the rainfall–runoff model was adapted and calibrated through three flood events without over‐bank flow. With the selected set of parameters, the hydrological model was found reasonably accurate in simulating the discharge hydrograph of the three events, whereas the flood event of October 1998 was simulated poorly, with an error in peak discharge and time to peak of −58% and 20%, respectively. This discrepancy was ascribed to the combined effect of the rainfall spatial variability and a partial obstruction of the bridge located in Pieve S. Stefano. In fact, taking account of the last hypothesis, the hydraulic model reproduced with a fair accuracy the observed flooded urban area. Moreover, incorporating into the hydrological model the flow resulting from a sudden cleaning of the obstruction, which was simulated by a ‘shock‐capturing’ one‐dimensional hydraulic model, the discharge hydrograph at the basin outlet was well represented if the rainfall was supposed to have occurred in the region near the main channel. This was simulated by reducing considerably the dynamic parameter, the lag time, of the instantaneous unit hydrograph for each homogeneous element into which the basin is divided. The error in peak discharge and time to peak decreased by a few percent. A sensitivity analysis of both the flooding volume involved in the shock wave and the lag time showed that this latter parameter requires a careful evaluation. Moreover, the analysis of the hydrograph peak prediction due to error in rainfall input showed that the error in peak discharge was lower than that of the same input error quantity. Therefore, the obtained results allowed us to support the hypothesis on the causes which triggered the complex event occurring in October 1998, and pointed out that the proposed procedure can be conveniently adopted for flood risk evaluation in ungauged river basins where a built‐up area is located. The need for a more detailed analysis regarding the processes of runoff generation and flood routing is also highlighted. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulating the surface waters of the Amazon River basin: impacts of new river geomorphic and flow parameterizations

This paper describes the impacts of new river geomorphic and flow parameterizations on the simulated surface waters dynamics of the Amazon River basin. Three major improvements to a hydrologic model are presented: (1) the river flow velocity equation is expanded to be dependent on river sinuosity and friction in addition to gradient forces; (2) equations defining the morphological characteristics of the river, such as river height, width and bankfull volume, are derived from 31 622 measurements of river morphology and applied within the model; (3) 1 km resolution topographic data from the Shuttle Radar Topography Mission (SRTM) are used to provide physically based fractional flooding of grid cells from a statistical representation of sub‐grid‐scale floodplain morphology. The discharge and floodplain inundation of the Amazon River is simulated for the period 1968–1998, validated against observations, and compared with results from a previous version of the model. These modifications result in considerable improvement in the simulations of the hydrological features of the Amazon River system. The major impact is that the average wet‐season flooded area on the Amazon mainstem for the period 1983–1988 is now within 5% of satellite‐derived estimates of flooded area, whereas the previous model overestimates the flooded area by about 80%. The improvements are a consequence of the new empirical river geomorphologic functions and the SRTM topography. The new formulation of the flow velocity equation results in increased river velocity on the mainstem and major tributaries and a better correlation between the mean monthly simulated and observed discharge. Copyright © 2007 John Wiley & Sons, Ltd.     

The mass transport velocity induced by free oscillations at a single frequency

Abstract

Second order effects due to the presence of a first order free oscillation at a single frequency in a variable depth rotating ocean are examined. It is found that the second order Lagrangian mean velocity (mass transport velocity) satisfies the linearized equations for unforced steady geostrophic motion. This implies that if the ocean basin is laterally bounded and contains no closed geostrophic contours, the second order Lagrangian mean velocity vanishes everywhere.     

Two-dimensional simulations of katabatic layers observed during the GIMEX experiment

The hydrostatic model SALSA is used to simulate a particular event observed during the Greenland Ice Margin EXperiment “GIMEX” (on July 12th, 1991). The time evolution of the large-scale flow was incorporated in the model through time dependent boundary conditions which were updated using the closest upwind sounding. A turbulent scheme for the stable boundary layer and an appropriate parametrization of the surface fluxes implemented in the same model, are used for this study. The simulation results are discussed and compared to the available observations. The computed turbulent fluxes are correctly estimated. The model predicts a mixing zone of about 1500 m high which is in good agreement with tundra site observations. Over the ice cap, the katabatic layer is correctly simulated by the model. Its height of 80–300 m is well estimated. The comparison between the simulation and observations taken at ice cap sites is reasonably valid. The ablation computed along the ice cap corresponds well to the values reconstructed of observations at sites 4 and 9. Finally, a sensibility study to a specified westward geostrophic wind (2 ms⁻¹) shows that the consideration of this latter improves the simulated tundra wind evolution.     

Spatial scale effects on sediment concentration in runoff during flood events for hilly areas of the Loess Plateau,China

The spatial scale effect on sediment concentration in runoff has received little attention despite numerous studies on sediment yield or sediment delivery ratio in the context of multiple spatial scales. We have addressed this issue for hilly areas of the Loess Plateau, north China where fluvial processes are mainly dominated by hyperconcentrated flows. The data on 717 flow events observed at 17 gauging stations and two runoff experimental plots, all located in the 3906 km² Dalihe watershed, are presented. The combination of the downstream scour of hyperconcentrated flows and the downstream dilution, which is mainly caused by the base flow and is strengthened as a result of the strong patchy storms, determines the spatial change of sediment concentration in runoff during flood events. At the watershed scale, the scouring effect takes predominance first but is subordinate to the downstream dilution with a further increase in spatial scale. As a result, the event mean sediment concentration first increases following a power function with drainage basin area and then declines at the drainage basin area of about 700 km². The power function in combination with the proportional model of the runoff‐sediment yield relationship we proposed before was used to establish the sediment‐yield model, which is neither the physical‐based model nor the regression model. This model, with only two variables (runoff depth and drainage basin area) and two parameters, can provide fairly accurate prediction of event sediment yield with model efficiency over 0·95 if small events with runoff depth lower than 1 mm are excluded. Copyright © 2011 John Wiley & Sons, Ltd.     

Quasi-geostrophic motions in a galactic flow

Summary We continue the theoretical investigation of geostrophically balanced motions of a star gas in a rotating galaxy. The method of small perturbations is employed to derive a quasi geostrophic vorticity equation suitable for solving initial value problems for slow galactic eddying motions. The results derived for this simple theoretical model are then generalized to obtain approximate geostrophic forecast equations suitable for study of nonlinear galactic motions over a whole disk galaxy. A simple discussion of the stability of our approximate vorticity equation is presented. It is expected that in the absence of rapid gravitational instabilities, actual eddy motions observed in spiral galaxies will be at least approximately of the type of motion discussed in this paper. The Rossby adjustment problem is investigated in this context.     

Detecting the flow relationships between deep and shallow aquifers in an exploited groundwater system,using long‐term monitoring data and quantitative hydrogeology: the Acque Albule basin case (Rome,Italy)

Five years of hydrogeological monitoring and field activities performed in the complex hydrogeological system of the Acque Albule basin (AAB) were conducted to define the hydrogeological setting, the relationship between deep and shallow aquifers and a conceptual groundwater flow model of this exploited area using conventional quantitative techniques. The basin, which is located close to Rome (Italy) on the west side of the Apennine chain and just north of the Colli Albani volcano, subsided after development of a north–south fault system (about 115 000 y bp). The AAB experiences intense hydrothermal activity, which has produced a large travertine deposit (80‐m thick). The travertine deposit constitutes a fractured aquifer that is the final destination of more than 5 m³ s^‐1 of water and is strongly dewatered by quarry activities. The complex hydrogeology of this basin was investigated, revealing two main hydraulically connected aquifers, one thermalised and partly confined into the limestone bedrock and one unconfined in the travertine. The two aquifers are separated by a non‐continuous clayey aquiclude. The hydrogeological survey and geological characterisation contributed to the development of the groundwater flow conceptual model. Analysis and comparison of the monitored levels highlighted the pattern of flow between the deep and shallow parts of the flow system. Copyright © 2012 John Wiley & Sons, Ltd.     

Lateral dynamic analysis and classification of barotropic tidal inlets

The dynamical balances at shallow tidal inlets are highly nonlinear, and can vary substantially over sub-kilometer scales. In this study, barotropic dynamics are examined with numerical experiments on a series of idealized inlets with differing inlet widths and lengths. Circulation and elevation fields obtained from fully nonlinear depth-integrated circulation models are used to reconstruct the contribution of each term in the momentum equations. Momentum terms are rotated into a streamline coordinate system to simplify interpretation of the dynamics. Spatial patterns in momentum reveal that the lateral balances at inlets can vary from nearly geostrophic to strongly cyclostrophic. Marked dynamical differences are seen between inlets with different lengths and widths. Inlet regions of geostrophic or cyclostrophic balances can be predicted using two dimensionless parameters, the dynamic length L^* and dynamic width W^*. A classification scheme is proposed using L^* and W^* to compare the idealized inlets analyzed here with inlets from 20 previous studies. Four distinct inlet types are identified and discussed.     

Drainage area and the variation of channel geometry downstream

Recent geomorphological studies tend to deal with small basins. The understanding of small basin dynamics provides important information for the understanding of large basin dynamics assuming that the extrapolation of small basin data to larger basins is valid. This work tests the validity of this extrapolation of data with reference to channel geometry. An analysis of the variation of channel width downstream reveals that the value b =0.5 (W = aQ^b) is a ‘good’ average. However, the use of a one-line model consisting of a simple power function incurs a loss of a considerable amount of relevant information concerning the channel form and hence the channel processes. It has been shown that the –b– value for small basins and very big basins is lower than the one for the intermediate basins.