Flash floods and their effects on the development in El-Qaá plain area in South Sinai, Egypt

The arid regions are dominated to a much larger degree than humid regions by major catastrophic events. Although most of Egypt lies within the great hot desert belt; it experiences especially in the north some torrential rainfall, which causes flash floods all over Sinai Peninsula. Flash floods in hot deserts are characterized by high velocity and low duration with a sharp discharge peak. Large sediment loads may be carried by floods threatening fields and settlements in the wadis and even people who are living there. The extreme spottiness of rare heavy rainfall, well known to desert people everywhere, precludes any efficient forecasting. Thus, although the limitation of data still reflects pre-satellite methods, chances of developing a warning system for floods in the desert seem remote. The relatively short flood-to-peak interval, a characteristic of desert floods, presents an additional impediment to the efficient use of warning systems. The present thesis contains introduction and five chapters, chapter one points out the physical settings of the study area. There are the geological settings such as outcrop lithology of the study area and the deposits. The alluvial deposits of Wadi Moreikh had been analyzed using OSL dating to know deposits and palaeoclimatic conditions. The chapter points out as well the stratigraphy and the structure geology containing main faults and folds. In addition, it manifests the pesent climate conditions such as temperature, humidity, wind and evaporation. Besides, it presents type of soils and natural vegetation cover of the study area using unsupervised classification for ETM+ images. Chapter two points out the morphometric analysis of the main basins and their drainage network in the study area. It is divided into three parts: The first part manifests the morphometric analysis of the drainage networks which had been extracted from two main sources, topographic maps and DEM images. Basins and drainage networks are considered as major influencing factors on the flash floods; Most of elements were studied which affect the network such as stream order, bifurcation ratio, stream lengths, stream frequency, drainage density, and drainage patterns. The second part of this chapter shows the morphometric analysis of basins such as area, dimensions, shape and surface. Whereas, the third part points the morphometric analysis of alluvial fans which form most of El-Qaá plain. Chapter three manifests the surface runoff through rainfall and losses analysis. The main subject in this chapter is rainfall which has been studied in detail; it is the main reason for runoff. Therefore, all rainfall characteristics are regarded here such as rainfall types, distribution, rainfall intensity, duration, frequency, and the relationship between rainfall and runoff. While the second part of this chapter concerns with water losses estimation by evaporation and infiltration which are together the main losses with direct effect on the high of runoff. Finally, chapter three points out the factors influencing desert runoff and runoff generation mechanism. Chapter four is concerned with assessment of flood hazard, it is important to estimate runoff and tocreate a map of affected areas. Therefore, the chapter consists of four main parts; first part manifests the runoff estimation, the different methods to estimate runoff and its variables such as runoff coefficient lag time, time of concentration, runoff volume, and frequency analysis of flash flood. While the second part points out the extreme event analysis. The third part shows the map of affected areas for every basin and the flash floods degrees. In this point, it has been depending on the DEM to extract the drainage networks and to determine the main streams which are normally more dangerous than others. Finally, part four presets the risk zone map of total study area which is of high inerest for planning activities. Chapter five as the last chapter concerns with flash flood Hazard mitigation. It consists of three main parts. First flood prediction and the method which can be used to predict and forecast the flood. The second part aims to determine the best methods which can be helpful to mitigate flood hazard in the arid zone and especially the study area. Whereas, the third part points out the development perspective for the study area indicating the suitable places in El-Qaá plain for using in economic activities.

Headward retreat of streams in the Late Oligocene to Early Miocene Swiss Alps

