In this article, computation for the purpose of spatial visualization is presented in the context of understanding the variability in global environmental processes. Here, we generate synthetic but realistic global data sets and input them into computational algorithms that have a visualization capability; we call this a simulation-visualization system. Visualization is key here, because the algorithms which we are evaluating must respect the spatial structure of the input. We modify, augment, and integrate four existing component technologies: statistical conditional simulation, Discrete Global Grids (DGGs), Array Set Addressing, and a visualization platform for displaying our results on a globe. The internal representation of the data to be visualized is built around the need for efficient storage and computation as well as the need to move up and downresolutions in a mutually consistent way. In effect, we have constructed a Geographic Information System that is based on a DGG and has desirable data storage, computation, and visualization capabilities. We provide an example of how our simulation-visualization system may be used, by evaluating a computational algorithm called Spatial Statistical Data Fusion that was developed for use on big, remote-sensing data sets.

Climate change is expected to impact the severity of harmful algal blooms in lakes and reservoirs through a number of mechanisms related to the influence of warming temperatures and changes to precipitation patterns. Evidence on the prevalence of individual mechanisms is lacking, however, with knowledge of many mechanisms restricted to studies of individual or small subsets of lakes. Here, we leverage over twelve hundred summertime lake observations from across the continental U.S. to explore evidence for the hypothesized risks from climate change attributable to specific mechanisms. Using a statistical model selection approach, we examine associations between temperature and precipitation variables and indicators of total phytoplankton abundance, species dominance, and toxicity. We find evidence in support of the hypotheses that summer temperatures drive total abundance, that the length of the summer drives cyanobacterial abundance, and that increased temperatures may reduce the observed toxicity of blooms in some cases. We find that nutrient concentrations are also likely to be impacted by lake warming, as increased temperatures are robustly associated with increased total phosphorus concentrations. Evidence for the impact of precipitation is mixed, however, as there is evidence to support that increased nutrient runoff from precipitation could support blooms but also that nutrient concentrations could be reduced through greater flushing due to precipitation. While statistical associations are not definitive evidence of formal mechanistic links, the geographic scale of the results is useful for identifying hypothesized mechanisms that are widespread across the continental U.S., and therefore for informing understanding of the influence of climate change.

Freshwater blooms of phytoplankton affect public health and ecosystem services globally(1,2). Harmful effects of such blooms occur when the intensity of a bloom is too high, or when toxin-producing phytoplankton species are present. Freshwater blooms result in economic losses of more than US$4 billion annually in the United States alone, primarily from harm to aquatic food production, recreation and tourism, and drinkingwater supplies(3). Studies that document bloom conditions in lakes have either focused only on individual or regional subsets of lakes(4-6), or have been limited by a lack of longterm observations(7-9). Here we use three decades of high-resolution Landsat 5 satellite imagery to investigate long-term trends in intense summertime near-surface phytoplankton blooms for 71 large lakes globally. We find that peak summertime bloom intensity has increased in most (68 per cent) of the lakes studied, revealing a global exacerbation of bloom conditions. Lakes that have experienced a significant (P < 0.1) decrease in bloom intensity are rare (8 per cent). The reason behind the increase in phytoplankton bloom intensity remains unclear, however, as temporal trends do not track consistently with temperature, precipitation, fertilizer-use trends or other previously hypothesized drivers. We do find, however, that lakes with a decrease in bloom intensity warmed less compared to other lakes, suggesting that lake warming may already be counteracting management efforts to ameliorate eutrophication(10,11). Our findings support calls for water quality management efforts to better account for the interactions between climate change and local hydrological conditions(12,13).

Terrestrial vegetation removes CO2 from the atmosphere; an important climate regulation service that slows global warming. This 119 Pg C per annum transfer of CO2 into plants-gross primary productivity (GPP)-is the largest land carbon flux globally. While understanding past and anticipated future GPP changes is necessary to support carbon management, the factors driving long-term changes in GPP are largely unknown. Here we show that 1901 to 2010 changes in GPP have been dominated by anthropogenic activity. Our dual constraint attribution approach provides three insights into the spatiotemporal patterns of GPP change. First, anthropogenic controls on GPP change have increased from 57% (1901 decade) to 94% (2001 decade) of the vegetated land surface. Second, CO2 fertilization and nitro gen deposition are the most important drivers of change, 19.8 and 11.1 Pg C per annum (2001 decade) respectively, especially in the tropics and industrialized areas since the 1970's. Third, changes in climate have functioned as fertilization to enhance GPP (1.4 Pg C per annum in the 2001 decade). These findings suggest that, from a land carbon balance perspective, the Anthropocene began over 100 years ago and that global change drivers have allowed GPP uptake to keep pace with anthropogenic emissions.

