Rice paddy soils are a quantitatively important carbon store according to a global synthesis

Rice paddies account for ~9% or the world’s cropland area and are characterized by environmental conditions promoting soil organic carbon storage, methane emissions and to a lesser extent nitrous oxide emissions. Here, we synthesize data from 612 sites across 51 countries to estimate global carbon stocks in paddy soils and determine the main factors affecting paddy soil carbon storage. Paddy soils (0–100 cm) contain 18 Pg carbon worldwide. Paddy soil carbon stocks decrease with increasing mean annual temperature and soil pH, whereas mean annual precipitation and clay content had minor impacts. Meta-analysis shows that paddy soil carbon stocks can be increased through several management practices. However, greenhouse gas mitigation through paddy soil carbon storage is generally outweighed by increases in methane and nitrous oxide emissions. Our results emphasize the key role of paddies in the global carbon cycle, and the importance of paddy management in minimizing anthropogenic greenhouse gas emissions. Rice paddies represent approximately 1.2% of the global soil organic carbon pool and contain 20% more soil organic carbon per hectare than croplands on average, according to a global synthesis.

S oils contain the largest reservoir of terrestrial organic carbon (C) and they are a main natural source of atmospheric carbon dioxide (CO 2 ) 1 . Soil organic carbon (SOC) is widely recognized as a key element of soil fertility, and croplands with high SOC contents have better structure and lower risks of erosion 2 . Insights into the global distribution of SOC stocks and the effects of environmental variables will thus improve estimates of C-climate feedbacks, and may contribute to agricultural policies designed to improve soil quality [3][4][5] . Over the last decade, SOC stocks have been increasingly estimated at global and regional scales for numerous ecosystems, including croplands 6 , grasslands 7 , wetlands 8,9 , and forests 10,11 . However, even though rice paddies cover~9% of the global cropland area and provide staple food for roughly half the world's population 12 , a global assessment of SOC stocks in rice paddies is still lacking.
Paddy soils are anthropogenic soils (Anthrosols) for cultivation of rice, which are intentionally flooded and puddled, i.e., tilled under water saturated conditions. Paddy soils are widely distributed from temperate to tropical climates on all continents, but mainly in Asia. Rice paddies can be established on various natural and previously agriculturally used soil types, and on various parent materials, but are highly modified by management practices during rice-paddy cultivation 13 . Because rice paddies are frequently flooded and puddled, their properties differ substantially from those of all other arable upland soil. Anaerobic conditions induced by flooding slow down organic matter decomposition, and thus beneficial to SOC accumulation 14 . At the same time, these anaerobic conditions promote CH 4 production by methanogens, making rice paddies a main source of anthropogenic CH 4 emissions 15 .
The development of efficient irrigation techniques led to expansion of the global paddy area by >30% since the 1960s 16 . During this period, rising levels of mineral fertilizer application and subsequent increased straw return to soil stimulated SOC storage in paddy soils around the world 17,18 . For example, the topsoil layer (0-30 cm) of rice paddies in China store~30% more SOC (45 Mg ha −1 ) than corresponding upland soils (35 Mg ha −1 ) 19 . Therefore, changes in the C pool size of paddies could strongly affect atmospheric CO 2 concentrations. However, the size of global rice-paddy SOC pool is still unclear. On a global scale, SOC stocks in upland soils increase with precipitation and clay content and decrease with temperature 20 , but the main environmental and management factors affecting paddy SOC stocks at different climates have not yet been determined. This information could help to optimize agronomical management designed to enhance SOC sequestration, inform agricultural policy measures designed to improve soil quality, and predict the potential impacts of climate changes on SOC stocks.
Several recent studies have reported paddy SOC stocks in regions that were previously underrepresented in rice-paddy research, such as South America and Africa (e.g., ref. 21 ). With SOC inventories now being available for most of the world's ricegrowing areas, a data synthesis may reduce the uncertainty regarding paddy soil C stocks and identify practices and areas with high potential for soil C storage. We thus conducted a global synthesis of SOC stocks in the topsoil (0-30 cm) and subsoil (30-100 cm) of rice paddies, including data from 612 sites around the world ( Fig. 1a; see Methods and Supplementary Data 1). Our objectives were (1) to determine climatic factors, soil properties, and management practices that affect SOC stocks of paddy topsoils on a global scale; (2) to compare paddy SOC storage between the main rice-producing countries and their contribution to the global paddy SOC pool; and (3) to determine the contribution of SOC storage in paddy soils to the global terrestrial and agricultural SOC pool.
We found that paddy soils (0-100 cm) contain 18 Pg SOC worldwide,~1.2% of the global SOC pool, corresponding to 14% of the total SOC pool in croplands. Paddy SOC stocks decrease with increasing mean annual temperature and soil pH, but mean annual precipitation and clay content had minor impacts. Metaanalysis further indicates that paddy SOC stocks (0-30 cm) increase with fertilization (9-32%), straw return (13%), and conservation tillage (8-10%). However, climate benefits of SOC storage in paddies are generally negated by increases in CH 4 and N 2 O emissions.

