Preserving a system’s viability in the presence of diversity erosion is critical if the goal is to sustainably support biodiversity. Reduction in population heterogeneity, whether inter- or intraspecies, may increase population fragility, either decreasing its ability to adapt effectively to environmental changes or facilitating the survival and success of ordinarily rare phenotypes. The latter may result in over-representation of individuals who may participate in resource utilization patterns that can lead to over-exploitation, exhaustion, and, ultimately, collapse of both the resource and the population that depends on it. Here, we aim to identify regimes that can signal whether a consumer–resource system is capable of supporting viable degrees of heterogeneity. The framework used here is an expansion of a previously introduced consumer–resource type system of a population of individuals classified by their resource consumption. Application of the Reduction Theorem to the system enables us to evaluate the health of the system through tracking both the mean value of the parameter of resource (over)consumption, and the population variance, as both change over time. The article concludes with a discussion that highlights applicability of the proposed system to investigation of systems that are affected by particularly devastating overly adapted populations, namely cancerous cells. Potential intervention approaches for system management are discussed in the context of cancer therapies.

In the field of infectious disease epidemiology, the assessment of model robustness outcomes plays a significant role in the identification, reformulation, and evaluation of preparedness strategies aimed at limiting the impact of catastrophic events (pandemics or the deliberate release of biological agents) or used in the management of disease prevention strategies, or employed in the identification and evaluation of control or mitigation measures. The research work in this dissertation focuses on: The comparison and assessment of the role of exponentially distributed waiting times versus the use of generalized non-exponential parametric distributed waiting times of infectious periods on the quantitative and qualitative outcomes generated by Susceptible-Infectious-Removed (SIR) models. Specifically, Gamma distributed infectious periods are considered in the three research projects developed following the applications found in (Bailey 1964, Anderson 1980, Wearing 2005, Feng 2007, Feng 2007, Yan 2008, lloyd 2009, Vergu 2010). i) The first project focuses on the influence of input model parameters, such as the transmission rate, mean and variance of Gamma distributed infectious periods, on disease prevalence, the peak epidemic size and its timing, final epidemic size, epidemic duration and basic reproduction number. Global uncertainty and sensitivity analyses are carried out using a deterministic Susceptible-Infectious-Recovered (SIR) model. The quantitative effect and qualitative relation between input model parameters and outcome variables are established using Latin Hypercube Sampling (LHS) and Partial rank correlation coefficient (PRCC) and Spearman rank correlation coefficient (RCC) sensitivity indices. We learnt that: For relatively low (R0 close to one) to high (mean of R0 equals 15) transmissibility, the variance of the Gamma distribution for the infectious period, input parameter of the deterministic age-of-infection SIR model, is key (statistically significant) on the predictability of the epidemiological variables such as the epidemic duration and the peak size and timing of the prevalence of infectious individuals and therefore, for the predictability these variables, it is preferable to utilize a nonlinear system of Volterra integral equations, rather than a nonlinear system of ordinary differential equations. The predictability of epidemiological variables such as the final epidemic size and the basic reproduction number are unaffected by (or independent of) the variance of the Gamma distribution for the infectious period and therefore for the choice on which type of nonlinear system for the description of the SIR model (VIE's or ODE's) is irrelevant. Although, for practical proposes, with the aim of lowering the complexity and number operations in the numerical methods, a nonlinear system of ordinary differential equations is preferred. The main contribution lies in the development of a model based decision-tool that helps determine when SIR models given in terms of Volterra integral equations are equivalent or better suited than SIR models that only consider exponentially distributed infectious periods. ii) The second project addresses the question of whether or not there is sufficient evidence to conclude that two empirical distributions for a single epidemiological outcome, one generated using a stochastic SIR model under exponentially distributed infectious periods and the other under the non-exponentially distributed infectious period, are statistically dissimilar. The stochastic formulations are modeled via a continuous time Markov chain model. The statistical hypothesis test is conducted using the non-parametric Kolmogorov-Smirnov test. We found evidence that shows that for low to moderate transmissibility, all empirical distribution pairs (generated from exponential and non-exponential distributions) for each of the epidemiological quantities considered are statistically dissimilar. The research in this project helps determine whether the weakening exponential distribution assumption must be considered in the estimation of probability of events defined from the empirical distribution of specific random variables. iii) The third project involves the assessment of the effect of exponentially distributed infectious periods on estimates of input parameter and the associated outcome variable predictions. Quantities unaffected by the use of exponentially distributed infectious period within low transmissibility scenarios include, the prevalence peak time, final epidemic size, epidemic duration and basic reproduction number and for high transmissibility scenarios only the prevalence peak time and final epidemic size. An application designed to determine from incidence data whether there is sufficient statistical evidence to conclude that the infectious period distribution should not be modeled by an exponential distribution is developed. A method for estimating explicitly specified non-exponential parametric probability density functions for the infectious period from epidemiological data is developed. The methodologies presented in this dissertation may be applicable to models where waiting times are used to model transitions between stages, a process that is common in the study of life-history dynamics of many ecological systems.