Divergence in land carbon cycle simulation is persistent and widespread. Regardless of model intercomparison project, results from individual models diverge significantly from each other and, in consequence, from reference datasets. Here we link model spread to structure using a 15-member ensemble of land surface models from the Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) as a test case. Our analysis uses functional benchmarks and model structure as predicted by model skill in a machine learning framework to isolate discrete aspects of model structure associated with divergence. Wealso quantify how initial conditions prejudice present-day model outcomes after centennial-scale transient simulations. Overall, the functional benchmark and machine learning exercises emphasize the importance of ecosystem structure in correctly simulating carbon and water cycling, highlight uncertainties in the structure of carbon pools, and advise against hard parametric limits on ecosystem function. Wealso find that initial conditions explain 90% of variation in global satellite-era values-initial conditions largely predetermine transient endpoints, historical environmental change not withstanding. As MsTMIP prescribes forcing data and spin-up protocol, the range in initial conditions and high levels of predetermination are also structural. Our results suggest that methodological tools linking divergence to discrete aspects of model structure would complement current community best practices in model development.

Given the magnitude of soil carbon stocks in northern ecosystems, and the vulnerability of these stocks to climate warming, land surface models must accurately represent soil carbon dynamics in these regions. We evaluate soil carbon stocks and turnover rates, and the relationship between soil carbon loss with soil temperature and moisture, from an ensemble of eleven global land surface models. We focus on the region of NASA's Arctic-Boreal vulnerability experiment (ABoVE) in North America to inform data collection and model development efforts. Models exhibit an order of magnitude difference in estimates of current total soil carbon stocks, generally under- or overestimating the size of current soil carbon stocks by greater than 50 PgC. We find that a model's soil carbon stock at steady-state in 1901 is the prime driver of its soil carbon stock a hundred years later-overwhelming the effect of environmental forcing factors like climate. The greatest divergence between modeled and observed soil carbon stocks is in regions dominated by peat and permafrost soils, suggesting that models are failing to capture the frozen soil carbon dynamics of permafrost regions. Using a set of functional benchmarks to test the simulated relationship of soil respiration to both soil temperature and moisture, we find that although models capture the observed shape of the soil moisture response of respiration, almost half of the models examined show temperature sensitivities, or Q10 values, that are half of observed. Significantly, models that perform better against observational constraints of respiration or carbon stock size do not necessarily perform well in terms of their functional response to key climatic factors like changing temperature. This suggests that models may be arriving at the right result, but for the wrong reason. The results of this work can help to bridge the gap between data and models by both pointing to the need to constrain initial carbon pool sizes, as well as highlighting the importance of incorporating functional benchmarks into ongoing, mechanistic modeling activities such as those included in ABoVE.

While substantial attention has been paid to the effects of both global climate oscillations and local meteorological conditions on the interannual variability of ecosystem carbon exchange, the relationship between the interannual variability of synoptic meteorology and ecosystem carbon exchange has not been well studied. Here we use a clustering algorithm to identify a summertime cyclonic precipitation system northwest of the Great Lakes to determine (a) the association at a daily scale between the occurrence of this system and the local meteorology and net ecosystem exchange at three Great Lakes region forested eddy covariance sites and (b) the association between the seasonal prevalence of this system and the summertime net ecosystem exchange of these sites. We find that temperature, in addition to precipitation and cloud cover, is an important explanatory factor for the suppression of net ecosystem productivity that occurs during these cyclonic events in this region. In addition, the prevalence of this cyclonic system can explain a significant proportion of the interannual variability in summertime forest ecosystem exchange in this region. This explanatory power is not due to a simple accumulation of low-productivity days that cooccur with this meteorological event, but rather a broader association between the frequency of these events and several aspects of prevailing seasonal conditions. This work demonstrates the usefulness of conceptualizing meteorology in terms of synoptic systems for explaining the interannual variability of regional carbon fluxes.

Large uncertainties in North American terrestrial carbon fluxes hinder regional climate projections. Terrestrial biosphere models (TBMs), the essential tools for understanding continental-scale carbon cycle, diverge on whether temperate forests or croplands dominate carbon uptake in North America. Evidence from novel photosynthetic proxies, such as those based on chlorophyll fluorescence, has cast doubt on the "weak cropland, strong forest" carbon uptake patterns simulated by most TBMs. However, no systematic evaluation of TBMs has yet been attempted to pin down space-time patterns that are most consistent with regional CO2 observational constraints. Here, we leverage atmospheric CO2 observations and satellite-observed photosynthetic proxies to understand emergent space-time patterns in North American carbon fluxes from a large suite of TBMs and data-driven models. To do so, we evaluate how well the atmospheric signals resulting from carbon flux estimates reproduce the space-time variability in atmospheric CO2, as is observed by a network of continuous-monitoring towers over North America. Models with gross or net carbon fluxes that are consistent with the observed CO2 variability share a salient feature of growing-season carbon uptake in Midwest US croplands. Conversely, the remaining models place most growing-season uptake in boreal or temperate forests. Differences in model explanatory power depend mainly on the simulated annual cycles of cropland uptake-especially, the timing of peak uptake-rather than the distribution of annual mean fluxes across biomes. Our results suggest that improved model representation of cropland phenology is crucial to robust, policy-relevant estimation of North American carbon exchange.