Adjoint free four-dimensional variational data assimilation for a storm surge model of the German North Sea

In this paper, a data assimilation scheme based on the adjoint free Four-Dimensional Variational(4DVar) method is applied to an existing storm surge model of the German North Sea. To avoid the need of an adjoint model, an ensemble-like method to explicitly represent the linear tangent equation is adopted. Results of twin experiments have shown that the method is able to recover the contaminated low dimension model parameters to their true values. The data assimilation scheme was applied to a severe storm surge event which occurred in the North Sea in December 5, 2013. By adjusting wind drag coefficient, the predictive ability of the model increased significantly. Preliminary experiments have shown that an increase in the predictive ability is attained by narrowing the data assimilation time window.     

Runoff and sediment yield from land under various uses in a Mediterranean mountain area: long‐term results from an experimental station

In this study we analyzed runoff and sediment yield from land under various traditional and current land uses in Mediterranean mountain areas, using long‐term data from an experimental station in the Aísa Valley, Central Spanish Pyrenees. Monitoring at this station has provided 20 years of data that can help explain the hydrological and geomorphological changes that have been observed at larger spatial scales, and also the changes that have occurred to some of the most characteristic landscapes of the Mediterranean middle mountains. In spite of the problems associated with the use of small experimental plots, the results obtained are consistent with other studies in the Mediterranean region, and confirm the strong influence of land use changes on runoff generation and sediment yield. The results indicate that: (i) cereal cultivation on steep slopes (both alternating cereal cultivation and fallow on sloping fields and shifting agriculture on the steepest slopes) represents a major problem for soil conservation. This explains the occurrence throughout the Mediterranean mountains of many degraded hillslopes, which show evidence of sheet wash erosion, rilling, gullying and shallow landsliding; (ii) farmland abandonment has led to a marked reduction in runoff and sediment yield as a consequence of rapid plant recolonization, particularly by dense shrubs; (iii) the natural transformation of abandoned fields into grazing meadows has reduced runoff and sediment yield. Land use trends in the Mediterranean mountains are mainly characterized by generalized farmland abandonment and a decrease in livestock pressure. From a hydrological and geomorphological point of view the main consequences have been a reduction in overland flow from the hillslopes, and a reduction in sediment sources, with differences up to one order of magnitude in sediment yield from dense shrub cover and grazing meadow areas compared with areas under shifting agriculture. Copyright © 2012 John Wiley & Sons, Ltd.     

Climatology of extratropical cyclones over the South American–southern oceans sector

A climatology of extratropical cyclones is presented. Extratropical cyclones, their main characteristics and their predominant tracks, as well as their interannual variability, affect weather in South America. For that purpose, a storm track database has been compiled by applying a cyclone tracking scheme to six-hourly sea level pressure fields, available from the National Center for Environmental Prediction–National Center for Atmospheric Research reanalyses II for the 1979–2003 period. The spatial distribution of the cyclogenesis frequency shows two main centers: one around Northern Argentina, Uruguay, and Southern Brazil in all seasons and the other near to the North Antarctic Peninsula. The lifetime of extratropical cyclones in the South American sector exhibits small seasonality, being typically of the order of 3.0 days during most of the year and slightly higher (3.5 days) in austral summer. The distance travelled by the cyclones formed in the South American sector tends to be smaller than the total paths found in other areas of the Southern Hemisphere. A k-mean clustering technique is used to summarize the analysis of the 25-year climatology of cyclone tracks. Three clusters were found: one storm-track cluster in Northeast Argentina; a second one west of the Andes Cordillera; and a third cluster located to the north of the Antarctic Peninsula (around the Weddell Sea). The influence of the Antarctic Oscillation (AAO) in the variability of extratropical cyclones is explored, and some signals of the impacts of the variability of the AAO can be observed in the position of the extratropical cyclones around 40°S, while the impacts on the intensity is detected around 55°S.     

Sorption-desorption kinetics and toxic cell concentration in marine phytoplankton microalgae exposed to Linear Alkylbenzene Sulfonate