Extraordinary medical advances have led to significant reductions in the burden of infectious diseases in humans. However, infectious diseases still account for more than 13 million annual deaths. This large burden is partly due to some pathogens having found suitable conditions to emerge and spread in denser and more connected host populations, and others having evolved to escape the pressures imposed by the rampant use of antimicrobials. It is then critical to improve our understanding of how diseases spread in these modern landscapes, characterized by new host population structures and socio-economic environments, as well as containment measures such as the deployment of drugs. Thus, the motivation of this dissertation is two-fold. First, we study, using both data-driven and modeling approaches, the the spread of infectious diseases in urban areas. As a case study, we use confirmed-cases data on sexually transmitted diseases (STDs) in the United States to assess the conduciveness of population size of urban areas and their socio-economic characteristics as predictors of STD incidence. We find that the scaling of STD incidence in cities is superlinear, and that the percent of African-Americans residing in cities largely determines these statistical patterns. Since disparities in access to health care are often exacerbated in urban areas, within this project we also develop two modeling frameworks to study the effect of health care disparities on epidemic outcomes. Discrepant results between the two approaches indicate that knowledge of the shape of the recovery period distribution, not just its mean and variance, is key for assessing the epidemiological impact of inequalities. The second project proposes to study, from a modeling perspective, the spread of drug resistance in human populations featuring vital dynamics, stochasticity and contact structure. We derive effective treatment regimes that minimize both the overall disease burden and the spread of resistance. Additionally, targeted treatment in structured host populations may lead to higher levels of drug resistance, and if drug-resistant strains are compensated, they can spread widely even when the wild-type strain is below its epidemic threshold.

Mathematical modeling of infectious diseases can help public health officials to make decisions related to the mitigation of epidemic outbreaks. However, over or under estimations of the morbidity of any infectious disease can be problematic. Therefore, public health officials can always make use of better models to study the potential implication of their decisions and strategies prior to their implementation. Previous work focuses on the mechanisms underlying the different epidemic waves observed in Mexico during the novel swine origin influenza H1N1 pandemic of 2009 and showed extensions of classical models in epidemiology by adding temporal variations in different parameters that are likely to change during the time course of an epidemic, such as, the influence of media, social distancing, school closures, and how vaccination policies may affect different aspects of the dynamics of an epidemic. This current work further examines the influence of different factors considering the randomness of events by adding stochastic processes to meta-population models. I present three different approaches to compare different stochastic methods by considering discrete and continuous time. For the continuous time stochastic modeling approach I consider the continuous-time Markov chain process using forward Kolmogorov equations, for the discrete time stochastic modeling I consider stochastic differential equations using Wiener's increment and Poisson point increments, and also I consider the discrete-time Markov chain process. These first two stochastic modeling approaches will be presented in a one city and two city epidemic models using, as a base, our deterministic model. The last one will be discussed briefly on a one city SIS and SIR-type model.

