Wet and dry spell characteristics of global tropical rainfall

In this study, we analyze satellite-based daily rainfall observations to compare and contrast the wet and dry spell characteristics of tropical rainfall. Defining a wet (dry) spell as the number of consecutive rainy (nonrainy) days, we find that the distributions of wet spells appear to exhibit universality in the following sense. While both ocean and land regions with high seasonal rainfall accumulation (humid regions; e. g., India, Amazon, Pacific Ocean) show a predominance of 2-4 day wet spells, those regions with low seasonal rainfall accumulation (arid regions; e. g., South Atlantic, South Australia) exhibit a wet spell duration distribution that is essentially exponential in nature, with a peak at 1 day. The behavior that we observed for wet spells is reversed for the dry spell characteristics. In other words, the main contribution to the dry part of the season, in terms of the number of nonrainy days, appears to come from 3-4 day dry spells in the arid regions, as opposed to 1 day dry spells in the humid regions. The total rainfall accumulated in each wet spell has also been analyzed, and we find that the major contribution to seasonal rainfall for arid regions comes from 1-5 day wet spells; however, for humid regions, this contribution comes from wet spells of duration as long as 30 days. We also explore the role of chance as well as the influence of organized convection in determining some of the observed features.

Aptamers as Theranostic Agents: Modifications, Serum Stability and Functionalisation

Aptamers, and the selection process known as Systematic Evolution of Ligands by Exponential Enrichment (SELEX) used to generate them, were first described more than twenty years ago. Since then, there have been numerous modifications to the selection procedures. This review discusses the use of modified bases as a means of enhancing serum stability and producing effective therapeutic tools, as well as functionalising these nucleic acids to be used as potential diagnostic agents.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Perturbation of the stability of amyloid fibrils through alteration of electrostatic interactions.

The self-assembly of proteins and peptides into polymeric amyloid fibrils is a process that has important implications ranging from the understanding of protein misfolding disorders to the discovery of novel nanobiomaterials. In this study, we probe the stability of fibrils prepared at pH 2.0 and composed of the protein insulin by manipulating electrostatic interactions within the fibril architecture. We demonstrate that strong electrostatic repulsion is sufficient to disrupt the hydrogen-bonded, cross-β network that links insulin molecules and ultimately results in fibril dissociation. The extent of this dissociation correlates well with predictions for colloidal models considering the net global charge of the polypeptide chain, although the kinetics of the process is regulated by the charge state of a single amino acid. We found the fibrils to be maximally stable under their formation conditions. Partial disruption of the cross-β network under conditions where the fibrils remain intact leads to a reduction in their stability. Together, these results support the contention that a major determinant of amyloid stability stems from the interactions in the structured core, and show how the control of electrostatic interactions can be used to characterize the factors that modulate fibril stability.

Análise global da estabilidade termodinâmica de misturas: um estudo com o método do conjunto gerador

O cálculo do equilíbrio de fases é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Mas para resolvê-lo é aconselhável que se estude a priori a estabilidade termodinâmica do sistema, a qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. Tal problema pode ser abordado como um problema de otimização, conhecido como a minimização da função distância do plano tangente à energia livre de Gibbs molar, onde modelos termodinâmicos, de natureza não convexa e não linear, são utilizados para descrevê-lo. Esse fato tem motivado um grande interesse em técnicas de otimização robustas e eficientes para a resolução de problemas relacionados com a termodinâmica do equilíbrio de fases. Como tem sido ressaltado na literatura, para proporcionar uma completa predição do equilíbrio de fases, faz-se necessário não apenas a determinação do minimizador global da função objetivo do teste de estabilidade, mas também a obtenção de todos os seus pontos estacionários. Assim, o desenvolvimento de metodologias para essa tarefa desafiadora tem se tornado uma nova área de pesquisa da otimização global aplicada à termodinâmica do equilíbrio, com interesses comuns na engenharia química e na engenharia do petróleo. O foco do presente trabalho é uma nova metodologia para resolver o problema do teste de estabilidade. Para isso, usa-se o chamado método do conjunto gerador para realizar buscas do tipo local em uma rede de pontos previamente gerada por buscas globais efetuadas com uma metaheurística populacional, no caso o método do enxame de partículas.Para se obter mais de um ponto estacionário, minimizam-se funções de mérito polarizadas, cujos pólos são os pontos previamente encontrados. A metodologia proposta foi testada na análise de quatorze misturas polares previamente consideradas na literatura. Os resultados mostraram que o método proposto é robusto e eficiente a ponto de encontrar, além do minimizador global, todos os pontos estacionários apontados previamente na literatura, sendo também capaz de detectar, em duas misturas ternárias estudadas, pontos estacionários não obtidos pelo chamado método de análise intervalar, uma técnica confiável e muito difundida na literatura. A análise do teste de estabilidade pela simples utilização do método do enxame de partículas associado à técnica de polarização mencionada acima, para a obtenção de mais de um ponto estacionário (sem a busca local feita pelo método do conjunto gerador em uma dada rede de pontos), constitui outra metodologia para a resolução do problema de interesse. Essa utilização é uma novidade secundária deste trabalho. Tal metodologia simplificada exibiu também uma grande robustez, sendo capaz de encontrar todos os pontos estacionários pesquisados. No entanto, quando comparada com a abordagem mais geral proposta aqui, observou-se que tal simplificação pode, em alguns casos onde a função de mérito apresenta uma geometria mais complexa, consumir um tempo de máquina relativamente grande, dessa forma é menos eficiente.