This study uses the widths, the spacing and the grain-size pattern of Oligo/Miocene alluvial fan conglomerates in the central segment of the Swiss Alpine foreland to reconstruct the topographic development of the Alps. These data are analysed with models of longitudinal stream profile development, to propose that the Alpine topography evolved from an early transient state where streams adjusted to rock uplift by headward retreat, to a mature phase where any changes in rock uplift were accommodated by vertical incision. The first stage comprises the time interval between ca 31 Ma and 22 Ma, when the Alpine streams deposited many small fans with a lateral spacing of <30 km in the north Alpine foreland. As the range evolved, the streams joined and the fans coalesced into a few large depositional systems with a lateral spacing of ca 80 to 100 km at 22 Ma. The models used here suggest that the overall elevation of the Alps increased rapidly within <5 Myr. The variability in pebble size increased either due to variations in sediment supply, enhanced orographic effects, or preferentially due to a change towards a stormier palaeoclimate. By 22 Ma, only two large rivers carried material into the foreland fans, suggesting that the major Alpine streams had established themselves. This second phase of stable drainage network was maintained until ca 5 Ma, when the uplift and erosion of the Molasse started and streams were redirected both in the Alps and in the foreland. This study illustrates that sedimentological archives of foreland basins can be used to reconstruct the chronology of the topographic development of mountain belts. It is suggested that the finite elevation of mountainous landscapes is reached early during orogeny and can be maintained for millions of years, provided that erosion is efficient.

Geomorphology and evolution of the late Pleistocene to Holocene fluvial system in the south-eastern Llanos de Moxos, Bolivian Amazon

In the Bolivian Amazon several paleochannel generations are preserved. Their wide spectrum of morphologies clearly provides crucial information on the type and magnitude of geomorphic and hydrological changes within the drainage network of the Andean foreland. Therefore, in this study we mapped geomorphological characteristics of paleochannels, and applied radiocarbon and optically stimulated luminescence dating. Seven paleochannel generations are identified. Significant changes in sinuosity, channel widths and river pattern are observed for the successive paleochannel generations. Our results clearly reflect at least three different geomorphic and hydrological periods in the evolution of the fluvial system since the late Pleistocene. Changes in discharge and sediment load may be controlled by combinations of two interrelated mechanisms: (i) spatial changes and re-organizations of the drainage network in the upper catchment, and/or (ii) climate changes with their associated local to catchment-scale modifications in vegetation cover, and changes in discharge, inundation frequencies and magnitudes, which have likely affected the evolution of the fluvial system in the Llanos de Moxos. In summary, our study has revealed the enormous potential which geomorphic mapping and analysis combined with luminescence based chronologies hold for the reconstruction of the late Pleistocene to recent fluvial system in a large portion of Amazonia.

A high-resolution map of direct and indirect connectivity of erosion risk areas to surface waters in Switzerland: A risk assessment tool for planning and policy-making

Off-site effects of soil erosion are becoming increasingly important, particularly the pollution of surface waters. In order to develop environmentally efficient and cost effective mitigation options it is essential to identify areas that bear both a high erosion risk and high connectivity to surface waters. This paper introduces a simple risk assessment tool that allows the delineation of potential critical source areas (CSA) of sediment input into surface waters concerning the agricultural areas of Switzerland. The basis are the erosion risk map with a 2 m resolution (ERM2) and the drainage network, which is extended by drained roads, farm tracks, and slope depressions. The probability of hydrological and sedimentological connectivity is assessed by combining soil erosion risk and extended drainage network with flow distance calculation. A GIS-environment with multiple-flow accumulation algorithms is used for routing runoff generation and flow pathways. The result is a high resolution connectivity map of the agricultural area of Switzerland (888,050 ha). Fifty-five percent of the computed agricultural area is potentially connected with surface waters, 45% is not connected. Surprisingly, the larger part of 34% (62% of the connected area) is indirectly connected with surface waters through drained roads, and only 21% are directly connected. The reason is the topographic complexity and patchiness of the landscape due to a dense road and drainage network. A total of 24% of the connected area and 13% of the computed agricultural area, respectively, are rated with a high connectivity probability. On these CSA an adapted land use is recommended, supported by vegetated buffer strips preventing sediment load. Even areas that are far away from open water bodies can be indirectly connected and need to be included in planning of mitigation measures. Thus, the connectivity map presented is an important decision-making tool for policy-makers and extension services. The map is published on the web and thus available for application.

Calidad del agua en el área regadía del Río Mendoza (Argentina)

