Effluent from on-site sewage treatment and disposal systems (OSTDS) is generally known to impact groundwaters and surface waters with nitrogen (N) and phosphorus (P) and other contaminants. Little research has quantified this problem along the Indian River Lagoon (IRL), especially in Indian River County (IRC) where there are 26,660 active systems. This study assessed the effects of OSTDS on contamination of surface and groundwaters along three urbanized canals and the St. Sebastian River in IRC, all of which flow into the Central IRL. Multiple lines of evidence were used to define the source of the nutrient loadings including the novel approach of using the artificial sweetener, sucralose, as an indicator of human sewage impact. Results indicate that areas with high densities of OSTDS are contributing N to surface waters and elevating N:P ratios through submarine groundwater discharge and promoting eutrophication in the Central IRL. / Includes bibliography. / Thesis (M.S.)--Florida Atlantic University, 2014. / FAU Electronic Theses and Dissertations Collection

L’identification des efforts vibratoires agissant sur les structures et la détection des défauts à partir des mesures opérationnelles sont des sujets importants dans des projets académiques et industriels. Le choix de l’outil ou de la méthode utilisée dépend de la bande de fréquences d'étude puisqu’il existe des approches appropriées pour chaque domaine fréquentiel. Une méthode énergétique appelée la méthode énergétique simplifiée (MES) est utilisée pour prédire la répartition de la densité d'énergie vibroacoustique en moyennes et hautes fréquences. L'objectif de ce mémoire est d'étendre cette méthode pour résoudre les problèmes vibro-acoustiques inverses pour identifier d'une part les sources de vibration et d'autre part pour détecter les défauts. La formulation MES inverse (IMES) est numériquement validée pour des systèmes continus basés sur le couplage tel que le cas d’un système comportant plusieurs plaques couplées et celui d’un système composé d’une cavité acoustique couplée avec une plaque. En plus, une nouvelle méthodologie numérique est proposée, pour étendre cet outil d'identification IMES pour la détection des défauts. Une analyse paramétrique est effectuée pour le cas d’un modèle présentant des défauts afin de tester la robustesse et l’efficacité de cette approche. Finalement, une étude expérimentale est effectuée pour valider la technique IMES à fin d'identifier et localiser les charges exercées pour plusieurs cas, et détecter les défauts. / The identification of inputs forces acting on structures and the detection of defects from operating measurement have been important topics in both academic and industrial projects. The choice of the used tool or method depends on the frequency band of study since there are appropriate approaches for each frequency domain. An energetic method so called the simplified energy method (MES) is used to predict the distribution of the vibro-acoustic energy density in the medium and high frequency band. The objective of this thesis is to extend this energy method to solve inverse vibro-acoustic problems and to identify the sources of vibrations on one hand and to detect the defects on the other hand. The inverse MES formulation (IMES) is numerically validated for continuous coupling-based systems such as the case of a system composed with several coupled plates and the case of a system composed of an acoustic cavity coupled with a plate. In addition, a new numerical methodology is proposed to extend this IMES identification tool for the detection of defects. A parametric analysis is performed in the case of plate with defects in order to test the robustness and the efficiency of this approach. Finally, an experimental study is carried out to validate the IMES technique to identify and locate the input loads for several scenarios, and detecting the defects.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Em engenharia, a modelagem computacional desempenha um papel importante na concepção de produtos e no desenvolvimento de técnicas de atenuação de ruído. Nesse contexto, esta tese investiga a intensidade acústica gerada pela radiação sonora de superfícies vibrantes. De modo específico, a pesquisa enfoca a identificação das regiões de uma fonte sonora que contribuem efetivamente para potência sonora radiada para o campo afastado, quando a frequência de excitação ocorre abaixo da frequência crítica de coincidência. São descritas as fundamentações teóricas de duas diferentes abordagens. A primeira delas, denominada intensidade supersônica (analítica) é calculada via transformadas de Fourier para fontes sonoras com geometrias separáveis. A segunda, denominada intensidade útil (numérica) é calculada através do método dos elementos de contorno clássico para fontes com geometrias arbitrárias. Em ambas, a identificação das regiões é feita pela filtragem das ondas não propagantes (evanescentes). O trabalho está centrado em duas propostas, a saber. A primeira delas, é a apresentação implementação e análise de uma nova técnica numérica para o cálculo da grandeza intensidade útil. Essa técnica constitui uma variante do método dos elementos de contorno (MEC), tendo como base o fato de as aproximações para as variáveis acústicas pressão e velocidade normal serem tomadas como constantes em cada elemento. E também no modo peculiar de obter a velocidade constante através da média de um certo número de velocidades interiores a cada elemento. Por esse motivo, a técnica recebe o nome de método de elemento de contorno com velocidade média (AVBEMAverage Velocity Boundary Element Method). A segunda, é a obtenção da solução forma fechada do campo de velocidade normal para placas retangulares com oito diferentes combinações de condições contorno clássicas. Então, a intensidade supersônica é estimada e comparada à intensidade acústica. Nos ensaios numéricos, a comparação da intensidade útil obtida via MEC clássico e via AVBEM é mostrada para ilustrar a eficiência computacional da técnica aqui proposta, que traz como benefício adicional o fato de poder ser utilizada uma malha menos refinada para as simulações e, consequentemente, economia significativa de recursos computacionais. / In engineering computational modeling plays an important role in product design and development techniques for noise attenuation. In such context, this thesis investigates the acoustic intensity generated by sound radiation from vibrating surfaces. Specifically, the research focuses on the identification that effectively contribute to the radiated sound power into the far-field, especially when the driven frequency occurs below the critical coincidence frequency. The theoretical formulations of two different approaches are described. The first one, called supersonic intensity (analytical) is calculated via Fourier transforms for noise sources with separable geometries. The second, called useful intensity (numerical) is calculated by the method of boundary elements to classic fonts with arbitrary geometries. In both, the identification of regions is done by filtering the non propagating waves (evanescent). The work is focused on the presentation, implementation and analysis of a new numerical technique for calculating the magnitude useful intensity. This technique is a variant of the method-boundary element (BEM), based on the fact that the approximations for the acoustic pressure and normal speed variables are taken as constant in each element. Also in particular way to obtain constant speed by mean of a number of speeds interior of each element. For this reason, the technique is called the Average Velocity Boundary Element Method . The closed form solution of the normal velocity field for rectangular plates in eight cases with distinct combinations of classical boundary conditions is also obtained. Then, the supersonic intensity is estimated and compared to the acoustic intensity. Numerical experiments are performed comparing to the useful intensity via the conventional BEM with theAVBEM in order to illustrate the positive computational features of the method. The numerical results indicate that the proposed method is much more computationally effcient than its standard BEM counterpart, it enabling the use of a coarser mesh.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Em engenharia, a modelagem computacional desempenha um papel importante na concepção de produtos e no desenvolvimento de técnicas de atenuação de ruído. Nesse contexto, esta tese investiga a intensidade acústica gerada pela radiação sonora de superfícies vibrantes. De modo específico, a pesquisa enfoca a identificação das regiões de uma fonte sonora que contribuem efetivamente para potência sonora radiada para o campo afastado, quando a frequência de excitação ocorre abaixo da frequência crítica de coincidência. São descritas as fundamentações teóricas de duas diferentes abordagens. A primeira delas, denominada intensidade supersônica (analítica) é calculada via transformadas de Fourier para fontes sonoras com geometrias separáveis. A segunda, denominada intensidade útil (numérica) é calculada através do método dos elementos de contorno clássico para fontes com geometrias arbitrárias. Em ambas, a identificação das regiões é feita pela filtragem das ondas não propagantes (evanescentes). O trabalho está centrado em duas propostas, a saber. A primeira delas, é a apresentação implementação e análise de uma nova técnica numérica para o cálculo da grandeza intensidade útil. Essa técnica constitui uma variante do método dos elementos de contorno (MEC), tendo como base o fato de as aproximações para as variáveis acústicas pressão e velocidade normal serem tomadas como constantes em cada elemento. E também no modo peculiar de obter a velocidade constante através da média de um certo número de velocidades interiores a cada elemento. Por esse motivo, a técnica recebe o nome de método de elemento de contorno com velocidade média (AVBEMAverage Velocity Boundary Element Method). A segunda, é a obtenção da solução forma fechada do campo de velocidade normal para placas retangulares com oito diferentes combinações de condições contorno clássicas. Então, a intensidade supersônica é estimada e comparada à intensidade acústica. Nos ensaios numéricos, a comparação da intensidade útil obtida via MEC clássico e via AVBEM é mostrada para ilustrar a eficiência computacional da técnica aqui proposta, que traz como benefício adicional o fato de poder ser utilizada uma malha menos refinada para as simulações e, consequentemente, economia significativa de recursos computacionais. / In engineering computational modeling plays an important role in product design and development techniques for noise attenuation. In such context, this thesis investigates the acoustic intensity generated by sound radiation from vibrating surfaces. Specifically, the research focuses on the identification that effectively contribute to the radiated sound power into the far-field, especially when the driven frequency occurs below the critical coincidence frequency. The theoretical formulations of two different approaches are described. The first one, called supersonic intensity (analytical) is calculated via Fourier transforms for noise sources with separable geometries. The second, called useful intensity (numerical) is calculated by the method of boundary elements to classic fonts with arbitrary geometries. In both, the identification of regions is done by filtering the non propagating waves (evanescent). The work is focused on the presentation, implementation and analysis of a new numerical technique for calculating the magnitude useful intensity. This technique is a variant of the method-boundary element (BEM), based on the fact that the approximations for the acoustic pressure and normal speed variables are taken as constant in each element. Also in particular way to obtain constant speed by mean of a number of speeds interior of each element. For this reason, the technique is called the Average Velocity Boundary Element Method . The closed form solution of the normal velocity field for rectangular plates in eight cases with distinct combinations of classical boundary conditions is also obtained. Then, the supersonic intensity is estimated and compared to the acoustic intensity. Numerical experiments are performed comparing to the useful intensity via the conventional BEM with theAVBEM in order to illustrate the positive computational features of the method. The numerical results indicate that the proposed method is much more computationally effcient than its standard BEM counterpart, it enabling the use of a coarser mesh.