Uma dedução dos critérios de multicriticalidade e o cálculo de pontos tricríticos via otimização global

Uma dedução dos critérios de multicriticalidade para o cálculo de pontos críticos de qualquer ordem representa a formalização de ideias utilizadas para calcular pontos críticos e tricríticos e ainda amplia tais ideias. De posse desta dedução pode-se compreender os critérios de tricriticalidade e, com isso, através de uma abordagem via problema de otimização global pode-se fazer o cálculo de pontos tricríticos utilizando um método numérico adequado de otimização global. Para evitar um excesso de custo computacional com rotinas numéricas utilizou-se aproximações na forma de diferenças finitas dos termos que compõem a função objetivo. Para simular a relação P v - T optou-se pela equação de estado cúbica de Peng-Robinson e pela regra clássica de fluidos de van der Vaals, para modelagem do problema também se calculou os tensores de ordem 2, 3, 4 e 5 da função do teste de estabilidade. Os resultados obtidos foram comparados com dados experimentais e por resultados obtidos com outros autores que utilizaram métodos numéricos, equação de estado ou abordagem diferente das utilizadas neste trabalho.

Stability analysis and observer design for discrete-time SEIR epidemic models

This paper applies Micken's discretization method to obtain a discrete-time SEIR epidemic model. The positivity of the model along with the existence and stability of equilibrium points is discussed for the discrete-time case. Afterwards, the design of a state observer for this discrete-time SEIR epidemic model is tackled. The analysis of the model along with the observer design is faced in an implicit way instead of obtaining first an explicit formulation of the system which is the novelty of the presented approach. Moreover, some sufficient conditions to ensure the asymptotic stability of the observer are provided in terms of a matrix inequality that can be cast in the form of a LMI. The feasibility of the matrix inequality is proved, while some simulation examples show the operation and usefulness of the observer.

On the Stability and Equilibrium Points of Multistaged SI (n) R Epidemic Models

This paper relies on the concept of next generation matrix defined ad hoc for a new proposed extended SEIR model referred to as SI(n)R-model to study its stability. The model includes n successive stages of infectious subpopulations, each one acting at the exposed subpopulation of the next infectious stage in a cascade global disposal where each infectious population acts as the exposed subpopulation of the next infectious stage. The model also has internal delays which characterize the time intervals of the coupling of the susceptible dynamics with the infectious populations of the various cascade infectious stages. Since the susceptible subpopulation is common, and then unique, to all the infectious stages, its coupled dynamic action on each of those stages is modeled with an increasing delay as the infectious stage index increases from 1 to n. The physical interpretation of the model is that the dynamics of the disease exhibits different stages in which the infectivity and the mortality rates vary as the individual numbers go through the process of recovery, each stage with a characteristic average time.

Formulação teórica dos fundamentos da otimização global topográfica com análise de desempenho e aplicações à estabilidade de fases de misturas termodinâmicas