Results
Our database included information about rice paddies between 48°N and 38°S and between 147°E and 90°W. The distribution of sites was skewed towards low elevations, with most sites located below 200 m a.s.l. (Fig. 1b). The SOC content in the topsoil of most sites (>70%) ranged from 7 to 16 g kg −1 , with a mean of 13.8 g kg −1 . The bulk density (BD) of the topsoil at most sites (>70%) ranged between 1.2 and 1.6 g cm −3 , with a mean of 1.3 g cm −3 ( Supplementary Fig. 1).
The estimated global average SOC stock of rice paddies is 108 Mg ha −1 for the 0-100 cm layer,~10% higher than the global average for all soils (Table 1). Average SOC stocks in rice paddies are lower than for mangroves, forests, and wetlands, but substantially higher than for grasslands and croplands (Table 1). Totaled across the globe, the upper 1 m of paddy soils contains 18 Pg (95% CI: 17.2-18.9) organic C. This amounts to~1.2% of the global SOC pool, or~14.2% of the total SOC pool in croplands worldwide (Table 2).
Topsoil paddy SOC stocks ranged between 7 and 330 Mg ha −1 (Fig. 2a). Mean SOC stocks increased with latitude (p < 0.01), from 50 Mg ha −1 in the tropics to 62 Mg ha −1 in temperate regions ( Supplementary Fig. 2). Topsoil SOC stocks in rice paddies differed more than three fold between main rice-producing countries ( Fig. 2b): paddies in Indonesia and Vietnam had the highest SOC stocks (~78 Mg ha −1 ), whereas paddies in Pakistan, Cambodia, Africa, and Central and South America contained less than 30 Mg C ha −1 . Paddies in China, India, and Indonesia together accounted for~56% of the global paddy SOC pool (Fig. 2c).
Correlation analyses indicated that paddy SOC stocks are mainly determined by soil pH and mean annual temperature (MAT), and to a much smaller extent by mean annual precipitation (MAP) and clay content (Table 3). SOC stocks decreased with increases in pH and MAT, and slightly increased with increasing MAP and clay content.
Our meta-analysis indicates that N fertilization increased SOC stocks by 9% on average, whereas combined NPK application doubled the increase in SOC stocks compared to sole N fertilization application (Fig. 3). Organic fertilizer application alone and combined with NPK increased SOC stocks by 19% and 32%, respectively. Returning straw to the soil increased C stocks by 13%. Compared to conventional tillage practices, no-till and reduced tillage increased SOC stocks by 10% and 8%, respectively.

