Assessing Spatial Representativeness of Global Flux Tower Eddy-Covariance Measurements Using Data from FLUXNET2015.

Large datasets of carbon dioxide, energy, and water fluxes were measured with the eddy-covariance (EC) technique, such as FLUXNET2015. These datasets are widely used to validate remote-sensing products and benchmark models. One of the major challenges in utilizing EC-flux data is determining the spatial extent to which measurements taken at individual EC towers reflect model-grid or remote sensing pixels. To minimize the potential biases caused by the footprint-to-target area mismatch, it is important to use flux datasets with awareness of the footprint. This study analyze the spatial representativeness of global EC measurements based on the open-source FLUXNET2015 data, using the published flux footprint model (SAFE-f). The calculated annual cumulative footprint climatology (ACFC) was overlaid on land cover and vegetation index maps to create a spatial representativeness dataset of global flux towers. The dataset includes the following components: (1) the ACFC contour (ACFCC) data and areas representing 50%, 60%, 70%, and 80% ACFCC of each site, (2) the proportion of each land cover type weighted by the 80% ACFC (ACFCW), (3) the semivariogram calculated using Normalized Difference Vegetation Index (NDVI) considering the 80% ACFCW, and (4) the sensor location bias (SLB) between the 80% ACFCW and designated areas (e.g. 80% ACFCC and window sizes) proxied by NDVI. Finally, we conducted a comprehensive evaluation of the representativeness of each site from three aspects: (1) the underlying surface cover, (2) the semivariogram, and (3) the SLB between 80% ACFCW and 80% ACFCC, and categorized them into 3 levels. The goal of creating this dataset is to provide data quality guidance for international researchers to effectively utilize the FLUXNET2015 dataset in the future.

Towards a better diagnosis of mouth breathing: validity and reliability of a protocol for assessing the awake breathing pattern in preschool children.

PURPOSE: The Awake Breathing Pattern Assessment (ABPA) is a prototypical clinical grid recently designed through an international consensus of Speech and Language Pathologists (SLPs) to categorize the awake and habitual breathing pattern during the orofacial myofunctional assessment. This cross-sectional study aims to explore the psychometric properties of the ABPA in a preschool population. METHODS: 133 children from 2;11 to 6 years old were assessed with the ABPA. The percentage of time spent breathing through the mouth was objectively measured by a CO2 sensor and used as a baseline measurement. We first performed a multivariate Latent Profile Analysis based on the CO2 measurement and a parental questionnaire to define the number of categories that best characterize the breathing pattern. Subsequently, we assessed the intra- and inter-rater reliability, internal consistency criterion validity, construct validity and sensitivity and specificity. RESULTS: The awake breathing pattern can best be described by two groups: nasal and mouth breathing. The ABPA, initially designed in three groups, was adjusted accordingly. This final version showed excellent intra-rater and inter-rater reliability. There was a significant correlation between the ABPA and the CO2 measurement. The ABPA showed a fair sensitivity and a good specificity. CONCLUSION: The reference tool based on CO2 data was used in children for the first time and was found to be reliable. The ABPA is a suitable tool for SLPs to confirm the diagnosis of mouth breathing in preschool children if more sensitive screening tools, like parental questionnaires, are used beforehand.

Cover crops support the climate change mitigation potential of agroecosystems.

Cover crops have the potential to mitigate climate change by reducing negative impacts of agriculture on ecosystems. This study is first to quantify the net climate change mitigation impact of cover crops including land-use effects. A systematic literature and data review was conducted to identify major drivers for climate benefits and costs of cover crops in maize (Zea maize L.) production systems. The results indicate that cover crops lead to a net climate change mitigation impact (NCCMI) of 3.30 Mg CO2e ha-1 a-1. We created four scenarios with different impact weights of the drivers and all of them showing a positive NCCMI. Carbon land benefit, the carbon opportunity costs based on maize yield gains following cover crops, is the major contributor to the NCCMI (34.5% of all benefits). Carbon sequestration is the second largest contributor (33.8%). The climate costs of cover crops are mainly dominated by emissions from their seed production and foregone benefits due to land use for cover crops seeds. However, these two costs account for only 15.8% of the benefits. Extrapolating these results, planting cover crops before all maize acreage in the EU results in a climate change mitigation of 49.80 million Mg CO2e a-1, which is equivalent to 13.0% of the EU's agricultural emissions. This study highlights the importance of incorporating cover crops into sustainable cropping systems to minimize the agricultural impact to climate change.

