Carbon Dioxide Fluxes Research Articles

River network‐scale drying impacts the spatiotemporal dynamics of greenhouse gas fluxes

AbstractRivers significantly contribute to global biogeochemical cycles; however, we have a limited understanding of how drying may influence these cycles. Drying fragments river networks, thereby influencing important ecosystem functions such as the processing of carbon and nitrogen, and associated fluxes of greenhouse gases (GHGs) both locally, and at the river network scale. Our objective was to assess, using a network‐scale approach, the lateral, longitudinal, and temporal dynamics of GHG fluxes in a river network naturally fragmented by drying. We used a closed‐loop chamber with automated analyzers to measure carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) fluxes from dry sediments, flowing waters, isolated pools, and riparian soils, along with a suite of environmental variables, over 9 months at 20 sites across a non‐perennial river network in France. Network‐scale drying had a spatial and temporal legacy effect on GHG fluxes. On average, CO2 fluxes were up to 29 times higher from perennial than non‐perennial sites under flowing conditions. At non‐perennial sites, CO2 and N2O fluxes positively covaried with time since rewetting. In addition, CO2 and N2O fluxes at perennial sites positively covaried with the percent of non‐perennial reaches upstream, indicating a spatial effect of drying. GHG fluxes from riparian soil and dry riverbed sediments had markedly different magnitudes and covariates. This research demonstrates that drying not only has a local‐scale impact but also influences GHG fluxes at the network scale, contributing valuable insights for upscaling global riverine GHG estimates.

Lawns and meadows in urban green space – a comparison from perspectives of greenhouse gases, drought resilience and plant functional types

Abstract. Today, city planners design urban futures by considering environmental degradation and climate mitigation. Here, we studied the greenhouse gas fluxes of urban lawns and meadows and linked the observations with plant functional types and soil properties. In eight lawns and eight meadows in the Helsinki metropolitan area, Finland, carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O) fluxes were measured using manual chambers, and plant functional types were recorded. Four of these sites, i.e. an irrigated lawn, an old mesic meadow, a non-irrigated lawn and a young dry meadow, were more intensively studied in 2021–2022. The process-based ecosystem model JSBACH was utilized together with the momentary observations collected approximately every second week on CO2 exchange to quantify the annual carbon (C) balance of these sites. On the remaining sites, we studied the initial dynamics of conversion from lawns to meadows by transforming parts of lawns to meadows in late 2020 and conducting measurements from 2020 to 2022. The mean photosynthetic production (GPP) of the irrigated lawn and mesic meadow was the highest in this study, whereas the dry meadow had the lowest GPP. The studied lawns were stronger C sinks compared to the meadows. However, the net exchange values were uncertain as the soils were not in equilibrium with the vegetation at all sites, which is common for urban habitats, and modelling the heterotrophic emissions was therefore challenging. The conversion from a lawn to a meadow did not affect the fluxes of CH4 and N2O. Moreover, the mesic meadow was more resistant to drought events than the non-irrigated lawn. Lastly, the proportion of herbaceous flowering plants other than grasses was higher in meadows than in lawns. Even though social and economic aspects also steer urban development, these results can guide planning when considering environmentally friendlier green spaces and carbon smartness.

Seasonal dynamics of soil CO2 emissions from different semi-arid land-use systems

ABSTRACT Soil carbon dioxide (CO2) fluxes are a critical component in understanding carbon sequestration. In sub-Saharan Africa, empirically measured CO2 emissions data from diverse land-use systems is limited. Soil CO2 emission rates were measured in the Limpopo Province, South Africa for 12 months at two-week intervals in natural systems (forest and shrubland) and commercially managed orchards (avocado and citrus) to establish seasonal dynamics of soil CO2 emissions across these land-use systems. The results showed a variation in emission rates with the variation depending on the season. In the spring and winter, soil CO2 emission rates in citrus were four times higher than in the shrubland due to higher moisture levels. However, in the summer season, the forest emission rates were 40% higher than in citrus due to higher soil organic carbon content. Organic carbon stocks were higher in the forest (1.19 kg/m2) compared to the other land uses. This study revealed differences in soil CO2 emission rates among land-use systems, with the cumulative amount of CO2 emitted over a 12-month period following the order: forest (39.3 tons/ha) > citrus (36.1 tons/ha) > shrubland (28.1 tons/ha) > avocado (26.9 tons/ha). Thus, understanding the emission patterns from various ecosystems can inform strategies for mitigating greenhouse gas emissions.