El río Mendoza conforma el oasis norte que es el más importante de la provincia. El crecimiento urbano ha avanzado sobre áreas originalmente agrícolas, rodeando la red de canales y desagües, que también recibe los desagües pluviales urbanos, producto de tormentas convectivas. La actividad antropogénica utiliza el recurso para bebida, saneamiento, riego, recreación, etc., y vuelca sus excedentes a la red, contaminándola. Para conocer la calidad del agua de esta cuenca se seleccionaron, estratégicamente, 15 sitios de muestreo: 3 a lo largo del río y a partir del dique derivador Cipolletti (R_I a R_III), 5 en la red de canales (C_I a C_V) y 7 ubicados en los colectores de drenaje (D_I a D_VII). Se realizaron los siguientes análisis físico-químicos y microbiológicos; en el río y en la red de canales: conductividad eléctrica, temperatura, pH, aniones y cationes (cálculo de RAS), oxígeno disuelto (OD), sólidos sedimentables, demanda química de oxígeno (DQO), bacterias aerobias mesófilas (BAM), coliformes totales y fecales y metales pesados. En la red de drenaje sólo se realizaron los cuatro primeros. Los resultados de los análisis, se incorporaron a una base de datos y se sometieron a un análisis estadístico descriptivo e inferencial. Este último consistió en la aplicación de diversas pruebas en busca de posibles diferencias entre los sitios de muestreo, para cada variable respuesta, a un α = 0.05. Se realizó el análisis de la varianza de efectos fijos y de efectos aleatorios y se probaron los supuestos de homocedasticidad y de normalidad de los errores. En el caso de violación de los supuestos, se utilizó la prueba de Kruskal- Wallis. Se compararon los siguientes sitios de muestreo entre sí: ríos, R_I-canales y drenajes. Se concluyó que hay un aumento significativo de la salinidad y la sodicidad en R_II, que los cambios de calidad ocurridos entre R_II y R_III podrían deberse al aporte de otras aguas. Con respecto a la comparación de los parámetros entre la cabeza del sistema (R_I) y la red de canales se puede decir que los aportes realizados por los escurrimientos urbanos ubicados hacia el oeste del canal Cacique Guaymallén, sumados a los vuelcos de Campo Espejo (detectados en C_II), incrementan significativamente la salinidad (+55 %) y sodicidad del agua (+95 %) respecto del punto R_I, aunque el valor de sodicidad sigue siendo bajo. También se han encontrado incrementos de salinidad (+80 %), de DQO (+1159 %) y BAM (+2873 %) con lógica disminución de OD (-58 %) en el punto C_V (canal Auxiliar Tulumaya) respecto del punto R_I, ocasionados por aportes urbanos (Gran Mendoza) sumados a la carga contaminante del canal Pescara. Los metales pesados no presentan grandes diferencias entre sitios de muestreo.

Water balance measurements on Samoylov Island during campaign Lena2011

The summer water balance of a typical Siberian polygonal tundra catchment is investigated in order to identify the spatial and temporal dynamics of its main hydrological processes. The results show that, besides precipitation and evapotranspiration, lateral flow considerably influences the site-specific hydrological conditions. The prominent microtopography of the polygonal tundra strongly controls lateral flow and storage behaviour of the investigated catchment. Intact rims of low-centred polygons build hydrological barriers, which release storage water later in summer than polygons with degraded rims and troughs above degraded ice wedges. The barrier function of rims is strongly controlled by soil thaw, which opens new subsurface flow paths and increases subsurface hydrological connectivity. Therefore, soil thaw dynamics determine the magnitude and timing of subsurface outflow and the redistribution of storage within the catchment. Hydraulic conductivities in the elevated polygonal rims sharply decrease with the transition from organic to mineral layers. This interface causes a rapid shallow subsurface drainage of rainwater towards the depressed polygon centres and troughs. The re-release of storage water from the centres through deeper and less conductive layers helps maintain a high water table in the surface drainage network of troughs throughout the summer.

Late Weichselian fluvial evolution of the southern Kara Sea

