Urbanization and edge effects interact to drive mutualism breakdown and the rise of unstable pathogenic communities in forest soil.

Temperate forests are threatened by urbanization and fragmentation, with over 20% (118,300 km2) of U.S. forest land projected to be subsumed by urban land development. We leveraged a unique, well-characterized urban-to-rural and forest edge-to-interior gradient to identify the combined impact of these two land use changes-urbanization and forest edge creation-on the soil microbial community in native remnant forests. We found evidence of mutualism breakdown between trees and their fungal root mutualists [ectomycorrhizal (ECM) fungi] with urbanization, where ECM fungi colonized fewer tree roots and had less connectivity in soil microbiome networks in urban forests compared to rural forests. However, urbanization did not reduce the relative abundance of ECM fungi in forest soils; instead, forest edges alone led to strong reductions in ECM fungal abundance. At forest edges, ECM fungi were replaced by plant and animal pathogens, as well as copiotrophic, xenobiotic-degrading, and nitrogen-cycling bacteria, including nitrifiers and denitrifiers. Urbanization and forest edges interacted to generate new "suites" of microbes, with urban interior forests harboring highly homogenized microbiomes, while edge forest microbiomes were more heterogeneous and less stable, showing increased vulnerability to low soil moisture. When scaled to the regional level, we found that forest soils are projected to harbor high abundances of fungal pathogens and denitrifying bacteria, even in rural areas, due to the widespread existence of forest edges. Our results highlight the potential for soil microbiome dysfunction-including increased greenhouse gas production-in temperate forest regions that are subsumed by urban expansion, both now and in the future.

Urbanization exacerbates climate sensitivity of eastern United States broadleaf trees.

Tree growth is a key mechanism driving carbon sequestration in forest ecosystems. Environmental conditions are important regulators of tree growth that can vary considerably between nearby urban and rural forests. For example, trees growing in cities often experience hotter and drier conditions than their rural counterparts while also being exposed to higher levels of light, pollution, and nutrient inputs. However, the extent to which these intrinsic differences in the growing conditions of trees in urban versus rural forests influence tree growth response to climate is not well known. In this study, we tested for differences in the climate sensitivity of tree growth between urban and rural forests along a latitudinal transect in the eastern United States that included Boston, Massachusetts, New York City, New York, and Baltimore, Maryland. Using dendrochronology analyses of tree cores from 55 white oak trees (Quercus alba), 55 red maple trees (Acer rubrum), and 41 red oak trees (Quercus rubra) we investigated the impacts of heat stress and water stress on the radial growth of individual trees. Across our three-city study, we found that tree growth was more closely correlated with climate stress in the cooler climate cities of Boston and New York than in Baltimore. Furthermore, heat stress was a significant hindrance to tree growth in higher latitudes while the impacts of water stress appeared to be more evenly distributed across latitudes. We also found that the growth of oak trees, but not red maple trees, in the urban sites of Boston and New York City was more adversely impacted by heat stress than their rural counterparts, but we did not see these urban-rural differences in Maryland. Trees provide a wide range of important ecosystem services and increasing tree canopy cover was typically an important component of urban sustainability strategies. In light of our findings that urbanization can influence how tree growth responds to a warming climate, we suggest that municipalities consider these interactions when developing their tree-planting palettes and when estimating the capacity of urban forests to contribute to broader sustainability goals in the future.

Underestimation of Thermogenic Methane Emissions in New York City.

Recent studies have shown that methane emissions are underestimated by inventories in many US urban areas. This has important implications for climate change mitigation policy at the city, state, and national levels. Uncertainty in both the spatial distribution and sectoral allocation of urban emissions can limit the ability of policy makers to develop appropriately focused emission reduction strategies. Top-down emission estimates based on atmospheric greenhouse gas measurements can help to improve inventories and inform policy decisions. This study presents a new high-resolution (0.02 × 0.02°) methane emission inventory for New York City and its surrounding area, constructed using the latest activity data, emission factors, and spatial proxies. The new high-resolution inventory estimates of methane emissions for the New York-Newark urban area are 1.3 times larger than those for the gridded Environmental Protection Agency inventory. We used aircraft mole fraction measurements from nine research flights to optimize the high-resolution inventory emissions within a Bayesian inversion. These sectorally optimized emissions show that the high-resolution inventory still significantly underestimates methane emissions within the New York-Newark urban area, primarily because it underestimates emissions from thermogenic sources (by a factor of 2.3). This suggests that there remains a gap in our process-based understanding of urban methane emissions.