Herpes zoster in older adults: Impact on carbon footprint in the United States.

We provide estimates for (I) annual herpes zoster (HZ) cases, (II) carbon costs related to healthcare utilization, and (III) annual carbon emissions due to HZ among ≥50 years of age (YOA) United States (US) population. We estimated the annual number of HZ cases in the US based on available incidence data and demographic data of individuals ≥50 YOA. Both the healthcare resource utilization (HCRU) associated with HZ cases and the unit carbon dioxide equivalent (i.e. CO2e) costs associated with each type of HCRU in the US were estimated based on literature and studies available online. The carbon footprint associated with HZ annually among US adults ≥50 YOA was estimated by multiplying the unit carbon estimates by the HCRU. In the US population aged ≥50 YOA in 2020 (i.e. approximately 118 million), approximately 1.1 million cases of HZ occur annually assuming no vaccination. Based on 2 sources of HCRU the average kgCO2e per HZ patient ranged from 61.0 to 97.6 kgCO2e, with values by age group ranging from 40.9 kgCO2e in patients aged 50-59 to 195.9 kgCO2e in patients ≥80 YOA. The total annual HZ associated carbon ranged between 67,000 and 107,000 tons of CO2e in the US population aged ≥50 YOA. The impact of HZ on carbon footprint in the US results in considerable greenhouse gas (GHG)emissions. Assuming no vaccination, the burden of HZ is projected to rise over the coming years with the aging populations consequently worsening its impact on GHG emissions. (Figure 1).

Impact of vegetation composition and seasonality on sensitivity of modelled CO2 exchange in temperate raised bogs.

Encroachment of vascular plants (VP) in temperate raised bogs, as a consequence of altered hydrological conditions and nutrient input, is widely observed. Effects of such vegetation shift on water and carbon cycles are, however, largely unknown and identification of responsible plant physiological traits is challenging. Process-based modelling offers the opportunity of gaining insights into ecosystem functioning beyond observations, and to infer decisive trait shifts of plant functional groups. We adapted the Soil-Vegetation-Atmosphere Transfer model pyAPES to a temperate raised bog site by calibration against measured peat temperature, water table and surface CO2 fluxes. We identified the most important traits determining CO2 fluxes by conducting Morris sensitivity analysis (MSA) under changing conditions throughout the year and simulated VP encroachment. We further investigated transferability of results to other sites by extending MSA to parameter ranges derived from literature review. We found highly variable intra-annual plant traits importance determining ecosystem CO2 fluxes, but only a partial shift of importance of photosynthetic processes from moss to VP during encroachment. Ecosystem respiration was dominated by peat respiration. Overall, carboxylation rate, base respiration rate and temperature sensitivity (Q10) were most important for determining bog CO2 balance and parameter ranking was robust even under the extended MSA.

Individual and combined impacts of carbon dioxide enrichment, heatwaves, flow velocity variability, and fine sediment deposition on stream invertebrate communities.

Climate change and land-use change are widely altering freshwater ecosystem functioning and there is an urgent need to understand how these broad stressor categories may interact in future. While much research has focused on mean temperature increases, climate change also involves increasing variability of both water temperature and flow regimes and increasing concentrations of atmospheric CO2, all with potential to alter stream invertebrate communities. Deposited fine sediment is a pervasive land-use stressor with widespread impacts on stream invertebrates. Sedimentation may be managed at the catchment scale; thus, uncovering interactions with these three key climate stressors may assist mitigation of future threats. This is the first experiment to investigate the individual and combined effects of enriched CO2, heatwaves, flow velocity variability, and fine sediment on realistic stream invertebrate communities. Using 128 mesocosms simulating small stony-bottomed streams in a 7-week experiment, we manipulated dissolved CO2 (ambient; enriched), fine sediment (no sediment; 300 g dry sediment), temperature (ambient; two 7-day heatwaves), and flow velocity (constant; variable). All treatments changed community composition. CO2 enrichment reduced abundances of Orthocladiinae and Chironominae and increased Copepoda abundance. Variable flow velocity had only positive effects on invertebrate abundances (7 of 13 common taxa and total abundance), in contrast to previous experiments showing negative impacts of reduced velocity. CO2 was implicated in most stressor interactions found, with CO2 × sediment interactions being most common. Communities forming under enriched CO2 conditions in sediment-impacted mesocosms had ~20% fewer total invertebrates than those with either treatment alone. Copepoda abundances doubled in CO2-enriched mesocosms without sediment, whereas no CO2 effect occurred in mesocosms with sediment. Our findings provide new insights into potential future impacts of climate change and land use in running freshwaters, in particular highlighting the potential for elevated CO2 to interact with fine sediment deposition in unpredictable ways.