Linear Alkylbenzene Sulfonates (LAS) are ubiquitous surfactants. Traces can be found in coastal environments. Sorption and toxicity of C₁₂-LAS congeners were studied in controlled conditions (2-3500 μg C₁₂LAS/L) in five marine phytoplanktonic species, using standardized methods. IC₅₀ values ranged from 0.5 to 2 mg LAS/L. Sorption of ¹⁴C₁₂-6 LAS isomer was measured at environmentally relevant trace levels (4 μg/L) using liquid scintillation counting. Steady-state sorption on algae was reached within 5 h in the order dinoflagellate > diatoms > green algae. The sorption data, fitted a L-type Freundlich isotherm, indicating saturation. Desorption was rapid but a low LAS fraction was still sorbed after 24 h. Toxic cell concentration was 0.38 ± 0.09 mg/g for the studied species. LAS toxicity results from sorption on biological membranes leading to non-specific disturbance of algal growth. Results indicate that LAS concentrations in coastal environments do not represent a risk for these organisms.     

An open source GIS‐based Planning Support System: Application to the land use plan of La Troncal,Ecuador

Planning Support Systems (PSS) comprise a wide variety of geo‐technological tools related to GIS and spatial modeling aimed at addressing land planning processes. This article describes the OpenRules system, a PSS based on a previous system called RULES. Among OpenRules new features are its architecture, based exclusively on free and open source software, and its applicability to all land use types, including rural and urban uses. In addition, OpenRules incorporates an unlimited number of land evaluation factors and a new objective in land use spatial allocation. OpenRules has been programmed in Java and implemented as a module of the free GIS software gvSIG, with full integration between the GIS and the decision support tools. Decision support tools include multicriteria evaluation, multiobjective linear programming and heuristic techniques, which support three basic stages of land use planning processes, namely land suitability evaluation, land use area optimization and land use spatial allocation. The application of OpenRules to the region of La Troncal, Ecuador, demonstrates its capability to generate alternative and coherent solutions through a scientific and justified procedure at low cost in terms of time and resources.     

Comparison of P‐wave attenuation models of wave‐induced flow

Wave‐induced oscillatory fluid flow in the vicinity of inclusions embedded in porous rocks is one of the main causes for P‐wave dispersion and attenuation at seismic frequencies. Hence, the P‐wave velocity depends on wave frequency, porosity, saturation, and other rock parameters. Several analytical models quantify this wave‐induced flow attenuation and result in characteristic velocity–saturation relations. Here, we compare some of these models by analyzing their low‐ and high‐frequency asymptotic behaviours and by applying them to measured velocity–saturation relations. Specifically, the Biot–Rayleigh model considering spherical inclusions embedded in an isotropic rock matrix is compared with White's and Johnson's models of patchy saturation. The modeling of laboratory data for tight sandstone and limestone indicates that, by selecting appropriate inclusion size, the Biot‐Rayleigh predictions are close to the measured values, particularly for intermediate and high water saturations.     

Groundwater behaviour in Madeira, volcanic island (Portugal)

Madeira Island is a hot-spot originating from a mantle plume. K-Ar age determinations indicate that the emerged part of the island was generated during Post-Miocene times 6000–7000 years B.P. Groundwater occurs in perched-water bodies, spring discharge from them is high, about 3,650 l/s; in dike-impounded water and basal groundwater. Basal groundwater is exploited by tunnels (1,100 l/s) and wells (1,100 l/s). Hydraulic gradients range from 10^–4 to 10^–2 and transmissivity ranges from 1.16×10^–2 to 2.89×10^–1 m²/s, indicating the heterogeneity of the volcanic aquifers. Water mineralisation is variable, and electrical conductivity ranges from 50 to 3,300 S/cm. There is a difference between groundwater discharging from perched-water bodies (43–201 S/cm) and from basal groundwater (109–3,300 S/cm). Groundwater average pH is 7.37, with waters acid to slightly alkaline (6.13–8.4), and generally cold. Nevertheless, a few samples associated with fault zones can be classified as thermal waters.

---

Resumen La Isla de Madeira tuvo origen en un punto caliente proveniente de una pluma del manto. Las dataciones K – Ar, indican que la parte emergida de la isla se formó durante en Post—Mioceno, hace unos 6000–7000 años. El agua subterránea se encuentra en acuíferos colgados, cuya descarga a través de manantiales es alta, alrededor de 3650 l/s; también como agua retenida en diques o bien como agua subterránea propiamente dicha. Esta última es explotada por túneles (1100 l/s) y pozos (1100 l/s). Los gradientes hidráulicos varían desde 10–4 hasta 10–2 y la transmisividad varía desde 1.16×10^–2 hasta 2.89×10^–1 m²/s, indicando la heterogeneidad de los acuíferos volcánicos. La mineralización de agua es variable y la conductividad eléctrica varía desde 50 mS/cm hasta 3300 mS/cm. Hay una diferencia entre el agua subterránea proveniente de acuíferos colgados (43–201 mS/cm) y aquella del acuífero principal (109–3300 mS/cm). El pH promedio para el agua subterránea es 7.37; con variación desde ácida hasta ligeramente alcalina (6.13–8.4), y es generalmente fría. Sin embargo unos pocos ejemplos asociados con zonas de falla, pueden ser clasificados como agua termal.

