Influence of soil stiffness on the dynamic response of an offshore wind turbine

La industria de la energía eólica marina ha crecido de forma significativa durante los últimos 15 años, y se espera que siga creciendo durante los siguientes. La construcción de torres en aguas cada vez más profundas y el aumento en potencia y tamaño de las turbinas han creado la necesidad de diseñar estructuras de soporte cada vez más fiables y optimizadas, lo que requiere un profundo conocimiento de su comportamiento. Este trabajo se centra en la respuesta dinámica de una turbina marina con cimentación tipo monopilote y sobre la que actúa la fuerza del viento. Se han realizado cálculos con distintas propiedades del suelo para cubrir un rango de rigideces que va desde una arena muy suelta a una muy densa. De este modo se ha analizado la influencia que tiene la rigidez del suelo en el comportamiento de la estructura. Se han llevado a cabo análisis estáticos y dinámicos en un modelo de elementos finitos implementado en Abaqus. El desplazamiento en la cabeza de la torre y la tensión en su base se han obtenido en función de la rigidez del suelo, y con ellos se ha calculado la amplificación dinámica producida cuando la frecuencia natural del sistema suelo‐cimentación torre se aproxima a la frecuencia de la carga. Dos diferentes enfoques a la hora de modelizar el suelo se han comparado: uno utilizando elementos continuos y otro utilizando muelles elásticos no lineales. Por último, un análisis de fiabilidad se ha llevado a cabo con un modelo analítico para calcular la probabilidad de resonancia del sistema, en el que se han considerado las propiedades de rigidez del suelo como variables aleatorias. Offshore wind energy industry has experienced a significant growth over the past 15 years, and it is expected to continue its growth in the coming years. The expansion to increasingly deep waters and the rise in power and size of the turbines have led to a need for more reliable and optimized support designs, which requires an extensive knowledge of the behaviour of these structures. This work focuses on the dynamic response of an offshore wind turbine founded on a monopile and subjected to wind loading. Different soil properties have been considered in order to cover the range of stiffness from a very loose to a very dense sand. In this way, the influence of stiffness on the structure behaviour has been assessed. Static and dynamic analyses have been carried out by means of a finite element model implemented in Abaqus. Head displacement and stress at the tower base have been obtained as functions of soil stiffness, and they have been used to calculate the dynamic amplification that is produced when the natural frequency of the system soil‐foundation‐tower approaches the load frequency. Two different approaches of soil modelling have been compared: soil modelled as a continuum and soil simulated with non linear elastic springs. Finally, a reliability analysis to assess the probability of resonance has been performed with an analytical model, in which soil stiffness properties are considered as stochastic variables.

Real-Time Dynamic PGD Calculation of Non-Linear Soil Behavior