Majority of US urban natural gas emissions unaccounted for in inventories.

Across many cities, estimates of methane emissions from natural gas (NG) distribution and end use based on atmospheric measurements have generally been more than double bottom-up estimates. We present a top-down study of NG methane emissions from the Boston urban region spanning 8 y (2012 to 2020) to assess total emissions, their seasonality, and trends. We used methane and ethane observations from five sites in and around Boston, combined with a high-resolution transport model, to calculate methane emissions of 76 ± 18 Gg/yr, with 49 ± 9 Gg/yr attributed to NG losses. We found no significant trend in the NG loss rate over 8 y, despite efforts from the city and state to increase the rate of repairing NG pipeline leaks. We estimate that 2.5 ± 0.5% of the gas entering the urban region is lost, approximately three times higher than bottom-up estimates. We saw a strong correlation between top-down NG emissions and NG consumed on a seasonal basis. This suggests that consumption-driven losses, such as in transmission or end-use, may be a large component of emissions that is missing from inventories, and require future policy action. We also compared top-down NG emission estimates from six US cities, all of which indicate significant missing sources in bottom-up inventories. Across these cities, we estimate NG losses from distribution and end use amount to 20 to 36% of all losses from the US NG supply chain, with a total loss rate of 3.3 to 4.7% of NG from well pad to urban consumer, notably larger than the current Environmental Protection Agency estimate of 1.4% [R. A. Alvarez et al., Science 361, 186-188 (2018)].

Urbanization and fragmentation have opposing effects on soil nitrogen availability in temperate forest ecosystems.

Nitrogen (N) availability relative to plant demand has been declining in recent years in terrestrial ecosystems throughout the world, a phenomenon known as N oligotrophication. The temperate forests of the northeastern U.S. have experienced a particularly steep decline in bioavailable N, which is expected to be exacerbated by climate change. This region has also experienced rapid urban expansion in recent decades that leads to forest fragmentation, and it is unknown whether and how these changes affect N availability and uptake by forest trees. Many studies have examined the impact of either urbanization or forest fragmentation on nitrogen (N) cycling, but none to our knowledge have focused on the combined effects of these co-occurring environmental changes. We examined the effects of urbanization and fragmentation on oak-dominated (Quercus spp.) forests along an urban to rural gradient from Boston to central Massachusetts (MA). At eight study sites along the urbanization gradient, plant and soil measurements were made along a 90 m transect from a developed edge to an intact forest interior. Rates of net ammonification, net mineralization, and foliar N concentrations were significantly higher in urban than rural sites, while net nitrification and foliar C:N were not different between urban and rural forests. At urban sites, foliar N and net ammonification and mineralization were higher at forest interiors compared to edges, while net nitrification and foliar C:N were higher at rural forest edges than interiors. These results indicate that urban forests in the northeastern U.S. have greater soil N availability and N uptake by trees compared to rural forests, counteracting the trend for widespread N oligotrophication in temperate forests around the globe. Such increases in available N are diminished at forest edges, however, demonstrating that forest fragmentation has the opposite effect of urbanization on coupled N availability and demand by trees.

Diverging patterns at the forest edge: Soil respiration dynamics of fragmented forests in urban and rural areas.

As urbanization and forest fragmentation increase around the globe, it is critical to understand how rates of respiration and carbon losses from soil carbon pools are affected by these processes. This study characterizes soils in fragmented forests along an urban to rural gradient, evaluating the sensitivity of soil respiration to changes in soil temperature and moisture near the forest edge. While previous studies found elevated rates of soil respiration at temperate forest edges in rural areas compared to the forest interior, we find that soil respiration is suppressed at the forest edge in urban areas. At urban sites, respiration rates are 25% lower at the forest edge relative to the interior, likely due to high temperature and aridity conditions near urban edges. While rural soils continue to respire with increasing temperatures, urban soil respiration rates asymptote as temperatures climb and soils dry. Soil temperature- and moisture-sensitivity modeling shows that respiration rates in urban soils are less sensitive to rising temperatures than those in rural soils. Scaling these results to Massachusetts (MA), which encompasses 0.25 Mha of the urban forest, we find that failure to account for decreases in soil respiration rates near urban forest edges leads to an overestimate of growing-season soil carbon fluxes of >350,000 Mg C. This difference is almost 2.5 times that for rural soils in the analogous comparison (underestimate of <143,000 Mg C), even though rural forest area is more than four times greater than urban forest area in MA. While a changing climate may stimulate carbon losses from rural forest edge soils, urban forests may experience enhanced soil carbon sequestration near the forest edge. These findings highlight the need to capture the effects of forest fragmentation and land use context when making projections about soil behavior and carbon cycling in a warming and increasingly urbanized world.