Glaciations had a profound impact on the global sea-level and particularly on the Arctic environments. One of the key questions related to this topic is, how did the discharge of the Siberian Ob and Yenisei rivers interact with a proximal ice sheet? In order to answer this question high-resolution (1-12 kHz), shallow-penetration seismic profiles were collected on the passive continental margin of the Kara Sea Shelf to study the paleo-drainage pattern of the Ob and Yenisei rivers. Both rivers incised into the recent shelf, leaving filled and unfilled river channels and river canyons/valleys connecting to a complex paleo-drainage network. These channels have been subaerially formed during a regressive phase of the global sea-level during the Last Glacial Maximum. Beyond recent shelf depths of 120 m particle transport is manifested in submarine channel-levee complexes acting as conveyor for fluvial-derived fines. In the NE area, uniform draping sediments are observed. Major morphology determining factors are (1) sea-level fluctuations and (2) LGM ice sheet influence. Most individual channels show geometries typical for meandering rivers and appear to be an order of magnitude larger than recent channel profiles of gauge stations on land. The Yenisei paleo-channels have larger dimensions than the Ob examples and could be originated by additional water release during the melt of LGM Putoran ice masses. Asymmetrical submarine channel-levee complexes with channel depths of 60 m and more developed, in some places bordered by glacially dominated morphology, implying deflection by the LGM ice masses. A total of more than 12,000 km of acoustic profiles reveal no evidence for an ice-dammed lake of greater areal extent postulated by several workers. Furthermore, the existence of the channel-levee complexes is indicative of unhindered sediment flow to the north. Channels situated on the shelf above 120-m water depth exhibit no phases of ponding and or infill during sea-level lowstand. These findings denote the non-existence of an ice sheet on large areas of the Kara Sea shelf.

A GIS-based methodology to quantitatively define an Adjacent Protected Area in a shallow karst cavity: The case of Altamira cave

Different types of land use are usually present in the areas adjacent to many shallow karst cavities. Over time, the increasing amount of potentially harmful matter and energy, of mainly anthropic origin or influence, that reaches the interior of a shallow karst cavity can modify the hypogeal ecosystem and increase the risk of damage to the Palaeolithic rock art often preserved within the cavity. This study proposes a new Protected Area status based on the geological processes that control these matter and energy fluxes into the Altamira cave karst system. Analysis of the geological characteristics of the shallow karst system shows that direct and lateral infiltration, internal water circulation, ventilation, gas exchange and transmission of vibrations are the processes that control these matter and energy fluxes into the cave. This study applies a comprehensive methodological approach based on Geographic Information Systems (GIS) to establish the area of influence of each transfer process. The stratigraphic and structural characteristics of the interior of the cave were determined using 3D Laser Scanning topography combined with classical field work, data gathering, cartography and a porosity–permeability analysis of host rock samples. As a result, it was possible to determine the hydrogeological behavior of the cave. In addition, by mapping and modeling the surface parameters it was possible to identify the main features restricting hydrological behavior and hence direct and lateral infiltration into the cave. These surface parameters included the shape of the drainage network and a geomorphological and structural characterization via digital terrain models. Geological and geomorphological maps and models integrated into the GIS environment defined the areas involved in gas exchange and ventilation processes. Likewise, areas that could potentially transmit vibrations directly into the cave were identified. This study shows that it is possible to define a Protected Area by quantifying the area of influence related to each transfer process. The combined maximum area of influence of all the processes will result in the new Protected Area. This area will thus encompass all the processes that account for most of the matter and energy carried into the cave and will fulfill the criteria used to define the Protected Area. This methodology is based on the spatial quantification of processes and entities of geological origin and can therefore be applied to any shallow karst system that requires protection.

Age determinaton and age model of sediments from the Bay of Biscay

Continuous high-resolution mass accumulation rates (MAR) and X-ray fluorescence (XRF) measurements from marine sediment records in the Bay of Biscay (NE Atlantic) have allowed the determination of the timing and the amplitude of the 'Fleuve Manche' (Channel River) discharges during glacial stages MIS 10, MIS 8, MIS 6 and MIS 4-2. These results have yielded detailed insight into the Middle and Late Pleistocene glaciations in Europe and the drainage network of the western and central European rivers over the last 350 kyr. This study provides clear evidence that the 'Fleuve Manche' connected the southern North Sea basin with the Bay of Biscay during each glacial period and reveals that 'Fleuve Manche' activity during the glaciations MIS 10 and MIS 8 was significantly less than during MIS 6 and MIS 2. We correlate the significant 'Fleuve Manche' activity, detected during MIS 6 and MIS 2, with the extensive Saalian (Drenthe Substage) and the Weichselian glaciations, respectively, confirming that the major Elsterian glaciation precedes the glacial MIS 10. In detail, massive 'Fleuve Manche' discharges occurred at ca 155 ka (mid-MIS 6) and during Termination I, while no significant discharges are found during Termination II. It is assumed that a substantial retreat of the European ice sheet at ca 155 kyr, followed by the formation of ice-free conditions between the British Isles and Scandinavia until Termination II, allowed meltwater to flow northwards through the North Sea basin during the second part of the MIS 6. We assume that this glacial pattern corresponds to the Warthe Substage glacial maximum, therefore indicating that the data presented here equates to the Drenthe and the Warthe glacial advances at ca 175-160 ka and ca 150-140 ka, respectively. Finally, the correlation of our records with ODP site 980 reveals that massive 'Fleuve Manche' discharges, related to partial or complete melting of the European ice masses, were synchronous with strong decreases in both the rate of deep-water formation and the strength of the Atlantic thermohaline circulation. 'Fleuve Manche' discharges over the last 350 kyr probably participated, with other meltwater sources, in the collapse of the thermohaline circulation by freshening the northern Atlantic surface water.

