Modeling GHG emissions, N and C dynamics in Spanish agricultural soils.

To date, only few initiatives have been carried out in Spain in order to use mathematical models (e.g. DNDC, DayCent, FASSET y SIMSNIC) to estimate nitrogen (N) and carbon (C) dynamics as well as greenhouse gases (GHG) in Spanish agrosystems. Modeling at this level may allow to gain insight on both the complex relationships between biological and physicochemical processes, controlling the processes leading to GHG production and consumption in soils (e.g. nitrification, denitrification, decomposing, etc.), and the interactions between C and N cycles within the different components of the continuum plant-soil-environment. Additionally, these models can simulate the processes behind production, consumition and transport of GHG (e.g. nitrous oxide, N2O, and carbon dioxide, CO2) in the short and medium term and at different scales. Other sources of potential pollution from soils can be identified and quantified using these process-based models (e.g. NO3 y NH3).

Modeling the dynamics of human hair cycles by a follicular automaton

The hair follicle cycle successively goes through the anagen, catagen, telogen, and latency phases, which correspond, respectively, to hair growth, arrest, shedding, and absence before a new anagen phase is initiated. Experimental observations collected over a period of 14 years in a group of 10 male volunteers, alopecic and nonalopecic, allowed us to determine the characteristics of scalp hair follicle cycles. On the basis of these observations, we propose a follicular automaton model to simulate the dynamics of human hair cycles. The automaton model is defined by a set of rules that govern the stochastic transitions of each follicle between the successive states anagen, telogen, and latency, and the subsequent return to anagen. The transitions occur independently for each follicle, after time intervals given stochastically by a distribution characterized by a mean and a variance. The follicular automaton model accounts both for the dynamical transitions observed in a single follicle and for the behavior of an ensemble of independently cycling follicles. Thus, the model successfully reproduces the evolution of the fractions of follicle populations in each of the three phases, which fluctuate around steady-state or slowly drifting values. We apply the follicular automaton model to the study of spatial patterns of follicular growth that result from a spatially heterogeneous distribution of parameters such as the mean duration of anagen phase. When considering that follicles die or miniaturize after going through a critical number of successive cycles, the model can reproduce the evolution to hair patterns similar to well known types of diffuse or androgenetic alopecia.

Modeling contact formation between atomic-sized gold tips via molecular dynamics

The formation and rupture of atomic-sized contacts is modelled by means of molecular dynamics simulations. Such nano-contacts are realized in scanning tunnelling microscope and mechanically controlled break junction experiments. These instruments routinely measure the conductance across the nano-sized electrodes as they are brought into contact and separated, permitting conductance traces to be recorded that are plots of conductance versus the distance between the electrodes. One interesting feature of the conductance traces is that for some metals and geometric configurations a jump in the value of the conductance is observed right before contact between the electrodes, a phenomenon known as jump-to-contact. This paper considers, from a computational point of view, the dynamics of contact between two gold nano-electrodes. Repeated indentation of the two surfaces on each other is performed in two crystallographic orientations of face-centred cubic gold, namely (001) and (111). Ultimately, the intention is to identify the structures at the atomic level at the moment of first contact between the surfaces, since the value of the conductance is related to the minimum cross-section in the contact region. Conductance values obtained in this way are determined using first principles electronic transport calculations, with atomic configurations taken from the molecular dynamics simulations serving as input structures.

Modeling Depression Progression Dynamics from Electronic Health Record