Monitoring the ventilation of living spaces to assess the risk of airborne transmission of infection using a novel Pocket CO2 Logger to track carbon dioxide concentrations in Tokyo.

We employed carbon dioxide (CO2) concentration monitoring using mobile devices to identify location-specific risks for airborne infection transmission. We lent a newly developed, portable Pocket CO2 Logger to 10 participants, to be carried at all times, for an average of 8 days. The participants recorded their location at any given time as cinema, gym, hall, home, hospital, other indoors, other outgoings, pub, restaurant, university, store, transportation, or workplace. Generalized linear mixed model was used for statistical analysis, with the objective variable set to the logarithm of CO2 concentration. Analysis was performed by assigning participant identification as the random effect and location as the fixed effect. The data were collected per participant (seven males, four females), resulting in a total of 12,253 records. Statistical analysis identified three relatively poorly ventilated locations (median values > 1,000 ppm) that contributed significantly (p < 0.0001) to CO2 concentrations: homes (1,316 ppm), halls (1,173 ppm), and gyms (1005ppm). In contrast, two locations were identified to contribute significantly (p < 0.0001) to CO2 concentrations but had relatively low average values (<1,000 ppm): workplaces (705 ppm) and stores (620 ppm). The Pocket CO2 Logger can be used to visualize airborne infectious transmission risk by location to help guide recommendation regarding infectious disease policies, such as restrictions on human flow and ventilation measures and guidelines. In the future, large-scale surveys are expected to utilize the global positioning system, Wi-Fi, or Bluetooth of an individual's smartphone to improve ease and accuracy.

Estimated greenhouse gas emissions in the Mbalmayo thermal power plant between 2016-2020 using the genetic-Gaussian algorithm coupling.

The greenhouse gas (GHG) emissions inventories in our context result from the production of electricity from fuel oil at the Mbalmayo thermal power plant between 2016 and 2020. Our study area is located in the Central Cameroon region. The empirical method of the second level of industrialisation was applied to estimate GHG emissions and the application of the genetic algorithm-Gaussian (GA-Gaussian) coupling method was used to optimise the estimation of GHG emissions. Our work is of an experimental nature and aims to estimate the quantities of GHG produced by the Mbalmayo thermal power plant during its operation. The search for the best objective function using genetic algorithms is designed to bring us closer to the best concentration, and the Gaussian model is used to estimate the concentration level. The results obtained show that the average monthly emissions in kilograms (kg) of GHGs from the Mbalmayo thermal power plant are: 526 kg for carbon dioxide (CO2), 971.41 kg for methane (CH4) and 309.41 kg for nitrous oxide (N2O), for an average monthly production of 6058.12 kWh of energy. Evaluation of the stack height shows that increasing the stack height helps to reduce local GHG concentrations. According to the Cameroonian standards published in 2021, the limit concentrations of GHGs remain below 30 mg/m3 for CO2 and 200 µg/m3 for N2O, while for CH4 we reach the limit value of 60 µg/m3. These results will enable the authorities to take appropriate measures to reduce GHG concentrations.

Carbon footprint of a laser unit: a study of two centres in the UK.

PURPOSE: Climate change has serious consequences for our wellbeing. Healthcare systems themselves contribute significantly to our total carbon footprint, of which emissions from surgical practice are a major component. The primary sources of emissions identified are anaesthetic gases, disposal of single-use equipment, energy usage, and travel to and from clinical areas. We sought to quantify the waste generated by laser surgery which, to our knowledge, has not been previously reported. METHODS: The carbon footprint of two laser centres operating within the United Kingdom were measured. The internationally recognised Greenhouse Gas Protocol was used as a guiding framework to classify sources of waste and conversion factors issued by the UK government were used to quantify emissions. RESULTS: The total carbon footprints per day at each unit were 299.181 carbon dioxide equivalents (kgCo2eq) and 121.512 kgCO2eq, respectively. We found the carbon footprint of individual laser treatments to be approximately 15 kgCO2eq per procedure. The biggest overall contributor to the carbon footprint was found to be the emissions generated from staff, patient and visitor travel. This was followed by electricity usage, and indirect emissions from physical waste and laundry. CONCLUSIONS: The carbon footprint of laser procedures was considerably less than the average surgical operation in the UK. This initial study measures the carbon footprint of a laser center in a clinical setting and allows us to identify where improvements can be made to eventually achieve a net carbon zero health care system.