El estudio sísmico en los últimos 50 años y el análisis del comportamiento dinámico del suelo revelan que el comportamiento del suelo es altamente no lineal e histéretico incluso para pequeñas deformaciones. El comportamiento no lineal del suelo durante un evento sísmico tiene un papel predominante en el análisis de la respuesta de sitio. Los análisis unidimensionales de la respuesta sísmica del suelo son a menudo realizados utilizando procedimientos lineales equivalentes, que requieren generalmente pocos parámetros conocidos. Los análisis de respuesta de sitio no lineal tienen el potencial para simular con mayor precisión el comportamiento del suelo, pero su aplicación en la práctica se ha visto limitada debido a la selección de parámetros poco documentadas y poco claras, así como una inadecuada documentación de los beneficios del modelado no lineal en relación al modelado lineal equivalente. En el análisis del suelo, el comportamiento del suelo es aproximado como un sólido Kelvin-Voigt con un módulo de corte elástico y amortiguamiento viscoso. En el análisis lineal y no lineal del suelo se están considerando geometrías y modelos reológicos más complejos. El primero está siendo dirigido por considerar parametrizaciones más ricas del comportamiento linealizado y el segundo mediante el uso de multi-modo de los elementos de resorte-amortiguador con un eventual amortiguador fraccional. El uso del cálculo fraccional está motivado en gran parte por el hecho de que se requieren menos parámetros para lograr la aproximación exacta a los datos experimentales. Basándose en el modelo de Kelvin-Voigt, la viscoelasticidad es revisada desde su formulación más estándar a algunas descripciones más avanzada que implica la amortiguación dependiente de la frecuencia (o viscosidad), analizando los efectos de considerar derivados fraccionarios para representar esas contribuciones viscosas. Vamos a demostrar que tal elección se traduce en modelos más ricos que pueden adaptarse a diferentes limitaciones relacionadas con la potencia disipada, amplitud de la respuesta y el ángulo de fase. Por otra parte, el uso de derivados fraccionarios permite acomodar en paralelo, dentro de un análogo de Kelvin-Voigt generalizado, muchos amortiguadores que contribuyen a aumentar la flexibilidad del modelado para la descripción de los resultados experimentales. Obviamente estos modelos ricos implican muchos parámetros, los asociados con el comportamiento y los relacionados con los derivados fraccionarios. El análisis paramétrico de estos modelos requiere técnicas numéricas eficientemente capaces de simular comportamientos complejos. El método de la Descomposición Propia Generalizada (PGD) es el candidato perfecto para la construcción de este tipo de soluciones paramétricas. Podemos calcular off-line la solución paramétrica para el depósito de suelo, para todos los parámetros del modelo, tan pronto como tales soluciones paramétricas están disponibles, el problema puede ser resuelto en tiempo real, porque no se necesita ningún nuevo cálculo, el solucionador sólo necesita particularizar on-line la solución paramétrica calculada off-line, que aliviará significativamente el procedimiento de solución. En el marco de la PGD, parámetros de los materiales y los diferentes poderes de derivación podrían introducirse como extra-coordenadas en el procedimiento de solución. El cálculo fraccional y el nuevo método de reducción modelo llamado Descomposición Propia Generalizada han sido aplicado en esta tesis tanto al análisis lineal como al análisis no lineal de la respuesta del suelo utilizando un método lineal equivalente. ABSTRACT Studies of earthquakes over the last 50 years and the examination of dynamic soil behavior reveal that soil behavior is highly nonlinear and hysteretic even at small strains. Nonlinear behavior of soils during a seismic event has a predominant role in current site response analysis. One-dimensional seismic ground response analysis are often performed using equivalent-linear procedures, which require few, generally well-known parameters. Nonlinear analyses have the potential to more accurately simulate soil behavior, but their implementation in practice has been limited because of poorly documented and unclear parameter selection, as well as inadequate documentation of the benefits of nonlinear modeling relative to equivalent linear modeling. In soil analysis, soil behaviour is approximated as a Kelvin-Voigt solid with a elastic shear modulus and viscous damping. In linear and nonlinear analysis more complex geometries and more complex rheological models are being considered. The first is being addressed by considering richer parametrizations of the linearized behavior and the second by using multi-mode spring-dashpot elements with eventual fractional damping. The use of fractional calculus is motivated in large part by the fact that fewer parameters are required to achieve accurate approximation of experimental data. Based in Kelvin-Voigt model the viscoelastodynamics is revisited from its most standard formulation to some more advanced description involving frequency-dependent damping (or viscosity), analyzing the effects of considering fractional derivatives for representing such viscous contributions. We will prove that such a choice results in richer models that can accommodate different constraints related to the dissipated power, response amplitude and phase angle. Moreover, the use of fractional derivatives allows to accommodate in parallel, within a generalized Kelvin-Voigt analog, many dashpots that contribute to increase the modeling flexibility for describing experimental findings. Obviously these rich models involve many parameters, the ones associated with the behavior and the ones related to the fractional derivatives. The parametric analysis of all these models require efficient numerical techniques able to simulate complex behaviors. The Proper Generalized Decomposition (PGD) is the perfect candidate for producing such kind of parametric solutions. We can compute off-line the parametric solution for the soil deposit, for all parameter of the model, as soon as such parametric solutions are available, the problem can be solved in real time because no new calculation is needed, the solver only needs particularize on-line the parametric solution calculated off-line, which will alleviate significantly the solution procedure. Within the PGD framework material parameters and the different derivation powers could be introduced as extra-coordinates in the solution procedure. Fractional calculus and the new model reduction method called Proper Generalized Decomposition has been applied in this thesis to the linear analysis and nonlinear soil response analysis using a equivalent linear method.

