skip to main content


Title: Probability-based mechanisms in biological networks with parameter uncertainty
Mathematical models of biomolecular networks are commonly used to study mechanisms of cellular processes, but their usefulness is often questioned due to parameter uncertainty. Here, we employ Bayesian parameter inference and dynamic network analysis to study dominant reaction fluxes in models of extrinsic apoptosis. Although a simplified model yields thousands of parameter vectors with equally good fits to data, execution modes based on reaction fluxes clusters to three dominant execution modes. A larger model with increased parameter uncertainty shows that signal flow is constrained to eleven execution modes that use 53 out of 2067 possible signal subnetworks. Each execution mode exhibits different behaviors to in silico perturbations, due to different signal execution mechanisms. Machine learning identifies informative parameters to guide experimental validation. Our work introduces a probability-based paradigm of signaling mechanisms, highlights systems-level interactions that modulate signal flow, and provides a methodology to understand mechanistic model predictions with uncertain parameters.  more » « less
Award ID(s):
1942255
NSF-PAR ID:
10302778
Author(s) / Creator(s):
Date Published:
Journal Name:
bioRxiv
ISSN:
2692-8205
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change. 
    more » « less
  2. Abstract

    Secondary forest regrowth shapes community succession and biogeochemistry for decades, including in the Upper Great Lakes region. Vegetation models encapsulate our understanding of forest function, and whether models can reproduce multi‐decadal succession patterns is an indication of our ability to predict forest responses to future change. We test the ability of a vegetation model to simulate C cycling and community composition during 100 years of forest regrowth following stand‐replacing disturbance, asking (a) Which processes and parameters are most important to accurately model Upper Midwest forest succession? (b) What is the relative importance of model structure versus parameter values to these predictions? We ran ensembles of the Ecosystem Demography model v2.2 with different representations of processes important to competition for light. We compared the magnitude of structural and parameter uncertainty and assessed which sub‐model–parameter combinations best reproduced observed C fluxes and community composition. On average, our simulations underestimated observed net primary productivity (NPP) and leaf area index (LAI) after 100 years and predicted complete dominance by a single plant functional type (PFT). Out of 4,000 simulations, only nine fell within the observed range of both NPP and LAI, but these predicted unrealistically complete dominance by either early hardwood or pine PFTs. A different set of seven simulations were ecologically plausible but under‐predicted observed NPP and LAI. Parameter uncertainty was large; NPP and LAI ranged from ~0% to >200% of their mean value, and any PFT could become dominant. The two parameters that contributed most to uncertainty in predicted NPP were plant–soil water conductance and growth respiration, both unobservable empirical coefficients. We conclude that (a) parameter uncertainty is more important than structural uncertainty, at least for ED‐2.2 in Upper Midwest forests and (b) simulating both productivity and plant community composition accurately without physically unrealistic parameters remains challenging for demographic vegetation models.

     
    more » « less
  3. Abstract

    Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance‐based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain‐scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. These results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

     
    more » « less
  4. Turbulence is a major source of momentum, heat, moisture, and aerosol transport in the atmosphere. Hence, it is crucial to understand and accurately characterize turbulence mechanisms in atmospheric flows. Many complex factors in the atmosphere influence the turbulence structures including stratification and background shear. However, our understanding of the interacting effects of these factors on coherent turbulence structure evolutions is still limited. In this talk, we aim to bridge this knowledge gap by using mode decomposition techniques and a wide range of large-eddy simulation (LES) data. By developing a data-driven technique, we will characterize unique features of atmospheric boundary layer (ABL) turbulence under different forcing scenarios. We will present 3D LES wind speed snapshots of different ABL flows that will be used as dynamic mode decomposition (DMD) input data. Then, the obtained modes and eigenvalues will be employed to gain insights into coherent turbulence structures in ABLs. We will explain the physical meaning of dominant modes and how each mode relates to the physical cause of turbulence structures. The dominant modes, which are selected based on the mode amplitude, contain the most important spatial and temporal characteristics of the flow. We will evaluate the accuracy of the performance of this method by reconstructing the flow field with only a small number of modes, and then calculate the mean average error between the real flow and the reconstructed flow fields. We will present different data frequencies, wind speeds, and surface heat fluxes. This allows us to elucidate the modes and determine the conditions in which the mode decomposition provides more accurate results for the ABL flows. Our findings can be used to identify the major causes of turbulence in real atmospheric flows and could provide a deeper insight into the dynamics of turbulence in ABLs. Our results will also be useful for developing reduced-order models that can rapidly predict the turbulent ABL flow fields. 
    more » « less
  5. Abstract

    The mechanisms that control the export of freshwater from the East Greenland Current, in both liquid and solid form, are explored using an idealized numerical model and scaling theory. A regional, coupled ocean–sea ice model is applied to a series of calculations in which key parameters are varied and the scaling theory is used to interpret the model results. The offshore ice flux, occurring in late winter, is driven primarily by internal stresses and is most sensitive to the thickness of sea ice on the shelf coming out of Fram Strait and the strength of alongshore winds over the shelf. The offshore liquid freshwater flux is achieved by eddy fluxes in late summer while there is an onshore liquid freshwater flux in winter due to the ice–ocean stress, resulting in only weak annual mean flux. The scaling theory identifies the key nondimensional parameters that control the behavior and reproduces the general parameter dependence found in the numerical model. Climate models predict that winds will increase and ice export from the Arctic will decrease in the future, both of which will lead to a decrease in the offshore flux of sea ice, while the influence on liquid freshwater may increase or decrease, depending on the relative changes in the onshore Ekman transport and offshore eddy fluxes. Additional processes that have not been considered here, such as more complex topography and synoptic wind events, may also contribute to cross-shelf exchange.

    Significance Statement

    The purpose of this study is to provide a basic understanding of what controls the flux of sea ice and low-salinity water from the East Greenland shelf into the interior of the Greenland and Iceland Seas. This is a potentially important process since it has been shown that sufficient freshening of the surface waters in the interior of the Nordic seas can inhibit deep convection and the associated air–sea heat flux and water mass transformation. A combination of idealized computer models and basic theory indicates that the fluxes of liquid and solid freshwater are controlled by different mechanisms and occur at different times of the year. Accurate representation in climate models will require representation of small-scale processes such as mesoscale eddies and gradients of ice thickness across the shelf.

     
    more » « less