Characteristics of PM2.5 bounded carbonaceous aerosols, carbon dioxide and its stable carbon isotopes (δ13C) in rural households in northwest China: Effect of different fuel combustion.

In order to fully understand the carbon emission from different fuels in rural villages of China, especially in the typical atmospheric pollution areas. The characteristics of carbonaceous aerosols and carbon dioxide (CO2) with its stable carbon isotope (δ13C) were investigated in six households, which two households used coal, two households used wood as well as two households used biogas and liquefied petroleum gas (LPG), from two rural villages in Fenwei Plain from March to April 2021. It showed that the fine particulate matter (PM2.5) emitted from biogas and LPG couldn't be as lower as expected in this area. However, the clean fuels could relatively reduce the emissions of organic carbon (OC) and element carbon (EC) in PM2.5 compare to the solid fuels. The pyrolyzed carbon (OP) accounted more total carbon (TC) in coal than the other fuels use households, indicating that more water-soluble OC existed, and it still had the highest secondary organic carbon (SOC) than the other fuels. Meantime, the coal combustions in the two villages had the highest CO2 concentration of 527.6 ppm and 1120.6 ppm, respectively, while the clean fuels could effectively reduce it. The average δ13C values (-26.9‰) was much lighter than almost all the outdoor monitoring and similar to the δ13C values for coal combustion and vehicle emission, showing that they might be the main contributors of the regional atmospheric aerosol in this area. During the sandstorm, the indoor PM2.5 mass and CO2 were increasing obviously. The indoor cancer risk of PAHs for adults and children were greater than 1 × 10-6, exert a potential carcinogenic risk to human of solid fuels combustion in rural northern China. It is important to continue concern the solid fuel combustion and its health impact in rural areas.

Financial regulations and sustainability: The role of energy price and climate policy uncertainty.

Russia ranks among the top five countries worldwide in terms of carbon emissions, with the energy, transportation, and manufacturing sectors as the major contributors. This poses a significant threat to both current and future generations. Russia faces challenges in achieving Sustainable Development Goal 13, necessitating the implementation of more innovative policies to promote environmental sustainability. Considering this alarming situation, this study investigates the role of financial regulations, energy price uncertainty, and climate policy uncertainty in reshaping sectoral CO2 emissions in Russia. This study utilizes a time-varying bootstrap rolling-window causality (BRW) approach using quarterly data from 1990 to 2021. The stability test for parameters indicates instability, suggesting that the full sample causality test may yield incorrect inferences. Thus, the BRW approach is employed for valid inferences. Our findings confirm the time-varying negative impact of financial regulations on CO2 emissions from energy, manufacturing, and transportation sectors. Additionally, findings confirm time-varying positive impact of energy prices and climate policy uncertainty on CO2 emissions from the energy, manufacturing, and transportation sectors. Strong financial regulations and stable energy and climate policies are crucial for achieving sustainability, highlighting significant policy implications for policymakers and stakeholders.

Greenhouse gas mitigation and soil carbon stabilization potential of forest biochar varied with biochar type and characteristics.