Métodos de otimização que utilizam condições de otimalidade de primeira e/ou segunda ordem são conhecidos por serem eficientes. Comumente, esses métodos iterativos são desenvolvidos e analisados à luz da análise matemática do espaço euclidiano n-dimensional, cuja natureza é de caráter local. Consequentemente, esses métodos levam a algoritmos iterativos que executam apenas as buscas locais. Assim, a aplicação de tais algoritmos para o cálculo de minimizadores globais de uma função não linear,especialmente não-convexas e multimodais, depende fortemente da localização dos pontos de partida. O método de Otimização Global Topográfico é um algoritmo de agrupamento, que utiliza uma abordagem baseada em conceitos elementares da teoria dos grafos, a fim de gerar bons pontos de partida para os métodos de busca local, a partir de pontos distribuídos de modo uniforme no interior da região viável. Este trabalho tem dois objetivos. O primeiro é realizar uma nova abordagem sobre método de Otimização Global Topográfica, onde, pela primeira vez, seus fundamentos são formalmente descritos e suas propriedades básicas são matematicamente comprovadas. Neste contexto, propõe-se uma fórmula semi-empírica para calcular o parâmetro chave deste algoritmo de agrupamento, e, usando um método robusto e eficiente de direções viáveis por pontos-interiores, estendemos o uso do método de Otimização Global Topográfica a problemas com restrições de desigualdade. O segundo objetivo é a aplicação deste método para a análise de estabilidade de fase em misturas termodinâmicas,o qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. A solução deste problema de otimização global é necessária para o cálculo do equilíbrio de fases, que é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Além disso, afim de ter uma avaliação inicial do potencial dessa técnica, primeiro vamos resolver 70 problemas testes, e então comparar o desempenho do método proposto aqui com o solver MIDACO, um poderoso software recentemente introduzido no campo da otimização global.

Global nonlinear dynamics of thin aerofoil wakes

In the present investigation of thin aerofoil wakes we compare the global nonlinear dynamics, obtained by direct numerical simulations, to the associated local instability features, derived from linear stability analyses. A given configuration depends on two control parameters: the Reynolds number Re and the adverse pressure gradient m (with m < 0) prevailing at the aerofoil trailing edge. Global instability is found to occur for large enough Re and |m|; the naturally selected frequency is determined by the local absolute frequency prevailing at the trailing edge. © 2010 Springer Science+Business Media B.V.

Surface tension-induced global instability of planar jets and wakes

The effect of surface tension on global stability of co-flow jets and wakes at a moderate Reynolds number is studied. The linear temporal two-dimensional global modes are computed without approximations. All but one of the flow cases under study are globally stable without surface tension. It is found that surface tension can cause the flow to be globally unstable if the inlet shear (or equivalently, the inlet velocity ratio) is strong enough. For even stronger surface tension, the flow is re-stabilized. As long as there is no change of the most unstable mode, increasing surface tension decreases the oscillation frequency. Short waves appear in the high-shear region close to the nozzle, and their wavelength increases with increasing surface tension. The critical shear (the weakest inlet shear at which a global instability is found) gives rise to antisymmetric disturbances for the wakes and symmetric disturbances for the jets. However, at stronger shear, the opposite symmetry can be the most unstable one, in particular for wakes at high surface tension. The results show strong effects of surface tension that should be possible to reproduce experimentally as well as numerically.

Feedback mechanisms for global oscillations in Lure systems

The paper presents two mechanisms for global oscillations in feedback systems, based on bifurcations in absolutely stable systems. The external characterization of the oscillators provides the basis for a (energy-based) dissipativity theory for oscillators, thereby opening new possibilities for rigorous stability analysis of high-dimensional systems and interconnected oscillators. © 2004 Elsevier B.V. All rights reserved.

Global modes with multiple saddle points

Significant progress has been made towards understanding the global stability of slowly-developing shear flows. The WKBJ theory developed by Patrick Huerre and his co-authors has proved absolutely central, with the result that both the linear and the nonlinear stability of a wide range of flows can now be understood in terms of their local absolute/convective instability properties. In many situations, the local absolute frequency possesses a single dominant saddle point in complex X-space (where X is the slow streamwise coordinate of the base flow), which then acts as a single wavemaker driving the entire global linear dynamics. In this paper we consider the more complicated case in which multiple saddles may act as the wavemaker for different values of some control parameter. We derive a frequency selection criterion in the general case, which is then validated against numerical results for the linearized third-order Ginzburg-Landau equation (which possesses two saddle points). We believe that this theory may be relevant to a number of flows, including the boundary layer on a rotating disk and the eccentric Taylor-Couette-Poiseuille flow. © 2014 Elsevier Masson SAS. All rights reserved.