Impacts of forestry drainage on surface peat stoichiometry and physical properties in boreal peatlands in Finland

Management of drained peatlands may pose a risk or a solution on the way towards climate change mitigation, which creates a need to evaluate the current state of forestry-drained peatlands, the magnitude of degradation processes and indicators for carbon (C) loss. Using a large dataset (778 profiles, 891 peat samples, collected between 1977 and 2017) from peatlands having different fertility classes across Finland, we investigate whether the surface peat profiles of undrained and forestry-drained peatlands differ in C:N, von Post and dry bulk density. The utility of element ratios (C:N:H stoichiometry) as site indicators for degradation were further analyzed from a subsample of 16 undrained and 30 drained peat profiles. This subsample of drained sites had carbon dioxide (CO2) and methane (CH4) fluxes measured allowing us to link peat element ratios to annual C gas effluxes. Element ratios H:C, O:C and C:N and degree of unsaturation (combining C, N, H changes) were found widely valid: they captured both differences in the botanical origin of peat as well as its potential decomposition pathway (C lost via a combination of dissolved organic C and C gas loss and/or the gaseous loss predominantly as CO2). Of the stoichiometric indexes, peat H:C ratio seemed to be the best proxy for degradation following drainage, it indicated not only long-term degradation but also explained 48% of the variation in annual CO2 emission. The O:C ratio positively correlated with annual CH4 flux, presumably because high O:C in peat reflected the availability of easily degradable substrate for methanogenesis. The differences in C:N ratio indicated notable decomposition processes for Sphagnum-dominated peatlands but not in Carex-dominated peatlands. Degree of unsaturation showed potential for an integrative proxy for drainage-induced lowering water table and post-drainage changes in peat substrate quality.

Temporally dynamic carbon dioxide and methane emission factors for rewetted peatlands

Rewetting drained peatlands is recognized as a leading and effective natural solution to curb greenhouse gas emissions. However, rewetting creates novel ecosystems whose emission behaviors are not adequately captured by currently used emission factors. These emission factors are applied immediately after rewetting, thus do not reflect the temporal dynamics of greenhouse gas emissions during the period wherein there is a transition to a rewetted steady-state. Here, we provide long-term data showing a mismatch between actual emissions and default emission factors and revealing the temporal patterns of annual carbon dioxide and methane fluxes in a rewetted peatland site in northeastern Germany. We show that site-level annual emissions of carbon dioxide and methane approach the IPCC default emission factors and those suggested for the German national inventory report only between 13 to 16 years after rewetting. Over the entire study period, we observed a source-to-sink transition of annual carbon dioxide fluxes with a decreasing trend of −0.36 t CO2-C ha−1 yr−1 and a decrease in annual methane emissions of −23.6 kg CH4 ha−1 yr−1. Our results indicate that emission factors should represent the temporally dynamic nature of peatlands post-rewetting and consider the effect of site characteristics to better estimate associated annual emissions.

Inverse problem for retrieving greenhouse gas fluxes at the non-uniform underlying surface from measurements of their concentrations at several levels

The study focuses on the formulation, analysis, and solution of the remote sensing inverse problem to retrieve surface carbon dioxide (CO2) fluxes from measurements of CO2 concentrations at different levels within the atmospheric boundary layer. A three-dimensional hydrodynamic model of turbulent greenhouse gas (GHG) transport was used as a forward model to link the surface GHG fluxes to the drone observations of GHG concentrations. The 3D model provides a GHG concentration distribution by solving the diffusion-advection equation using information on wind speed, its direction, and turbulent exchange coefficients. The surface GHG fluxes are considered as a boundary condition. The spatial distributions of wind speed and turbulence coefficient “for a moment in time” are computed from the relaxation problem for the averaged Navier-Stokes and continuity equations, using a 1.5 order closure scheme (E-ω model). The inverse problem is to retrieve a surface GHG flux by minimizing the difference between the measured and modelled concentrations at several levels. The algorithm was applied to estimate CO2 fluxes over a non-uniform forest canopy at the Roshny-Chu experimental site in the foothills of the Greater Caucasus (Chechen Republic). To test the forward numerical problem, data on surface topography, vegetation height and density, spatial distribution of photosynthetically active solar radiation, as well as data on plant photosynthesis and soil CO2 fluxes were used.

