Complexity, contingency, and criticality.

Complexity originates from the tendency of large dynamical systems to organize themselves into a critical state, with avalanches or "punctuations" of all sizes. In the critical state, events which would otherwise be uncoupled become correlated. The apparent, historical contingency in many sciences, including geology, biology, and economics, finds a natural interpretation as a self-organized critical phenomenon. These ideas are discussed in the context of simple mathematical models of sandpiles and biological evolution. Insights are gained not only from numerical simulations but also from rigorous mathematical analysis.

Estudo da expansão da Leishmaniose Visceral no estado de São Paulo e de fatores intervenientes

INTRODUÇÃO: Leishmaniose visceral (LV) é uma doença negligenciada que afeta milhões de pessoas no mundo e que constitui um grave problema de saúde pública. OBJETIVOS: Descrever no tempo e no espaço, a dispersão de Lutzomyia longipalpis e a expansão da LV no estado de São Paulo (SP); identificar fatores associados a estes processos. MÉTODOS: Foram realizados estudos descritivos, ecológicos e de análise de sobrevida. Informações sobre o vetor e os casos foram obtidas na Superintendência de Controle de Endemias e no Sistema de Informações de Agravos de Notificação para o período de 1997 a 2014. A área de estudo foi composta pelos 645 municípios de SP. Foram produzidos mapas temáticos e de fluxo e calcularam-se incidência, mortalidade e letalidade por LV em humanos (LVH). Utilizou-se a técnica de análise de sobrevida (Curvas de Kaplan-Meier e Regressão de Cox) para a identificação de fatores associados à dispersão do vetor e expansão da LV. RESULTADOS: A partir da detecção de Lu. longipalpis em Araçatuba em 1997, deram-se a ocorrência do primeiro caso canino (LVC) (1998) e o primeiro caso humano (LVH) autóctones (1999) em SP. Até 2014, foi detectada a presença do vetor em 173 (26,8 por cento ) municípios, LVC em 108 (16,7 por cento ) e LVH em 84 (13 por cento ). A expansão dos três fenômenos ocorreu no sentido noroeste para sudeste e se deram a velocidades constantes. Na região de São José do Rio Preto, a dispersão do vetor deu-se por vizinhança com municípios anteriormente infestados, a expansão da LV relacionou-se com os municípios sede das microrregiões e a doença ocorreu com maior intensidade nas áreas periféricas dos municípios. A presença da Rodovia Marechal Rondon e a divisa com o Mato Grosso do Sul foram fatores associados à ocorrência dos três eventos, assim como a presença da Rodovia Euclides da Cunha para presença do vetor e casos caninos, e, presença de presídios para casos humanos. CONCLUSÕES: A dispersão do vetor e da LV em SP iniciou-se, a partir de 1997, próximo à divisa com o estado do Mato Grosso do Sul, avançou no sentido noroeste para sudeste, na trajetória da rodovia Marechal Rondon, e ocorreu em progressão aritmética, com as sedes das microrregiões de SP tendo papel preponderante neste processo. A ocorrência autóctone de LVC e LVH iniciou-se na sequência da detecção de Lu. longipalpis em Araçatuba e de seu espalhamento por SP e não a partir dos locais onde anteriormente ele já estava presente. O uso da análise de sobrevida permitiu identificar fatores associados à dispersão do vetor e a expansão da LV. Os resultados deste estudo podem ser úteis para aprimorar as atividades de vigilância e controle da LV, no sentido de retardar sua expansão e/ou de mitigar seus efeitos, quando de sua ocorrência.

Modelo numérico mecanobiológico para a obtenção da matriz de rigidez estrutural e da densidade mineral na remodelagem óssea