Discussion
Our data synthesis and meta-analysis reveal the importance of rice paddies for the global C cycle. Per unit area, paddy soils contain more SOC than upland agricultural soils 14,22 . Whereas rice paddies occupy less than 9% of the global cropland area, they harbor more than 14% of its SOC stocks ( Table 2). These large SOC stores can be explained by anaerobic conditions of rice paddies after flooding, slowing down decomposition rates and thus, increasing soil C accumulation compared to other cropland types 13,23-25 .
SOC stocks were best predicted by soil pH (Table 3). This supports previous reports that soil acidity strongly affects ecosystem C balances (e.g., ref. 26 ). pH regulates several soil properties and processes that play key roles in determining C stocks. For example, the solubility of organic matter decreases under low pH by formation of organic matter complex with polyvalent metal ions such as iron and aluminum 26,27 . Consequently, leaching of dissolved organic matter will be reduced. Furthermore, low soil pH values slow down litter decomposition by reducing enzyme activity 28,29 and by changing the composition of microbial communities 30,31 .
The MAP was less important than MAT in determining C stocks because of the regular flooding of paddy fields. High temperatures and rainfall in the tropics typically stimulate plant productivity 32 , but accelerated SOM decomposition negate or even override the effects of increased C inputs from plant production on C stocks 33 . Similarly, the slowing down of SOC decomposition rates with decreasing temperature 34 explains the increase in C stocks with latitude.
Since paddy soils can develop from different parent materials, their initial mineralogy, texture, and fertility can vary considerably 13 . However, prolonged rice cultivation masks initial soil characteristics and minimizes the influence of parent material on pedogenic features (e.g., ref. 35 ). This likely explains why clay content explained less of the variation in SOC stock than did other environmental factors (e.g., pH and MAT, Table 3). Upland soils with high clay contents generally store more C than sandy soils 36 , because clay minerals provide binding surfaces for organic matter and creates anoxic microsites within aggregates. Soil aggregation has minor impact on C dynamics by regular puddling 35 , and strongly variable redox conditions may reduce the formation and stability of organic matter-clay complexes 37 . Other recent SOC inventories also suggested that clay contents accounted for a small amount of variation in SOC stocks in rice paddies 38 . Rather, SOC stabilization in paddies is largely regulated by thermodynamic constraints of organic matter decomposition under anaerobic conditions 39 .
The impact of environmental factors on paddy soil C stocks can explain some of the main differences between countries. For instance, average paddy SOC stocks were 1.8-fold larger for China (65 Mg ha −1 ) than India (36 Mg ha −1 ). The difference between  countries partly reflects the specifics of climate; Chinese ricegrowing regions predominantly have subtropical climates, whereas Indian rice-growing regions have predominantly hot tropical climates 40 . Compared to other climate zones, the rate of paddy SOC decomposition in tropical climates is fast 40 .
Management practices also strongly affect paddy SOC stocks, which probably explains the relatively low amount of variation explained by environmental factors (Table 3) compared to nonagricultural upland soils (e.g., ref. 19 ). Our global meta-analysis corroborates previous national syntheses of paddy SOC dynamics under various management practices 41,42 . C gains in fertilized soils are explained by N and other nutrients stimulating plant growth and rhizodeposition, thereby increasing soil C input rates 43 . Fertilizer N addition can also stimulate soil C storage by slowing down the decomposition of plant litter and SOM (e.g., refs. 44,45 ). Specifically, N additions might reduce so-called "microbial N mining", whereby nutrient-poor conditions (e.g., low N) stimulate recalcitrant SOC decomposition by N-acquiring microbes [46][47][48] . Furthermore, organic fertilizers are additional C input into soil. Finally, organic fertilizers stimulate the succession of microbial communities favorable to SOC accumulation 49 .
However, in addition to previous findings (e.g., ref. 42 ), significant differences in SOC stocks were recorded between organic fertilizers (OF) and mineral plus organic fertilizers (NPK + OF). High N and other nutrient levels increase microbial growth on the available C pools, and so more necromass will be produced, which is a main component of SOM 48 .
Our estimates of no-till effects on paddy SOC stocks are quantitatively similar to no-till effects observed in upland soils 50 and indicate potential for SOC storage. Conservation tillage increases soil C storage by reducing aeration and the oxidative decomposition of SOC at periods of paddy soil drainage 40 and by increasing the physical and chemical protection of C from microbial attacks through organo-mineral associations 51 . Furthermore, conservation tillage reduces C losses associated with erosion 40 . Many paddies, especially in Asia, have been continuously tilled for hundreds to thousands of years. These practices likely depleted soil C stocks 27 , suggesting further potential for C storage under no-till. However, due to data paucity, our meta-analysis only considered C stocks in the top 30 cm. SOC gains in upper soil layers under no-till can be partly offset by losses at lower depths 52 , thereby reducing the SOC storage potential. Thus, to improve estimates of paddy SOC storage potential under no-till, whole profile analyses are still needed.
Our meta-analysis also explains some of the differences in SOC stocks between countries. For instance, rice paddies in eastern Asia (China, South Korea, and Japan), western Indonesian islands and Madagascar contained more SOC per area unit than paddies in western Africa, southern Asia, and South and Central America (Fig. 2a). These differences can be partly explained by management practices: farmers in southern/southeastern Asia and Africa often cannot afford sufficient mineral fertilizers to improve crop yield and support soil C storage 53 . Rice straw is also often removed for fodder and other uses in these regions 54 , thereby reducing soil C input rates and C stocks even further (Fig. 3). In contrast, high fertilizer application rates and high levels of crop residue return in China 55-57 contribute to high soil C stocks in Chinese rice paddies.
Another factor that might explain low soil C stocks in Africa might be the age of rice paddies. To feed a growing world population and to accommodate changing diets, global croplands have expanded by an average of 4 million hectares per year in recent decades 58 . Paddy expansion rates during this time differed strongly between continents, with half of the new global paddy area being located in Africa 16 . African rice paddies contain relatively low amounts of initial SOC compared to other continents 8 , suggesting a high potential of C sequestration. Our results also indicate considerable potential for paddy soil C sequestration in southern Asia, the central Indochina Peninsula and eastern South America (Fig. 2a). Realizing this potential requires adoption of recommended management practices such as crop residue incorporation, conservation tillage associated with seldom (once per 10-15 years) deep tillage, crop rotations with inclusion of grasses, legumes and deep rooting crops, and integrated nutrient management with pH adjustment through liming 2,59 . Promoting these practices will require new environmental and economic policies. For instance, aggregating small cropland patches can facilitate efficient fertilizer application, whereas farmer subsidies could provide an incentive for straw incorporation and rotations with deep rooting crops 56 .
Even though rice paddies store more SOC than the global average, this does not necessarily mean that the recent expansion in paddy area equates to a net climate benefit. First, new rice paddies are often established in ecosystems with relatively high soil C stocks, such as wetlands 60,61 . Second, rice paddies require a considerable amount of global irrigation water, accounting for 20% of total freshwater withdrawals by crops 62 . Pumping this water requires energy, which in turn causes ancillary CO 2 emissions. In gravity-fed irrigation systems or when pumping water from shallow aquifers this energy requirement can be minimal, but it can be high with diesel-based groundwater extraction systems or when using electricity not generated by hydropower 63 . Finally, and most importantly, CH 4 emissions from rice paddies are substantially higher than for other staple crops 62,64 , and rice paddies also produce considerable amounts of N 2 O (e.g., ref. 15 ). Thus, any benefits in terms of soil C sequestration with ricepaddy establishment need to be considered against a backdrop of increased greenhouse gas (GHG) emissions.
Soil C storage and GHG emissions can be compared directly by expressing them in CO 2 equivalents, using the global warming potential (GWP) values over a 100 year time horizon relative to CO 2 , i.e., 34 for CH 4 and 298 for N 2 O 65 . Global average GHG emissions from rice paddies have previously been estimated as 6300 kg CO 2 -eq ha −1 yr -1 for CH 4 16 and 280 kg CO 2 -eq ha −1 yr −1 Table 3 Correlation coefficients (R) and stepwise regression model between SOC stocks, climatic variables (MAT, n = 612, and MAP, n = 612), and soil properties (pH, n = 519, and clay content, n = 440). for N 2 O 15 (Supplementary Table 1). Average annual soil C storage in rice paddies can be estimated from average rice yields 16 and previously reported conversion factors 66,67 , and amounts to roughly 314 kg CO 2 -eq ha −1 yr −1 , i.e., an order of magnitude less than the combined emissions of CH 4 Table 1). Our estimates are corroborated by field studies showing that even after accounting for soil C storage, rice paddies remain a large net source of GHGs (e.g., refs. [68][69][70] . Supplementary Table 2). Management practices that increase SOC sequestration in rice paddies need to account for increased CH 4 emissions as well. Even though rice straw incorporation stimulates soil C storage, it more than doubles CH 4 emissions from rice paddies on average 71 . Previous syntheses suggest that the net effect of these two responses is negative, i.e., straw incorporation constitutes a net source of GHG emissions 72 . Whereas reduced till and no-till practices generally increase paddy soil C stocks, their effect on CH 4 emissions remains uncertain, with recent syntheses suggesting either increases 73 or decreases 74 in CH 4 emissions with no-till. Higher surface SOC with no-till may stimulate CH 4 production by increasing the availability of organic substrates 75 . On the other hand, increased soil macroporosity and soil pore continuity with no-till may accelerate gas diffusion and increase CH 4 oxidation 76 . The net effect of these opposing mechanisms is still unclear, and further research is needed to determine which of these mechanisms dominates under which conditions. Fertilizer N addition not only stimulates paddy SOC storage; it also stimulates N 2 O emissions from rice paddies 77 . The effect of fertilizer N on CH 4 emissions depends on application rates, with positive effects at low and medium rates, but negative effects at very high rates 77 . The increase in GHG emissions with fertilizer addition generally outweighs the climate benefit of soil C storage 68,69 . Moreover, the manufacturing and distribution of fertilizer requires energy and thus produces ancillary CO 2 emissions, possibly negating climate benefits 78 . In addition, excessive fertilizer N application in rice paddies causes a range of other environmental problems 79 .

