Exploratory ensemble designs for environmental models using k-extended Latin Hypercubes

Vector bionomics and vectorial capacity as emergent properties of mosquito behaviors and ecology

Mosquitoes are important vectors for pathogens that infect humans and other vertebrate animals. Some aspects of adult mosquito behavior and mosquito ecology play an important role in determining the capacity of vector populations to transmit pathogens. Here, we re-examine factors affecting the transmission of pathogens by mosquitoes using a new approach. Unlike most previous models, this framework considers the behavioral states and state transitions of adult mosquitoes through a sequence of activity bouts. We developed a new framework for individual-based simulation models called MBITES (Mosquito Bout-based and Individual-based Transmission Ecology Simulator). In MBITES, it is possible to build models that simulate the behavior and ecology of adult mosquitoes in exquisite detail on complex resource landscapes generated by spatial point processes. We also developed an ordinary differential equation model which is the Kolmogorov forward equations for models developed in MBITES under a specific set of simplifying assumptions. While mosquito infection and pathogen development are one possible part of a mosquito's state, that is not our main focus. Using extensive simulation using some models developed in MBITES, we show that vectorial capacity can be understood as an emergent property of simple behavioral algorithms interacting with complex resource landscapes, and that relative density or sparsity of resources and the need to search can have profound consequences for mosquito populations' capacity to transmit pathogens.

Rainfall From Resolved Rather Than Parameterized Processes Better Represents the Present-Day and Climate Change Response of Moderate Rates in the Community Atmosphere Model

Deficiencies in the parameterizations of convection used in global climate models often lead to a distorted representation of the simulated rainfall intensity distribution (i.e., too much rainfall from weak rain rates). While encouraging improvements in high percentile rainfall intensity have been found as the horizontal resolution of the Community Atmosphere Model is increased to ∼25 km, we demonstrate no corresponding improvement in the moderate rain rates that generate the majority of accumulated rainfall. Using a statistical framework designed to emphasize links between precipitation intensity and accumulated rainfall beyond just the frequency distribution, we show that CAM cannot realistically simulate moderate rain rates, and cannot capture their intensification with climate change, even as resolution is increased. However, by separating the parameterized convective and large-scale resolved contributions to total rainfall, we find that the intensity, geographic pattern, and climate change response of CAM's large-scale rain rates are more consistent with observations (TRMM 3B42), superparameterization, and theoretical expectations, despite issues with parameterized convection. Increasing CAM's horizontal resolution does improve the representation of total rainfall intensity, but not due to changes in the intensity of large-scale rain rates, which are surprisingly insensitive to horizontal resolution. Rather, improvements occur through an increase in the relative contribution of the large-scale component to the total amount of accumulated rainfall. Analysis of sensitivities to convective timescale and entrainment rate confirm the importance of these parameters in the possible development of scale-aware parameterizations, but also reveal unrecognized trade-offs from the entanglement of precipitation frequency and total amount.

A comparison of statistical emulation methodologies for multi-wave calibration of environmental models

Expensive computer codes, particularly those used for simulating environmental or geological processes, such as climate models, require calibration (sometimes called tuning). When calibrating expensive simulators using uncertainty quantification methods, it is usually necessary to use a statistical model called an emulator in place of the computer code when running the calibration algorithm. Though emulators based on Gaussian processes are typically many orders of magnitude faster to evaluate than the simulator they mimic, many applications have sought to speed up the computations by using regression-only emulators within the calculations instead, arguing that the extra sophistication brought using the Gaussian process is not worth the extra computational power. This was the case for the analysis that produced the UK climate projections in 2009. In this paper, we compare the effectiveness of both emulation approaches upon a multi-wave calibration framework that is becoming popular in the climate modeling community called "history matching." We find that Gaussian processes offer significant benefits to the reduction of parametric uncertainty over regression-only approaches. We find that in a multi-wave experiment, a combination of regression-only emulators initially, followed by Gaussian process emulators for refocussing experiments can be nearly as effective as using Gaussian processes throughout for a fraction of the computational cost. We also discover a number of design and emulator-dependent features of the multi-wave history matching approach that can cause apparent, yet premature, convergence of our estimates of parametric uncertainty. We compare these approaches to calibration in idealized examples and apply it to a well-known geological reservoir model.