Neste trabalho é proposto um modelo mecanobiológico de remodelagem óssea para a estimativa de variações, provocadas por perturbações mecânicas ou biológicas, na matriz de rigidez estrutural da escala macroscópica e na densidade mineral em uma região do osso. Na cooperação entre as áreas da saúde e da engenharia, como nos estudos estruturais de biomecânica no sistema esquelético, as propriedades mecânicas dos materiais devem ser conhecidas, entretanto os ossos possuem uma constituição material altamente complexa, dinâmica e variante entre indivíduos. Sua dinâmica decorre dos ciclos de absorção e deposição de matriz óssea na remodelagem óssea, a qual ocorre para manter a integridade estrutural do esqueleto e adaptá-lo aos estímulos do ambiente, sejam eles biológicos, químicos ou mecânicos. Como a remodelagem óssea pode provocar alterações no material do osso, espera-se que suas propriedades mecânicas também sejam alteradas. Na literatura científica há modelos matemáticos que preveem a variação da matriz de rigidez estrutural a partir do estímulo mecânico, porém somente os modelos mais recentes incluíram explicitamente processos biológicos e químicos da remodelagem óssea. A densidade mineral óssea é um importante parâmetro utilizado no diagnóstico de doenças ósseas na área médica. Desse modo, para a obtenção da variação da rigidez estrutural e da densidade mineral óssea, propõe-se um modelo numérico mecanobiológico composto por cinco submodelos: da dinâmica da população de células ósseas, da resposta das células ao estímulo mecânico, da porosidade óssea, da densidade mineral óssea e, baseado na Lei de Voigt para materiais compósitos, da rigidez estrutural. Os valores das constantes das equações dos submodelos foram obtidos de literatura. Para a solução das equações do modelo, propõe-se uma implementação numérica e computacional escrita em linguagem C. O método de Runge-Kutta-Dorman-Prince, cuja vantagem consiste no uso de um passo de solução variável, é utilizado no modelo para controlar o erro numérico do resultado do sistema de equações diferenciais. Foi realizada uma avaliação comparativa entre os resultados obtidos com o modelo proposto e os da literatura dos modelos de remodelagem óssea recentes. Conclui-se que o modelo e a implementação propostos são capazes de obter variações da matriz de rigidez estrutural macroscópica e da densidade mineral óssea decorrentes da perturbação nos parâmetros mecânicos ou biológicos do processo de remodelagem óssea.

Identificação de sistemas para o estudo de controle motor.

Qualquer tarefa motora ativa se dá pela ativação de uma população de unidades motoras. Porém, devido a diversas dificuldades, tanto técnicas quanto éticas, não é possível medir a entrada sináptica dos motoneurônios em humanos. Por essas razões, o uso de modelos computacionais realistas de um núcleo de motoneurônios e as suas respectivas fibras musculares tem um importante papel no estudo do controle humano dos músculos. Entretanto, tais modelos são complexos e uma análise matemática é difícil. Neste texto é apresentada uma abordagem baseada em identificação de sistemas de um modelo realista de um núcleo de unidades motoras, com o objetivo de obter um modelo mais simples capaz de representar a transdução das entradas do núcleo de unidades motoras na força do músculo associado ao núcleo. A identificação de sistemas foi baseada em um algoritmo de mínimos quadrados ortogonal para achar um modelo NARMAX, sendo que a entrada considerada foi a condutância sináptica excitatória dendrítica total dos motoneurônios e a saída foi a força dos músculos produzida pelo núcleo de unidades motoras. O modelo identificado reproduziu o comportamento médio da saída do modelo computacional realista, mesmo para pares de sinal de entrada-saída não usados durante o processo de identificação do modelo, como sinais de força muscular modulados senoidalmente. Funções de resposta em frequência generalizada do núcleo de motoneurônios foram obtidas do modelo NARMAX, e levaram a que se inferisse que oscilações corticais na banda-beta (20 Hz) podem influenciar no controle da geração de força pela medula espinhal, comportamento do núcleo de motoneurônios até então desconhecido.

Modelo preditivo fuzzy da mobilização do solo no sistema de plantio direto.

As culturas do milho e da soja respondem pela maior parte da produção nacional de grãos, predominando o sistema de plantio direto. Para uma semeadura direta de qualidade, o bom aterramento do sulco é indispensável, pois garante um ambiente adequado às sementes. Neste sentido, é importante estimar a mobilização de solo promovida por uma haste sulcadora estreita durante esta operação. O modelo analítico existente visa representar a mobilização do solo no sistema de plantio convencional. Como consequência, há situações em que este não pode se aplicado, como no caso de hastes sulcadoras estreitas utilizadas em semeadoras de plantio direto. Nestas situações, o mecanismo de falha do solo pode se alterar, assumindo um comportamento não modelado na literatura. Essa pesquisa propõe um modelo fuzzy capaz de representar estas situações, aproveitando conhecimento da teoria de mecânica dos solos e da análise de resultados experimentais. No modelo proposto, parte das regras descrevem situações não abrangidas pelo modelo analítico, as quais foram formuladas a partir da estimativa das prováveis áreas de solo mobilizado. O modelo fuzzy foi testado com dados de experimentos conduzidos durante a pesquisa, em duas condições de granulometria de solo (arenoso e argiloso). O modelo proposto reproduziu as tendências observadas nos dados experimentais, mas superestimou os valores de área observados, sendo esse efeito bem mais intenso para os dados do experimento em solo arenoso. A superestimativa ocorreu devido à soma de diversos fatores. Um deles é a diferença entre as leituras experimentais, as quais consideram apenas o solo realmente movimentado, e a premissa do modelo analítico, que considera toda a área de solo incluindo aquela cisalhada, porém não mobilizada. Outro fator foi devido ao efeito do disco de corte da palha, que pré-cisalha o solo à frente da ferramenta. No ensaio em solo arenoso os valores observados de área de solo mobilizado foram menores que os esperados, intensificando o efeito de superestimativa do modelo fuzzy, sendo que este efeito não representa uma deficiência deste modelo.

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.