Le but de cette thèse est de modéliser la propagation d’un contaminant au sein d’un réseau d’eau potable muni de capteurs temps réel. Elle comporte les trois axes de développement suivant: la résolution des équations de transport, celle du problème d’identification des sources de contamination et le placement des capteurs.Le transport d’un produit chimique est modélisé dans un réseau d’eau potable par l’équation de transport réaction 1-D avec l’hypothèse de mélange parfait aux noeuds. Il est proposé d’améliorer sa prédiction par l’ajout d’un modèle de mélange imparfait aux jonctions double T et d’un modèle de dispersion prenant en compte un profil de vitesse 3-D et la diffusion radiale. Le premier modèle est créé à l’aide d’un plan d’expériences avec triangulation de Delaunay, de simulations CFD, et de la méthode d’interpolation krigeage. Le second utilise les équations adjointes du problème de transport avec l’ajout de particules évoluant à l’aide d’une marche aléatoire, cette dernière modélisant la diffusion radiale dans la surface droite du tuyau.Le problème d’identification des sources consiste, à l’aide de réponses positives ou négatives à la contamination des noeuds capteurs, à trouver l’origine, le temps d’injection et la durée de la contamination. La résolution de ce problème inverse est faite par la résolution des équations de transport adjointes par formulation backtracking. La méthode donne la liste des sources potentielles ainsi que le classement de celles-ci selon leur probabilité d’être la vraie source de contamination. Elle s’exprime en fonction de combien, en pourcentage, cette source potentielle peut expliquer les réponses positives aux capteurs.Le placement des capteurs est optimisé pour l’identification des sources. L’objectif est la maximisation du potentiel de détection de la véritable source de contamination. Deux résolutions sont testées. La première utilise un algorithme glouton combiné à une méthode de Monte Carlo.La seconde utilise une méthode de recherche locale sur graphe.Finalement les méthodes sont appliquées à un cas test réel avec dans l’ordre : le placement des capteurs, l’identification de la source de contamination et l’estimation de sa propagation. / The aim of this thesis is to model the propagation of a contaminant inside a water distribution network equipped with real time sensors. There are three research directions: the solving of the transport equations, the source identification and the sensor placement. Classical model for transport of a chemical product in a water distribution network isusing 1D-advection-reaction equations with the hypothesis of perfect mixing at junctions. It isproposed to improve the predictions by adding a model of imperfect mixing at double T-junctions and by considering dispersion effect in pipes which takes into account a 3-D velocity profile. The first enhancement is created with the help of a design of experiment based on the Delaunay triangulation, CFD simulations and the interpolation method Kriging. The second one uses the adjoint formulation of the transport equations applied with an algorithm of particle backtracking and a random walk, which models the radial diffusion in the cross-section of a pipe.The source identification problem consists in finding the contamination origin, itsinjection time and its duration from positive and negative responses given by the sensors. The solution to this inverse problem is computed by solving the adjoint transport equations with a backtracking formulation. The method gives a list of potential sources and the ranking of thosemore likely to be the real sources of contamination. It is function of how much, in percentage, they can explain the positive responses of the sensors.The sensor placement is chosen in order to maximize the ranking of the real source of contamination among the potential sources. Two solutions are proposed. The first one uses agreedy algorithm combined with a Monte Carlo method. The second one uses a local search method on graphs. Finally the methods are applied to a real test case in the following order: the sensor placement, the source identification and the estimation of the contamination propagation.