As Escadarias de Terraços do Ródão à Chamusca (Baixo Tejo) - Caracterização e Interpretação de Dados Sedimentares, Tectónicos, Climáticos e do Paleolítico

As unidades estratigráficas que resultaram da evolução do rio Tejo em Portugal, aqui analisadas em pormenor entre Vila Velha de Ródão e Chamusca, possuem distintas características sedimentares e indústrias líticas: uma unidade culminante do enchimento sedimentar (o ancestral Tejo, antes do início da etapa de incisão fluvial) - SLD13 (+142 a 262 m acima do leito actual; com provável idade 3,6 a 1,8 Ma), sem indústrias identificadas; terraço T1 (+76 a 180 m; ca. 1000? a 900 ka), sem indústrias; terraço T2 (+57 a 150 m; idade estimada em ca. 600 ka), sem indústrias; terraço T3 (+36 a 113 m; ca. 460 a 360? ka), sem indústrias; terraço T4 (+26 a 55 m; ca. 335 a 155 ka), Paleolítico Inferior (Acheulense) em níveis da base e intermédios mas Paleolítico Médio inicial em níveis do topo; terraço T5 (+5 a 34 m; 135 a 73 ka), Paleolítico Médio (com talhe Mustierense, Levallois); terraço T6 (+3 a 14 m; 62 a 32 ka), Paleolítico Médio final (Mustierense final); Areias da Carregueira (areias eólicas) e coluviões (+3 a ca. 100 m; 32 a 12 ka), Paleolítico Superior a Epipaleolítico; enchimento da planície aluvial (+0 a 8 m; ca. 12 ka a actual), Mesolítico e indústrias mais recentes. As diferenças na elevação (a.r.b.) das escadarias de terraços resultam de soerguimento diferencial, devido a falhas ativas. Numa dada escadaria datada, a projeção da elevação da superfície de cada terraço (a.r.b.) versus a sua idade permitiu estimar a idade do topo do terraço T2 (ca. 600 ka) e a provável idade do início da etapa de incisão (ca. 1,8 Ma). Obteve-se a duração da fase de agradação dos terraços baixos e médios: T6 – 30 ka; T5 – 62 ka; T4 – ca. 180 ka; T3 – ca. 100? ka. Conclui-se que durante o Plistocénico médio e final, as fases de incisão e alargamento do vale foram curtas (ca. 11-25 ka) e ocorreram durante períodos de nível do mar muito baixo, alternando com mais longas fases de inundação e agradação do vale durante níveis do mar mais altos. Estas oscilações eustáticas de causa climática estão sobrepostas a um contexto de soerguimento de longo termo, controlando o desenvolvimento das escadarias. Calculou-se que para os últimos ca. 155 ka as taxas de incisão de curto-termo apresentam valores (0,09 a 0,41 m/ka), aproximadamente, duplos dos calculados para o intervalo ca. 155 a 900 ka (0,04 a 0,28 m/ka). Este aumento na taxa de incisão deve estar relacionado com um aumento na taxa de soerguimento por intensificação da compressão devido à convergência entre as placas Africana e Eurasiática. Abstract: The terrace staircases of the Lower Tagus River (Ródão to Chamusca) – characterization and interpretation of the sedimentary, tectonic, climatic and Palaeolithic data The stratigraphic units that record the evolution of the Tagus River in Portugal (study area between Vila Velha de Ródão and Chamusca villages) have different sedimentary characteristics and lithic industries: a culminant sedimentary unit (the ancestral Tagus, before the drainage network entrenchment) – SLD13 (+142 to 262 m above river bed – a.r.b.; with probable age 3.6 to 1.8 Ma), without artefacts; T1 terrace (+76 to 180 m; ca. 1000? to 900 ka), without artefacts; T2 terrace (+57 to 150 m; top deposits with a probable age ca. 600 ka), without artefacts; T3 terrace (+36 to 113 m; ca. 460 to 360? ka), without artefacts; T4 terrace (+26 to 55 m; ca. 335 a 155 ka), Lower Paleolithic (Acheulian) at basal and middle levels but early Middle Paleolithic at top levels; T5 terrace (+5 to 34 m; 135 to 73 ka), Middle Paleolithic (Mousterian; Levallois technique); T6 terrace (+3 to 14 m; 62 to 32 ka), late Middle Paleolithic (late Mousterian); Carregueira Sands (aeolian sands) and colluvium (+3 a ca. 100 m; 32 to 12 ka), Upper Paleolithic to Epipaleolithic; alluvial plain (+0 to 8 m; ca. 12 ka to present), Mesolithic and more recent industries. The differences in elevation (a.r.b.) of the several terrace staircases results from differential uplift due to active faults. The age interval for each aggradation phase of T3 to T6 terraces was obtained: T3 – ca. 100? ka; T4 – ca. 180 ka; T5 – 62 ka; T6 – 30 ka. The intervals of river down-cutting and widening of the valley floor were short (ca. 11-25 ka) and coincided with periods of very low sea-level. The plotting of the elevation (a.r.b.) versus the age of each terrace surface allows to estimate the age of the T2 terrace (ca. 600 ka) and the probable age of the beginning of the incision stage (ca. 1.8 Ma). So, the high amplitude sea-level changes that characterized the Middle and Late Pleistocene strongly determined the episodic down-cutting phases of the river during the low stands sea levels that alternated with the flooding and aggradation phases of the incised valley during highstand sea levels. These climate related eustatic oscillations are superimposed onto a long term uplift pattern, controlling the river terrace staircase development. During the last ca. 155 ka, the short-term incision rates (0.09 a 0.41 m/ka) were twice the values determined for the interval 155 to 900 ka (0.04 to 0.28 m/ka). This increase in incision rate should be related with an increase in uplift rate resulting from an intensification of compression due to the convergence between African - Eurasian plates.