Water Deficit and Spatial Pattern of Leaf Development. Variability in Responses Can Be Simulated Using a Simple Model of Leaf Development1

We analyzed the effect of short-term water deficits at different periods of sunflower (Helianthus annuus L.) leaf development on the spatial and temporal patterns of tissue expansion and epidermal cell division. Six water-deficit periods were imposed with similar and constant values of soil water content, predawn leaf water potential and [ABA] in the xylem sap, and with negligible reduction of the rate of photosynthesis. Water deficit did not affect the duration of expansion and division. Regardless of their timing, deficits reduced relative expansion rate by 36% and relative cell division rate by 39% (cells blocked at the G0-G1 phase) in all positions within the leaf. However, reductions in final leaf area and cell number in a given zone of the leaf largely differed with the timing of deficit, with a maximum effect for earliest deficits. Individual cell area was only affected during the periods when division slowed down. These behaviors could be simulated in all leaf zones and for all timings by assuming that water deficit affects relative cell division rate and relative expansion rate independently, and that leaf development in each zone follows a stable three-phase pattern in which duration of each phase is stable if expressed in thermal time (C. Granier and F. Tardieu [1998b] Plant Cell Environ 21: 695–703).

Biological impact on mineral dissolution: Application of the lichen model to understanding mineral weathering in the rhizosphere

Microorganisms modify rates and mechanisms of chemical and physical weathering and clay growth, thus playing fundamental roles in soil and sediment formation. Because processes in soils are inherently complex and difficult to study, we employ a model based on the lichen–mineral system to identify the fundamental interactions. Fixed carbon released by the photosynthetic symbiont stimulates growth of fungi and other microorganisms. These microorganisms directly or indirectly induce mineral disaggregation, hydration, dissolution, and secondary mineral formation. Model polysaccharides were used to investigate direct mediation of mineral surface reactions by extracellular polymers. Polysaccharides can suppress or enhance rates of chemical weathering by up to three orders of magnitude, depending on the pH, mineral surface structure and composition, and organic functional groups. Mg, Mn, Fe, Al, and Si are redistributed into clays that strongly adsorb ions. Microbes contribute to dissolution of insoluble secondary phosphates, possibly via release of organic acids. These reactions significantly impact soil fertility. Below fungi–mineral interfaces, mineral surfaces are exposed to dissolved metabolic byproducts. Through this indirect process, microorganisms can accelerate mineral dissolution, leading to enhanced porosity and permeability and colonization by microbial communities.

USLE and WEPP Soil Erosion Models: Their Applicability and Use in Generating Cover Values to Stabilize Soil in the Los Alamos, New Mexico Area

Soil erosion is a naturally occurring process that involves the detachment, transport, and deposition of soil particles. Disturbances such as thinning and wildfire can reduce cover greatly and increase erosion rates. Forest managers may use erosion prediction tools, such as the Universal Soil Loss Equation (USLE) and Water Erosion Prediction Project (WEPP) to estimate erosion rates and develop techniques to manage erosion. However, it is important to understand the differences and the applications of each model. Erosion rates were generated by each model and the model most applicable to the study site, Los Alamos, New Mexico was determined. It was also used to find the amount of cover needed to stabilize soil. The USLE is a simpler model and less complicated than a computer model like WEPP, and thus easier to manipulate to estimate cover values. Predicted cover values were compared to field cover values. Cover is necessary to establish effective erosion control guidelines.

Feedbacks between vegetation pattern and resource loss dramatically decrease ecosystem resilience and restoration potential in a simple dryland model

