Skip to main content
Download PDF

The impact of human and livestock respiration on CO2 emissions from 14 global cities

Carbon Balance and Management volume 17, Article number: 17 (2022) Cite this article

Abstract

Background

The CO2 released by humans and livestock through digestion and decomposition is an important part of the urban carbon cycle, but is rarely considered in studies of city carbon budgets since its annual magnitude is usually much lower than that of fossil fuel emissions within the boundaries of cities. However, human and livestock respiration may be substantial compared to fossil fuel emissions in areas with high population density such as Manhattan or Beijing. High-resolution datasets of CO2 released from respiration also have rarely been reported on a global scale or in cities globally. Here, we estimate the CO2 released by human and livestock respiration at global and city scales and then compare it with the carbon emissions inventory from fossil fuels in 14 cities worldwide.

Results

The results show that the total magnitude of human and livestock respiration emissions is 38.2% of the fossil fuel emissions in Sao Paulo, highest amongst the 14 cities considered here. The proportion is larger than 10% in cities of Delhi, Cape Town and Tokyo. In other cities, it is relatively small with a proportion around 5%. In addition, almost 90% of respiratory carbon comes from urban areas in most of the cities, while up to one-third comes from suburban areas in Beijing on account of the siginificant livestock production.

Conclution

The results suggest that the respiration of human and livestock represents a significant CO2 source in some cities and is nonnegligible for city carbon budget analysis and carbon monitoring.

Background

Currently, approximately 55% of the world’s population [1] occupies only 0.37% of the global land surface [2, 3]. From the viewpoint of the process of urban material metabolism, the total harvested carbon transported into an urban system releases CO2 through the processes of human and livestock metabolism and constitutes a part of the global CO2 cycle (Fig. 1) [4,5,6,7].

Fig. 1
figure 1

The process of urban material metabolism with crop production, the digestion of food and feedstuff by humans and livestock, and finally the release of carbon to the atmosphere via respiration

Full size image

However, CO2 from human and livestock respiration is often neglected due to its perceived small magnitude compared to fossil fuel emissions (FFE) from the burning of fuels for electricity, heating and industrial purposes, other industrial processes and ground transportation within city boundaries [8, 9]. Some researchers have considered human respiration as a significant contribution only in street level or residential areas, where CO2 emissions from power plants and industry can be ruled out [10,11,12]. Others have regarded biogenic emissions to be a nonnegligible contribution at night or even during winter [13,14,15,16] but have not included CO2 release from humans and livestock respiration (HLR).

However, human and livestock respiration could be important in cities with high demands for food and feedstuff consumption. It has been estimated that, human and livestock respiration could represent 1.2 − 30% of FFE in densely populated regions (Table 1), such as some typical cities of Beijing, Greater Paris and Mexico City [17,18,19,20,21,22,23]. When comparing CO2 emissions on regional-scales from both bottom-up inventories and top-down approaches, human and livestock repiration is a significant component in the reconciliation of the differences [23]. Thus, accurately estimating the amount of HLR could improve the results of atmospheric CO2 flux inversion approach for estimating FFE and for comparison with the bottom-up technique [21, 23, 24]. However, the published studies (in Table 1) only focus on individual cities, and the HLR has been clearly calculated in only a few cities. Moreover, only the studies of Zhao et al. and Gurney et al. include livestock respiration [22, 23], while the CO2 release from respiration in other studies is only from humans. Thus, in areas with high population density, detailed estimation of the HLR is necessary for CO2 monitoring and CO2 flux inversion. However, a high-resolution dataset of HLR has rarely been reported on a global scale or in cities globally.

Table 1 Published results for CO2 releases from human and livestock respiration compared with total fossil fuel emissions
Full size table

The purpose of this study is to establish high-resolution datasets of global human and livestock carbon production and to compare with CO2 from FFE within large cities/metropolitan areas around the world. Excluding food loss/waste and livestock feed from local sources, the harvest carbon from crop production should correspond to the total human and livestock consumption of carbon. The fourth section discusses these carbon budgets and presents an uncertainty analysis.

Methods

Study area and in-boundary FFE

The 14 reported global major cities according to the research of Chen et al. [25] are selected as the study area and include Bangkok, Beijing, Shanghai, Delhi, Cape Town, Sao Paulo, Tokyo, Greater Paris, Greater London, Los Angeles, Manhattan, New York City, Washington D.C., and Greater Toronto (see Table 2 and Additional file 1: Figure S1). The definitions of the 14 cities ranges from ‘district’ to ‘metropolitan’ (see also Table 2) [26]. The CO2 from human and livestock respiration is directly emitted within city boundaries, which is belongs to the scope 1 for direct emissions, that is, produced in city boundaries mainly from fossil fuel combustion, transportation, industrial processes and production, land use and waste [27]. To compare the CO2 emissions from human and livestock respiration within these global cities, we retrieved the in-boundary anthropogenic FFE from Chen et al. [25], who estimated the total FFE directly within the city boundaries of these 14 cities and metropolitan areas around the world.

Table 2 The HRh of 6 age-sex groups in 7 global regions (unit: kg C yr− 1)
Full size table

