Impact of hydrographic data assimilation on the modelled Atlantic meridional overturning circulation

Here we make an initial step toward the development of an ocean assimilation system that can constrain the modelled Atlantic Meridional Overturning Circulation (AMOC) to support climate predictions. A detailed comparison is presented of 1° and 1/4° resolution global model simulations with and without sequential data assimilation, to the observations and transport estimates from the RAPID mooring array across 26.5° N in the Atlantic. Comparisons of modelled water properties with the observations from the merged RAPID boundary arrays demonstrate the ability of in situ data assimilation to accurately constrain the east-west density gradient between these mooring arrays. However, the presence of an unconstrained "western boundary wedge" between Abaco Island and the RAPID mooring site WB2 (16 km offshore) leads to the intensification of an erroneous southwards flow in this region when in situ data are assimilated. The result is an overly intense southward upper mid-ocean transport (0–1100 m) as compared to the estimates derived from the RAPID array. Correction of upper layer zonal density gradients is found to compensate mostly for a weak subtropical gyre circulation in the free model run (i.e. with no assimilation). Despite the important changes to the density structure and transports in the upper layer imposed by the assimilation, very little change is found in the amplitude and sub-seasonal variability of the AMOC. This shows that assimilation of upper layer density information projects mainly on the gyre circulation with little effect on the AMOC at 26° N due to the absence of corrections to density gradients below 2000 m (the maximum depth of Argo). The sensitivity to initial conditions was explored through two additional experiments using a climatological initial condition. These experiments showed that the weak bias in gyre intensity in the control simulation (without data assimilation) develops over a period of about 6 months, but does so independently from the overturning, with no change to the AMOC. However, differences in the properties and volume transport of North Atlantic Deep Water (NADW) persisted throughout the 3 year simulations resulting in a difference of 3 Sv in AMOC intensity. The persistence of these dense water anomalies and their influence on the AMOC is promising for the development of decadal forecasting capabilities. The results suggest that the deeper waters must be accurately reproduced in order to constrain the AMOC.

The soil-dwelling earthworm Allolobophora chlorotica modifies its burrowing behaviour in response to carbendazim applications

Carbendazim-amended soil was placed above or below unamended soil. Control tests comprised two layers of unamended soil. Allolobophora chlorotica earthworms were added to either the upper or the unamended soil. After 72 h vertical distributions of earthworms were compared between control and carbendazim-amended experiments. Earthworm distributions in the carbendazim-amended test containers differed significantly from the ‘normal’ distribution observed in the control tests. In the majority of the experiments, earthworms significantly altered their burrowing behaviour to avoid carbendazim. However, when earthworms were added to an upper layer of carbendazim-amended soil they remained in this layer. This non-avoidance is attributed to (1) the earthworms’ inability to sense the lower layer of unamended soil and (2) the toxic effect of carbendazim inhibiting burrowing. Earthworms modified their burrowing behaviour in response to carbendazim in the soil. This may explain anomalous results observed in pesticide field trials when carbendazim is used as a control substance.

Viscous coupling of shear-free turbulence across nearly flat fluid interfaces