Conceptual frameworks of dryland degradation commonly include ecohydrological feedbacks between landscape spatial organization and resource loss, so that decreasing cover and size of vegetation patches result in higher water and soil losses, which lead to further vegetation loss. However, the impacts of these feedbacks on dryland dynamics in response to external stress have barely been tested. Using a spatially-explicit model, we represented feedbacks between vegetation pattern and landscape resource loss by establishing a negative dependence of plant establishment on the connectivity of runoff-source areas (e.g., bare soils). We assessed the impact of various feedback strengths on the response of dryland ecosystems to changing external conditions. In general, for a given external pressure, these connectivity-mediated feedbacks decrease vegetation cover at equilibrium, which indicates a decrease in ecosystem resistance. Along a gradient of gradual increase of environmental pressure (e.g., aridity), the connectivity-mediated feedbacks decrease the amount of pressure required to cause a critical shift to a degraded state (ecosystem resilience). If environmental conditions improve, these feedbacks increase the pressure release needed to achieve the ecosystem recovery (restoration potential). The impact of these feedbacks on dryland response to external stress is markedly non-linear, which relies on the non-linear negative relationship between bare-soil connectivity and vegetation cover. Modelling studies on dryland vegetation dynamics not accounting for the connectivity-mediated feedbacks studied here may overestimate the resistance, resilience and restoration potential of drylands in response to environmental and human pressures. Our results also suggest that changes in vegetation pattern and associated hydrological connectivity may be more informative early-warning indicators of dryland degradation than changes in vegetation cover.

Comparison among monitoring strategies to assess water flow dynamic and soil hydraulic properties in agricultural soils

Irrigated agriculture is usually performed in semi-arid regions despite scarcity of water resources. Therefore, optimal irrigation management by monitoring the soil is essential, and assessing soil hydraulic properties and water flow dynamics is presented as a first measure. For this purpose, the control of volumetric water content, θ, and pressure head, h, is required. This study adopted two types of monitoring strategies in the same experimental plot to control θ and h in the vadose zone: i) non-automatic and more time-consuming; ii) automatic connected to a datalogger. Water flux was modelled with Hydrus-1D using the data collected from both acquisition strategies independently (3820 daily values for the automatic; less than 1000 for the non-automatic). Goodness-of-fit results reported a better adjustment in case of automatic sensors. Both model outputs adequately predicted the general trend of θ and h, but with slight differences in computed annual drainage (711 mm and 774 mm). Soil hydraulic properties were inversely estimated from both data acquisition systems. Major differences were obtained in the saturated volumetric water content, θs, and the n and α van Genuchten model shape parameters. Saturated hydraulic conductivity, Ks, shown lower variability with a coefficient of variation range from 0.13 to 0.24 for the soil layers defined. Soil hydraulic properties were better assessed through automatic data acquisition as data variability was lower and accuracy was higher.

Fate and Transport of Naproxen in a Sandy Aquifer Material: Saturated Column Studies and Model Evaluation

Naproxen-C14H14O3 is a nonsteroidal anti-inflammatory drug which has been found at detectable concentrations in wastewater, surface water, and groundwater. Naproxen is relatively hydrophilic and is in anionic form at pH between 6 and 8. In this study, column experiments were performed using an unconsolidated aquifer material from an area near Barcelona (Spain) to assess transport and reaction mechanisms of Naproxen in the aquifer matrix under different pore water fluxes. Results were evaluated using HYDRUS-1D, which was used to estimate transport parameters. Batch sorption isotherms for Naproxen conformed with the linear model with a sorption coefficient of 0.42 (cm3 g−1), suggesting a low sorption affinity. Naproxen breakthrough curves (BTCs) measured in soil columns under steady-state, saturated water flow conditions displayed similar behavior, with no apparent hysteresis in sorption or dependence of retardation (R, 3.85-4.24) on pore water velocities. Soil sorption did not show any significant decrease for increasing flow rates, as observed from Naproxen recovery in the effluent. Sorption parameters estimated by the model suggest that Naproxen has a low sorption affinity to aquifer matrix. Most sorption of Naproxen occurred on the instantaneous sorption sites, with the kinetic sorption sites representing only about 10 to 40% of total sorption.

Routine histological technique applied to a biological model for environmental studies