Olfaction is an important sensory modality for behavior since odors inform animals of the presence of food, potential mates, and predators. The fruit fly, Drosophila melanogaster, is a favorable model organism for the investigation of the biophysical mechanisms that contribute to olfaction because its olfactory system is anatomically similar to but simpler than that of vertebrates. In the Drosophila olfactory system, sensory transduction takes place in olfactory receptor neurons housed in the antennae and maxillary palps on the front of the head. The first stage of olfactory processing resides in the antennal lobe, where the structural unit is the glomerulus. There are at least three classes of neurons in the antennal lobe - excitatory projection neurons, excitatory local neurons, and inhibitory local neurons. The arborizations of the local neurons are confined to the antennal lobe, and output from the antennal lobe is carried by projection neurons to higher regions of the brain. Different views exist of how circuits of the Drosophila antennal lobe translate input from the olfactory receptor neurons into projection neuron output. We construct a conductance based neuronal network model of the Drosophila antennal lobe with the aim of understanding possible mechanisms within the antennal lobe that account for the variety of projection neuron activity observed in experimental data. We explore possible outputs obtained from olfactory receptor neuron input that mimic experimental recordings under different connectivity paradigms. First, we develop realistic minimal cell models for the excitatory local neurons, inhibitory local neurons, and projections neurons based on experimental data for Drosophila channel kinetics, and explore the firing characteristics and mathematical structure of these models. We then investigate possible interglomerular and intraglomerular connectivity patterns in the Drosophila antennal lobe, where olfactory receptor neuron input to the antennal lobe is modeled with Poisson spike trains, and synaptic connections within the antennal lobe are mediated by chemical synapses and gap junctions as described in the Drosophila antennal lobe literature. Our simulation results show that inhibitory local neurons spread inhibition among all glomeruli, where projection neuron responses are decreased relatively uniformly for connections of synaptic strengths that are homogeneous. Also, in the case of homogeneous excitatory synaptic connections, the excitatory local neuron network facilitates odor detection in the presence of weak stimuli. Excitatory local neurons can spread excitation from projection neurons that receive more input from olfactory receptor neurons to projection neurons that receive less input from olfactory receptor neurons. For the parameter values for the network models associated with these results, eLNs decrease the ability of the network to discriminate among single odors.

Solution methods for certain linear and nonlinear evolution equations are presented in this dissertation. Emphasis is placed mainly on the analytical treatment of nonautonomous differential equations, which are challenging to solve despite the existent numerical and symbolic computational software programs available. Ideas from the transformation theory are adopted allowing one to solve the problems under consideration from a non-traditional perspective. First, the Cauchy initial value problem is considered for a class of nonautonomous and inhomogeneous linear diffusion-type equation on the entire real line. Explicit transformations are used to reduce the equations under study to their corresponding standard forms emphasizing on natural relations with certain Riccati(and/or Ermakov)-type systems. These relations give solvability results for the Cauchy problem of the parabolic equation considered. The superposition principle allows to solve formally this problem from an unconventional point of view. An eigenfunction expansion approach is also considered for this general evolution equation. Examples considered to corroborate the efficacy of the proposed solution methods include the Fokker-Planck equation, the Black-Scholes model and the one-factor Gaussian Hull-White model. The results obtained in the first part are used to solve the Cauchy initial value problem for certain inhomogeneous Burgers-type equation. The connection between linear (the Diffusion-type) and nonlinear (Burgers-type) parabolic equations is stress in order to establish a strong commutative relation. Traveling wave solutions of a nonautonomous Burgers equation are also investigated. Finally, it is constructed explicitly the minimum-uncertainty squeezed states for quantum harmonic oscillators. They are derived by the action of corresponding maximal kinematical invariance group on the standard ground state solution. It is shown that the product of the variances attains the required minimum value only at the instances that one variance is a minimum and the other is a maximum, when the squeezing of one of the variances occurs. Such explicit construction is possible due to the relation between the diffusion-type equation studied in the first part and the time-dependent Schrodinger equation. A modication of the radiation field operators for squeezed photons in a perfect cavity is also suggested with the help of a nonstandard solution of Heisenberg's equation of motion.