and N 2 O (Supplementary
Although rice agriculture represents a large net source of GHGs compared to other staple crops, it also shows large potential for GHG mitigation through management 64 . For instance, mid-season drainage and intermittent irrigation can prevent the development of strong anaerobic conditions, thereby reducing CH 4 emissions by 53% 80 . While these practices stimulate N 2 O emissions, their net effect on GHG emissions is still negative 80 . Combining intermittent irrigation with several other management practices, the System of Rice Intensification may reduce both GHG emissions and the use of irrigation water (e.g., refs. 63,81 ). Applying rice straw off season rather than in season may reduce global CH 4 emissions by 4.1 Tg year -182 . Moreover, selecting high-yielding cultivars can simultaneously reduce CH 4 emissions and increase crop yields 83 , and so, the C input into soil. A full accounting of the mitigation potential of these measures calls for long-term experiments under real-world conditions that account for changes in soil C stocks, interactions between management practices, as well as direct and ancillary GHG emissions. In addition, models such as DNDC-rice (e.g., ref. 84 ) may be used to evaluate trade-offs between soil C sequestration and CH 4 and N 2 O emissions under a range of management practices 85 .
In summary, we present the first global assessment of paddy soil C stocks. Our results identify paddy soils as an important C pool, containing~20% more SOC per hectare than croplands on average. Our analysis underlines the role of both natural factors and agronomical management in determining paddy SOC stocks; fertilization, straw incorporation, and no-till practices all increased paddy SOC storage, whereas SOC levels decreased with MAT and soil pH. However, the climate benefit of SOC storage in rice paddies is generally outweighed by increases in GHG emissions. These data underline the importance of paddies in the terrestrial C cycle, and should be used to improve global C inventories and to inform policy advice related to land use.

Methods
Data collection. We used Web of Science, Google Scholar, and China National Knowledge Infrastructure to search for studies published between 1999 and 2019, applying the search terms "paddy AND soil organic carbon" and "rice AND soil organic carbon". We only considered studies reporting contents for soil organic carbon (SOC) or soil organic matter (SOM) and sampling depth for the quantitative determination of paddy SOC storage (excluding upland rice). For each study we tabulated SOC contents and sampling depth; we also tabulated soil bulk density (BD) data when these were reported. We selected a total of 239 publications reporting 2234 sets of raw data from rice paddies around the world for analyses (see Supplementary Data 1). To avoid data duplication, we checked the latitude and longitude of all the sites included in our dataset, and we eliminated duplicated sites. Finally, our dataset included 612 sampling locations in 50 countries (out of 118 rice-producing countries, Fig. 1a) encompassing 95% of the global paddy area and 98% of the global rice production.
We also tabulated the following information for each study: (1) geographical location of sites (latitude, longitude, elevation above sea level, and country), (2) climatic conditions (mean annual temperature, MAT; and mean annual precipitation, MAP), and (3) properties of the paddy soil (pH and clay content). Missing data for latitude, longitude, and elevation were estimated using Google Maps (https://maps.google.com/). Missing data for MAT and MAP were obtained from https://en.climate-data.org/.

Data processing
Estimating missing data for SOC content and BD. Paddy SOC stocks were calculated separately for the topsoil (0-30 cm) and subsoil (30-100 cm) to facilitate a comparison between our results and global SOC storage data in the Harmonized World Soil Database 86 . Original data reported as SOM content were converted to SOC content using the conventional "van Bemmelen factor" of 1.724 87 .
Most SOC data in our database were for soil layers to a maximum depth of 30 cm or less. To extrapolate these data to lower depths, we used a subset of 42 studies that reported SOC contents for 409 profiles ≥50 cm. We assume that soil compaction was similar across all the profiles and could be extrapolated to other SOC data in the deep soil. The relative ratio of SOC content (RR SOC , see ref. 88 ) was first calculated as: where SOC surface is the SOC content (g kg −1 ) of surface soil and SOC below is the SOC content below the surface soil for various depths in the profile. Combing the data for all studies in the subset, the relationship between RR SOC and soil depth could be described by a logarithmic curve ( Supplementary Fig. 3, R 2 = 0.63, n = 1227, see ref. 89 ). SOC content at depth i (cm) was then estimated as: The depth gradient for BD was less pronounced than for SOC content ( Supplementary Fig. 4a). The relationship between BD and SOC content of the topsoil in our dataset could be described by a negative power function ( Fig. S5; R 2 = 0.49, n = 1370): The availability of subsoil BD data was insufficient to perform a regression analysis with SOC content. Because the ratio of subsoil over topsoil bulk density averaged 1.18 across our dataset (Supplementary Fig. 4b; standard error = 0.01, n = 376), we estimated missing subsoil BDs by multiplying the topsoil BDs by 1.18.
Estimating paddy SOC stocks in the topsoil and subsoil. For each soil layer at each sampling location in our dataset, total SOC stock (SOCT, Mg ha −1 ) was calculated according to ref. 21 where SOC and BD are SOC content (g kg −1 ) and bulk density (g cm −3 ), respectively, H is soil thickness (cm) and δ 2mm is the fraction (%) of fragments >2 mm in the soil. Since the paddy soils were mostly derived from deposits in flat areas, the >2 mm fraction of the total mass is usually negligible 21 .
When SOC content and BD data were available for the entire 0-100 cm profile, we calculated SOC stock (Mg ha −1 ) in the 0-30 cm topsoil layer (SOCT 30 ) and the entire 0-100 cm profile (SOCT 100 ) by adding SOCT of all soil layers within the 0-30 cm and 0-100 cm range, respectively.
Estimating national and global paddy SOC stocks. National paddy SOC stocks (SOCT N , Pg) for any country in our dataset were estimated as: where SOCTmean (Mg ha −1 ) is the mean SOC stock across all sampling locations in that country and HA is the rice harvest area (ha) in that country. HA data for all countries in our analysis were derived from FAO 15 . The global SOC stock (SOCT G , Pg) was estimated as: where SOCTNa is paddy SOC stock (Pg) in country a, and SOCT Nn is paddy SOC stock (Pg) in country n. No SOC data were available for some countries where HA was small (e.g., Congo, Mali, and Peru). For these countries, we estimated SOCT N using the SOC N data from neighboring countries with the closest climate based on Köppen-Geiger climate classification. These estimates did not substantially affect the estimates of the global SOC stocks because of the small areas of these countries (~5% of total HA). Global mean SOC stocks per unit area was calculated based on global SOC stocks (SOCT G ) divided by global rice-paddy area.
Data analysis. The importance of the environmental variables was estimated using Pearson's and partial correlation coefficients, which are commonly used to measure the association between variables, implemented in SPSS 20.0 (SPSS, Chicago, USA). In these correlations, the p value defines whether two variables are statistically correlated. p values below 0.05 were accepted as significant correlations. Optimized model for SOC with environmental variables was determined by stepwise regression using forward selection criteria (p of 0.05 for entering and 0.1 for removal). ArcGIS 10.3 (Esri, Redlands, USA) was used to analyze and visualize the spatial distribution of SOC stocks. The 95% confidence interval (CI) of SOC stocks (SOCT N ) in each country was calculated by bootstrapping, using 4999 iterations 90 . The uncertainty (U) of the total SOC stock in each country was then calculated as: where x is the SOC stock in a country, and CI is the 95% confidence interval of x.
The total uncertainty (U total ) at the global scale was calculated as: where U a and U n are the uncertainties associated with x a and x n in country a and country n. Because some countries were represented by only one site, the CI of that country could not be calculated. In these cases, the coefficient of variation was conservatively set to 50% 91 . Global SOC stocks were then estimated as described by Eq. 8, and the 95% confidence CI of this estimate was calculated using Eq. 9.
Meta-analysis. We assessed the effects of management practices (fertilization, return of straw, and tillage) on paddy SOC stocks by creating subsets of experiments that included side-by-side comparisons between management practices. Studies had to meet specific criteria to be included in the dataset. First, growing conditions in the control and treatment plots had to be identical (except for the management practice being studied). Second, mean SOC stock and the number of field replicates had to be reported for both control and treatment plots. Studies were incorporated into seven datasets, based on seven types of management practices: (1) addition of mineral nitrogen (N) fertilizer, (2) addition of mineral nitrogen-phosphorus-potassium (NPK) fertilizer, (3) addition of organic fertilizer (e.g., green/farmyard manure, compost), (4) addition of mineral NPK and organic fertilizer, (5) no tillage, (6) reduced tillage, and (7) the return of straw. The effects of these management practices on SOC stocks were quantified as the natural log of the response ratio (lnRR), a metric commonly used in metaanalyses: 92 lnRR ¼ lnðSOCT t =SOCT c Þ ð 11Þ where SOCTt and SOCTc represent the mean SOC stock in the 0-30 cm layer of the treatment and control groups, respectively. SOC stocks for the 0-30 cm layer were calculated as described above. Treatments included (i) the application of mineral fertilizers (e.g., N or NPK), organic fertilizer (OF), and mineral plus organic fertilizers (NPK + OF) versus no application, (ii) straw return versus no return, and (iii) no or reduced tillage versus conventional tillage.
Most studies in our analysis did not report the standard deviations of the means. We therefore adopted a replication-based weighting method: 93,94 Weight ¼ ðn t n c Þ=ðn t þn c Þ ð 12Þ where nt and nc are the numbers of replicates of the treatment and control, respectively. Mean effect sizes and 95% confidence intervals (CIs) were generated by bootstrapping with 4999 iterations using MetaWin 2.1 90 . Effects of paddy management were considered significant if the 95% CIs did not overlap with zero. To ease interpretation, results were back-transformed to percent change ((RR − 1) × 100) in SOC stocks. Positive and negative changes indicate increases and decreases due to the management practices, respectively.