Introduction and aim Earthworms are an important test organisms used in several standard ecotoxicological tests (e.g. ISO, 2008, 2012; OECD, 2004, 1094) and they are also model organisms to test soil structure and permeability, as well as for important soil services as the degradation organic matter. Although histopathological changes have been demonstrated to be good biomarkers to assess the exposure of these animals to different physical and chemical stress agents, studies with clear and high quality images describing normal tissue conditions are scarce in the literature, so the aim of this work was to better characterize this biological model. Material and Methods Eight adult earthworms exposed to an artificial standard soil (OECD, 2004) for 28 days, were extracted and placed in a plastic box to depurate their gut content, fixed in 10% neutral-buffered formalin and processed for routine histopathological diagnosis. Results and discussion Satisfactory histological sections were obtained. Some difficulties were faced related with microtome sectioning, resulting in artefacts, namely lines across sections produced by a nick in the cutting edge of the microtome knife cutting tear, motivated by the presence of sand and other solid particles that persisted in the gut of earthworms. Nevertheless, it was possible to obtain representative figures from different earthworm sections. Conclusion Routine histological technique was effective for obtaining satisfactory histological sections and the knowledge of the histology of earthworms could be very useful for future application in environmental studies, using this biological model.

Characterization of histopathological lesions in a biological model for environmental studies: Eisenia andrei Bouché. Preliminar results

Introduction The strong expansion of the world production of plastics caused a severe accumulation of plastic debris in the environment, which makes them one of the most important contaminants, growing as a global environmental problem. Although the production in Europe has been relatively constant in the last 10 years, world plastic production continues to increase, affecting soil biota and their functions. Objectives Thus, in order to evaluate the effects of MP in soil-dwelling organisms, earthworms (Eisenia andrei Bouché), were exposed to standard artificial soil mixed with MPs and the authors documented, using microscopic figures, the pathological lesions found in this biological model. Material and Methods Eight adult earthworms extracted from soils contaminated with different concentrations of MP (mg/kgdw) with sizes ranging between 250-1000 m, were fixed in 10% neutral-buffered formalin and processed for routine histopathological diagnosis. Results and discussion Contrary to what would be expected, MP were not found throughout the GI tube of earthworms but several lesions were found in the individuals extracted from the soils with high MP concentrations, when compared with control group, namely epithelial intestinal atrophy and evidences of inflammatory responses to this stress agent. Conclusion Earthworms have probably avoided the consumption of the biggest MPs. However, evidences point for lesions that were likely caused by the smallest MPs that were likely egested during the depuration phase.

A daily rainfall erosivity model for Western Amazonia

Rainfall erosivities as defined by the R factor from the universal soil loss equation were determined for all events during a two-year period at the station La Cuenca in western Amazonia. Three methods based on a power relationship between rainfall amount and erosivity were then applied to estimate event and daily rainfall erosivities from the respective rainfall amounts. A test of the resulting regression equations against an independent data set proved all three methods equally adequate in predicting rainfall erosivity from daily rainfall amount. We recommend the Richardson model for testing in the Amazon Basin, and its use with the coefficient from La Cuenca in western Amazonia.

Global assessment of soil phosphorus retention potential

Limited availability of P in soils to crops may be due to deficiency and/or severe P retention. Earlier studies that drew on large soil profile databases have indicated that it is not (yet) feasible to present meaningful values for "plant-available" soil P, obtained according to comparable analytical methods, that may be linked to soil geographical databases derived from 1:5 million scale FAO Digital Soil Map of the World, such as the 5 x 5 arc-minute version of the ISRIC-WISE database. Therefore, an alternative solution for studying possible crop responses to fertilizer-P applied to soils, at a broad scale, was sought. The approach described in this report considers the inherent capacity of soils to retain phosphorus (P retention), in various forms. Main controlling factors of P retention processes, at the broad scale under consideration, are considered to be pH, soil mineralogy, and clay content. First, derived values for these properties were used to rate the inferred capacity for P retention of the component soil units of each map unit (or grid cell) using four classes (i.e., Low, Moderate, High, and Very High). Subsequently, the overall soil phosphorus retention potential was assessed for each mapping unit, taking into account the P-ratings and relative proportion of each component soil unit. Each P retention class has been assigned to a likely fertilizer P recovery fraction, derived from the literature, thereby permitting spatially more detailed, integrated model-based studies of environmental sustainability and agricultural production at the global and continental level (< 1:5 million). Nonetheless, uncertainties remain high; the present analysis provides an approximation of world soil phosphorus retention potential.

