We assess the detectability of city emissions via a tower-based greenhouse gas (GHG) network, as part of the Indianapolis Flux (INFLUX) experiment. By examining afternoon-averaged results from a network of carbon dioxide (CO₂), methane (CH₄), and carbon monoxide (CO) mole fraction measurements in Indianapolis, Indiana for 2011–2013, we quantify spatial and temporal patterns in urban atmospheric GHG dry mole fractions. The platform for these measurements is twelve communications towers spread across the metropolitan region, ranging in height from 39 to 136 m above ground level, and instrumented with cavity ring-down spectrometers. Nine of the sites were deployed as of January 2013 and data from these sites are the focus of this paper. A background site, chosen such that it is on the predominantly upwind side of the city, is utilized to quantify enhancements caused by urban emissions. Afternoon averaged mole fractions are studied because this is the time of day during which the height of the boundary layer is most steady in time and the area that influences the tower measurements is likely to be largest. Additionally, atmospheric transport models have better performance in simulating the daytime convective boundary layer compared to the nighttime boundary layer. Averaged from January through April of 2013, the mean urban dormant-season enhancements range from 0.3 ppm CO₂ at the site 24 km typically downwind of the edge of the city (Site 09) to 1.4 ppm at the site at the downwind edge of the city (Site 02) to 2.9 ppm at the downtown site (Site 03). When the wind is aligned such that the sites are downwind of the urban area, the enhancements are increased, to 1.6 ppm at Site 09, and 3.3 ppm at Site 02. Differences in sampling height affect the reported urban enhancement by up to 50%, but the overall spatial pattern remains similar. The time interval over which the afternoon data are averaged alters the calculated urban enhancement by an average of 0.4 ppm. The CO₂ observations are compared to CO₂ mole fractions simulated using a mesoscale atmospheric model and an emissions inventory for Indianapolis. The observed and modeled CO₂ enhancements are highly correlated (r² = 0.94), but the modeled enhancements prior to inversion average 53% of those measured at the towers. Following the inversion, the enhancements follow the observations closely, as expected. The CH₄ urban enhancement ranges from 5 ppb at the site 10 km predominantly downwind of the city (Site 13) to 21 ppb at the site near the landfill (Site 10), and for CO ranges from 6 ppb at the site 24 km downwind of the edge of the city (Site 09) to 29 ppb at the downtown site (Site 03). Overall, these observations show that a dense network of urban GHG measurements yield a detectable urban signal, well-suited as input to an urban inversion system given appropriate attention to sampling time, sampling altitude and quantification of background conditions.

The ‘Hestia Project’ uses a bottom-up approach to quantify fossil fuel CO₂(FFCO₂) emissions spatially at the building/street level and temporally at the hourly level. Hestia FFCO₂ emissions are provided in the form of a group of sector-specific vector layers with point, line, and polygon sources to support carbon cycle science and climate policy. Application to carbon cycle science, in particular, requires regular gridded data in order to link surface carbon fluxes to atmospheric transport models. However, the heterogeneity and complexity of FFCO₂ sources within regular grids is sensitive to spatial resolution. From the perspective of a data provider, we need to find a balance between resolution and data volume so that the gridded data product retains the maximum amount of information content while maintaining an efficient data volume.

The Shannon entropy determines the minimum bits that are needed to encode an information source and can serve as a metric for the effective information content. In this paper, we present an analysis of the Shannon entropy of gridded FFCO₂ emissions with varying resolutions in four Hestia study areas, and find: (1) the Shannon entropy increases with smaller grid resolution until it reaches a maximum value (the max-entropy resolution); (2) total emissions (the sum of several sector-specific emission fields) show a finer max-entropy resolution than each of the sector-specific fields; (3) the residential emissions show a finer max-entropy resolution than the commercial emissions; (4) the max-entropy resolution of the onroad emissions grid is closely correlated to the density of the road network. These findings suggest that the Shannon entropy can detect the information effectiveness of the spatial resolution of gridded FFCO₂ emissions. Hence, the resolution-entropy relationship can be used to assist in determining an appropriate spatial resolution for urban CO₂ flux studies. We conclude that the optimal spatial resolution for providing Hestia total FFCO₂ emissions products is centered around 100 m, at which the FFCO₂ emissions data can not only fully meet the requirement of urban flux integration, but also be effectively used in understanding the relationships between FFCO₂ emissions and various social-economic variables at the U.S. census block group level.