Hybrid numerical methods for exponential models of growth

A series of novel numerical methods for the exponential models of growth are proposed. Based on these methods, hybrid predictor-corrector methods are constructed. The hybrid numerical methods can increase the accuracy and the computing speed obviously, as well as enlarge the stability domain greatly. (c) 2005 Published by Elsevier Inc.

川西北高寒草甸土壤呼吸研究：季节动态和放牧效应

除植被冠层的光合作用之外，土壤的呼吸作用是陆地生态系统碳收支中最大的通量。土壤呼吸即使发生较小的变化也能显著地减缓或加剧大气中CO2浓度的增加，从而明显影响到全球气候变化。土壤呼吸速率变化与否以及变化的方向可以反映生态系统对环境变化的敏感程度和响应模式。尽管如此，土壤呼吸仍是一个为人们了解不多的生态系统过程。 草地生态系统是陆地生态系统的一个重要组成部分。针对草地土壤呼吸进行野外实验研究和相应方法论的探讨将对区域乃至全球碳源汇性质的准确估算具有重要的科学意义。然而，近几年来关于草地土壤呼吸的主要研究工作都集中在温带草原和部分热带草原，而针对高寒草甸生态系统土壤呼吸的研究报道还很少。 2008年4月至2009年4月期间，我分别在2008年6、8、10、12月和2009年2月和4月分6次对川西北的典型高寒草甸群落的土壤呼吸进行观测，分析了不同类型高寒草甸群落土壤呼吸的季节变化特征以及环境因子和放牧模式对其影响。主要研究结果如下： 1）该地区高寒草甸生态系统在生长季（6月～8月）土壤呼吸速率较大（6.07～9.30μmolCO2¡m-2¡s-1 ） ， 在非生长季（ 12 月～ 2 月） 较小（ 0.16 ～0.49μmolCO2¡m-2¡s-1 ） 。土壤CO2 年累积最大释放量为3963 ～ 5730gCO2¡m-2¡yr-1，其中，生长季土壤CO2的释放量占年总释放量的85%～90%。非生长季占10%～15%。非生长季所占比例略小于冬季积雪覆盖地区的冬季土壤呼吸占年土壤呼吸量的比例（14%～30%）。温度，尤其地温，是影响该地区高寒草甸生态系统土壤呼吸速率的最主要环境因子。土壤呼吸速率与地上生物量和土壤水分之间没有显著相关性，但是土壤含水量过大会导致土壤呼吸速率下降。 2）在观测期内，草丘区的土壤呼吸显著高于对照区的土壤呼吸，其最大土壤呼吸速率为16.77μmolCO2¡m-2¡s-1，土壤CO2 年累积最大释放量为8145gCO2¡m-2¡yr-1，是对照区的近2 倍。由于草丘在高寒草甸中占有较大的面积比例（近30%），因此，它将对高寒草甸生态系统的碳循环起着重要的作用。 3）放牧模式不仅可以影响高寒草甸群落的土壤CO2 排放，而且还可以改变土壤呼吸的温度敏感性（Q10）。本研究表明，在生长季有长期放牧活动干扰时将会增加土壤向大气中释放二氧化碳的速度，促使土壤碳库中碳的流失。禁牧样地的土壤呼吸速率在刚禁牧时先迅速增大，随着禁牧时间的延长土壤呼吸速率将会下降。此外，与其它放牧模式相比，冬季放牧将高寒草甸群落土壤呼吸速率在生长季达到最大值的时间明显向后推迟。不同放牧模式下高寒草甸群落土壤呼吸的Q10 值大小顺序为：禁牧一年群落>冬季放牧群落>禁牧三年群落>夏季放牧群落>自由放牧群落。 4）基于呼吸室技术的观测方法中，测量前的剪草处理可以明显改变该地区高寒草甸群落的土壤温度和土壤呼吸速率。在生长季，剪草处理将使土壤呼吸速率的瞬时响应增加90%左右。由于剪草处理明显增加了剪草样方白天的土壤温度，而土壤温度与土壤呼吸之间存在着极显著的指数相关关系，因而剪草处理导致土壤呼吸速率迅速增加。因此，在高寒地区基于呼吸室技术观测的土壤呼吸应当进行校正。 综上所述，川西北高寒草甸生态系统土壤呼吸速率在生长季较高，而在非生长季较低。土壤温度是影响该地区土壤呼吸的最主要环境因子。在实验观测期，草丘区土壤呼吸速率显著高于对照区的，是对照区土壤呼吸速率的近2倍。由于测量前的剪草处理可以明显改变待测点的土壤呼吸速率，因此，应对在高寒地区基于呼吸室技术观测的土壤呼吸进行校正。 Soil respiration is the second largest component (less than plant phtotosynthesis) of carbon dioxide flux between terrestrial ecosystems and the atmosphere. A minor change in soil respiration rate can significantly slow down or accelerate the increase of atmospheric CO2 concentration that is closely related to global climatic change. In turn, the change in the flux direction and rate of soil respiration may indicate the elasticity and stability of ecosystems to global changes and human disturbance. However, soil respiration is still an ecosystem process that has been poorly understood. Grassland ecosystem is an important component of the terrestrial ecosystem. Accurately estimating the CO2 flux from soil to atmosphere in situ is the key to evaluating the carbon resource and sink regionally or globally. Despite of extensive studies on the temperate and tropic grasslands, the soil respiration of alpine meadows has not substantially been measured. In the current study, soil respiration was measured for an annual cycle from April, 2008 to April, 2009 for the alpine meadow in northwestern Sichuan Province of China to determine the seasonal variation of soil respiration for the typical plant communities. The results are shown as follows: 1) Large seasonal variation of soil respiration was observed in the alpine meadow. The rate of soil respiration was the greatest (6.07～9.30μmolCO2¡m-2¡s-1) in June and the smallest (0.16 ～ 0.49μmolCO2¡m-2¡s-1) occurred from December to February in the non-growing season. The total emission of soil CO2 was estimated as 3963～5730 gCO2¡m-2¡yr-1, 85%~90% of which was released during the growing season, and 10%~15% during the non-growing season which was slightly less than the ratio of winter and annual CO2 flux from soil. Temperature, particularly the soil temperature, was the major environmental factor regulating the soil respiration. Significant and positive relationships were not found between soil respiration and soil moisture and between soil respiration and plant above-ground biomass, but excessive soil water content would decrease in the rate of soil respiration. 2) The rate of soil respiration in grass hummock communities was up to 16.77μmolCO2¡m-2¡s-1, which was about twice as great as in the controls (communities located in low and even sites). Considering the large proportion (about 30% on average) of hummock area in the meadow, it can be concluded that the hummocks played an important role in the carbon cycling of the study ecosystem. 3) Grazing patterns affected the flux of CO2 emission and the temperature sensitivity of soil respiration (Q10) in the alpine meadow. Grazing during growing season increased the rate of soil respiration. The rate of soil respiration increased significantly immediately after the alpine meadow being fenced, but thereafter decreased. In addition, grazing in winter delayed the peak respiration rate relative to the non-grazing mode. The Q10 value was the largest in the non-grazed area for one year, and next came the area with grazing in winter, followed by the non-grazed area for three years, the area with grazing in summer, and the non-limited grazed area. 4) In the chamber-based techniques, clipping manipulation before each measurement increased the transient rate of soil respiration by about 90% in the summer of the alpine meadow. As increase in soil temperature at daytime in the clipped plots by clipping and the exponential relationship between soil respiration and temperature, clipping manipulation led to increase in the rate of soil respiration. This suggested that a correction should be done for the techniques if employed in alpine and cold regions. In summary, the rate of soil respiration in the alpine meadow was the greatest in June and the smallest occurred from ecember to February in the non-growing season. Soil temperature was the major environmental factor regulating the soil respiration. The rate of soil respiration in grass hummock communities was up to 16.77μmolCO2¡m-2¡s-1, which was about twice as great as in the controls. A correction should be done for the techniques if employed in alpine and cold regions, because of the effect of clipping manipulation on soil temperature and respiration.