As climate continues to change, scientists are left to analyze the effects these changes will have on the public. In this article, a flexible class of distributed lag models are used to analyze the effects of heat on mortality in four major metropolitan areas in the U.S. (Chicago, Dallas, Los Angeles, and New York). Specifically, the proposed methodology uses Gaussian processes as a prior model for the distributed lag function. Gaussian processes are adequately flexible to capture a wide variety of distributed lag functions while ensuring smoothness properties of process realizations. Additionally, the proposed framework allows for probabilistic inference of the maximum lag. Applying the proposed methodology revealed that mortality displacement (or, harvesting) was present for most age groups and cities analyzed suggesting that heat advanced death in some individuals. Additionally, the estimated shape of the DL functions gave evidence that prolonged heat exposure and highly variable temperatures pose a threat to public health.