With the current proliferation of cyber attacks, safeguarding internet facing assets from network intrusions, is becoming a vital task in our increasingly digitalised economies. Although recent successes of machine learning (ML) models bode the dawn of a new generation of intrusion detection systems (IDS); current solutions struggle to implement these in an efficient manner, leaving many IDSs to rely on rule-based techniques. In this paper we begin by reviewing the different approaches to feature construction and attack source identification employed in such applications. We refer to these steps as the framework within which models are implemented, and use it as a prism through which we can identify the challenges different solutions face, when applied in modern network traffic conditions. Specifically, we discuss how the most popular framework -- the so called flow-based approach -- suffers from significant overhead being introduced by its resource heavy pre-processing step. To address these issues, we propose the Information Theoretic Framework for Network Anomaly Detection (ITF-NAD); whose purpose is to facilitate online application of statistical learning models onto high-speed network links, as well as provide a method of identifying the sources of traffic anomalies. Its development was inspired by previous work on information theoretic-based anomaly and outlier detection, and employs modern techniques of entropy estimation over data streams. Furthermore, a case study of the framework's detection performance over 5 different types of Denial of Service (DoS) attacks is undertaken, in order to illustrate its potential use for intrusion detection and mitigation. The case study resulted in state-of-the-art performance for time-anomaly detection of single source as well as distributed attacks, and show promising results regarding its ability to identify underlying sources. / I takt med att antalet cyberattacker växer snabbt blir det alltmer viktigt för våra digitaliserade ekonomier att skydda uppkopplade verksamheter från nätverksintrång. Maskininlärning (ML) porträtteras som ett kraftfullt alternativ till konventionella regelbaserade lösningar och dess anmärkningsvärda framgångar bådar för en ny generation detekteringssytem mot intrång (IDS). Trots denna utveckling, bygger många IDS:er fortfarande på signaturbaserade metoder, vilket förklaras av de stora svagheter som präglar många ML-baserade lösningar. I detta arbete utgår vi från en granskning av nuvarande forskning kring tillämpningen av ML för intrångsdetektering, med fokus på de nödvändiga steg som omger modellernas implementation inom IDS. Genom att sätta upp ett ramverk för hur variabler konstrueras och identifiering av attackkällor (ASI) utförs i olika lösningar, kan vi identifiera de flaskhalsar och begränsningar som förhindrar deras praktiska implementation. Särskild vikt läggs vid analysen av de populära flödesbaserade modellerna, vars resurskrävande bearbetning av rådata leder till signifikant tidsfördröjning, vilket omöjliggör deras användning i realtidssystem. För att bemöta dessa svagheter föreslår vi ett nytt ramverk -- det informationsteoretiska ramverket för detektering av nätverksanomalier (ITF-NAD) -- vars syfte är att möjliggöra direktanslutning av ML-modeller över nätverkslänkar med höghastighetstrafik, samt tillhandahåller en metod för identifiering av de bakomliggande källorna till attacken. Ramverket bygger på modern entropiestimeringsteknik, designad för att tillämpas över dataströmmar, samt en ASI-metod inspirerad av entropibaserad detektering av avvikande punkter i kategoriska rum. Utöver detta presenteras en studie av ramverkets prestanda över verklig internettrafik, vilken innehåller 5 olika typer av överbelastningsattacker (DoS) genererad från populära DDoS-verktyg, vilket i sin tur illustrerar ramverkets användning med en enkel semi-övervakad ML-modell. Resultaten visar på hög nivå av noggrannhet för detektion av samtliga attacktyper samt lovande prestanda gällande ramverkets förmåga att identifiera de bakomliggande aktörerna.