Thesis (Master's)--University of Washington, 2016-06

Functionality, Robustness and Control of Nonlinear Network Dynamics: Modeling and Understanding the C. elegans Connectome

Thesis (Ph.D.)--University of Washington, 2016-06

Effector dynamics of rhythmic wrist activity and its implications for (modeling) bimanual coordination

To examine the role of the effector dynamics of the wrist in the production of rhythmic motor activity, we estimated the phase shifts between the EMG and the task-related output for a rhythmic isometric torque production task and an oscillatory movement, and found a substantial difference (45-52degrees) between the two. For both tasks, the relation between EMG and task-related output (torque or displacement) was adequately reproduced with a physiologically motivated musculoskeletal model. The model simulations demonstrated the importance of the contribution of passive structures to the overall dynamics and provided an account for the observed phase shifts in the dynamic task. Additional simulations of the musculoskeletal model with added load suggested that particular changes in the phase relation between EMG and movement may follow largely from the intrinsic muscle dynamics, rather than being the result of adaptations in the neural control of joint stiffness. The implications of these results are discussed in relation to (models of) interlimb coordination in rhythmic tasks. (C) 2004 Elsevier B.V. All rights reserved.

Applied Aspects of Mathematical Modeling and Optimization of Dynamics of Charged Beams

The complex of questions connected with the analysis, estimation and structural-parametrical optimization of dynamic system is considered in this article. Connection of such problems with tasks of control by beams of trajectories is emphasized. The special attention is concentrated on the review and analysis of spent scientific researches, the attention is stressed to their constructability and applied directedness. Efficiency of the developed algorithmic and software is demonstrated on the tasks of modeling and optimization of output beam characteristics in linear resonance accelerators.

Modeling of CW Yb-doped fiber lasers with highly nonlinear cavity dynamics

We develop a theoretical framework for modeling of continuous wave Yb-doped fiber lasers with highly nonlinear cavity dynamics. The developed approach has shown good agreement between theoretical predictions and experimental results for particular scheme of Yb-doped laser with large spectral broadening during single round trip. The model is capable to accurately describe main features of the experimentally measured laser outputs such as power efficiency slope, power leakage through fibre Bragg gratings, spectral broadening and spectral shape of generated radiation. © 2011 Optical Society of America.

Mode-locking intracavity dynamics:modeling nonlinear power and energy fluctuations

Modern high-power, pulsed lasers are driven by strong intracavity fluctuations. Critical in driving the intracavity dynamics is the nontrivial phase profiles generated and their periodic modification from either nonlinear mode-coupling, spectral filtering or dispersion management. Understanding the theoretical origins of the intracavity fluctuations helps guide the design, optimization and construction of efficient, high-power and high-energy pulsed laser cavities. Three specific mode-locking component are presented for enhancing laser energy: waveguide arrays, spectral filtering and dispersion management. Each component drives a strong intracavity dynamics that is captured through various modeling and analytic techniques.

Simulation modeling of phytoplankton dynamics in a large eutrophic river, Hungary – Danubian Phytoplankton Growth Model (DPGM)

Ecological models have often been used in order to answer questions that are in the limelight of recent researches such as the possible effects of climate change. The methodology of tactical models is a very useful tool comparison to those complex models requiring relatively large set of input parameters. In this study, a theoretical strategic model (TEGM ) was adapted to the field data on the basis of a 24-year long monitoring database of phytoplankton in the Danube River at the station of G¨od, Hungary (at 1669 river kilometer – hereafter referred to as “rkm”). The Danubian Phytoplankton Growth Model (DPGM) is able to describe the seasonal dynamics of phytoplankton biomass (mg L−1) based on daily temperature, but takes the availability of light into consideration as well. In order to improve fitting, the 24-year long database was split in two parts in accordance with environmental sustainability. The period of 1979–1990 has a higher level of nutrient excess compared with that of the 1991–2002. The authors assume that, in the above-mentioned periods, phytoplankton responded to temperature in two different ways, thus two submodels were developed, DPGM-sA and DPGMsB. Observed and simulated data correlated quite well. Findings suggest that linear temperature rise brings drastic change to phytoplankton only in case of high nutrient load and it is mostly realized through the increase of yearly total biomass.

Modeling seasonal dynamics of small fish cohorts in fluctuating freshwater marsh landscapes

Small-bodied fishes constitute an important assemblage in many wetlands. In wetlands that dry periodically except for small permanent waterbodies, these fishes are quick to respond to change and can undergo large fluctuations in numbers and biomasses. An important aspect of landscapes that are mixtures of marsh and permanent waterbodies is that high rates of biomass production occur in the marshes during flooding phases, while the permanent waterbodies serve as refuges for many biotic components during the dry phases. The temporal and spatial dynamics of the small fishes are ecologically important, as these fishes provide a crucial food base for higher trophic levels, such as wading birds. We develop a simple model that is analytically tractable, describing the main processes of the spatio-temporal dynamics of a population of small-bodied fish in a seasonal wetland environment, consisting of marsh and permanent waterbodies. The population expands into newly flooded areas during the wet season and contracts during declining water levels in the dry season. If the marsh dries completely during these times (a drydown), the fish need refuge in permanent waterbodies. At least three new and general conclusions arise from the model: (1) there is an optimal rate at which fish should expand into a newly flooding area to maximize population production; (2) there is also a fluctuation amplitude of water level that maximizes fish production, and (3) there is an upper limit on the number of fish that can reach a permanent waterbody during a drydown, no matter how large the marsh surface area is that drains into the waterbody. Because water levels can be manipulated in many wetlands, it is useful to have an understanding of the role of these fluctuations.

SYSTEM DYNAMICS MODELING AS A QUANTITATIVE-QUALITATIVE FRAMEWORK FOR SUSTAINABLE WATER RESOURCES MANAGEMENT: INSIGHTS FOR WATER QUALITY POLICY IN THE GREAT LAKES REGION

Early water resources modeling efforts were aimed mostly at representing hydrologic processes, but the need for interdisciplinary studies has led to increasing complexity and integration of environmental, social, and economic functions. The gradual shift from merely employing engineering-based simulation models to applying more holistic frameworks is an indicator of promising changes in the traditional paradigm for the application of water resources models, supporting more sustainable management decisions. This dissertation contributes to application of a quantitative-qualitative framework for sustainable water resources management using system dynamics simulation, as well as environmental systems analysis techniques to provide insights for water quality management in the Great Lakes basin. The traditional linear thinking paradigm lacks the mental and organizational framework for sustainable development trajectories, and may lead to quick-fix solutions that fail to address key drivers of water resources problems. To facilitate holistic analysis of water resources systems, systems thinking seeks to understand interactions among the subsystems. System dynamics provides a suitable framework for operationalizing systems thinking and its application to water resources problems by offering useful qualitative tools such as causal loop diagrams (CLD), stock-and-flow diagrams (SFD), and system archetypes. The approach provides a high-level quantitative-qualitative modeling framework for "big-picture" understanding of water resources systems, stakeholder participation, policy analysis, and strategic decision making. While quantitative modeling using extensive computer simulations and optimization is still very important and needed for policy screening, qualitative system dynamics models can improve understanding of general trends and the root causes of problems, and thus promote sustainable water resources decision making. Within the system dynamics framework, a growth and underinvestment (G&U) system archetype governing Lake Allegan's eutrophication problem was hypothesized to explain the system's problematic behavior and identify policy leverage points for mitigation. A system dynamics simulation model was developed to characterize the lake's recovery from its hypereutrophic state and assess a number of proposed total maximum daily load (TMDL) reduction policies, including phosphorus load reductions from point sources (PS) and non-point sources (NPS). It was shown that, for a TMDL plan to be effective, it should be considered a component of a continuous sustainability process, which considers the functionality of dynamic feedback relationships between socio-economic growth, land use change, and environmental conditions. Furthermore, a high-level simulation-optimization framework was developed to guide watershed scale BMP implementation in the Kalamazoo watershed. Agricultural BMPs should be given priority in the watershed in order to facilitate cost-efficient attainment of the Lake Allegan's TP concentration target. However, without adequate support policies, agricultural BMP implementation may adversely affect the agricultural producers. Results from a case study of the Maumee River basin show that coordinated BMP implementation across upstream and downstream watersheds can significantly improve cost efficiency of TP load abatement.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.