Comparación de las aplicaciones de Google y Yahoo para la geocodificación de direcciones postales con fines epidemiológicos

Fundamentos: Geocodificar es asignar coordenadas geográficas a puntos del espacio, frecuentemente direcciones postales. El error cometido al aplicar este proceso puede introducir un sesgo en las estimaciones de modelos espacio-temporales en estudios epidemiológicos. No se han encontrado estudios que midan este error en ciudades españolas. El objetivo es evaluar los errores en magnitud y direccionalidad de dos recursos gratuitos (Google y Yahoo) respecto a GPS en dos ciudades de España. Método: Se geocodificaron 30 direcciones aleatorias con los dos recursos citados y con GPS en Santa Pola (Alicante) y en Alicante. Se calculó la mediana y su IC95% del error en metros entre los recursos y GPS, para el total y por el status reportado. Se evaluó la direccionalidad del error calculando el cuadrante de localización y aplicando un test Chi-Cuadrado. Se evaluó el error del GPS midiendo 11 direcciones dos veces en un intervalo de 4 días. Resultados: La mediana del error total desde Google-GPS fue de 23,2 metros (16,0-32,2) para Santa Pola y 21,4 metros (14,9-31,1) en Alicante. Para Yahoo fue de 136,0 (19,2-318,5) para Santa Pola y 23,8 (13,6-29,2) para Alicante. Por status, se geocodificó entre un 73% y 90% como ‘exactas o interpoladas’ (menor error), tanto Google como Yahoo tuvieron una mediana del error de entre 19 y 22 metros en las dos ciudades. El error del GPS fue de 13,8 (6,7-17,8) metros. No se detectó direccionalidad. Conclusiones: El error de Google es asumible y estable en las dos ciudades, siendo un recurso fiable para geocodificar direcciones postales en España en estudios epidemiológicos.

Integration of different models in the design of chemical processes: Application to the design of a power plant

With advances in the synthesis and design of chemical processes there is an increasing need for more complex mathematical models with which to screen the alternatives that constitute accurate and reliable process models. Despite the wide availability of sophisticated tools for simulation, optimization and synthesis of chemical processes, the user is frequently interested in using the ‘best available model’. However, in practice, these models are usually little more than a black box with a rigid input–output structure. In this paper we propose to tackle all these models using generalized disjunctive programming to capture the numerical characteristics of each model (in equation form, modular, noisy, etc.) and to deal with each of them according to their individual characteristics. The result is a hybrid modular–equation based approach that allows synthesizing complex processes using different models in a robust and reliable way. The capabilities of the proposed approach are discussed with a case study: the design of a utility system power plant that has been decomposed into its constitutive elements, each treated differently numerically. And finally, numerical results and conclusions are presented.

A methodology for evacuation route planning inside buildings using geospatial technology

Evacuation route planning is a fundamental task for building engineering projects. Safety regulations are established so that all occupants are driven on time out of a building to a secure place when faced with an emergency situation. As an example, Spanish building code requires the planning of evacuation routes on large and, usually, public buildings. Engineers often plan these routes on single building projects, repeatedly assigning clusters of rooms to each emergency exit in a trial-and-error process. But problems may arise for a building complex where distribution and use changes make visual analysis cumbersome and sometimes unfeasible. This problem could be solved by using well-known spatial analysis techniques, implemented as a specialized software able to partially emulate engineer reasoning. In this paper we propose and test an easily reproducible methodology that makes use of free and open source software components for solving a case study. We ran a complete test on a building floor at the University of Alicante (Spain). This institution offers a web service (WFS) that allows retrieval of 2D geometries from any building within its campus. We demonstrate how geospatial technologies and computational geometry algorithms can be used for automating the creation and optimization of evacuation routes. In our case study, the engineers’ task is to verify that the load capacity of each emergency exit does not exceed the standards specified by Spain’s current regulations. Using Dijkstra’s algorithm, we obtain the shortest paths from every room to the most appropriate emergency exit. Once these paths are calculated, engineers can run simulations and validate, based on path statistics, different cluster configurations. Techniques and tools applied in this research would be helpful in the design and risk management phases of any complex building project.