A sequence of models is developed to describe urban population growth in the context of the embedded physical, social and economic environments and an urban disease are developed. This set of models is focused on urban growth and the relationship between the desire to move and the utility derived from city life. This utility is measured in terms of the economic opportunities in the city, the level of human constructed amenity, and the level of amenity caused by the natural environment. The set of urban disease models is focused on examining prospects of eliminating a disease for which a vaccine does not exist. It is inspired by an outbreak of the vector-borne disease dengue fever in Peru, during 2000-2001.

In complex consumer-resource type systems, where diverse individuals are interconnected and interdependent, one can often anticipate what has become known as the tragedy of the commons, i.e., a situation, when overly efficient consumers exhaust the common resource, causing collapse of the entire population. In this dissertation I use mathematical modeling to explore different variations on the consumer-resource type systems, identifying some possible transitional regimes that can precede the tragedy of the commons. I then reformulate it as a game of a multi-player prisoner's dilemma and study two possible approaches for preventing it, namely direct modification of players' payoffs through punishment/reward and modification of the environment in which the interactions occur. I also investigate the questions of whether the strategy of resource allocation for reproduction or competition would yield higher fitness in an evolving consumer-resource type system and demonstrate that the direction in which the system will evolve will depend not only on the state of the environment but largely on the initial composition of the population. I then apply the developed framework to modeling cancer as an evolving ecological system and draw conclusions about some alternative approaches to cancer treatment.

Dopamine (DA) is a neurotransmitter involved in attention, goal oriented behavior, movement, reward learning, and short term and working memory. For the past four decades, mathematical and computational modeling approaches have been useful in DA research, and although every modeling approach has limitations, a model is an efficient way to generate and explore hypotheses. This work develops a model of DA dynamics in a representative, single DA neuron by integrating previous experimental, theoretical and computational research. The model consists of three compartments: the cytosol, the vesicles, and the extracellular space and forms the basis of a new mathematical paradigm for examining the dynamics of DA synthesis, storage, release and reuptake. The model can be driven by action potentials generated by any model of excitable membrane potential or even from experimentally induced depolarization voltage recordings. Here the model is forced by a previously published model of the excitable membrane of a mesencephalic DA neuron in order to study the biochemical processes involved in extracellular DA production. After demonstrating that the model exhibits realistic dynamics resembling those observed experimentally, the model is used to examine the functional changes in presynaptic mechanisms due to application of cocaine. Sensitivity analysis and numerical studies that focus on various possible mechanisms for the inhibition of DAT by cocaine provide insight for the complex interactions involved in DA dynamics. In particular, comparing numerical results for a mixed inhibition mechanism to those for competitive, non-competitive and uncompetitive inhibition mechanisms reveals many behavioral similarities for these different types of inhibition that depend on inhibition parameters and levels of cocaine. Placing experimental results within this context of mixed inhibition provides a possible explanation for the conflicting views of uptake inhibition mechanisms found in experimental neuroscience literature.

There have been many studies on the dynamics of infectious diseases considering the age structure of the population. This study analyzes the dynamics when the population is stratified by size. This kind of models are useful in the spread of a disease in fisheries where size matters, for microorganism populations or even human diseases that are driven by weight. A simple size structured SIR model is introduced for which a threshold condition, R0, equilibria and stability are established in special cases. Hethcote's approach is used to derive, from first principles, a parallel ODE size-structure system involving n-size classes.The specific case of n = 2 is partially analyzed. Constant effort harvesting is added to this model with the purpose of exploring the role of controls and harvesting. Different harvesting policies are proposed and analyzed through simulations.