The interactions between shear-free turbulence in two regions (denoted as + and − on either side of a nearly flat horizontal interface are shown here to be controlled by several mechanisms, which depend on the magnitudes of the ratios of the densities, ρ+/ρ−, and kinematic viscosities of the fluids, μ+/μ−, and the root mean square (r.m.s.) velocities of the turbulence, u0+/u0−, above and below the interface. This study focuses on gas–liquid interfaces so that ρ+/ρ− ≪ 1 and also on where turbulence is generated either above or below the interface so that u0+/u0− is either very large or very small. It is assumed that vertical buoyancy forces across the interface are much larger than internal forces so that the interface is nearly flat, and coupling between turbulence on either side of the interface is determined by viscous stresses. A formal linearized rapid-distortion analysis with viscous effects is developed by extending the previous study by Hunt & Graham (J. Fluid Mech., vol. 84, 1978, pp. 209–235) of shear-free turbulence near rigid plane boundaries. The physical processes accounted for in our model include both the blocking effect of the interface on normal components of the turbulence and the viscous coupling of the horizontal field across thin interfacial viscous boundary layers. The horizontal divergence in the perturbation velocity field in the viscous layer drives weak inviscid irrotational velocity fluctuations outside the viscous boundary layers in a mechanism analogous to Ekman pumping. The analysis shows the following. (i) The blocking effects are similar to those near rigid boundaries on each side of the interface, but through the action of the thin viscous layers above and below the interface, the horizontal and vertical velocity components differ from those near a rigid surface and are correlated or anti-correlated respectively. (ii) Because of the growth of the viscous layers on either side of the interface, the ratio uI/u0, where uI is the r.m.s. of the interfacial velocity fluctuations and u0 the r.m.s. of the homogeneous turbulence far from the interface, does not vary with time. If the turbulence is driven in the lower layer with ρ+/ρ− ≪ 1 and u0+/u0− ≪ 1, then uI/u0− ~ 1 when Re (=u0−L−/ν−) ≫ 1 and R = (ρ−/ρ+)(v−/v+)1/2 ≫ 1. If the turbulence is driven in the upper layer with ρ+/ρ− ≪ 1 and u0+/u0− ≫ 1, then uI/u0+ ~ 1/(1 + R). (iii) Nonlinear effects become significant over periods greater than Lagrangian time scales. When turbulence is generated in the lower layer, and the Reynolds number is high enough, motions in the upper viscous layer are turbulent. The horizontal vorticity tends to decrease, and the vertical vorticity of the eddies dominates their asymptotic structure. When turbulence is generated in the upper layer, and the Reynolds number is less than about 106–107, the fluctuations in the viscous layer do not become turbulent. Nonlinear processes at the interface increase the ratio uI/u0+ for sheared or shear-free turbulence in the gas above its linear value of uI/u0+ ~ 1/(1 + R) to (ρ+/ρ−)1/2 ~ 1/30 for air–water interfaces. This estimate agrees with the direct numerical simulation results from Lombardi, De Angelis & Bannerjee (Phys. Fluids, vol. 8, no. 6, 1996, pp. 1643–1665). Because the linear viscous–inertial coupling mechanism is still significant, the eddy motions on either side of the interface have a similar horizontal structure, although their vertical structure differs.

The link between the Barents Sea and ENSO events simulated by NEMO model

An analysis of observational data in the Barents Sea along a meridian at 33°30' E between 70°30' and 72°30' N has reported a negative correlation between El Niño/La Niña Southern Oscillation (ENSO) events and water temperature in the top 200 m: the temperature drops about 0.5 °C during warm ENSO events while during cold ENSO events the top 200 m layer of the Barents Sea is warmer. Results from 1 and 1/4-degree global NEMO models show a similar response for the whole Barents Sea. During the strong warm ENSO event in 1997–1998 an anomalous anticyclonic atmospheric circulation over the Barents Sea enhances heat loses, as well as substantially influencing the Barents Sea inflow from the North Atlantic, via changes in ocean currents. Under normal conditions along the Scandinavian peninsula there is a warm current entering the Barents Sea from the North Atlantic, however after the 1997–1998 event this current is weakened. During 1997–1998 the model annual mean temperature in the Barents Sea is decreased by about 0.8 °C, also resulting in a higher sea ice volume. In contrast during the cold ENSO events in 1999–2000 and 2007–2008, the model shows a lower sea ice volume, and higher annual mean temperatures in the upper layer of the Barents Sea of about 0.7 °C. An analysis of model data shows that the strength of the Atlantic inflow in the Barents Sea is the main cause of heat content variability, and is forced by changing pressure and winds in the North Atlantic. However, surface heat-exchange with the atmosphere provides the means by which the Barents sea heat budget relaxes to normal in the subsequent year after the ENSO events.

Geochemistry of iron and manganese in soils and sediments of a mangrove system, Island of Pai Matos (Cananeia - SP, Brazil)

Five zones along a transect of 180 m were selected for study on the Island of Pai Matos (Sao Paulo, Brazil). Four of the zones are colonised by vascular plants (Spartina SP, Laguncularia LG, Avicennia AV and Rhizophora RH) and were denominated soils, and the other zone, which lacks vegetation, was denominated sediment (SD). The geochemical conditions differed significantly in soils and sediment and also at different depths. The soils were oxic (Eh > 350 mV) or suboxic (Eh: 350-100 mV) at the surface and anoxic (Eh < 100 mV) at depth, whereas in the sediment anoxic conditions prevailed at all depths, but with a lower concentration of sulphides in the pore water and pyrite in the solid fraction. Under these geochemical conditions Fe is retained in the soils, while the Mn tends to be mobilized and lost. The most abundant form of iron oxyhydroxide was lepidocrocite (mean concentration for all sites and depths, 45 +/- 19 mu mol g(-1)), followed by goethite (30 19 mu mol g(-1))and ferrihydrite (19 +/- 11 mu mol g(-1)),with significant differences among the mean concentrations. There was a significant decrease with depth in all the types of Fe oxyhydroxides measured, particularly the poorly crystalline forms. The pyrite fraction was an important component of the free Fe pool (non-silicate Fe) in all soils as well as in the sediment, especially below 20 cm depth (mean concentration for all sites and depths, 60 +/- 54 mu mol CI). Furthermore, the mean concentration of Fe-pyrite for all sites and depths was higher than that obtained for any of the three Fe oxyhydroxides measured. The Fe-AVS was a minor fraction, indicating that the high concentrations of dissolved Fe in the soils in the upper area of the transect result from the oxidation of Fe sulphides during low tide. Mossbauer spectroscopy also revealed that most of the Fe (III) was associated with silicates, in this case nontronite. The presence of crystals of pyrite associated with phyllosilicates in samples from the upper layer of the soils may indicate that pyritization of this form of Fe(III) is more rapid than usually reported for ocean bed sediments. The sequential extraction of Mn did not reveal any clearly dominant fraction, with the Mn-carbonate fraction being the most prevalent, followed by exchangeable Mn and oxides of Mn, whereas pyrite-Mn and Mn associated with crystalline Fe-oxides were present at significantly lower concentrations. The high concentration of dissolved Mn found in the soils in the lower part of the transect is consistent with the fact that the solubility is determined by the carbonate fraction. Unlike for Fe, in the soils in the higher zone, which are subject to intense drainage during low tide, there was loss of Mn, as reflected by the concentration of total Mn. (C) 2008 Elsevier B.V. All rights reserved.

Tide-topography interaction along the eastern Brazilian shelf

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

Efeitos da eliminação de neurônios infragranulares sobre a especificação de neurônios supragranulares do córtex cerebral

The cerebral cortex of mammals is histologically organized into different layers of excitatory neurons that have distinct patterns of connections with cortical or subcortical targets. During development, these cortical layers are established through an intricate combination of neuronal specification and migration in a radial pattern known as "insideout": deep-layer neurons are generated prior to upper-layer neurons. In the last few decades, several genes encoding transcription factors involved in the sequential specification of neurons destined to different cortical layers have been identified. However, the influence of early-generated neurons in the specification of subsequent neuronal cohorts remains unclear. To investigate this possible influence, we induced the selective death of cortical neurons from layer V and VI before the generation of layer II, III and IV neurons. Thus, we can evaluate the effects of ablation of early born neurons on the phenotype of late born neurons. Our data shows that one-day after ablation, layer VI neurons expressing the transcription factor TBR1 are newly generated while virtually no neuron expressing TBR1 was generated in the same age in control animals. This suggests that progenitors involved in the generation of neurons destined for superficial layers suffer interference from the selective death of neurons in deep layers, changing their specification. We also observed that while TBR1-positive neurons are located exclusively in deep cortical layers of control animals, many TBR1-positive neurons are misplaced in superficial layers of ablated animals, suggesting that the migration of cortical neurons could be controlled independently of neuronal phenotypes. Furthermore, we observed an increase in layer V neurons expressing CTIP2 and neurons expressing SATB2 and that these cells have changed their distributions. As a conclusion, our data indicate the existence of a mechanism of control exercised by the early-generated neurons in the cerebral cortex on the fate of the progenitors involved in the generation of the following cortical neurons. This mechanism could help to control the number of neurons in different layers and contribute to the establishment of different cortical areas

Differentiation of epithelial cells in the urinary tract

Uroplakins, cytokeratins and the apical plasma membrane were studied in the epithelia of mouse urinary tract. In the simple epithelium covering the inner medulla of the renal pelvis, no uroplakins or cytokeratin 20 were detected and cells had microvilli on their apical surface. The epithelium covering the inner band of the outer medulla became pseudostratified, with the upper layer consisting of large cells with stalks connecting them to the basal lamina. Uroplakins and cytokeratin 20 were not expressed in these cells. However, some superficial cells appeared without connections to the basal lamina; these cells expressed uroplakins Ia, Ib, II and III and cytokeratin 20, they contained sparse small uroplakin-positive cytoplasmic vesicles and their apical surface showed both microvilli and ridges. Cytokeratin 20 was seen as dots in the cytoplasm. This epithelium therefore showed partial urothelial differentiation. The epithelium covering the outer band of the outer medulla gradually changed from a two-layered to a three-layered urothelium with typical umbrella cells that contained all four uroplakins. Cytokeratin 20 was organized into a complex network. The epithelium possessed an asymmetric unit membrane at the apical cell surface and fusiform vesicles. Umbrella cells were also observed in the ureter and urinary bladder. In males and females, the urothelium ended in the bladder neck and was continued by a non-keratinized stratified epithelium in the urethra in which no urothelial cell differentiation markers were detected. We thus show here the expression, distribution and organization of specific proteins associated with the various cell types in the urinary tract epithelium.

Dry matter production of shoots and root density of two cultivars of Lablab purpureus (L.) Sweet

This experiment was conducted in green house conditions to evaluate the DM accumulation in the shoots and in the roots of two cultivars of Lablab purpureus (L.) Sweet. A 2x3 factorial (two cultivars and three evaluation dates) was conducted according to a randomized complete block design with four replications, being the cultivars Highworth and Rongai evaluated at 42, 56, and 70 days after seedling emergence (DASE). The results indicated that the cvs. Highworth and Rongai have the same pattern of DM accumulation in the shoots. In the upper layer of the soil (0-0.20 in) it was found 38.83% and 43.64% of the DM accumulated in the roots down to 2.00 in depth, in the cvs. Highworth and Rongai, respectively. In the deepest layer (1.80-2.00 in) it was found 3.02% and 1.5% of the DM accumulated in the roots of the cvs. Highworth and Rongai, respectively. The root density showed a striking decrease upper layer from the soil (0-0.2 m) down to the depth of 0.60 0.80 in (from 10.83 to 1.75 cm.cm(-3) in the cv. Highworth and from 10.76 to 1.28 cm.cm(-3) in the cv. Rongai). At the bottom layer (1.80-2.00 in) the root density values were 0.98 cm.cm(-3) and 0.59 cm.cm(-3), respectively for the cvs. Highworth and Rongai. The root/shoot ratios were similar in both cvs. and decreased from 42 to 70 DASE showing that the cvs. evaluated had the same dynamics of DM accumulation.

Potential use of the veris apparent EC sensor to predict soil texture under the semi-arid conditions of central Arizona

The operational details of the apparent electrical conductivity (ECa) sensor manufactured by Veris Technologies have been extensively documented in literature reports, but the geographical distribution of these research studies indicate a strong regional concentration in the US Mid-west and Southern states. The agricultural lands of these states diverge significantly to the soil conditions and water regime of irrigated land in the US South-western states such as Arizona where there is no previous research reports of the use of this particular sensor. The objectives of the present study were to analyze the performance of this sensor under the conditions of typical soils in irrigated farms of Central Arizona. We tested under static conditions the performance of the sensor on three soils of contrasting texture. Observations were collected as time series data as soil moisture changed from saturation to permanent wilting point. Observations were repeated at the hours of lowest and highest temperatures. In addition, this study included soil penetration resistance and salinity determinations. Preliminary results indicate that soil temperature of the upper layer caused the most dynamic change in the sensor output. The ECa curves of the three soil textures tested had well defined distinctive characteristics. Final multivariate analysis is pending.

Aspects of Mechanical Behavior and Modeling of a Tropical Unsaturated Soil

The research studies the applicability of two elastoplastic models for the collapse prediction of the lateritic soil profile from Southeastern Brazil. These tropical soils have peculiar geotechnical behavior, due to their mineralogical composition and porous structure coming from intense process of formation. Two elastoplastic models were analyzed: the Barcelona Basic Model (BBM) and another one based on BBM, however developed for tropical soils. Oedometric tests with suction control were performed at three distinct depths of the soil profile. The BBM was not suitable for the upper layer of the soil profile, because BBM considers the compressible behavior of the soil in function of the reduction of the elastoplastic compressibility index with the increase of the matric suction. The model developed for tropical soils showed better suited to the compressible behavior of the soil profile, resulting in good prediction of the collapse potential, mainly by accepting increasing values of the elastoplastic compressibility index of the soil profile with the matric suction rise. © 2013 Springer Science+Business Media Dordrecht.

Lixiviação de sulfentrazone e amicarbazone com adição de óleo em resposta à precipitação e emergência de Ipomea spp. em função da profundidade de semeadura e cobertura com palha

Pós-graduação em Agronomia (Produção Vegetal) - FCAV

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

Estudo da viabilidade de implantação de plantas para conversão de energia térmica do oceano (OTEC) no Brasil

Pós-graduação em Engenharia Mecânica - FEIS

Constructed Wetland for Treating Effluent from Subtropical Aquaculture Farm

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)