Atmospheric CO₂ inversions estimate surface carbon fluxes from an optimal fit to atmospheric CO₂ measurements, usually including prior constraints on the flux estimates. Eleven sets of carbon flux estimates are compared, generated by different inversions systems that vary in their inversions methods, choice of atmospheric data, transport model and prior information. The inversions were run for at least 5 yr in the period between 1990 and 2010. Mean fluxes for 2001–2004, seasonal cycles, interannual variability and trends are compared for the tropics and northern and southern extra-tropics, and separately for land and ocean. Some continental/basin-scale subdivisions are also considered where the atmospheric network is denser. Four-year mean fluxes are reasonably consistent across inversions at global/latitudinal scale, with a large total (land plus ocean) carbon uptake in the north (−3.4 Pg C yr^-1 (±0.5 Pg C yr^-1 standard deviation), with slightly more uptake over land than over ocean), a significant although more variable source over the tropics (1.6 ± 0.9 Pg C yr^-1) and a compensatory sink of similar magnitude in the south (−1.4 ± 0.5 Pg C yr^-1) corresponding mainly to an ocean sink. Largest differences across inversions occur in the balance between tropical land sources and southern land sinks. Interannual variability (IAV) in carbon fluxes is larger for land than ocean regions (standard deviation around 1.06 versus 0.33 Pg C yr[superscript −1] for the 1996–2007 period), with much higher consistency among the inversions for the land. While the tropical land explains most of the IAV (standard deviation ~ 0.65 Pg C yr^-1), the northern and southern land also contribute (standard deviation ~ 0.39 Pg C yr^-1). Most inversions tend to indicate an increase of the northern land carbon uptake from late 1990s to 2008 (around 0.1 Pg C yr^-1, predominantly in North Asia. The mean seasonal cycle appears to be well constrained by the atmospheric data over the northern land (at the continental scale), but still highly dependent on the prior flux seasonality over the ocean. Finally we provide recommendations to interpret the regional fluxes, along with the uncertainty estimates.

Quantifying greenhouse gas (GHG) emissions from cities is a key challenge towards effective emissions management. An inversion analysis from the INdianapolis FLUX experiment (INFLUX) project, as the first of its kind, has achieved a top-down emission estimate for a single city using CO₂ data collected by the dense tower network deployed across the city. However, city-level emission data, used as a priori emissions, are also a key component in the atmospheric inversion framework. Currently, fine-grained emission inventories (EIs) able to resolve GHG city emissions at high spatial resolution, are only available for few major cities across the globe. Following the INFLUX inversion case with a global 1.1 km ODIAC fossil fuel CO₂ emission dataset, we further improved the ODIAC emission field and examined its utility as a prior for the city scale inversion. We disaggregated the 1.1 km ODIAC non-point source emissions using geospatial datasets such as the global road network data and satellite-data driven surface imperviousness data to a 30 . 30 m resolution. We assessed the impact of the improved emission field on the inversion result, relative to priors in previous studies (Hestia and ODIAC). The posterior total emission estimate (5.1 MtC/yr) remains statistically similar to the previous estimate with ODIAC (5.3 MtC/yr). However, the distribution of the flux corrections was very close to those of Hestia inversion and the model-observation mismatches were significantly reduced both in forward and inverse runs, even without hourly temporal changes in emissions. EIs reported by cities often do not have estimates of spatial extents. Thus, emission disaggregation is a required step when verifying those reported emissions using atmospheric models. Our approach offers gridded emission estimates for global cities that could serves as a prior for inversion, even without locally reported EIs in a systematic way to support city-level Measuring, Reporting and Verification (MRV) practice implementation.

This paper presents an analysis of methane emissions from the Los Angeles Basin at monthly timescales across a 4-year time period – from September 2011 to August 2015. Using observations acquired by a ground-based near-infrared remote sensing instrument on Mount Wilson, California, combined with atmospheric CH₄–CO₂ tracer–tracer correlations, we observed −18 to +22 % monthly variability in CH₄ : CO₂ from the annual mean in the Los Angeles Basin. Top-down estimates of methane emissions for the basin also exhibit significant monthly variability (−19 to +31 % from annual mean and a maximum month-to-month change of 47 %). During this period, methane emissions consistently peaked in the late summer/early fall and winter. The estimated annual methane emissions did not show a statistically significant trend over the 2011 to 2015 time period.

In spite of the recent interest and advances in linear controllability of complex networks, controlling nonlinear network dynamics remains an outstanding problem. Here we develop an experimentally feasible control framework for nonlinear dynamical networks that exhibit multistability. The control objective is to apply parameter perturbation to drive the system from one attractor to another, assuming that the former is undesired and the latter is desired. To make our framework practically meaningful, we consider restricted parameter perturbation by imposing two constraints: it must be experimentally realizable and applied only temporarily. We introduce the concept of attractor network, which allows us to formulate a quantifiable controllability framework for nonlinear dynamical networks: a network is more controllable if the attractor network is more strongly connected. We test our control framework using examples from various models of experimental gene regulatory networks and demonstrate the beneficial role of noise in facilitating control.

We investigate the emergence of extreme events in interdependent networks. We introduce an inter-layer traffic resource competing mechanism to account for the limited capacity associated with distinct network layers. A striking finding is that, when the number of network layers and/or the overlap among the layers are increased, extreme events can emerge in a cascading manner on a global scale. Asymptotically, there are two stable absorption states: a state free of extreme events and a state of full of extreme events, and the transition between them is abrupt. Our results indicate that internal interactions in the multiplex system can yield qualitatively distinct phenomena associated with extreme events that do not occur for independent network layers. An implication is that, e.g., public resource competitions among different service providers can lead to a higher resource requirement than naively expected. We derive an analytical theory to understand the emergence of global-scale extreme events based on the concept of effective betweenness. We also articulate a cost-effective control scheme through increasing the capacity of very few hubs to suppress the cascading process of extreme events so as to protect the entire multi-layer infrastructure against global-scale breakdown.

A remarkable phenomenon in spatiotemporal dynamical systems is chimera state, where the structurally and dynamically identical oscillators in a coupled networked system spontaneously break into two groups, one exhibiting coherent motion and another incoherent. This phenomenon was typically studied in the setting of non-local coupling configurations. We ask what can happen to chimera states under systematic changes to the network structure when links are removed from the network in an orderly fashion but the local coupling topology remains invariant with respect to an index shift. We find the emergence of multicluster chimera states. Remarkably, as a parameter characterizing the amount of link removal is increased, chimera states of distinct numbers of clusters emerge and persist in different parameter regions. We develop a phenomenological theory, based on enhanced or reduced interactions among oscillators in different spatial groups, to explain why chimera states of certain numbers of clusters occur in certain parameter regions. The theoretical prediction agrees well with numerics.

Successful identification of directed dynamical influence in complex systems is relevant to significant problems of current interest. Traditional methods based on Granger causality and transfer entropy have issues such as difficulty with nonlinearity and large data requirement. Recently a framework based on nonlinear dynamical analysis was proposed to overcome these difficulties. We find, surprisingly, that noise can counterintuitively enhance the detectability of directed dynamical influence. In fact, intentionally injecting a proper amount of asymmetric noise into the available time series has the unexpected benefit of dramatically increasing confidence in ascertaining the directed dynamical influence in the underlying system. This result is established based on both real data and model time series from nonlinear ecosystems. We develop a physical understanding of the beneficial role of noise in enhancing detection of directed dynamical influence.

Attributing observed CO2 variations to human or natural cause is critical to deducing and tracking emissions from observations. We have used in situ CO2, CO, and planetary boundary layer height (PBLH) measurements recorded during the CalNex-LA (CARB et al., 2008) ground campaign of 15 May-15 June 2010, in Pasadena, CA, to deduce the diurnally varying anthropogenic component of observed CO2 in the megacity of Los Angeles (LA). This affordable and simple technique, validated by carbon isotope observations and WRF-STILT (Weather Research and Forecasting model - Stochastic Time-Inverted Lagrangian Transport model) predictions, is shown to robustly attribute observed CO2 variation to anthropogenic or biogenic origin over the entire diurnal cycle. During CalNex-LA, local fossil fuel combustion contributed up to similar to 50% of the observed CO2 enhancement overnight, and similar to 100% of the enhancement near midday. This suggests that sufficiently accurate total column CO2 observations recorded near midday, such as those from the GOSAT or OCO-2 satellites, can potentially be used to track anthropogenic emissions from the LA megacity.