Anthropogenic and biogenic CO2 fluxes in the Boston urban region.

With the pending withdrawal of the United States from the Paris Climate Accord, cities are now leading US actions toward reducing greenhouse gas emissions. Implementing effective mitigation strategies requires the ability to measure and track emissions over time and at various scales. We report CO2 emissions in the Boston, MA, urban region from September 2013 to December 2014 based on atmospheric observations in an inverse model framework. Continuous atmospheric measurements of CO2 from five sites in and around Boston were combined with a high-resolution bottom-up CO2 emission inventory and a Lagrangian particle dispersion model to determine regional emissions. Our model-measurement framework incorporates emissions estimates from submodels for both anthropogenic and biological CO2 fluxes, and development of a CO2 concentration curtain at the boundary of the study region based on a combination of tower measurements and modeled vertical concentration gradients. We demonstrate that an emission inventory with high spatial and temporal resolution and the inclusion of urban biological fluxes are both essential to accurately modeling annual CO2 fluxes using surface measurement networks. We calculated annual average emissions in the Boston region of 0.92 kg C·m-2·y-1 (95% confidence interval: 0.79 to 1.06), which is 14% higher than the Anthropogenic Carbon Emissions System inventory. Based on the capability of the model-measurement approach demonstrated here, our framework should be able to detect changes in CO2 emissions of greater than 18%, providing stakeholders with critical information to assess mitigation efforts in Boston and surrounding areas.

Edge effects enhance carbon uptake and its vulnerability to climate change in temperate broadleaf forests.

Forest fragmentation is a ubiquitous, ongoing global phenomenon with profound impacts on the growing conditions of the world's remaining forest. The temperate broadleaf forest makes a large contribution to the global terrestrial carbon sink but is also the most heavily fragmented forest biome in the world. We use field measurements and geospatial analyses to characterize carbon dynamics in temperate broadleaf forest fragments. We show that forest growth and biomass increase by 89 ± 17% and 64 ± 12%, respectively, from the forest interior to edge, but ecosystem edge enhancements are not currently captured by models or approaches to quantifying regional C balance. To the extent that the findings from our research represent the forest of southern New England in the United States, we provide a preliminary estimate that edge growth enhancement could increase estimates of the region's carbon uptake and storage by 13 ± 3% and 10 ± 1%, respectively. However, we also find that forest growth near the edge declines three times faster than that in the interior in response to heat stress during the growing season. Using climate projections, we show that future heat stress could reduce the forest edge growth enhancement by one-third by the end of the century. These findings contrast studies of edge effects in the world's other major forest biomes and indicate that the strength of the temperate broadleaf forest carbon sink and its capacity to mitigate anthropogenic carbon emissions may be stronger, but also more sensitive to climate change than previous estimates suggest.

Cities, traffic, and CO2: A multidecadal assessment of trends, drivers, and scaling relationships.

Emissions of CO2 from road vehicles were 1.57 billion metric tons in 2012, accounting for 28% of US fossil fuel CO2 emissions, but the spatial distributions of these emissions are highly uncertain. We develop a new emissions inventory, the Database of Road Transportation Emissions (DARTE), which estimates CO2 emitted by US road transport at a resolution of 1 km annually for 1980-2012. DARTE reveals that urban areas are responsible for 80% of on-road emissions growth since 1980 and for 63% of total 2012 emissions. We observe nonlinearities between CO2 emissions and population density at broad spatial/temporal scales, with total on-road CO2 increasing nonlinearly with population density, rapidly up to 1,650 persons per square kilometer and slowly thereafter. Per capita emissions decline as density rises, but at markedly varying rates depending on existing densities. We make use of DARTE's bottom-up construction to highlight the biases associated with the common practice of using population as a linear proxy for disaggregating national- or state-scale emissions. Comparing DARTE with existing downscaled inventories, we find biases of 100% or more in the spatial distribution of urban and rural emissions, largely driven by mismatches between inventory downscaling proxies and the actual spatial patterns of vehicle activity at urban scales. Given cities' dual importance as sources of CO2 and an emerging nexus of climate mitigation initiatives, high-resolution estimates such as DARTE are critical both for accurately quantifying surface carbon fluxes and for verifying the effectiveness of emissions mitigation efforts at urban scales.