Our study also separated city areas into two subcategories, urban and suburban. The urban extent of each city is based on the 1:10 m urban areas shapefile from Nature Earth (https://www.naturalearthdata.com), which is derived from 2002 to 2003 MODIS satellite data at 1 km resolution [28]. Urban areas are defined as built-up areas with high population densities, high radiance levels in commercial/industrial areas and high-density residential land cover, instead of being based on impervious surfaces [29]. The suburban refers to the remaining area within the administrative boundary of the city that is not included in the built-up area.

Estimate methods for HLR

The CO2 released from respiration of per person (HRh) or per head of livestock (LRl) is obtained according to the basal metabolic rate (BMR). The BMR refers to the minimum level of energy required to sustain vital functions of organs at complete rest in a neutrally temperate environment and in a fasting state. It is measured by heat production or oxygen consumption and can be expressed as Cal m− 2 h− 1, Cal kg− 1 h− 1 or O2 g− 1 h− 1 for individuals [30, 31]. For various mammals, the oxygen consumption rate per body mass consistently decreases with increasing body size, while the rate of oxygen consumption for individuals against body mass tends to decrease along regression lines in logarithmic coordinates (poultry have a similar equation to mammals) [32]. Additionally, oxygen is combined with carbon according to the respiration reaction. Therefore, based on the BMR of each species, we can estimate the CO2 produced by respiration according to the oxygen consumption. What’s more, concerning the metabolic enhancement caused by exercise metabolism and other factors, the physical activity level (PAL) was defined in terms of three levels of physical activity [33]. For simplicity, we assume that the WHO recommended PAL = 1.55 could be used as an uniform parameter for global countries and for different gender and age groups for both human and livestock [34].

The method of HR h

In this study, the BMR is given by 6 age-sex groups and 7 global regions (Additional file 1: Table S1), which was obtained according to the body mass of each age-sex group in each region and the daily BMR predicted by the FAO for different age groups and for both sexes [35]. Then we convert the BMR in the unit of heat production (MJ day− 1) into oxygen consumption (L O2 day− 1, see Additional file 1: Table S2) by introducing the thermal equivalent of oxygen (20.2 kJ L− 1). Finally, the HR in age-sex group h (HRh) is convert as carbon release in kg C yr− 1 (Table 3). The equation of HRh could be:

$$H{R_h} = \frac{\begin{array}{l}BM{R_h} \times Body\_weigh{t_h} \times \frac{{M\left( {{O_2}} \right)}}{{{V_m}}}\\\times \frac{{12}}{{32}} \times 24 \times 365\end{array}}{{{{10}^6}}}$$
(1)
$${Body\_weight}_{h}={BMI}_{h}\times {{Body\_height}_{h}}^{2}$$
(2)
Table 3 The parameters of eight types of livestock
Full size table

where BMRh is the BMR with units of ml O2 g− 1 h− 1; \(M\left({O}_{2}\right)\) is the molecular mass of O2 in g mol− 1; and Vm is the molar volume of gas in 22.4 L mol− 1; The Body_weighth is the body mass estimated as the product of the body mass index (BMI) and the mean height for each age-sex group in each region as formulae (2). The mean body height and mean BMI of each age-sex group in each region are obtained from NCD Risk Factor Collaboration (NCD-RisC) [36, 37] (detail method of the BMR estimation is introduced in Additional file 1:Method). The ratio of carbon (C) and O2 is set to 12/32 according to the processes of respiration, which can be expressed by the following chemical equation [38]:

$${(CH}_{2}O)+{O}_{2}\to {H}_{2}O+{CO}_{2}$$

where (CH20) represents the composition of biological material.

The method of LR l

The BMR (ml O2 g− 1 h− 1) of mammalian livestock and chickens (Additional file 1: Table S4) are measured values from previous experimental results that controlled the environmental temperature, nutrition, age and activity level [39,40,41]. The amount of LR of species l (LRl) represents the total carbon released during the days when animals are alive (Table 4). Therefore, we also assumed that the life span of poultry is 42 days, that the life span of pigs is half a year [42, 43], and that all species except poultry and pigs live for more than one year. The LRl is estimated from the following equation:

$$L{R_l} = \frac{{\left( \begin{array}{l}BM{R_l} \times {\rm{Bod}}y\_weigh{t_l} \times \frac{{M\left( {{O_2}} \right)}}{{{V_m}}}\\\times \frac{{12}}{{32}} \times 24 \times 365\end{array} \right)}}{{{{10}^6}}}$$
(3)
Table 4 Sources of data on humans and livestock
Full size table

where BMRl is the BMR with units of ml O2 g− 1 h− 1; Body_weightl is the average of different breeding ages and genders of each species l; \(M\left({O}_{2}\right)\), Vm and the ratio of carbon (C) and O2 are the same with folume (1) for HRh.

The method and data source of HLR

In this study, the population and livestock production we use are reported as high-resolution datasets (see Sect. Datasets of humans and livestock). The total HLR is the sum of HR and LR, which are estimated by multiplying the CO2 emission of each individual (HRh and LRl) by the total population/livestock production in each grid within city boundaries. We assume that in the same region, the population in each grid have a unified fraction of each age-sex group. Thus, the HR in each grid is actually a weighted sum of the HRh in each age-sex group. The formulaes for HR and LR in grid i are given as:

$$\text{HR}=Population\left(i\right)\times \sum ({f}_{h}\times {HR}_{h})$$
(2)
$$\text{LR}=\sum \left({Livestock}_{l}{\times LR}_{l}\right)$$
(3)

where Population(i) is the total population in grid i; fh is the fraction of age-sex group h in the region where grid i belongs to; Livestockl is the production of species l in grid i. The annual fractions for 6 age-sex groups of the total population in 7 regions come from the World Bank [44].

Datasets of humans and livestock

The HLR in each city are extracted from high-resolution vector datasets (see Table 5). The Gridded Livestock of the World (GLW) datasets include global distributions of eight major livestock species (also see Additional file 1: Table S4). It should be noted that the total cattle and poultry production in Beijing from the high-resolution datasets is 17 times higher than the statistical data from the National Bureau of Statistics of China (NBS, http://data.stats.gov.cn/english/), while cattle production is consistent with census statistics for Shanghai, Delhi and Sao Paulo (Additional file 1: Table S5-S6). As the detailed cattle census statistics of Beijing for GLW were mined from the NBS, we consider the values from official source of China are more reliable. Thus, we first corrected the livestock production in each grid in Beijing according to the spatial distribution from GLW and the total livestock production from NBS.

Table 5 Definition, population and carbon emissions of the 14 cities and metropolitan regions in this study
Full size table

The city boundaries we used in this study come from the database of Global Administrative Areas (GADM) version 2.0 (http://gadm.org/). The shapefile with polygon features of 14 cities was first converted to a high-resolution vector-form dataset at a resolution of 30 arc-second and can be used as a region mask to extract values for population and livestock production within different cities.

To match the years with FFE in 14 cities from Chen et al. [25], we use linear interpolation to obtain the annual human respiration after extracting the total CO2 released in each city every 5 years based on GPWv3 and GPWv4. For livestock, the GLW only provides the high-resolution livestock production in 2010. Thus, we assume that the trend of livestock production in each city is the same as that of the country where the city belongs to (data source see Table 5). The total LR in countries was obtained according to the method described in Sect. 2.2. Then, the total annual livestock respiration within cities were estimated from the annual livestock production of country from FAO and the LR of cities in 2010.

Results

Global CO2 release contributed by human and livestock respiration

We now examine results for the global total CO2 contributed from human and livestock respiration. Globally, the annual total HR was increased from 0.52 Gt C in 1995 to 0.68 Gt C in 2015 (approximately 0.64 Gt C in 2010). In terms of global total LR, it was approximately 0.81 Gt C in 2010, with the majority contributed by cattle, followed by sheep and buffalo, accounting for 60.2%, 10.6% and 8.1% of the total livestock carbon, respectively (Additional file 1: Figure S2(b)).

The spatial distribution of total HLR is shown in Fig. 2. High values of HLR are found in East Asia, southern Asia and Europe and are mainly contributed by humans (Additional file 1: Figure S2(a)). Meanwhile, the high values found in Australia, New Zealand, South America and central North America are mostly related to livestock respiration (Additional file 1: Figure S2(b)). Moreover, the spatial distribution of human respiration shows a consistent spatial pattern with population because it is simply from population multiplied by the BMR.

Fig. 2
figure 2

The spatial distribution of total CO2 released from human and livestock respiration (HLR) (kg C m− 2 yr− 1)

Full size image

Carbon release from human and livestock respiration in 14 cities

Based on the high-resolution gridded data, we extracted the value of HLR in 14 large cities and metropolises. Figure 3 shows the variation of HLR per unit area in the 14 cities. The highest HLR per area in 2014 occurred in Manhattan, followed by Delhi and New York City, with values of 2235.7, 1055.0 and 875.4 g C m− 2 yr− 1, respectively. Greater Toronto featured the lowest HLR per area, with a value as 89.3 g C m− 2 yr− 1.

From the viewpoint of the variation trend, all 14 cities have an increasing trend in HLR over the studied 25 years. The CO2 emission increased by up to 107.8% with a value of 8.7 g C m− 2 yr− 1 over 25 year in Shanghai, while Delhi and Beijng increased by 97.8% and 94.0% with the values of 22.6 and 3.4 g C m− 2 yr− 1, respectively, while the CO2 emission only increased 0.7% in Washington DC with a value of 1.3 g C m− 2 yr− 1. Among the total HLR in 14 cities, human contribute much more than livestock, and the trend of HLR is dominated by the increasing population (see Fig. 3 and S3).

Fig. 3
figure 3

Variation in the total HLR in 14 cities. (a) Greater than 200 g C m− 2 yr− 1 in from 1990, (b) less than 200 g C m− 2 yr− 1 befor 2014

Full size image

For the values of livestock respiration, the differences among cities are noticeable. The total LR in Beijing (502.0 Kt C yr− 1 in 2014) is 4 times those of Delhi (122.0 Kt C yr− 1 in 2014) and Shanghai (147.8 Kt C yr− 1 in 2014) and might be related to cattle farms in the southern and northeastern portions of Beijing. Among the 14 cities, 9 cities had a descending trend from 1990 to 2014, including Greater Paris, Los Angeles, Cape Town, Tokyo, Bangkok, Greater London, Washington DC, New York City and Manhattan (Additional file 1: Figure S3(b)). Part of the reason for this decline is that, with the development of cities, the livestock industry has gradually shifted to the surrounding areas outside of the city.

Figure 4 shows the spatial distribution of HLR in the 14 cities. In most of the cities, the human respiration is more than one order of magnitude larger than that of livestock (Additional file 1: Figures S4 and S5); thus, the spatial distribution of total HLR is dominated by the spatial distribution of humans and is very similar to that of the population (Additional file 1: Figures S1 and S4). In most of the cities, the total livestock respiration in each grid is less than 10 t C yr− 1. The areas around Paris have somewhat higher values, with values of almost 10 t C yr− 1. Moreover, the livestock respiration in eastern Beijing were greater than 60 t C yr− 1 and might be contributed mainly by cattle.

Fig. 4
figure 4

The HLR per in 14 cities and metropolitan regions in 2010 at a resolution of 1 km (ton C yr− 1 per grid cell). The areas within the gray line show the urban extent. Three scale bars are included. Scale bar (1) corresponds to Beijing, Shanghai, Greater Paris, Greater Toronto and Los Angeles; scale bar (2) corresponds to New York City, Delhi, Tokyo, Sao Paulo, Bangkok, Cape Town and Greater London; scale bar (3) corresponds to Washington D.C. and Manhattan

Full size image

Comparison with FFE

Figure 5 shows the HLR compared with the FFE in cities with different areas and populations (the actual value of HLR is shown in Table 2). The HLR amounts to up to 38.2% of the FFE in Sao Paulo, almost all of which is contributed by humans. Delhi has the second highest proportion compared to FFE (nearly 35.6%), with human respiration amounting to almost 32.0% of the FFE. The high proportion to FFE in Delhi and Sao Paulo is mostly due to the relatively low amounts of FFE and their large populations. The contribution of humans and livestock in other cities is approximately 7.5% relative to higher FFE. Otherwise, the dataset of population distribution is an annual average state, which dose not reflect the diurnal variation of population in cities. In large cities, people usually live in the outskirts and commute from the suburbs to commercial areas or industrial parks. That is another reason for the lower ratio of HLR to FFE in cities such as Washington DC with a value of 3.1%.

Fig. 5
figure 5

The proportions of HLR to FFE (in %). The values beside the bar show the sum value of the proportion of human and livestock respiration

Full size image

We also compared the ratio of total respiration to in-boundary FFE associated with other sectors (Fig. 6), including heating and industrial fuels, industrial processes, ground transportation, and large point sources, according to the results from Chen et al. [25]. Based on the average of the 14 cities, the average contribution of human and livestock respiration is comparable to that of large point sources (in-boundary) and is greater than that of industrial processes.

Fig. 6
figure 6

Comparison of the ratio of each sector to in-boundary FFE. The blue point represents the average ratio for the 14 cities, and the gray line is the standard error bar for the 14 cities

Full size image

Moreover, we separated the source of HLR into urban and suburban. In the 14 cities overall, the HLR from urban areas is approximately 90% of the total HLR, and that from suburban areas varies from 0.4 to 32.6% (Additional file 1: Figure S6 and Table S7). In cities with relatively large proportions of HLR compared to FFE such as Delhi and Cape Town, the HLR from suburban accounts for approximately 3.3% and 3.1% of the total in-boundary FFE. As for Sao Paulo, although HLR is 36.2% of total FFE, the HLR from suburban is only 0.4% of total FFE. Furthermore, the HLR from suburban areas in Beijing accounts for up to one-third of the total HLR (2.5% of total FFE). Such a relatively high proportion is mainly contributed by the significant livestock production in suburban area, wchich accounts for 80% of the total area. Therefore, in studies with high-precision CO2 concentration measurements, the stations located in suburban areas of cities such as Beijing, Delhi and Cape Town should also be taken into consideration [69].

Discussion

Carbon monitoring in urban areas has shown that the partern of human metabolism can partially explain the diurnal pattern of CO2 flux, as well as differences in CO2 flux between working days and non-working days, in densely populated urban areas [58, 59]. The research of Ciais et al. [24] shows that, in total, humans and livestock contribute 5% to FFE in urban areas globally, and the CO2 releases from respiration is even larger than oil burning in India and larger than FFE from gas in Chinese cities. Also, to explain the difference in the central estimate of bottom-up and top-down approaches must consider all flux contributions including the important contribution of animal respiration[23]. Thus, the HLR is not a negligible contribution in some populated cities.

Uncertainty associated with the parameters

We compared the HLR with in-boundary FFE via a bottom-up method. The uncertainty in this study mainly comes from the uncertainty of the data source of population and livestock production, as well as the uncertainty in the assumptions of the parameters of BMR and PAL.

The BMR is affected by a variety of factors, such as age, gender, exercise, body temperature, nutritional status, or lactation [31, 60,61,62,63]. Although several equations have been developed to predict BMR [31, 40, 64] and often take into account weight, height, age, gender and other factors, considering the high-resolution data of population and livestock production used in this study, it is difficult to distinguish the above factors on this spatial resolution. Therefore, we assume that for each livestock species (mammalian livestock, poultry) individuals of different ages and genders and in different areas are assigned the same value of BMR. As for humans, we distinguished the BMR for each age-sex group in 7 global regions. In future research on more detailed CO2 emissions for specific urban areas, we will consider the physical differences of people in different countries and the natural environment.

In fact, BMR is only the lowest estimate of metabolic rate, and the actual respiratory intensity of each individual is related to climate and daily activity. The PAL recommended by the FAO/WHO/UNU Expert Committee in 1985 was defined in three levels: the minimum was set at 1.55 and 1.56 BMR for men and women, while the highest was defined as 2.10 and 1.82 BMR for men and women, respectively [33]. Research also showed that the PAL were not significantly different between age groups[65]. Due to lack of data, we roughly assume that the average daily activity is at the lowest level since most people do not have long-term high-intensity activity.

The carbon emissions from respiration per individual in this study are only a first-order approximation based on the above assumption. Through the weighted average of the HRh for each age-sex group in 7 regions, the global averaged HRh is assumed to be 89.90 kg C yr− 1 in this study. This result is comparable with that in other studies, which is varies over a range of 52.9–160 kg C yr− 1 (Additional file 1: Table S8) based on different methods or without considering different age groups and gender [7, 17, 19, 21,22,23, 66]. The value used by Huang et al. is much higher than other studies. As a result, their estimated annual averaged global carbon emissions from human respiration is approximately 1.2 Gt C yr− 1 [66] from 1990 to 2005, which is much greater than the value estimated in this study (0.54 Gt C yr− 1). Additionally, their estimation of livestock respiration is approximately 0.8 Gt C yr− 1 [66], which is comparable with this study (0.62 Gt C yr− 1).

The HLR within other 118 world urban areas

In addition to the 14 selected cities, we provide a supplementary dataset on the HLR in 2010 within other 118 world urban areas larger than 100 km2, which were selected on the basis of being capitals or areas with large populations (greater than 1 million) or that have carbon monitoring sites, as a reference for carbon monitoring and carbon emission research (Additional file 2: Table S9). Among the 118 selected urban areas, approximately 11 urban areas have the proportions of HLR to FFE greater than 50%, and 50 urban areas have the proportions greater than 10% (Additional file 1: Figure S7 and Additional file 2: Table S9). We can see that the cities where HLR to FFE is greater than 50% are almost all the cities in countries of low income in Sub-Saharan Africa ( see Additional file 2: Table S9). In these cities, the main reason of the high ratio of HLR to FFE is the lower levels of economic development and lower fossil fuel emissions. But for the cities in high income and upper middle income countries, the ratio of HLR to FFE varies from 0.2 − 22.5% with the mean and medians of 8.2% and 6.4%, respectively. There does not appear to be a clear correlation between population density/GDP per capital and the ratio of HLR to FFE. In cities with higher population density (greater than 10,000 per km2), regions with higher GDP per capital tend to have lower ratio of HLR to FFE, but this is not absolutely the case. For example, Istanbul, Turkey and Bangalore, India have similar population density. In 2010, GDP per capita in Turkey (10,742 US$ in 2010) was much higher than that of India (1358 US$ in 2010), and the ratio of HLR to FFE in Istanbul (20.3%) was also higher than that of Bangalore (8.5%). In this part, due to the lack of the statistics on the GDP per capita of each city, our analysis is based on the GDP per capita of the country where the city is located. Even so, we can see that the ratio of HLR to FFE is caused by factors such as population density, urban development level, and perhaps urban cleanliness. This study provides only a global overview through first-order approximation. It is not enough to judge the ratio of HLR to FFE simply by the population density or GDP per capital of a city, and to determine whether the impact of HLR needs to be considered in the city’s carbon monitoring.

The budget of harvested crop carbon

The processes that release CO2 (human and livestock respiration) and the processes of crop carbon harvest result in fluxes of CO2 to and from the atmosphere and constitute a part of the global CO2 cycle. In this section, we will briefly discuss the budget of the carbon cycle from the bottom-up approach and top-down approach.

Based on the results in the previous sections, humans and livestock together released 1.45 Gt C in 2010 (0.64 and 0.81 Gt C, respectively). Forage grass is a local feed resource of horses, sheep and goats, as well as a portion of the feed for poultry, pigs and housed cattle in developing countries [18]. We follow the assumption of Ciais et al. that only 20% of poultry, pigs and cattle in developing countries received grain-based feeds, whereas 100% of these livestock in developed countries received grain-based feeds. Other livestock received local feed resources in our assumption. Based on the livestock production in each country from the FAO [55], we find that the carbon emissions from livestock utilizing local feed sources was approximately 0.64 Gt C (approximately 79% of total livestock respiration) in 2010. Hence, the carbon released by livestock respiration comes from grain-based feeds is approximately 0.17 Gt C in 2010 (about 21% of total livestock respiration).

Here, we also considered food waste. According to the FAO, approximately one-third of food produced by humans is lost or wasted globally every year [67]. Crop production is estimated to have been approximately 1.50 Gt C in 2010 based on crop production from the FAO [68], which is slightly larger than the human and livestock consumption estimated in this study.

The budget of harvested crop carbon is constructed as follows:

$$P-W={R}_{human}+{R}_{livestock}\times 0.21$$
(4)

where P is crop carbon production via the top-down modeling method, W is the lost and wasted crop production, Rhuman and Rlivestock are the carbon releases from human and livestock respiration estimated by the bottom-up approach, and only 0.21 of Rlivestock comes from grain-based feeds. Based on the above discussion, we estimate that W is approximately 0.5 Gt C. Excluding food loss and waste, crop production used for metabolization is approximately 1.00 Gt C, which has the same magnitude and is comparable to human and livestock consumption (0.81 Gt C, excluding livestock feeding on local forage grass).

In addition, we should note that the carbon released by livestock respiration in 14 cities also includes local feed sources. In particular, in Delhi and Beijing, the percentages of livestock respiratory carbon are 3.4% and 1.5% relative to fossil fuel. Considering livestock fed with local feed, livestock respiration should have less influence on city carbon emissions than our estimation with in cities.

Conclusion

In this study, we used the global population and livestock production, as well as the parameter of BMR, to estimate the amount of carbon released by human and livestock respiration from 1990 to 2014. Then, we calculated the carbon emissions from human and livestock respiration in 14 of the world’s largest cities and compared them with the in-boundary FFE. The results showed that the proportion of total carbon released from humans and livestock is approximately 5–10% relative to in- boundary FFE in most of the 14 typical cities. In cities such as Delhi, Sao Paulo and Cape Town, humans and livestock contribute up to 38.2% relative to FFE. In studies monitoring FFE from ground stations or satellites, neglicting human/livestock emissions could overestimate the in-boundary FFE [21] or leave ambiguity in the evaluation of the FFE trend [24]. In addition, approximately 90% of the respiratory carbon is released in the urban areas of most cities, while the suburban HLR has a noteable contribution compared to FFE in Beijing, Delhi and Cape Town. This means that the setup of carbon monitoring sites should not neglect suburban areas in those cities. Further more, the results in suburban areas also helps to analyse the vertical distribution of CO2 in the boundary layer, and provide data for validating transport models[69]. We should also note that this study used global unified values for BMR for livestock and PAL for all species. To estimate the carbon from human and livestock respiration more accurately and to compare it with FFE, it is better to adopt parameters appropriate for regional or national climate and livestock feeding conditions.

Data availability

The datasets supporting the conclusions of this article are included within the article and its additional files. High-resolution data are available upon request to corresponding author.

References

  1. The World Bank. Urban population [Internet]. https://www.worldbank.org/. 2019. Available from: https://data.worldbank.org/indicator/SP.URB.TOTL.IN.ZS.

  2. Klein Goldewijk K, Beusen A, Janssen P. Long-term dynamic modeling of global population and built-up area in a spatially explicit way: HYDE 3.1. The Holocene. 2010;20:565–73.

    Article  Google Scholar 

  3. Klein Goldewijk K, Beusen A, Doelman J, Stehfest E. Anthropogenic land use estimates for the Holocene – HYDE 3.2. Earth Syst Sci Data. 2017;9:927–53.

    Article  Google Scholar 

  4. Pataki DE, Alig RJ, Fung AS, Golubiewski NE, Kennedy CA, Mcpherson EG, et al. Urban ecosystems and the North American carbon cycle. Glob Change Biol. 2006;12:2092–102.

    Article  Google Scholar 

  5. Luo T, Ouyang Z, Frostick LE. Food carbon consumption in Beijing urban households. Int J Sustainable Dev World Ecol. 2008;15:189–97.

    Article  Google Scholar 

  6. Churkina G. Modeling the carbon cycle of urban systems. Ecol Model. 2008;216:107–13.

    Article  Google Scholar 

  7. Churkina G, Brown DG, Keoleian G. Carbon stored in human settlements: the conterminous United States: Carbon in human settlements. Glob Change Biol. 2010;16:135–43.

    Article  Google Scholar 

  8. Helfter C, Famulari D, Phillips GJ, Barlow JF, Wood CR, Grimmond CSB, et al. Controls of carbon dioxide concentrations and fluxes above central London. Atmos Chem Phys. 2011;11:1913–28.

    Article  CAS  Google Scholar 

  9. Lietzke B, Vogt R. Variability of CO2 concentrations and fluxes in and above an urban street canyon. Atmos Environ. 2013;74:60–72.

    Article  CAS  Google Scholar 

  10. Kellett R, Christen A, Coops NC, van der Laan M, Crawford B, Tooke TR, et al. A systems approach to carbon cycling and emissions modeling at an urban neighborhood scale. Landsc Urban Plann. 2013;110:48–58.

    Article  Google Scholar 

  11. Moriwaki R, Kanda M. Seasonal and diurnal fluxes of radiation, heat, water vapor, and carbon dioxide over a suburban area. J Appl Meteorol. 2004;43:1700–10.

    Article  Google Scholar 

  12. Widory D, Javoy M. The carbon isotope composition of atmospheric CO2 in Paris. Earth Planet Sci Lett. 2003;215:289–98.

    Article  CAS  Google Scholar 

  13. Cecchi G, Wint W, Shaw A, Marletta A, Mattioli R, Robinson T. Geographic distribution and environmental characterization of livestock production systems in Eastern Africa. Agric Ecosyst Environ. 2010;135:98–110.

    Article  Google Scholar 

  14. Henninger S, Kuttler W. Near surface carbon dioxide within the urban area of Essen, Germany. 35. Parts A/B/C: Physics and Chemistry of the Earth; 2010. pp. 76–84.

    Google Scholar 

  15. Lauvaux T, Miles NL, Deng A, Richardson SJ, Cambaliza MO, Davis KJ, et al. High-resolution atmospheric inversion of urban CO2 emissions during the dormant season of the Indianapolis Flux Experiment (INFLUX): URBAN INVERSION. J Geophys Research: Atmos. 2016;121:5213–36.

    Article  CAS  Google Scholar 

  16. Wada R, Pearce JK, Nakayama T, Matsumi Y, Hiyama T, Inoue G, et al. Observation of carbon and oxygen isotopic compositions of CO2 at an urban site in Nagoya using Mid-IR laser absorption spectroscopy. Atmos Environ. 2011;45:1168–74.

    Article  CAS  Google Scholar 

  17. Koerner B, Klopatek J. Anthropogenic and natural CO2 emission sources in an arid urban environment. Environ Pollut. 2002;116:45–51.

    Article  Google Scholar 

  18. Ciais P, Bousquet P, Freibauer A, Naegler T. Horizontal displacement of carbon associated with agriculture and its impacts on atmospheric CO2: AGRICULTURE CARBON DISPLACEMENT. Glob Biogeochem Cycles. 2007;21:1–12.

    Article  Google Scholar 

  19. West TO, Marland G, Singh N, Bhaduri BL, Roddy AB. The human carbon budget: an estimate of the spatial distribution of metabolic carbon consumption and release in the United States. Biogeochemistry. 2009;94:29–41.

    Article  CAS  Google Scholar 

  20. Velasco E, Roth M. Cities as net sources of CO2: review of atmospheric CO2 exchange in urban environments measured by eddy covariance technique: urban CO2 flux measurements by eddy covariance. Geogr Compass. 2010;4:1238–59.

    Article  Google Scholar 

  21. Bréon FM, Broquet G, Puygrenier V, Chevallier F, Xueref-Rémy I, Ramonet M, et al. An attempt at estimating Paris area CO2 emissions from atmospheric concentration measurements. Atmospheric Chem Phys Discuss. 2014;14:9647–703.

    Google Scholar 

  22. Zhao R, Huang X, Liu Y, Zhong T, Ding M, Chuai X. Urban carbon footprint and carbon cycle pressure: The case study of Nanjing. J Geog Sci. 2014;24:159–76.

    Article  Google Scholar 

  23. Gurney KR, Liang J, Patarasuk R, O’Keeffe D, Huang J, Hutchins M, et al. Reconciling the differences between a bottom-up and inverse-estimated FFCO2 emissions estimate in a large US urban area. Elem Sci Anth. 2017;5:44.

    Article  Google Scholar 

  24. Ciais P, Wang Y, Andrew R, Bréon FM, Chevallier F, Broquet G, et al. Biofuel burning and human respiration bias on satellite estimates of fossil fuel CO2 emissions. Environ Res Lett. 2020;15:074036.

    Article  CAS  Google Scholar 

  25. Chen J, Zhao F, Zeng N, Oda T. Comparing a global high-resolution downscaled fossil fuel CO2 emission dataset to local inventory-based estimates over 14 global cities. Carbon Balance and Management [Internet]. 2020 [cited 2020 Jun 15];15. Available from: https://cbmjournal.biomedcentral.com/articles/https://doi.org/10.1186/s13021-020-00146-3.

  26. United Nations Department of Economic and Social Affairs. The World’s Cities in 2018 [Internet]. UN; 2016 [cited 2020 Feb 9]. Available from: https://www.un-ilibrary.org/population-and-demography/the-world-s-cities-in-2016_8519891f-en.

  27. Kennedy C, Steinberger J, Gasson B, Hansen Y, Hillman T, Havránek M, et al. Methodology for inventorying greenhouse gas emissions from global cities. Energy Policy. 2010;38:4828–37.

    Article  CAS  Google Scholar 

  28. Schneider A, Friedl MA, McIver DK, Woodcock CE. Mapping urban areas by fusing multiple sources of coarse resolution remotely sensed data. Photogram Eng Remote Sens. 2003;69:1377–86.

    Article  Google Scholar 

  29. Schneider A, Friedl MA, Potere D. A new map of global urban extent from MODIS satellite data. Environ Res Lett. 2009;4:044003.

    Article  Google Scholar 

  30. Gellman MD, Turner JR, editors. Encyclopedia of Behavioral Medicine [Internet]. New York, NY: Springer New York; 2013 [cited 2019 Nov 26]. Available from: http://link.springer.com/https://doi.org/10.1007/978-1-4419-1005-9.

  31. Henry C. Mechanisms of changes in basal metabolism during ageing. Eur J Clin Nutr. 2000;54:77–91.

    Article  Google Scholar 

  32. Schmidt-Nielsen K. Animal Physiology: Adaptation and Environment [Internet]. Second. Cambridge University Press; 1979 [cited 2020 Feb 16]. Available from: https://www.jstor.org/stable/2403350?origin=crossref.

  33. Alfonzo-González G, Doucet E, Alméras N, Bouchard C, Tremblay A. Estimation of daily energy needs with the FAO/WHO/UNU 1985 procedures in adults: comparison to whole-body indirect calorimetry measurements. Eur J Clin Nutr. 2004;58:1125–31.

    Article  Google Scholar 

  34. Black A. The sensitivity and specificity of the Goldberg cut-off for EI:BMR for identifying diet reports of poor validity. Eur J Clin Nutr. 2000;54:359–404.

    Article  Google Scholar 

  35. J.V.G.A.Durnin. Basal metabolic rate in man [Internet]. Rome: FAO; 1981 Oct. Available from: http://www.fao.org/3/M2845E/M2845E00.htm.

  36. Abarca-Gómez L, Abdeen ZA, Hamid ZA, Abu-Rmeileh NM, Acosta-Cazares B, Acuin C, et al. Worldwide trends in body-mass index, underweight, overweight, and obesity from 1975 to 2016: a pooled analysis of 2416 population-based measurement studies in 128·9 million children, adolescents, and adults. The Lancet. 2017;390:2627–42.

    Article  Google Scholar 

  37. Rodriguez-Martinez A, Zhou B, Sophiea MK, Bentham J, Paciorek CJ, Iurilli ML, et al. Height and body-mass index trajectories of school-aged children and adolescents from 1985 to 2019 in 200 countries and territories: a pooled analysis of 2181 population-based studies with 65 million participants. The Lancet. 2020;396:1511–24.

    Article  Google Scholar 

  38. Keeling RF. Measuring correlations between atmospheric oxygen and carbon dioxide mole fractions: A preliminary study in urban air. J Atmos Chem. 1988;7:153–76.

    Article  CAS  Google Scholar 

  39. Elgar MA, Harvey PH. Basal Metabolic Rates in Mammals: Allometry, Phylogeny and Ecology. Funct Ecol. 1987;1:25.

    Article  Google Scholar 

  40. Meltzer A. Thermoneutral zone and resting metabolic rate of broilers. Br Poult Sci. 1983;24:471–6.

    Article  CAS  Google Scholar 

  41. Freeman BM. Gaseous metabolism of the domestic chicken: IV. The effect of temperature on the resting metabolism of the fowl during the first month of life. Br Poult Sci. 1963;4:275–8.

    Article  Google Scholar 

  42. Bruce Hoar JA. Production Cycle of Swine [Internet]. Western Institute for Food Safety and Security at the University of California; 2015. Available from: https://www.wifss.ucdavis.edu/wp-content/uploads/2015/FDA/feed/animalclass_swine_FINAL.pdf.

  43. Food and Agriculture Organization of the United Nations, editor. Poultry development review [Internet]. 2013. Available from: http://www.fao.org/3/i3531e/i3531e.pdf.

  44. The World Bank: Featured indicators of Health. [Internet]. 2019. Available from: https://data.worldbank.org/indicator/.

  45. Center for International Earth Science Information Network - CIESIN - Columbia University, Centro Internacional de Agricultura Tropical - CIAT. Gridded Population of the World, Version 3 (GPWv3): Population Density Grid [Internet]. NASA Socioeconomic Data and Applications Center (SEDAC); 2005. Available from: https://doi.org/10.7927/H4XK8CG2.

  46. Center for International Earth Science Information Network - CIESIN - Columbia University. Gridded Population of the World, Version 4 (GPWv4): Population Count Adjusted to Match 2015 Revision of UN WPP Country Totals, Revision 11 [Internet]. NASA Socioeconomic Data and Applications Center (SEDAC); 2018. Available from: https://doi.org/10.7927/H4PN93PB.

  47. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint WGR, et al. Global sheep distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:10.7910/DVN/BLWPZN.

  48. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint WGR, et al. Global pigs distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:https://doi.org/10.7910/DVN/33N0JG.

  49. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint WGR, et al. Global horses distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:https://doi.org/10.7910/DVN/7Q52MV.

  50. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint WGR, et al. Global goats distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:https://doi.org/10.7910/DVN/OCPH42.

  51. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint WGR, et al. Global ducks distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:10.7910/DVN/ICHCBH.

  52. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint WGR, et al. Global cattle distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:https://doi.org/10.7910/DVN/GIVQ75.

  53. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint GRW, et al. Global chickens distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:10.7910/DVN/SUFASB.

  54. Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke S, Wint WGR, et al. Global buffaloes distribution in 2010 (5 minutes of arc) [Internet]. Harvard Dataverse; 2018 [cited 2020 Feb 11]. Available from: https://dataverse.harvard.edu/citation?persistentId=doi:https://doi.org/10.7910/DVN/5U8MWI.

  55. FAO. FAOSTAT. Production Live animals. (Accessed: 2019-08-15). [Internet]. 2019. Available from: http://www.fao.org/faostat/en/#home.

  56. Kennedy CA, Ramaswami A, Carney S, Dhakal S. Greenhouse Gas Emission Baselines for Global Cities and Metropolitan Regions. Cities and Climate Change [Internet]. The World Bank; 2011 [cited 2020 Feb 7]. p. 15–54. Available from: http://elibrary.worldbank.org/doi/abs/https://doi.org/10.1596/9780821384930_CH02.

  57. Kennedy CA, Ibrahim N, Hoornweg D. Low-carbon infrastructure strategies for cities. Nat Clim Change. 2014;4:343–6.

    Article  CAS  Google Scholar 

  58. Grimmond CSB. Flux and turbulence measurements at a densely built-up site in Marseille: Heat, mass (water and carbon dioxide), and momentum. Journal of Geophysical Research [Internet]. 2004 [cited 2020 Apr 28];109. Available from: https://doi.org/10.1029/2004JD004936.

  59. Ward HC, Kotthaus S, Grimmond CSB, Bjorkegren A, Wilkinson M, Morrison WTJ, et al. Effects of urban density on carbon dioxide exchanges: Observations of dense urban, suburban and woodland areas of southern England. Environ Pollut. 2015;198:186–200.

    Article  CAS  Google Scholar 

  60. Lafortuna CL, Proietti M, Agosti F, Sartorio A. The energy cost of cycling in young obese women. Eur J Appl Physiol. 2006;97:16–25.

    Article  Google Scholar 

  61. M’Kaouar H, Péronnet F, Massicotte D, Lavoie C. Gender difference in the metabolic response to prolonged exercise with [13 C]glucose ingestion. Eur J Appl Physiol. 2004;92:462–9.

    Article  Google Scholar 

  62. Owen OE, Holup JL, D’Alessio DA, Craig ES, Polansky M, Smalley KJ, et al. A reappraisal of the caloric requirements of men. Am J Clin Nutr. 1987;46:875–85.

    Article  CAS  Google Scholar 

  63. Speakman JR, Król E, Johnson MS. The functional significance of individual variation in basal metabolic rate. Physiol Biochem Zool. 2004;77:900–15.

    Article  Google Scholar 

  64. Hayssen V, Lacy RC. Basal metabolic rates in mammals: Taxonomic differences in the allometry of BMR and body mass. Comp Biochem Physiol Part A: Physiol. 1985;81:741–54.

    Article  CAS  Google Scholar 

  65. Piers LS, Soares MJ, McCormack LM, O’Dea K. Is there evidence for an age-related reduction in metabolic rate? J Appl Physiol. 1998;85:2196–204.

    Article  CAS  Google Scholar 

  66. Huang J, Huang J, Liu X, Li C, Ding L, Yu H. The global oxygen budget and its future projection. Sci Bull. 2018;63:1180–6.

    Article  CAS  Google Scholar 

  67. Gustavsson J, Cederberg C, Sonesson U. Global food losses and food waste: extent, causes and prevention; study conducted for the International Congress Save Food! at Interpack 2011, [16–17 May], Düsseldorf, Germany. Rome: Food and Agriculture Organization of the United Nations; 2011.

    Google Scholar 

  68. Han P, Zeng N, Zhao F, Lin X. Estimating global cropland production from 1961 to 2010. Earth Sys Dyn. 2017;8:875–87.

    Article  Google Scholar 

  69. Bao Z, Han P, Zeng N, Liu D, Cai Q, Wang Y, et al. Observation and modeling of vertical carbon dioxide distribution in a heavily polluted suburban environment. Atmospheric and Oceanic Science Letters. 2020;13:371–9.

    Article  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Key R&D Program of China (no. 2017YFB0504000) and the National Natural Science Foundation of China (no. 41975050).

Funding

Funding was provided by the National Key R&D Program of China (No. 2017YFB0504000).

Author information

Authors and Affiliations

  1. State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China

    Qixiang Cai, Pengfei Han & Di Liu

  2. Department of Atmospheric and Oceanic Science, and Earth System Science Interdisciplinary Center, University of Maryland, College Park, MD, USA

    Ning Zeng

  3. Key Laboratory of Geographic Information Science (Ministry of Education), School of Geographic Sciences, East China Normal University, Shanghai, China

    Fang Zhao

  4. State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry, Institute of Atmospheric Physics, Chinese Academ y of Sciences, Beijing, China

    Xiaohui Lin

  5. State Key Laboratory of Hydrology- Water Resources and Hydraulic Engineering, Hohai University, Nanjing, China

    Jingwen Chen

Authors
  1. Qixiang Cai
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ning Zeng
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Fang Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Pengfei Han
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Di Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xiaohui Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jingwen Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

QXC, PFH and NZ conceived and designed the study. QXC and PFH collected and analyzed the datasets. QXC, PFH, and FZ led the paper writing with contributions from all coauthors. All authors made suggestions and contributed to the revision of the final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Qixiang Cai.

Ethics declarations

Competing interests

The authors declare that they have no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

13021_2022_217_MOESM1_ESM.xlsx

Supplementary Material 1: Table S9 The HLR within 118 worldurban areas. Notes: The 118 world urban areas are selected onthe basis of area larger than 100 km2, being capitals, with populationsgreater than 1 million or have carbon monitoring sites.

Supplementary Material 2: Supplementary tables and figures.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cai, Q., Zeng, N., Zhao, F. et al. The impact of human and livestock respiration on CO2 emissions from 14 global cities. Carbon Balance Manage 17, 17 (2022). https://doi.org/10.1186/s13021-022-00217-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13021-022-00217-7

Keywords

Carbon Balance and Management

ISSN: 1750-0680

Contact us