---

Résumé Lîle de Madeira est un point chaud, provenant dun panache mantélique. Les datations par la méthode de K-Ar indiquent que la part émergée de lîle a été générée pendant le Post Miocène, il y a 6000–7000 ans. Les eaux souterraines sont cantonnées dans des aquifères perchées qui se déchargent par des sources à grands débits, avec des valeurs dapproximativement de 3500 l/s, ainsi que dans les dikes. Une importante quantité des eaux souterraine se trouve sur la forme deau de fonds, formant une nappe basale, qui est exploitée par des tunnels (1100 l/s) et par des puits (1100 l/s). Les gradients hydrauliques ont des valeurs dans lintervalle de 10^–2 à 10^–4, tandis que les transmissivités se rangent entre 1.16×10^–2 m²/s et 2.81×10^–1 m²/s, ce quindique la hétérogénéité de laquifère volcanique. La minéralisation des eaux est aussi variable avec des conductivités électriques qui se rangent ente 50 mS/cm et 3300 mS/cm. Il y a une différence entre la conductivité des eaux provenant des aquifères perchées (30–201 mS/cm) par rapport à la conductivité les eaux de la nappe basale (109–3300 mS/cm). Le pH moyen est de 7.37 avec des valeurs de 6.13 pour les eaux acides et de 8.4 pour les eaux à faible alcalinité. En général il sagit des eaux froides, mais quelques échantillons provenant des zones faillées peuvent être classifiés comme des eaux thermales.

     

Reconstructing fluvial bar surfaces from compound cross‐strata and the interpretation of bar accretion direction in large river deposits

The interpretation of fluvial styles from the rock record is based for a significant part on the identification of different types of fluvial bars, characterized by the geometric relationship between structures indicative of palaeocurrent and surfaces interpreted as indicative of bar form and bar accretion direction. These surfaces of bar accretion are the boundaries of flood‐related bar increment elements, which are typically less abundant in outcrops than what would be desirable, particularly in large river deposits in which each flood mobilizes large volumes of sediment, causing flood‐increment boundary surfaces to be widely spaced. Cross‐strata set boundaries, on the other hand, are abundant and indirectly reflect the process of unit bar accretion, inclined due to the combined effect of the unit bar surface inclination and the individual bedform climbing angle, in turn controlled by changes in flow structure caused by local bar‐scale morphology. This work presents a new method to deduce the geometry of unit bar surfaces from measured pairs of cross‐strata and cross‐strata set boundaries. The method can be used in the absence of abundant flood‐increment bounding surfaces; the study of real cases shows that, for both downstream and laterally accreting bars, the reconstructed planes are very similar to measured bar increment surfaces.     

LES of turbulent surface shear stress and pressure-gradient-driven flow on shallow continental shelves

Turbulent shear flows on shallow continental shelves (here shallow means that the interaction with the solid, no-slip bottom is important) are of great importance because tide- and wind-driven flows on the shelf are drivers of the transfer of momentum, heat, and mass (gas) across the air–sea interface. These turbulent flows play an important role because vertical mixing and current are vectors for the transport of sediment and bioactive material on continental shelves. Understanding the dynamics of this class of flows presents complications because of the presence of a free surface and also because the flow can be driven by a pressure gradient (a tidal current), a stress at the free surface (a wind-driven current), or a combination of both. In addition, the flow can be modified by the presence of a wave field that can induce Langmuir circulation (Langmuir, Science 87:119–123, 1938). Large eddy simulation is used to quantify the effects of pressure gradient and wind shear on the distinctive structures of the turbulent flow. From these computations, an understanding of the physics governing the turbulence of pressure-driven and wind-driven flows, how they can interact in a normal or a tangential direction, and the effect of wave forcing on these flows is obtained.