Methane emissions from natural gas infrastructure and use in the urban region of Boston, Massachusetts.

Methane emissions from natural gas delivery and end use must be quantified to evaluate the environmental impacts of natural gas and to develop and assess the efficacy of emission reduction strategies. We report natural gas emission rates for 1 y in the urban region of Boston, using a comprehensive atmospheric measurement and modeling framework. Continuous methane observations from four stations are combined with a high-resolution transport model to quantify the regional average emission flux, 18.5 ± 3.7 (95% confidence interval) g CH4 ⋅ m(-2) ⋅ y(-1). Simultaneous observations of atmospheric ethane, compared with the ethane-to-methane ratio in the pipeline gas delivered to the region, demonstrate that natural gas accounted for ∼ 60-100% of methane emissions, depending on season. Using government statistics and geospatial data on natural gas use, we find the average fractional loss rate to the atmosphere from all downstream components of the natural gas system, including transmission, distribution, and end use, was 2.7 ± 0.6% in the Boston urban region, with little seasonal variability. This fraction is notably higher than the 1.1% implied by the most closely comparable emission inventory.

Characterizing forest structure variations across an intact tropical peat dome using field samplings and airborne LiDAR.

Tropical peat swamp forests (PSF) are one of the most carbon dense ecosystems on the globe and are experiencing substantial natural and anthropogenic disturbances. In this study, we combined direct field sampling and airborne LiDAR to empirically quantify forest structure and aboveground live biomass (AGB) across a large, intact tropical peat dome in Northwestern Borneo. Moving up a 4 m elevational gradient, we observed increasing stem density but decreasing canopy height, crown area, and crown roughness. These findings were consistent with hypotheses that nutrient and hydrological dynamics co-influence forest structure and stature of the canopy individuals, leading to reduced productivity towards the dome interior. Gap frequency as a function of gap size followed a power law distribution with a shape factor (λ) of 1.76 ± 0.06. Ground-based and dome-wide estimates of AGB were 217.7 ± 28.3 Mg C/ha and 222.4 ± 24.4 Mg C/ha, respectively, which were higher than previously reported AGB for PSF and tropical forests in general. However, dome-wide AGB estimates were based on height statistics, and we found the coefficient of variation on canopy height was only 0.08, three times less than stem diameter measurements, suggesting LiDAR height metrics may not be a robust predictor of AGB in tall tropical forests with dense canopies. Our structural characterization of this ecosystem advances the understanding of the ecology of intact tropical peat domes and factors that influence biomass density and landscape-scale spatial variation. This ecological understanding is essential to improve estimates of forest carbon density and its spatial distribution in PSF and to effectively model the effects of disturbance and deforestation in these carbon dense ecosystems.

Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools.

Urban land-cover change threatens biodiversity and affects ecosystem productivity through loss of habitat, biomass, and carbon storage. However, despite projections that world urban populations will increase to nearly 5 billion by 2030, little is known about future locations, magnitudes, and rates of urban expansion. Here we develop spatially explicit probabilistic forecasts of global urban land-cover change and explore the direct impacts on biodiversity hotspots and tropical carbon biomass. If current trends in population density continue and all areas with high probabilities of urban expansion undergo change, then by 2030, urban land cover will increase by 1.2 million km(2), nearly tripling the global urban land area circa 2000. This increase would result in considerable loss of habitats in key biodiversity hotspots, with the highest rates of forecasted urban growth to take place in regions that were relatively undisturbed by urban development in 2000: the Eastern Afromontane, the Guinean Forests of West Africa, and the Western Ghats and Sri Lanka hotspots. Within the pan-tropics, loss in vegetation biomass from areas with high probability of urban expansion is estimated to be 1.38 PgC (0.05 PgC yr(-1)), equal to ∼5% of emissions from tropical deforestation and land-use change. Although urbanization is often considered a local issue, the aggregate global impacts of projected urban expansion will require significant policy changes to affect future growth trajectories to minimize global biodiversity and vegetation carbon losses.

Reduced impact logging minimally alters tropical rainforest carbon and energy exchange.