No description available.

Submitted by Biblioteca de Pós-Graduação em Geoquímica BGQ (bgq@ndc.uff.br) on 2016-05-02T17:58:54Z No. of bitstreams: 1 GALLOTTA_30_06_14.pdf: 62695484 bytes, checksum: 33335ee3affc2bcffb11fc10caf9ae3f (MD5) / Made available in DSpace on 2016-05-02T17:58:54Z (GMT). No. of bitstreams: 1 GALLOTTA_30_06_14.pdf: 62695484 bytes, checksum: 33335ee3affc2bcffb11fc10caf9ae3f (MD5) / Universidade Federal Fluminense. Instituto de Química. Programa de Pós-Graduação em Geociências-Geoquímica. Niterói, RJ / Os hidrocarbonetos presentes no ambiente consistem em misturas complexas de compostos derivados de múltiplas fontes. Os combustíveis fósseis representam a principal contribuição, devido à taxa e escala espacial, em que o petróleo tem sido usado como fonte de energia e matéria-prima para a indústria química. O objetivo deste estudo foi avaliar os níveis de concentração e identificar fontes de hidrocarbonetos na Bacia do Alto Iguaçu e, em especial, na área de influência da Refinaria Presidente Getúlio Vargas (REPAR). Além dos fatores de poluição crônica, a área de estudo foi alvo de um derrame acidental de petróleo em julho de 2000. Diversos indicadores em diferentes compartimentos ambientais (água superficial, sedimento, solo e água subterrânea) foram avaliados na fase pós-derrame e no monitoramento ambiental ao longo de várias campanhas por mais de uma década. Os esforços de avaliação foram concentrados nos seguintes indicadores: os nalcanos, os alcanos isoprenoides, os hidrocarbonetos policíclicos aromáticos (HPA), os biomarcadores de petróleo e ainda o total de hidrocarbonetos de petróleo (THP). Os resultados identificaram como principal aporte de matéria orgânica para as águas superficiais dos rios Barigüi e Iguaçu, os fluxos materiais originados em região a montante do acidente, refletindo a contribuição antropogênica crônica da cidade de Curitiba. Em 2007 e 2008, as concentrações de THP e HPA nas águas superficiais e sedimentos dos rios Barigüi e Iguaçu refletem uma expressiva melhoria nas condições desses rios em relação a 2000. Na área interna da refinaria (Ponto Zero), observou-se uma nítida diminuição (atenuação) natural das concentrações de THP no solo em todas as profundidades dos perfis amostrados ao longo do tempo nos Banhados 1 e 4 e, em particular, uma diminuição importante das concentrações entre 2004 e 2007. Para a identificação das fontes de hidrocarbonetos foram utilizadas razões diagnósticas e quimiometria. As razões diagnósticas calculadas a partir de concentrações de HPA sugeriram que, na maioria dos sedimentos dos Rios Barigüi e Iguaçu coletados nas campanhas de 2000 e 2001, a fonte petrogênica é a principal. Somente na estação a montante do acidente no Rio Barigüi, a fonte pirolítica predominou nestas duas campanhas. As razões diagnósticas que apresentaram maior eficiência na identificação de fontes de hidrocarbonetos nos sedimentos dos rios Barigüi e Iguaçu foram: ΣC1-Fenantrenos/Fenantreno; e (ΣHPA parentais de 3-6 anéis)/(Σ5 séries de HPA alquilados). A identificação de fontes através de razões diagnósticas calculadas a partir de áreas e alturas de picos cromatográficos demonstrou sua aplicabilidade verificando a relação entre os compostos encontrados em amostras de solo da área interna da refinaria com a amostra de petróleo derramado no acidente, após quase uma década da ocorrência do vazamento. A identificação de fontes através método de quimiométrico baseado na análise de componentes principais (ACP) de seções pré-processadas e combinadas dos Cromatogramas de Íons Selecionados (CIS) mostrou que as amostras mais contaminadas estão na área interna da refinaria. Essas amostras apresentam um padrão de distribuição petrogênica e diferentes graus de intemperismo. As amostras da área externa à refinaria (Guajuvira, General Lúcio e Balsa Nova) são menos ou não contaminadas e/ou contém uma mistura de contribuições diagenéticas, pirolíticas e petrogênicas onde predominam diferentes proporções. Os locais mais distantes da atividade industrial (Balsa Nova) contem, como esperado, os níveis mais baixos de contaminação por HPA. Os resultados de biomarcadores demonstraram que não há evidências para concluir que as amostras da área externa à refinaria e o óleo Cusiana vazado tenham a mesma origem. Os resultados ao longo dos rios Barigüi e Iguaçu e do Ponto Zero demonstraram que as ações de emergência para a contenção do óleo foram adequadas para os rios, e que a contaminação decorrente do derrame ficou predominantemente contida no Ponto Zero e diminuiu significativamente após uma década. / Hydrocarbons present in the environment consist of complex mixtures of compounds derived from multiple sources. The main contribution lies on fossil fuel inputs due to the rate and spatial scale by which petroleum has been used as an energy source and chemical feedstock. The aim of this study was to assess the concentration levels and identify sources of hydrocarbons in the Upper Iguaçu Watershed and, in particular, in the area of influence of the President Getulio Vargas Refinery (REPAR). In addition to the factors of chronic pollution, the study area was the scenario of an acute accidental oil spill in July 2000. Numerous indicators in different environmental compartments (surface water, sediment, soil and groundwater) were assessed in the post spill phase and during the environmental monitoring programs over the course of several campaigns for more than a decade. Assessment efforts were concentrated on the following indicators: n-alkanes, alkanes isoprenoids, polycyclic aromatic hydrocarbons (PAH), petroleum biomarkers and total petroleum hydrocarbons (TPH). The results identified as the main contribution of organic matter to surface waters of the Barigüi and Iguaçu Rivers the materials flows originated in the region upstream of the accident, reflecting chronic anthropogenic contribution of the city of Curitiba. In 2007 and 2008, the TPH and PAH concentrations in surface waters and sediments of the Barigüi and Iguaçu Rivers revealed a significant improvement in the conditions of these rivers when compared with 2000. Inside the refinery area (Point Zero), it was observed a clear natural decrease (attenuation) of the concentrations of TPH in the soil at all depths sampled over time in Marshes 1 and 4 and, in particular, an important decrease of concentrations between 2004 and 2007. Diagnostic ratios and chemometrics were used to identity hydrocarbon sources. The diagnostic ratios calculated from the concentrations of PAH suggested that, in the majority of sediments from the Barigüi and Iguaçu Rivers collected in 2000 and 2001 campaigns, the main source is petrogenic. Only in the station upstream the accident in the Barigüi River, the pyrolytic source predominated in these two campaigns. The diagnostic ratios that presented higher efficiency in identifying sources of hydrocarbons in sediments of the Barigüi and Iguaçu Rivers were: ΣC1- Phenantrenes/Phenanthrene; and (Σ3-6 rings parental PAH)/(Σ5 alkylated PAH series). The source identification through diagnostic ratios calculated from heights and areas of chromatographic peaks demonstrated its applicability establishing a relationship between the compounds found in soil samples of the internal area of the refinery and the sample of the oil spilled in the accident, after nearly a decade of occurrence of the spillage. The source identification through chemometric method based on principal components analysis (PCA) of preprocessed and combined sections of Selected Ion Chromatograms (SIC) showed that the most contaminated samples are inside the refinery area. These samples present a petrogenic pattern and different weathering degrees. Samples from outside the refinery area (Guajuvira, General Lúcio e Balsa Nova) are either less or not contaminated, and/or contain mixtures of diagenetic, pyrogenic and petrogenic inputs where different proportions predominate. The locations farthest away from industrial activity (Balsa Nova) contain, as expected, the lowest levels of PAH contamination. The biomarkers results do not show any evidences to conclude positive matches between the samples from outside the refinery area and the spilled Cusiana oil. The results along Barigüi and Iguaçu rivers and Point Zero demonstrated that emergency actions to contain the oil were appropriate for the rivers, and that the contamination resulting from the spill was mostly contained in the Point Zero and decreased significantly after one decade