On the decorrelation filtering of RL05 GRACE data for global applications

In autumn 2012, the new release 05 (RL05) of monthly geopotencial spherical harmonics Stokes coefficients (SC) from GRACE (Gravity Recovery and Climate Experiment) mission was published. This release reduces the noise in high degree and order SC, but they still need to be filtered. One of the most common filtering processing is the combination of decorrelation and Gaussian filters. Both of them are parameters dependent and must be tuned by the users. Previous studies have analyzed the parameters choice for the RL05 GRACE data for oceanic applications, and for RL04 data for global application. This study updates the latter for RL05 data extending the statistics analysis. The choice of the parameters of the decorrelation filter has been optimized to: (1) balance the noise reduction and the geophysical signal attenuation produced by the filtering process; (2) minimize the differences between GRACE and model-based data; (3) maximize the ratio of variability between continents and oceans. The Gaussian filter has been optimized following the latter criteria. Besides, an anisotropic filter, the fan filter, has been analyzed as an alternative to the Gauss filter, producing better statistics.

Saint Mathew Law and Bonini Paradox in Textual Theory of Complex Models

The mathematical models of the complex reality are texts belonging to a certain literature that is written in a semi-formal language, denominated L(MT) by the authors whose laws linguistic mathematics have been previously defined. This text possesses linguistic entropy that is the reflection of the physical entropy of the processes of real world that said text describes. Through the temperature of information defined by Mandelbrot, the authors begin a text-reality thermodynamic theory that drives to the existence of information attractors, or highly structured point, settling down a heterogeneity of the space text, the same one that of ontologic space, completing the well-known law of Saint Mathew, of the General Theory of Systems and formulated by Margalef saying: “To the one that has more he will be given, and to the one that doesn't have he will even be removed it little that it possesses.

Twenty-year advanced DInSAR analysis of severe land subsidence: The Alto Guadalentín Basin (Spain) case study

A twenty-year period of severe land subsidence evolution in the Alto Guadalentín Basin (southeast Spain) is monitored using multi-sensor SAR images, processed by advanced differential interferometric synthetic aperture radar (DInSAR) techniques. The SAR images used in this study consist of four datasets acquired by ERS-1/2, ENVISAT, ALOS and COSMO-SkyMed satellites between 1992 and 2012. The integration of ground surface displacement maps retrieved for different time periods allows us to quantify up to 2.50 m of cumulated displacements that occurred between 1992 and 2012 in the Alto Guadalentín Basin. DInSAR results were locally compared with global positioning system (GPS) data available for two continuous stations located in the study area, demonstrating the high consistency of local vertical motion measurements between the two different surveying techniques. An average absolute error of 4.6 ± 4 mm for the ALOS data and of 4.8 ± 3.5 mm for the COSMO-SkyMed data confirmed the reliability of the analysis. The spatial analysis of DInSAR ground surface displacement reveals a direct correlation with the thickness of the compressible alluvial deposits. Detected ground subsidence in the past 20 years is most likely a consequence of a 100–200 m groundwater level drop that has persisted since the 1970s due to the overexploitation of the Alto Guadalentín aquifer system. The negative gradient of the pore pressure is responsible for the extremely slow consolidation of a very thick (> 100 m) layer of fine-grained silt and clay layers with low vertical hydraulic permeability (approximately 50 mm/h) wherein the maximum settlement has still not been reached.

Logical Modeling and Dynamical Analysis of Cellular Networks

The logical (or logic) formalism is increasingly used to model regulatory and signaling networks. Complementing these applications, several groups contributed various methods and tools to support the definition and analysis of logical models. After an introduction to the logical modeling framework and to several of its variants, we review here a number of recent methodological advances to ease the analysis of large and intricate networks. In particular, we survey approaches to determine model attractors and their reachability properties, to assess the dynamical impact of variations of external signals, and to consistently reduce large models. To illustrate these developments, we further consider several published logical models for two important biological processes, namely the differentiation of T helper cells and the control of mammalian cell cycle.

Comparisons of models predicting ambient lake phosphorus concentrations /

Mode of access: Internet.

Comparison of numerical and physical hydraulic models, Masonboro Inlet, North Carolina /

"June 1977."