Biochar is increasingly used in climate-smart agriculture, yet its impact on greenhouse gas (GHG) emissions and soil carbon (C) sequestration remains poorly understood. This study examined biochar-mediated changes in soil properties and their contribution to C stabilization and GHG mitigation by evaluating four types of biochar. Soil carbon dioxide (CO2) and nitrous oxide (N2O) emissions, soil chemical and biological properties, and soil organic carbon (SOC) mineralization kinetics were monitored using greenhouse, laboratory, and modeling experiments. Three pine wood biochars pyrolyzed at 460 °C (PB-460), 500 °C (PB-500), 700 °C (PB-700), and one pine bark biochar from gasification at 760 °C (GB-760) were added into soil at 1 % w/w basis. Soils amended with biochar were used to cultivate sorghum for three months in a greenhouse, followed by three months of laboratory incubation. Data obtained from laboratory incubation was modeled using various statistical approaches. The PB-500 and PB-700 reduced cumulative N2O-N emissions by 68.5 % and 73.9 % and CO2 equivalent C emissions by 66.9 % and 72.4 %, respectively, compared to unamended control. The N2O emissions were positively associated with soil nitrate N, available P, and biochar ash content while negatively associated with SOC. The CO2 emission was negatively related to biochar C:N ratio and volatile matter content. Biochar amended soils had 49.2 % (PB-500) to 87.7 % (PB-700) greater SOC and 22.9 % (PB-700) to 48.1 % (GB-760) greater sorghum yield than the control. While PB-700 had more saprophytes than the control, the GB-760 yielded a greater yield than biochars prepared by pyrolysis. Microbial biomass C was 7.23 to 23.3 % greater in biochar amended soils than in control. The double exponential decay model best explained the dynamics of C mineralization, which was associated with initial soil nitrate N and available P positively and total fungi and protozoa biomass negatively. Biochar amendment could be a climate smart agricultural strategy. Pyrolysis pine wood biochar showed the greatest potential to reduce GHG emissions and enhance SOC storage and stability, and gasification biochar contributed more to SOC storage and increased crop yield.

Effects of extracorporeal CO2 removal on gas exchange and ventilator settings: a systematic review and meta-analysis.

PURPOSE: A systematic review and meta-analysis to evaluate the impact of extracorporeal carbon dioxide removal (ECCO2R) on gas exchange and respiratory settings in critically ill adults with respiratory failure. METHODS: We conducted a comprehensive database search, including observational studies and randomized controlled trials (RCTs) from January 2000 to March 2022, targeting adult ICU patients undergoing ECCO2R. Primary outcomes were changes in gas exchange and ventilator settings 24 h after ECCO2R initiation, estimated as mean of differences, or proportions for adverse events (AEs); with subgroup analyses for disease indication and technology. Across RCTs, we assessed mortality, length of stay, ventilation days, and AEs as mean differences or odds ratios. RESULTS: A total of 49 studies encompassing 1672 patients were included. ECCO2R was associated with a significant decrease in PaCO2, plateau pressure, and tidal volume and an increase in pH across all patient groups, at an overall 19% adverse event rate. In ARDS and lung transplant patients, the PaO2/FiO2 ratio increased significantly while ventilator settings were variable. "Higher extraction" systems reduced PaCO2 and respiratory rate more efficiently. The three available RCTs did not demonstrate an effect on mortality, but a significantly longer ICU and hospital stay associated with ECCO2R. CONCLUSIONS: ECCO2R effectively reduces PaCO2 and acidosis allowing for less invasive ventilation. "Higher extraction" systems may be more efficient to achieve this goal. However, as RCTs have not shown a mortality benefit but increase AEs, ECCO2R's effects on clinical outcome remain unclear. Future studies should target patient groups that may benefit from ECCO2R. PROSPERO Registration No: CRD 42020154110 (on January 24, 2021).

The role of environmental tax on the environmental quality in EU counties: evidence from panel vector autoregression approach.

This study intends to analyze the influence of environmental taxes on pollution in EU-27 nations. Furthermore, energy from renewable sources consumption and urbanization are employed to clarify CO2 emissions in this study that tests the EKC hypothesis. According to the findings, an increase in environmental taxes reduces CO2 emissions by 0.14%. Also, the data supported the validity of the EKC concept. The findings of the causality test demonstrated that there is a bidirectional causal link between CO2 emissions and environmental taxes. These results reflect that environmental tax revenues contribute to sustainability as an effective policy tool in EU countries. Policies regarding environmental tax enforcement come to the fore in terms of both keeping the balance in economic activities and serving sustainability.

Assessment of organizational carbon footprints in a denim-washing company: a systematic approach to indirect non-energy emissions.

As stated in the 2016 Paris Agreement, concerns about global climate change and carbon emissions have increased, and organizations, in particular, have embarked on an annual measurement process to estimate their contribution to global climate change. Carbon footprint, one of the measurement methods, is a widely applied tool to assess the environmental impact of organizations. This study presents a real case study of a denim-washing company's activities based on ISO standard calculation methods of greenhouse gas emissions. Accordingly, the annual carbon footprint of the denim-washing company was 2482.09 tCO2e for the year 2021 in total for the overall carbon footprint. Direct emission was calculated at 1575.75 tCO2e, indirect energy-related emission at 798.09 tCO2e, and indirect non-energy-related emission at 108.25 tCO2e. The highest CO2 emissions are related to heating from greenhouse gas direct emission sources, followed by purchased electricity consumption, and the lowest CO2 emissions are related to fire-CO2 tube storage. In conclusion, this study is particular in that it analyzes not only the specific processes of a denim-washing company but also the overall organizational carbon footprint calculation, assesses the importance of indirect non-energy in the total carbon footprint, and evaluates the calculation findings with sector-specific mitigation strategies.