We used eddy covariance and ecological measurements to investigate the effects of reduced impact logging (RIL) on an old-growth Amazonian forest. Logging caused small decreases in gross primary production, leaf production, and latent heat flux, which were roughly proportional to canopy loss, and increases in heterotrophic respiration, tree mortality, and wood production. The net effect of RIL was transient, and treatment effects were barely discernable after only 1 y. RIL appears to provide a strategy for managing tropical forest that minimizes the potential risks to climate associated with large changes in carbon and water exchange.

A bottom up approach to on-road CO2 emissions estimates: improved spatial accuracy and applications for regional planning.

On-road transportation is responsible for 28% of all U.S. fossil-fuel CO2 emissions. Mapping vehicle emissions at regional scales is challenging due to data limitations. Existing emission inventories use spatial proxies such as population and road density to downscale national or state-level data. Such procedures introduce errors where the proxy variables and actual emissions are weakly correlated, and limit analysis of the relationship between emissions and demographic trends at local scales. We develop an on-road emission inventory product for Massachusetts-based on roadway-level traffic data obtained from the Highway Performance Monitoring System (HPMS). We provide annual estimates of on-road CO2 emissions at a 1 × 1 km grid scale for the years 1980 through 2008. We compared our results with on-road emissions estimates from the Emissions Database for Global Atmospheric Research (EDGAR), with the Vulcan Product, and with estimates derived from state fuel consumption statistics reported by the Federal Highway Administration (FHWA). Our model differs from FHWA estimates by less than 8.5% on average, and is within 5% of Vulcan estimates. We found that EDGAR estimates systematically exceed FHWA by an average of 22.8%. Panel regression analysis of per-mile CO2 emissions on population density at the town scale shows a statistically significant correlation that varies systematically in sign and magnitude as population density increases. Population density has a positive correlation with per-mile CO2 emissions for densities below 2000 persons km(-2), above which increasing density correlates negatively with per-mile emissions.

Urban green space and albedo impacts on surface temperature across seven United States cities.

Extreme heat represents a growing threat to public health, especially across the densely populated, developed landscape of cities. Climate adaptation strategies that aim to manage urban microclimates through purposeful design can reduce the heat exposure of urban populations, however, it is unclear how the temperature impacts of urban green space and albedo vary across cities and background climate. This study quantifies the sensitivity of surface temperature to landcover characteristics tied to two widely used climate adaptation strategies, urban greening and albedo manipulation (e.g. white roofs), by combining long-term remote sensing observations of land surface temperature, albedo, and moisture with high-resolution landcover datasets in a spatial regression analysis at the census block scale across seven United States cities. We find tree cover to have an average cooling impact of -0.089 K per % cover, which is approximately four times stronger than the average grass cover cooling impact of -0.021 K per % cover. Variability in the magnitude of grass cover cooling impacts was primarily a function of vegetation moisture content, with the Land Surface Water Index (LSWI) explaining 89 % of the variability in grass cover cooling impacts across cities. Variability in tree cover cooling impacts was primarily a function of sunlight and vegetation moisture content, with solar irradiance and LSWI explaining 97 % of the cooling variability across cities. Albedo cooling impacts were consistent across cities with an average cooling impact of -0.187 K per increase of 0.01. While these interventions are broadly effective across cities, there are critical regional trade-offs between vegetation cooling efficiency, irrigation requirements, and the temporal duration and evolution of the cooling benefits. In warm, arid cities, high albedo surfaces offer multifaceted benefits such as cooling and water conservation, whereas temperate, mesic cities likely benefit from a combination of strategies, with greening efforts targeting highly paved neighborhoods.

Soils at the temperate forest edge: An investigation of soil characteristics and carbon dynamics.

Global proliferation of forest edges through anthropogenic land-use change and forest fragmentation is well documented, and while forest fragmentation has clear consequences for soil carbon (C) cycling, underlying drivers of belowground activity at the forest edge remain poorly understood. Increasing soil C losses via respiration have been observed at rural forest edges, but this process was suppressed at urban forest edges. We offer a comprehensive, coupled investigation of abiotic soil conditions and biotic soil activity from forest edge to interior at eight sites along an urbanization gradient to elucidate how environmental stressors are linked to soil C cycling at the forest edge. Despite significant diverging trends in edge soil C losses between urban and rural sites, we did not find comparable differences in soil % C or microbial enzyme activity, suggesting an unexpected decoupling of soil C fluxes and pools at forest edges. We demonstrate that across site types, soils at forest edges were less acidic than the forest interior (p < 0.0001), and soil pH was positively correlated with soil calcium, magnesium and sodium content (adj R2 = 0.37), which were also elevated at the edge. Compared to forest interior, forest edge soils exhibited a 17.8 % increase in sand content and elevated freeze-thaw frequency with probable downstream effects on root turnover and decomposition. Using these and other novel forest edge data, we demonstrate that significant variation in edge soil respiration (adj R2 = 0.46; p = 0.0002) and C content (adj R2 = 0.86; p < 0.0001) can be explained using soil parameters often mediated by human activity (e.g., soil pH, trace metal and cation concentrations, soil temperature), and we emphasize the complex influence of multiple, simultaneous global change drivers at forest edges. Forest edge soils reflect legacies of anthropogenic land-use and modern human management, and this must be accounted for to understand soil activity and C cycling across fragmented landscapes.