High-resolution spatial patterns and drivers of terrestrial ecosystem carbon dioxide, methane, and nitrous oxide fluxes in the tundra

Abstract. Arctic terrestrial greenhouse gas (GHG) fluxes of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) play an important role in the global GHG budget. However, these GHG fluxes are rarely studied simultaneously, and our understanding of the conditions controlling them across spatial gradients is limited. Here, we explore the magnitudes and drivers of GHG fluxes across fine-scale terrestrial gradients during the peak growing season (July) in sub-Arctic Finland. We measured chamber-derived GHG fluxes and soil temperature, soil moisture, soil organic carbon and nitrogen stocks, soil pH, soil carbon-to-nitrogen (C/N) ratio, soil dissolved organic carbon content, vascular plant biomass, and vegetation type from 101 plots scattered across a heterogeneous tundra landscape (5 km2). We used these field data together with high-resolution remote sensing data to develop machine learning models for predicting (i.e., upscaling) daytime GHG fluxes across the landscape at 2 m resolution. Our results show that this region was on average a daytime net GHG sink during the growing season. Although our results suggest that this sink was driven by CO2 uptake, it also revealed small but widespread CH4 uptake in upland vegetation types, almost surpassing the high wetland CH4 emissions at the landscape scale. Average N2O fluxes were negligible. CO2 fluxes were controlled primarily by annual average soil temperature and biomass (both increase net sink) and vegetation type, CH4 fluxes by soil moisture (increases net emissions) and vegetation type, and N2O fluxes by soil C/N (lower C/N increases net source). These results demonstrate the potential of high spatial resolution modeling of GHG fluxes in the Arctic. They also reveal the dominant role of CO2 fluxes across the tundra landscape but suggest that CH4 uptake in dry upland soils might play a significant role in the regional GHG budget.

The influence of the hydrometeorological factors on the CO2 fluxes from the oligotrophic bog surface.

The influence of the hydrometeorological factors on the CO2 fluxes from the oligotrophic bog surface.

A novel model for measuring carbon dioxide flux on soil surface using an accumulation chamber

Soil surface carbon dioxide (CO2) fluxes hold significant implications for comprehending carbon cycling, global climate change, and ecosystem functioning. Accurate understanding and quantification of this flux are crucial for assessing the carbon dynamics of ecosystems and predicting the impacts of climate change. In this study, we analyzed the variation pattern of CO2 concentration within the accumulated chamber from Fick's law, accounting for diffusion gradient effect and considering leakage effect. We developed a simple model to measure the soil surface CO2 flux, and used the model to calculate the values of CO2 fluxes on a specific soil surface. Comparing the fitting results of the proposed nonlinear model with traditional linear models for the measurement data, we found that the linear model underestimated by at least 18%, even up to 45%. Overall, our research indicates that the proposed model exhibits significant advancement and higher predictive reliability, successfully evaluating soil surface CO2 flux.

Changes in inundation drive carbon dioxide and methane fluxes in a temperate wetland

Wetlands cycle carbon by being net sinks for carbon dioxide (CO2) and net sources of methane (CH4). Daily and seasonal temporal patterns, dissolved oxygen (DO) availability, inundation status (flooded or dry/partially flooded), water depth, and vegetation can affect the magnitude of carbon uptake or emissions, but the extent and interactive effects of these variables on carbon gas fluxes are poorly understood. We characterized the linkages between carbon fluxes and these environmental and temporal drivers at the Old Woman Creek National Estuarine Research Reserve (OWC), OH. We measured diurnal gas flux patterns in an upstream side channel (called the cove) using chamber measurements at six sites (three vegetated and three non-vegetated). We sampled hourly from 7 AM to 7 PM and monthly from July to October 2022. DO concentrations and water levels were measured monthly. Water inundation status had the most influential effect on carbon fluxes with flooded conditions supporting higher CH4 fluxes (0.39 μmol CH4 m−2 s−1; −1.23 μmol CO2 m−2 s−1) and drier conditions supporting higher CO2 fluxes (0.03 μmol CH4 m−2 s−1; 0.86 μmol CO2 m−2 s−1). When flooded, the wetland was a net CO2 sink; however, it became a source for both CH4 and CO2 when water levels were low. We compared chamber-based gas fluxes from the cove in flooded (July) and dry (August) months to fluxes measured with an eddy covariance tower whose footprint covers flooded portions of the wetland. The diurnal pattern of carbon fluxes at the tower did not vary with changing water levels but remained a CO2 sink and a CH4 source even when the cove where we performed the chamber measurements dried out. These results emphasize the role of inundation status on wetland carbon cycling and highlight the importance of fluctuating hydrologic patterns, especially hydrologic drawdowns, under changing climatic conditions.