Compartimentação fisiográfica para análise ambiental do potencial erosivo a partir das propriedades da rede de drenagem

O presente trabalho teve como objetivo compartimentar a sub-bacia do Baixo Rio Piracicaba - SP, em unidades homogêneas quanto à potencialidade à erosão visando a subsidiar o gerenciamento ambiental. A importância desta pesquisa concentra-se no intenso desenvolvimento dos processos erosivos na área em foco e na sua importância socioeconômica em níveis estadual e nacional. O procedimento adotado para atingir o objetivo foi a análise da rede de drenagem e dos lineamentos obtidos pelas imagens TM/Landsat-5. O resultado obtido foi a divisão da área em quatro compartimentos quanto à potencialidade à erosão: muito alta, alta, média e baixa. Concluiu-se que a área estudada é heterogênea, com regiões sujeitas a diferentes intensidades de processos erosivos e que a sistemática adotada se mostrou eficiente para caracterizar as compartimentações fisiográficas da potencialidade erosiva, e essas compartimentações podem e devem ser utilizadas como ponto de partida para estudos ambientais e de utilização do território, em consonância com o desenvolvimento sustentável.

NEW DATA ON THE CHRONOLOGY OF THE VALE DO FORNO SEDIMENTARY SEQUENCE (LOWER TAGUS RIVER TERRACE STAIRCASE) AND ITS RELEVANCE AS FLUVIAL ARCHIVE OF THE MIDDLE PLEISTOCENE IN WESTERN IBERIA