School Greenness and Student-Level Academic Performance: Evidence From the Global South.

Greenspace in schools might enhance students' academic performance. However, the literature-dominated by ecological studies at the school level in countries from the Northern Hemisphere-presents mixed evidence of a beneficial association. We evaluated the association between school greenness and student-level academic performance in Santiago, Chile, a capital city of the Global South. This cross-sectional study included 281,695 fourth-grade students attending 1,498 public, charter, and private schools in Santiago city between 2014 and 2018. Student-level academic performance was assessed using standardized test scores and indicators of attainment of learning standards in mathematics and reading. School greenness was estimated using Normalized Difference Vegetation Index (NDVI). Linear and generalized linear mixed-effects models were fit to evaluate associations, adjusting for individual- and school-level sociodemographic factors. Analyses were stratified by school type. In fully adjusted models, a 0.1 increase in school greenness was associated with higher test scores in mathematics (36.9 points, 95% CI: 2.49; 4.88) and in reading (1.84 points, 95% CI: 0.73; 2.95); as well as with higher odds of attaining learning standards in mathematics (OR: 1.20, 95% CI: 1.12; 1.28) and reading (OR: 1.07, 95% CI: 1.02; 1.13). Stratified analysis showed differences by school type, with associations of greater magnitude and strength for students attending public schools. No significant associations were detected for students in private schools. Higher school greenness was associated with improved individual-level academic outcomes among elementary-aged students in a capital city in South America. Our results highlight the potential of greenness in the school environment to moderate educational and environmental inequalities in urban areas.

Inconsistent definitions of "urban" result in different conclusions about the size of urban carbon and nitrogen stocks.

There is conflicting evidence about the importance of urban soils and vegetation in regional C budgets that is caused, in part, by inconsistent definitions of "urban" land use. We quantified urban ecosystem contributions to C stocks in the Boston (Massachusetts, USA) Metropolitan Statistical Area (MSA) using several alternative urban definitions. Development altered aboveground and belowground C and N stocks, and the sign and magnitude of these changes varied by land use and development intensity. Aboveground biomass (live trees, dbh > or = 5 cm) for the MSA was 7.2 +/- 0.4 kg C/m2 (mean +/- SE), reflecting a high proportion of forest cover. Vegetation C was highest in forest (11.6 +/- 0.5 kg C/m2), followed by residential (4.6 +/- 0.5 kg C/m2), and then other developed (2.0 +/- 0.4 kg C/m2) land uses. Soil C (0-10 cm depth) followed the same pattern of decreasing C concentration from forest, to residential, to other developed land uses (4.1 +/- 0.1, 4.0 +/- 0.2, and 3.3 +/- 0.2 kg C/m2, respectively). Within a land use type, urban areas (which we defined as > 25% impervious surface area [ISA] within a 1-km(2) moving window) generally contained less vegetation C, but slightly more soil C, than nonurban areas. Soil N concentrations were higher in urban areas than nonurban areas of the same land use type, except for residential areas, which had similarly high soil N concentrations. When we compared our definition of urban to other commonly used urban extents (U.S. Census Bureau, Global Rural-Urban Mapping Project [GRUMP], and the MSA itself), we found that urban soil (1 m depth) and vegetation C stocks spanned a wide range, from 14.4 +/- 0.8 to 54.5 +/- 3.4 Tg C and from 4.2 +/- 0.4 to 27.3 +/- 3.2 Tg C, respectively. Conclusions about the importance of urban soils and vegetation to regional C and N stocks are very sensitive to the definition of urban used by the investigators. Urban areas, regardless of definition, are rapidly expanding in their extent; a systematic understanding of how our development patterns influence ecosystems is necessary to inform future development choices.