Contrasting impact of extreme soil and atmospheric dryness on the functioning of trees and forests

Recent studies indicate an increase in the frequency of extreme compound dryness days (days with both extreme soil AND air dryness) across central Europe in the future, with little information on their impact on the functioning of trees and forests. This study aims to quantify and assess the impact of extreme soil dryness, extreme air dryness, and extreme compound dryness on the functioning of trees and forests. For this, >15 years of ecosystem-level (carbon dioxide and water vapor fluxes) and 6–10 years of tree-level measurements (transpiration and growth) each from a montane mixed deciduous forest (CH-Lae) and a subalpine evergreen coniferous forest (CH-Dav) in Switzerland, is used. The results showed extreme air dryness limitation on CO2 fluxes and extreme soil dryness limitations on water vapor fluxes. Additionally, CH-Dav was mainly affected by extreme air dryness whereas CH-Lae was affected by both extreme soil dryness and extreme air dryness. The impact of extreme compound dryness on net CO2 uptake (about 75 % decrease) was more due to higher increased ecosystem respiration (40 % and 70 % increase at CH-Dav and CH-Lae, respectively) than decreased gross primary productivity (10 % and 40 % decrease at CH-Dav and CH-Lae, respectively). A significant negative impact on evapotranspiration and transpiration was only observed at CH-Lae during extreme soil and compound dryness (about 25 % decrease). Furthermore, with some differences, the tree-level impact on tree water deficit, transpiration, and growth were consistent with the ecosystem-level impact on carbon uptake and evapotranspiration. Finally, the impact of extreme dryness showed no significant relationship with tree allometry (diameter and height) but across different tree species. The projected future is likely to expose these forest areas to more extreme and frequent dryness conditions, thus compromising the functioning of trees and forests, thereby calling for management interventions to increase the adaptive capacity and resistance of these forests.

Total Soil CO2 Efflux from Drained Terric Histosols.

Histosols cover about 8-10% of Lithuania's territory and most of this area is covered with nutrient-rich organic soils (Terric Histosols). Greenhouse gas (GHG) emissions from drained Histosols contribute more than 25% of emissions from the Land Use, Land Use Change and Forestry (LULUCF) sector. In this study, as the first step of examining the carbon dioxide (CO2) fluxes in these soils, total soil CO2 efflux and several environmental parameters (temperature of air and topsoil, soil chemical composition, soil moisture, and water table level) were measured in drained Terric Histosols under three native forest stands and perennial grasslands in the growing seasons of 2020 and 2021. The drained nutrient-rich organic soils differed in terms of concentrations of soil organic carbon and total nitrogen, as well as soil organic carbon and total nitrogen ratio. The highest rate of total soil CO2 efflux was found in the summer months. Overall, the rate was statistically significant and strongly correlated only with soil and air temperature. A trend emerged that total soil CO2 efflux was 30% higher in perennial grassland than in forested land. Additional work is still needed to estimate the net CO2 balance of these soils.

Driving and limiting factors of CH4 and CO2 emissions from coastal brackish-water wetlands in temperate regions

Abstract. Coastal wetlands play a fundamental role in mitigating climate change thanks to their ability to store large amounts of organic carbon in the soil. However, degraded freshwater wetlands are also known to be the first natural emitter of methane (CH4). Salinity is known to inhibit CH4 production, but its effect in brackish ecosystems is still poorly understood. This study provides a contribution to understanding how environmental variables may affect greenhouse gas (GHG) emissions in coastal temperate wetlands. We present the results of over 1 year of measurements performed in four wetlands located along a salinity gradient on the northeast Adriatic coast near Ravenna, Italy. Soil properties were determined by coring soil samples, while carbon dioxide (CO2) and CH4 fluxes from soils and standing waters were monitored monthly by a portable gas flux meter. Additionally, water levels and surface and groundwater physical–chemical parameters (temperature, pH, electrical conductivity, and sulfate concentrations of water) were monitored monthly by multiparametric probes. We observed a substantial reduction in CH4 emissions when water depth exceeded the critical threshold of 50 cm. Regardless of the water salinity value, the mean CH4 flux was 5.04 gm-2d-1 in freshwater systems and 12.27 gm-2d-1 in brackish ones. In contrast, when water depth was shallower than 50 cm, CH4 fluxes reached an average of 196.98 gm-2d-1 in freshwater systems, while non-significant results are available for brackish/saline waters. Results obtained for CO2 fluxes showed the same behavior described for CH4 fluxes, even though they were statistically non-significant. Temperature and irradiance strongly influenced CH4 emissions from water and soil, resulting in higher rates during summer and spring.