Catchment characteristics and soil and water chemistry of the San Francisco basin and subbasins, Ecuador

In this study, the Mean Transit Time and Mixing Model Analysis methods are combined to unravel the runoff generation process of the San Francisco River basin (73.5 km**2) situated on the Amazonian side of the Cordillera Real in the southernmost Andes of Ecuador. The montane basin is covered with cloud forest, sub-páramo, pasture and ferns. Nested sampling was applied for the collection of streamwater samples and discharge measurements in the main tributaries and outlet of the basin, and for the collection of soil and rock water samples. Weekly to biweekly water grab samples were taken at all stations in the period April 2007-November 2008. Hydrometric data, Mean Transit Time and Mixing Model Analysis allowed preliminary evaluation of the processes controlling the runoff in the San Francisco River basin. Results suggest that flow during dry conditions mainly consists of lateral flow through the C-horizon and cracks in the top weathered bedrock layer, and that all subcatchments have an important contribution of this deep water to runoff, no matter whether pristine or deforested. During normal to low precipitation intensities, when antecedent soil moisture conditions favour water infiltration, vertical flow paths to deeper soil horizons with subsequent lateral subsurface flow contribute most to streamflow. Under wet conditions in forested catchments, streamflow is controlled by near surface lateral flow through the organic horizon. Exceptionally, saturation excess overland flow occurs. By absence of the litter layer in pasture, streamflow under wet conditions originates from the A horizon, and overland flow.

Characteristics of soil profiles from the Ouémé Upper Catchment (HVO), Térou, Benin

Soil degradation threatens agricultural production and food security in Sub-Saharan Africa. In the coming decades, soil degradation, in particular soil erosion, will become worse through the expansion of agriculture into savannah and forest and changes in climate. This study aims to improve the understanding of how land use and climate change affect the hydrological cycle and soil erosion rates at the catchment scale. We used the semi-distributed, time-continuous erosion model SWAT (Soil Water Assessment Tool) to quantify runoff processes and sheet and rill erosion in the Upper Ouémé River catchment (14500 km**2, Central Benin) for the period 1998-2005. We could then evaluate a range of land use and climate change scenarios with the SWAT model for the period 2001-2050 using spatial data from the land use model CLUE-S and the regional climate model REMO. Field investigations were performed to parameterise a soil map, to measure suspended sediment concentrations for model calibration and validation and to characterise erosion forms, degraded agricultural fields and soil conservation practices. Modelling results reveal current "hotspots" of soil erosion in the north-western, eastern and north-eastern parts of the Upper Ouémé catchment. As a consequence of rapid expansion of agricultural areas triggered by high population growth (partially caused by migration) and resulting increases in surface runoff and topsoil erosion, the mean sediment yield in the Upper Ouémé River outlet is expected to increase by 42 to 95% by 2025, depending on the land use scenario. In contrast, changes in climate variables led to decreases in sediment yield of 5 to 14% in 2001-2025 and 17 to 24% in 2026-2050. Combined scenarios showed the dominance of land use change leading to changes in mean sediment yield of -2 to +31% in 2001-2025. Scenario results vary considerably within the catchment. Current "hotspots" of soil erosion will aggravate, and a new "hotspot" will appear in the southern part of the catchment. Although only small parts of the Upper Ouémé catchment belong to the most degraded zones in the country, sustainable soil and plant management practices should be promoted in the entire catchment. The results of this study can support planning of soil conservation activities in Benin.

A numerical model for predicting shoreline changes /

A mathematical model for long-term, three-dimensional shoreline evolution is developed. The combined effects of variations of sea level; wave refraction and diffraction; loss of sand by density currents during storms, by rip currents, and by wind; bluff erosion and berm accretion; effects of manmade structures such as long groin or navigational structures; and beach nourishment are all taken into account. A computer program is developed with various subroutines which permit modification as the state-of-the-art progresses. The program is applied to a test case at Holland Harbor, Michigan. (Author).