NEW DATA ON THE CHRONOLOGY OF THE VALE DO FORNO SEDIMENTARY SEQUENCE (LOWER TAGUS RIVER TERRACE STAIRCASE) AND ITS RELEVANCE AS FLUVIAL ARCHIVE OF THE MIDDLE PLEISTOCENE IN WESTERN IBERIA Pedro P. Cunha 1, António A. Martins 2, Jan-Pieter Buylaert 3,4, Andrew S. Murray 4, Luis Raposo 5, Paolo Mozzi 6, Martin Stokes 7 1 MARE - Marine and Environmental Sciences Centre, Department of Earth Sciences, University of Coimbra, Portugal: pcunha@dct.uc.pt 2 MARE - Marine and Environmental Sciences Centre, Dep. Geociências, University of Évora, Portugal; aam@uevora.pt 3 Centre for Nuclear Technologies, Technical University of Denmark, Risø Campus, Denmark; jabu@dtu.dk 4 Nordic Laboratory for Luminescence Dating, Aarhus University, Risø DTU, Denmark; anmu@dtu.dk 5 Museu Nacional de Arqueologia, Lisboa, Portugal; 3raposos@sapo.pt 6 Department of Geosciences, University of Padova, Italy; paolo.mozzi@unipd.it 7 School of Geography, Earth and Environmental Sciences, University of Plymouth, UK; m.stokes@plymouth.ac.uk The stratigraphic units that record the evolution of the Tagus River in Portugal (study area between Vila Velha de Ródão and Porto Alto villages; Fig. 1) have different sedimentary characteristics and lithic industries (Cunha et al., 2012): - a culminant sedimentary unit (the ancestral Tagus, before the drainage network entrenchment) – SLD13 (+142 to 262 m above river bed – a.r.b.; with probable age ca. 3,6 to 1,8 Ma), without artefacts; - T1 terrace (+84 to 180 m; ca. 1000? to 900 ka), without artefacts; - T2 terrace (+57 to 150 m; top deposits with a probable age ca. 600 ka), without artefacts; - T3 terrace (+43 to 113 m; ca. 460 to 360? ka), without artefacts; - T4 terrace (+26 to 55 m; ca. 335 a 155 ka), Lower Paleolithic (Acheulian) at basal and middle levels but early Middle Paleolithic at top levels; - T5 terrace (+5 to 34 m; 135 to 73 ka), Middle Paleolithic (Mousterian; Levallois technique); - T6 terrace (+3 to 14 m; 62 to 32 ka), late Middle Paleolithic (late Mousterian); - Carregueira Sands (aeolian sands) and colluvium (+3 a ca. 100 m; 32 to 12 ka), Upper Paleolithic to Epipaleolithic; - alluvial plain (+0 to 8 m; ca. 12 ka to present), Mesolithic and more recent industries. The differences in elevation (a.r.b.) of the several terrace staircases results from differential uplift due to active faults. Longitudinal correlation with the terrace levels indicates that a graded profile ca. 200 km long was achieved during terrace formation periods and a strong control by sea base level was determinant for terrace formation. The Neogene sedimentary units constituted the main source of sediments for the fluvial terraces (Fig. 2). Geomorphological mapping, coupled with lithostratigraphy, sedimentology and luminescence dating (quartz-OSL and K-feldspar post-IRIR290) were used in this study focused on the T4 terrace, which comprises a Lower Gravels (LG) unit and an Upper Sand (US) unit. The thick, coarse and dominantly massive gravels of the LG unit indicate deposition by a coarse bed-load braided river, with strong sediment supply, high gradient and fluvial competence, during conditions of rapidly rising sea level. Luminescence dating only provided minimum ages but it is probable that the LG unit corresponds to the earlier part of the MIS9 (ca. 335 to 325 ka), immediately postdating the incision promoted by the very low sea level (reaching ca. -140 m) during MIS10 (362 to 337 ka), a period of relatively cold climate conditions with weak vegetation cover on slopes and low sea level. Fig. 1. Main Portuguese reaches in which the Tagus River can be divided (Lower Tagus Basin): I – from the Spanish border to Arneiro (a general E–W trend, mainly consisting of polygonal segments); II – from Arneiro to Gavião (NE–SW); III – from Gavião to Arripiado (E–W); IV – from Arripiado to Vila Franca de Xira (NNE-SSW); V – from Vila Franca de Xira to the Atlantic shoreline. The faults considered to be the limit of the referred fluvial sectors are: F1 – Ponsul-Arneiro fault (WSW-ENE); F2 – Gavião fault (NW-SE); F3 – Ortiga fault (NW-SE); F4 – Vila Nova da Barquinha fault (W-E); F5 – Arripiado-Chamusca fault (NNE-SSW). 1 – estuary; 2 – terraces; 3 – faults; 4 – Tagus main channel. The main Iberian drainage basins are also represented (inset). The lower and middle parts of the US unit, comprising an alternation of clayish silts with paleosols and minor sands to the east (flood-plain deposits) and sand deposits to the west (channel belt), have a probable age of ca. 325 to 200 ka. This points to formation during MIS9 to MIS7, under conditions of high to medium sea levels and warm to mild conditions. The upper part of the US unit, dominated by sand facies and with OSL ages of ca. 200 to 154 ka, correlates with the early part of the MIS6. During this period, progradation resulted from climate deterioration and relative depletion of vegetation that promoted enhanced sediment production in the catchment, coupled with initiation of sea-level lowering that increased the longitudinal slope. The Vale do Forno and Vale da Atela archaeological sites (Alpiarça, central Portugal) document the earliest human occupation in the Lower Tagus River, well established in geomorphological and environmental terms, within the Middle Pleistocene. The Lower Palaeolithic sites were found on the T4 terrace (+26 m, a.r.b.). The oldest artefacts previously found in the LG unit, display crude bifacial forms that can be attributed to the Acheulian, with a probable age of ca. 335 to 325 ka. The T4 US unit has archaeological sites stratigraphically documenting successive phases of an evolved Acheulian, that probably date ca. 325 to 300 ka. Notably, these Lower Palaeolithic artisans were able to produce tools with different sophistication levels, simply by applying different strategies: more elaborated reduction sequences in case of bifaces and simple reduction sequences to obtain cleavers. Fig. 2. . Simplified geologic map of the Lower Tagus Cenozoic basin, adapted from the Carta Geológica de Portugal, 1/500000, 1992). The study area (comprising the Vale do Forno and Vale de Atela sites) is located on the more upstream sector of the Lower Tagus River reach IV, between Arripiado and Chamusca villages. 1 – alluvium (Holocene); 2 – terraces (Pleistocene); 3 – sands, silts and gravels (Paleogene to Pliocene); 4 – Sintra Massif (Cretaceous); 5 – limestones, marls, silts and sandstones (Mesozoic); 6 – quartzites (Ordovician); 7 – basement (Proterozoic to Palaeozoic); 8 – main fault. The main Portuguese reaches of the Tagus River are identified (I to V). The VF3 site (Milharós), containing a Final Acheulian industry, with fine and elaborated bifaces) found in a stratigraphic level located between the T4 terrace deposits and a colluvium associated with Late Pleistocene aeolian sands (32 to 12 ka), has an age younger than ca. 154 ka but much older than 32 ka. In the study area, the sedimentary units of the T4 terrace seem to record the river response to sea-level changes and climatically-driven fluctuations in sediment supply. REFERENCES Cunha P. P., Almeida N. A. C., Aubry T., Martins A. A., Murray A. S., Buylaert J.-P., Sohbati R., Raposo L., Rocha L., 2012, Records of human occupation from Pleistocene river terrace and aeolian sediments in the Arneiro depression (Lower Tejo River, central eastern Portugal). Geomorphology, vol. 165-166, pp. 78-90.

A Network Model of Static Foam Drainage

On the basis of a more realistic tetrakaidecahedral structure of foam bubbles, a network model of static foam drainage has been developed. The model considers the foam to be made up of films and Plateau borders. The films drain into the adjacent Plateau borders, which in turn form a network through which the liquid moves from the foam to the liquid pool. From the structure, a unit flow cell was found, which constitutes the foam when stacked together both horizontally and vertically. Symmetry in the unit flow cell indicates that the flow analysis of a part of it can be employed to obtain the drainage for the whole foam. Material balance equations have been written for each segment of this subsection, ensuring connectivity, and solved with the appropriate boundary and initial conditions. The calculated rates of drainage, when compared with the available experimental results, indicate that the model predicts the experimental results well.