Impact of Marine Heatwaves on Air‐Sea CO2 Flux Along the US East Coast

AbstractMarine heatwaves (MHWs) are extremely warm ocean temperature events that significantly affect marine environments, but their effects on the coastal carbonate system are still uncertain. In this study, we systematically quantify MHWs' impacts on air‐sea carbon dioxide (CO2) flux anomalies (FCO2′) in the Mid‐Atlantic Bight (MAB) and South Atlantic Bight (SAB) from 1992 to 2020. During the longest MHW in both regions, oceanic CO2 uptake capabilities substantially decreased, primarily due to significant increases in the seawater partial pressure of CO2 (pCO2sea). For all cases, MHWs played a more significant role in driving pCO2sea changes in the MAB than the SAB, where non‐thermal drivers dominated pCO2sea variability. In the MAB, weakened wind speeds related to wintertime atmospheric perturbations increase ocean temperatures and pCO2sea, further reducing CO2 uptake during winter MHWs. This work is the first to connect extreme temperatures to coastal air‐sea CO2 fluxes. The reduction in CO2 absorption noted during MHWs in this study has important implications for coastal regions to act as continued sinks for excess CO2 emissions in the atmosphere.

Estimation Method of Covariance Matrix in Atmospheric Inversion of CO2 Emissions

Atmospheric inversion of CO2 Emissions is based on the correction of prior carbon dioxide flux inventories using concentration monitoring data and atmospheric transport models to obtain posterior carbon dioxide flux. In atmospheric inversion studies, fixed covariance functions are commonly used to generate covariance matrices, and the hyperparameters in the covariance functions are empirically estimated. In this study, we design and implement an ideal experiment based on meteorological data from the central urban area of Zhengzhou, using WRF-STILT to generate sensitivity matrices and construct real carbon emission inventories and prior inventories. Based on the real carbon emission inventories and sensitivity matrices of monitoring stations, simulated observation concentration values are generated. Firstly, based on the observed concentration values, sensitivity matrices of monitoring stations, prior inventories, and constructed covariance matrices, the values of hyperparameters are determined based on maximum marginal likelihood estimation. Then, the influence of different prior covariance functions on the inversion results is tested, and it is found that the prior covariance matrix generated by the balgovind covariance function is most suitable for the experimental data.

Measurements of Carbon Assimilation in Aquatic Systems.

The recent development of an infrared gas analyzer capable of making carbon dioxide flux measurements from aquatic samples has enabled a new sphere of photosynthesis research. This study details key photosynthesis measurements on four aquatic and hydrophytic species, diverse in their morphology, physiology, and habitat. This guide specifies the methods and procedures needed to make reliable and accurate gas exchange measurements, with examples of data correction and presentation.

Maize residue input rather than cover cropping influenced N2O emissions and soil–crop N dynamics during the intercrop and cash crop periods