Least-cost decarbonization pathways for electricity generation in Finland: A convex quantile regression approach.

This study investigates the least-cost decarbonization pathways in the Finnish electricity generation industry in order to achieve the national carbon neutrality goal by 2035. Various abatement measures, such as downscaling production, capital investment, and increasing labor and intermediate inputs, are considered. The marginal abatement costs (MACs) of greenhouse gas emissions are estimated using the convex quantile regression method and applied to unique register-based firm-level greenhouse gas emission data merged with financial statement data. We adjust the MAC estimates for the sample selection bias caused by zero-emission firms by applying the two-stage Heckman correction. Our empirical findings reveal that the median MAC ranges from 0.1 to 3.5 euros per tonne of CO2 equivalent. The projected economic cost of a 90% reduction in emissions is 62 million euros, while the estimated cost of achieving zero emissions is 83 million euros.

Optimizing bioelectromethanosynthesis of CO2 and membrane fouling mitigation in MECs via in-situ biogas recirculation.

The CO2 bioelectromethanosynthesis via two-chamber microbial electrolysis cell (MEC) holds tremendous potential to solve the energy crisis and mitigate the greenhouse gas emissions. However, the membrane fouling is still a big challenge for CO2 bioelectromethanosynthesis owing to the poor proton diffusion across membrane and high inter-resistance. In this study, a new MEC bioreactor with biogas recirculation unit was designed in the cathode chamber to enhance secondary-dissolution of CO2 while mitigating the contaminant adhesion on membrane surface. Biogas recirculation improved CO2 re-dissolution, reduced concentration polarization, and facilitated the proton transmembrane diffusion. This resulted in a remarkable increase in the cathodic methane production rate from 0.4 mL/L·d to 8.5 mL/L·d. A robust syntrophic relationship between anodic organic-degrading bacteria (Firmicutes 5.29%, Bacteroidetes 25.90%, and Proteobacteria 6.08%) and cathodic methane-producing archaea (Methanobacterium 65.58%) enabled simultaneous organic degradation, high CO2 bioelectromethanosynthesis, and renewable energy storage.

Contrasting impacts of fertilization on topsoil and subsoil greenhouse gas fluxes in a thinned Chinese fir plantation.

Globally, forest soils are considered as important sources and sinks of greenhouse gases (GHGs). However, most studies on forest soil GHG fluxes are confined to the topsoils (above 20 cm soil depths), with only very limited information being available regarding these fluxes in the subsoils (below 20 cm soil depths), especially in managed forests. This limits deeper understanding of the relative contributions of different soil depths to GHG fluxes and global warming potential (GWP). Here, we used a concentration gradient-based method to comprehensively investigate the effects of thinning intensity (15% vs. 35%) and nutrient addition (no fertilizer vs. NPK fertilizers) on soil GHG fluxes from the 0-40 cm soil layers at 10 cm depth intervals in a Chinese fir (Cunninghamia lanceolata) plantation. Results showed that forest soils were the sources of CO2 and N2O, but the sinks of CH4. Soil GHG fluxes decreased with increasing soil depth, with the 0-20 cm soil layers identified as the dominant producers of CO2 and N2O and consumers of CH4. Thinning intensity did not significantly affect soil GHG fluxes. However, fertilization significantly increased CO2 and N2O emissions and CH4 uptake at 0-20 cm soil layers, but decreased them at 20-40 cm soil layers. This is because fertilization alleviated microbial N limitation and decreased water filled pore space (WFPS) in topsoils, while it increased WFPS in subsoils, ultimately suggesting that soil WFPS and N availability (especially NH4+-N) were the predominant regulators of GHG fluxes along soil profiles. Generally, there were positive interactive effects of thinning and fertilization on soil GHG fluxes. Moreover, the 35% thinning intensity without fertilization had the lowest GWP among all treatments. Overall, our results suggest that fertilization may not only cause depth-dependent effects on GHG fluxes within soil profiles, but also impede efforts to mitigate climate change by promoting GHG emissions in managed forest plantations.