The well-known benefits of cover cropping for mitigating climate change and improving soil quality could be threatened by increased greenhouse gas emissions. Therefore, such practices could have negative consequences in the context of global climate change. The present year-long field experiment, carried out under irrigated semiarid conditions, evaluated the effect of returning maize crop residues at two topsoil input levels in combination with bare fallow, a cereal (barley, Hordeum vulgare L.) and a legume (vetch, Vicia sativa L.) as cover crops as part of an annual cover crop–cash crop (maize, Zea mays L.) rotation. In addition, control plots, without the addition of nitrogen fertilizer (either in the previous cropping season or during the experiment), were established following the same experimental design in order to assess the carry-over impact of residual nitrogen from prior crops. Nitrous oxide (N2O), methane (CH4) and carbon dioxide (CO2) fluxes were measured together with soil mineral nitrogen and crop yields and nitrogen use efficiency. The significant effect of maize residue input on N2O emissions during the cover crop phase (this was 43% greater for the lower maize residue input than for the higher) was offset after subsequent maize fertilization (there was a 170% increase following the higher maize residue input compared to the lower). Although no differences in soil mineral nitrogen content were reported between treatments before maize fertilization, N2O emissions were 72% higher on the fertilized plots than on the control plots. Before maize fertilization, N2O emissions were increased by the incorporation of cover crop residues and soil rewetting after irrigation, with the highest N2O fluxes being reported for the cereal crop. Cumulative CO2 fluxes were higher in the plots with cover crops (cereal and legume) than in bare fallow plots, while CH4 sink was not affected by cover cropping or maize residue input. The grain yield penalty for maize production (a 7% decrease) with the legume cover crop on the nitrogen fertilized plots (the opposite of what was observed in the control area) highlights the need for an accurate estimation of the nitrogen supply from the mineralization of legume residues. Overall, cover crop and maize residue input did not have any significant effect on annual cumulative N2O emissions, so the positive effects of these practices should be considered to recommend them under semiarid conditions for climate change mitigation and adaptation.

Results of daily measurements of carbon dioxide fluxes in an oligotrophic bog in southern Karelia

Results of daily measurements of carbon dioxide fluxes in an oligotrophic bog in southern Karelia

Long‐term changes in the daytime growing season carbon dioxide exchange following increased temperature and snow cover in arctic tundra

AbstractIncreasing temperatures and winter precipitation can influence the carbon (C) exchange rates in arctic ecosystems. Feedbacks can be both positive and negative, but the net effects are unclear and expected to vary strongly across the Arctic. There is a lack of understanding of the combined effects of increased summer warming and winter precipitation on the C balance in these ecosystems. Here we assess the short‐term (1–3 years) and long‐term (5–8 years) effects of increased snow depth (snow fences) (on average + 70 cm) and warming (open top chambers; 1–3°C increase) and the combination in a factorial design on all key components of the daytime carbon dioxide (CO2) fluxes in a wide‐spread heath tundra ecosystem in West Greenland. The warming treatment increased ecosystem respiration (ER) on a short‐ and long‐term basis, while gross ecosystem photosynthesis (GEP) was only increased in the long term. Despite the difference in the timing of responses of ER and GEP to the warming treatment, the net ecosystem exchange (NEE) of CO2 was unaffected in the short term and in the long term. Although the structural equation model (SEM) indicates a direct relationship between seasonal accumulated snow depth and ER and GEP, there were no significant effects of the snow addition treatment on ER or GEP measured over the summer period. The combination of warming and snow addition turned the plots into net daytime CO2 sources during the growing season. Interestingly, despite no significant changes in air temperature during the snow‐free time during the experiment, control plots as well as warming plots revealed significantly higher ER and GEP in the long term compared to the short term. This was in line with the satellite‐derived time‐integrated normalized difference vegetation index of the study area, suggesting that more factors than air temperature are drivers for changes in arctic tundra ecosystems.

Carbon Mineralization Kinetics and Soil Fertility Induced by Organic and Bio-Fertilization Through Different Incubation Periods

ABSTRACT Agricultural practices such as fertilization considerably influence soil carbon dioxide (CO2) fluxes that cause climate changes. The study aims to test the influence of various organic (vermicompost, VC, and rabbit manure, RM) and bio (Azotobacter and Azospirillum) fertilizers in addition to incubation times on CO2flux, carbon (C) mineralization kinetic, and some soil chemical properties. The dual exponential equation was the best one to describe carbon mineralization kinetics compared to the 1st and 2nd order equations. Carbon half-lives were resided by 18.98, 1.97, 2.00, 1.89, and 2.53 years for control, VC without and with bio, RM without and with bio treatments, respectively. The results exposed that applying VC and RM without or with bio-fertilizer reduced soil reaction (pH) that differ than control treatment after 2 months incubation. Applying VC and RM caused a significant increase in soil salinity (EC) compared to the control. Available nitrogen (N), phosphorus (P) and potassium (K) were significantly enhanced with adding different fertilizers compared to pre-treated soil. It is found that adding organic and bio-fertilizers improved fertility and soil quality.