CN107014951A - Forest ecosystem breathes Carbon flux assay method - Google Patents

Forest ecosystem breathes Carbon flux assay method Download PDF

Info

Publication number
CN107014951A
CN107014951A CN201710108896.9A CN201710108896A CN107014951A CN 107014951 A CN107014951 A CN 107014951A CN 201710108896 A CN201710108896 A CN 201710108896A CN 107014951 A CN107014951 A CN 107014951A
Authority
CN
China
Prior art keywords
measured
soil
forest
isotope
produced
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710108896.9A
Other languages
Chinese (zh)
Inventor
余新晓
李瀚之
贾国栋
樊登星
赵娜
孙佳美
路伟伟
刘自强
张永娥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Forestry University
Original Assignee
Beijing Forestry University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Forestry University filed Critical Beijing Forestry University
Priority to CN201710108896.9A priority Critical patent/CN107014951A/en
Publication of CN107014951A publication Critical patent/CN107014951A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The present invention provides a kind of forest ecosystem breathing Carbon flux assay method, including:Obtain the Atmospheric CO of different height in target forest layer in preset time period2Concentration and Atmospheric CO2In δ13C isotope values, and forest ecosystem breathing δ is obtained according to it13C isotope values;Obtain the CO that tree branches breathing to be measured is produced in target forest in preset time period2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13C isotope values;The CO produced is breathed according to the tree branches to be measured of acquisition2In δ13The CO that C isotope values and soil respiration to be measured are produced2In δ13C isotope values, determine the proportionate relationship of each breathing component in forest ecosystem;Obtain the Carbon flux of soil to be measured;The Carbon flux of the soil to be measured of acquisition is substituted into the proportionate relationship of each breathing component in forest ecosystem, forest ecosystem breathing Carbon flux is obtained.The present invention can more accurately determine the Carbon flux change between air and plant in the forest ecosystem of MODEL OVER COMPLEX TOPOGRAPHY.

Description

Forest ecosystem breathes Carbon flux assay method
Technical field
Breathed the present invention relates to forest ecosystem respiration monitoring technical field, more particularly to a kind of forest ecosystem Carbon flux assay method.
Background technology
Carbon cycle of forest ecosystem is most important in global carbon, has not in regulation atmospheric carbon dioxide Alternative effect.Forest ecosystem year carbon capacity is about 2.4 ± 0.4Pg, account for terrestrial ecosystems carbon capacity (2.6 ± 0.8Pg) 90%, or even than more than the carbon capacity (2.2 ± 0.2Pg) of the global ocean ecosystem.
During global carbon is studied, due to the complicated variety of terrestrial ecosystems, its Carbon flux is difficult to standard Really measure, therefore researcher generally converges the carbon of terrestrial ecosystems and handled as the remainder of Global carbon balance equation. And forest is due to the characteristics of its composition is various, tree crown is tall and big, geometric distortion is big, being all types of ecosystems in land again In it is the most complicated.Therefore, the direct assessment converged for Global Forests system carbon, is still the focus and difficult point studied at present.
At present, in the method for Forest Ecosystem Carbon flux measurement, mainly include:Case method and microclimate method.
Case method is the most commonly used method in current scientific research, and its operation principle is to be covered tested region with casing, Changes of concentrations of carbon dioxide in isolating inner air and outer air, time sight case, in the hope of carbon dioxide exchange amount in case.The letter of its method It is single, it is with low cost.But by being limited by box volume, case method is difficult to be applied to high megaphanerophyte, more at present to be used to measure soil The Carbon flux of breathing or short shrub is determined, and its time precision is relatively low, it is impossible to carry out long-time Continuous Observation.
Microclimate method is the gas turbulence and its changes of concentrations of carbon dioxide close to earth's surface by measurement, and then is derived Carbon dioxide flux between vegetation air.Eddy covariance technique is long-term continuous straight with it as one of representative method of microclimate method Obtain to obtain the advantage that Carbon flux changes under large scale, it has also become determine the standard method of carbon dioxide flux in the world at present.But That eddy covariance technique also has the restrictive condition of implementation, it is desirable to tested region atmospheric condition be in stable state, underlying surface it is uniform, There is no carbon exchange between underlying surface and instrument.And in actual applications, tested region tends not to be under ideal conditions, this Also the application to eddy covariance technique proposes problem.
In consideration of it, how a kind of reliable forest ecosystem breathing Carbon flux assay method is provided, can be intricately The Carbon flux change between air and plant is more accurately determined in the forest ecosystem of shape condition turns into what needs at present were solved Technical problem.
The content of the invention
To solve above-mentioned technical problem, the present invention provides a kind of forest ecosystem breathing Carbon flux assay method, energy Enough Carbon flux changes more accurately determined in the forest ecosystem of MODEL OVER COMPLEX TOPOGRAPHY between air and plant.
In a first aspect, the present invention provides a kind of forest ecosystem breathing Carbon flux assay method, including:
Obtain the Atmospheric CO of different height in target forest layer in preset time period2Concentration and Atmospheric CO2In δ13C Isotope value;
According to the Atmospheric CO of different height in the target forest of acquisition layer2Concentration and Atmospheric CO2In δ13C isotopes Value, obtains forest ecosystem breathing δ13C isotope values;
Obtain the CO that tree branches breathing to be measured is produced in target forest in preset time period2In δ13C isotope values and The CO that soil respiration to be measured is produced in target forest2In δ13C isotope values;
The CO produced is breathed according to the tree branches to be measured of acquisition2In δ13C isotope values and soil respiration to be measured are produced CO2In δ13C isotope values, determine the proportionate relationship of each breathing component in the forest ecosystem;
Obtain the Carbon flux of the soil to be measured;
The Carbon flux of the soil to be measured of acquisition is substituted into the proportionate relationship of each breathing component in the forest ecosystem, Obtain forest ecosystem breathing Carbon flux.
Alternatively, the Atmospheric CO for obtaining different height in target forest layer in preset time period2Concentration and air CO2In δ13C isotope values, including:
Using carbon isotope analyser coordinate Port Multiplier, in preset time period target forest layer in different height it is big Gas CO2Concentration and Atmospheric CO2In δ13C isotope values carry out layering measure;
Wherein, the carbon isotope analyser is connected by the Port Multiplier with multiple air sampling probes, described big Gas sampling probe is set in advance in different height in forest to be measured layer, the air sampling probe by gas collection pipeline with The Port Multiplier is connected.
Alternatively, Port Multiplier is coordinated to carry out layering survey to different height in target forest layer using carbon isotope analyser The fixed time includes:Buffer stabilization time and practical measurement time;
The buffering stabilization time, the practical measurement time, the buffering stabilization time was less than the reality rear preceding Border minute.
Alternatively, it is described to obtain the CO that tree branches breathing to be measured is produced in target forest in preset time period2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13C isotope values, including:
Coordinate Port Multiplier using carbon isotope analyser, in preset time period, respectively to being treated in target forest The branch air chamber and the soil air chamber on soil to be measured in target forest surveyed in tree branches are measured, and obtain pre- If the CO that tree branches breathing to be measured is produced in target forest in the period2In δ13It is to be measured in C isotope values and target forest The CO that soil respiration is produced2In δ13C isotope values;
Wherein, the carbon isotope analyser by the Port Multiplier respectively with the branch air chamber and the SOIL GAS Room is connected, and the branch air chamber and the soil air chamber are connected by gas collection pipeline with the Port Multiplier.
Alternatively, using carbon isotope analyser coordinate Port Multiplier branch air chamber and soil air chamber are measured when Between include:Buffer stabilization time and practical measurement time;
The buffering stabilization time, the practical measurement time, the buffering stabilization time was less than the reality rear preceding Border minute.
Alternatively, in the target forest layer according to acquisition different height Atmospheric CO2Concentration and Atmospheric CO2In δ13C isotope values, obtain forest ecosystem breathing δ13C isotope values, including:
According to the Atmospheric CO of different height in the target forest of acquisition layer2Concentration and Atmospheric CO2In δ13C isotopes Value, using keeling curve methods, draws Atmospheric CO2In δ13C isotope values and Atmospheric CO2Concentration scatterplot reciprocal Figure;
Using linear equation y=ax+b, the scatter diagram is fitted, forest ecosystem breathing δ is obtained13The same positions of C Plain value δ13Ceco
Wherein, y is δ13C isotope values, x is Atmospheric CO2Concentration inverse, b be forest ecosystem breathe δ13C is same The plain value δ in position13Ceco, a is equation parameter.
Alternatively, the CO that tree branches breathing to be measured is produced in target forest in the preset time period2In δ13C is same The CO that soil respiration to be measured is produced in position element value and target forest2In δ13C isotope values are by repeatedly measuring averaging What value was obtained.
Alternatively, it is described that the CO produced is breathed according to the tree branches to be measured of acquisition2In δ13C isotope values and to be measured The CO that soil respiration is produced2In δ13C isotope values, determine the proportionate relationship of each breathing component in the forest ecosystem, Including:
The CO produced is breathed according to the tree branches to be measured of acquisition2In δ13C isotope values and soil respiration to be measured are produced CO2In δ13C isotope values, by isotope conservation principle and the principle of mass conservation, obtain and are respectively exhaled in forest ecosystem Inhale the proportionate relationship of component;
Wherein, the isotope conservation principle is:
δ13Ceco×Reco13Ctree×Rtree13Csoil×Rsoil,
RecoCarbon flux, δ are breathed for forest ecosystem13Cecoδ is breathed for forest ecosystem13C isotope values, Rtree Carbon flux, δ are breathed for tree branches to be measured13CtreeThe CO produced is breathed for tree branches to be measured2In δ13C isotope values, RsoilFor soil respiration Carbon flux to be measured, δ13CsoilThe CO produced for soil respiration to be measured2In δ13C isotope values;
The principle of mass conservation is:
Reco=Rtree+Rsoil
The proportionate relationship f of each breathing component in the forest ecosystem of acquisitionsoilFor:
Alternatively, the Carbon flux for obtaining the soil to be measured, including:
Using the soil carbon flux measuring system being connected with the soil air chamber, the carbon for obtaining the soil to be measured leads to Amount.
As shown from the above technical solution, forest ecosystem of the invention breathing Carbon flux assay method, pre- by obtaining If in the period target forest layer in different height Atmospheric CO2Concentration and Atmospheric CO2In δ13C isotope values, and according to It obtains forest ecosystem breathing δ13C isotope values, obtain tree branches breathing to be measured in target forest in preset time period The CO of generation2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13C isotope values, and According to the proportionate relationship of each breathing component in its determination forest ecosystem, the Carbon flux of soil to be measured is obtained, by treating for acquisition The Carbon flux for surveying soil substitutes into the proportionate relationship of each breathing component in forest ecosystem, obtains forest ecosystem breathing carbon Flux, thus, can not only overcome case method is short, can not carry out Continuous Observation to forest ecosystem breathing Carbon flux to lack Point, additionally it is possible to overcome microclimate method to may not apply to the restrictive condition of complicated landform, can give birth in the forest of MODEL OVER COMPLEX TOPOGRAPHY The Carbon flux more accurately determined in state system between air and plant changes, and technology is provided for Global Scale carbon cycle process Support.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will to embodiment or The accompanying drawing used required in description of the prior art is briefly described, it should be apparent that, drawings in the following description are this hairs Some bright embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, can be with Other accompanying drawings are obtained according to these accompanying drawings.
A kind of flow for forest ecosystem breathing Carbon flux assay method that Fig. 1 provides for one embodiment of the invention is shown It is intended to;
Fig. 2 breathes the device that Carbon flux assay method is utilized by forest ecosystem provided in an embodiment of the present invention Structural representation.
In figure, 1 is air sampling probe, and 2 be branch air chamber, and 3 be soil air chamber, and 4 be soil carbon flux measuring system, 5 It is carbon isotope analyser for Port Multiplier, 6,7 be bioassay standard wood.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention In accompanying drawing, clear, complete description is carried out to the technical scheme in the embodiment of the present invention, it is clear that described embodiment is only Only it is a part of embodiment of the invention, rather than whole embodiments.Based on embodiments of the invention, ordinary skill The every other embodiment that personnel are obtained under the premise of creative work is not made, belongs to the model that the present invention is protected Enclose.
Fig. 1 shows that the flow for the forest ecosystem breathing Carbon flux assay method that one embodiment of the invention is provided is shown It is intended to, as shown in figure 1, the forest ecosystem breathing Carbon flux assay method of the present embodiment is as described below.
101st, the Atmospheric CO of different height in target forest layer in preset time period is obtained2Concentration and Atmospheric CO2In δ13C isotope values.
In a particular application, the step 101 can utilize the carbon isotope analyser 6 in Fig. 2 to coordinate Port Multiplier 5, right In preset time period target forest layer in different height Atmospheric CO2Concentration and Atmospheric CO2In δ13C isotope values are carried out Layering is determined;
Wherein, the carbon isotope analyser 6 is connected by the Port Multiplier 5 with multiple air sampling probes 1, institute State air sampling probe 1 and be set in advance in different height in forest layer to be measured, the air sampling probe 1 passes through gas collection Pipeline is connected with the Port Multiplier 5.
In a particular application, the carbon isotope analyser 6 can be LGR carbon isotope analysers.
Specifically, the time that different height in target forest layer be layered measure can be included:When buffering is stable Between and the practical measurement time;
The buffering stabilization time, the practical measurement time, the buffering stabilization time was less than the reality rear preceding Border minute.For example, the time to every layer of measure of different height in target forest layer can be 5 minutes, wherein, first 2 minutes It it is within latter 3 minutes the practical measurement time to buffer stabilization time.
102nd, according to the Atmospheric CO of different height in the target forest of acquisition layer2Concentration and Atmospheric CO2In δ13C is same Position element value, obtains forest ecosystem breathing δ13C isotope values.
In a particular application, the step 102 can be specifically included:
According to the Atmospheric CO of different height in the target forest of acquisition layer2Concentration and Atmospheric CO2In δ13C isotopes Value, using keeling curve methods, draws Atmospheric CO2In δ13C isotope values and Atmospheric CO2Concentration scatter diagram reciprocal (e.g., can be drawn in excel softwares);
Using linear equation y=ax+b, the scatter diagram is fitted, forest ecosystem breathing δ is obtained13The same positions of C Plain value δ13Ceco
Wherein, y is δ13C isotope values, x is Atmospheric CO2Concentration inverse, b be forest ecosystem breathe δ13C is same The plain value δ in position13Ceco, a is equation parameter.
103rd, trees (i.e. bioassay standard wood 7) branch breathing generation to be measured in target forest in preset time period is obtained CO2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13C isotope values.
In a particular application, the step 103 can utilize the carbon isotope analyser 6 in Fig. 2 to coordinate Port Multiplier 5, In preset time period, respectively to the branch air chamber 2 in target forest in tree branches to be measured and installed in target forest In soil air chamber 3 on soil to be measured be measured, obtain tree branches breathing production to be measured in target forest in preset time period Raw CO2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13C isotope values;
Wherein, the carbon isotope analyser 6 by the Port Multiplier 5 respectively with the branch air chamber 2 and the soil Air chamber 3 is connected, and the branch air chamber 2 and the soil air chamber 3 are connected by gas collection pipeline with the Port Multiplier 5 Connect.
Specifically, the time that branch air chamber 2 and soil air chamber 3 are measured can be included:Buffer stabilization time and The practical measurement time;
The buffering stabilization time, the practical measurement time, the buffering stabilization time was less than the reality rear preceding Border minute.For example, the time being measured to branch air chamber 2 and soil air chamber 3 can be 5 minutes, wherein, first 2 points Clock is buffering stabilization time, is within latter 3 minutes the practical measurement time.
Further, in order to obtain in more accurately measurement result, the preset time period trees to be measured in target forest The CO that branch breathing is produced2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13The same positions of C Plain value can be by repeatedly measuring acquisition of averaging.
104th, the CO produced is breathed according to the tree branches to be measured of acquisition2In δ13C isotope values and soil respiration to be measured The CO of generation2In δ13C isotope values, determine in the forest ecosystem each breathing component (including:Soil respiration and plant Thing breathe) proportionate relationship.
In a particular application, the step 104 can be specifically included:
The CO produced is breathed according to the tree branches to be measured of acquisition2In δ13C isotope values and soil respiration to be measured are produced CO2In δ13C isotope values, by isotope conservation principle and the principle of mass conservation, obtain and are respectively exhaled in forest ecosystem Inhale the proportionate relationship f of componentsoil
Wherein, the isotope conservation principle is:
δ13Ceco×Reco13Ctree×Rtree13Csoil×Rsoil (1)
RecoCarbon flux, δ are breathed for forest ecosystem13Cecoδ is breathed for forest ecosystem13C isotope values, Rtree Carbon flux, δ are breathed for tree branches to be measured13CtreeThe CO produced is breathed for tree branches to be measured2In δ13C isotope values, RsoilFor soil respiration Carbon flux to be measured, δ13CsoilThe CO produced for soil respiration to be measured2In δ13C isotope values;
The principle of mass conservation is:
Reco=Rtree+Rsoil (2)
The proportionate relationship f of each breathing component in the forest ecosystem of acquisitionsoilFor:
Wherein, the unit of each component breathing Carbon flux is μm olm in forest ecosystem-2·s-1, isotope abundance value Unit be ‰.
105th, the Carbon flux of the soil to be measured is obtained.
In a particular application, the step 105 can utilize the soil carbon being connected with the soil air chamber 3 in Fig. 2 Flux measurement system 4, obtains the Carbon flux of the soil to be measured.
For example, the soil carbon flux measuring system 4 can be LI-8100 soil carbon flux measuring systems.
106th, the ratio that the Carbon flux of the soil to be measured of acquisition is substituted into each breathing component in the forest ecosystem is closed System, obtains forest ecosystem breathing Carbon flux.
It is understood that according to the proportionate relationship of each breathing component in the forest ecosystem:Understand,AndThe carbon of soil to be measured leads to Measure RsoilObtain in step 105, forest ecosystem breathing Carbon flux R can be obtained by calculatingeco
The forest ecosystem breathing Carbon flux assay method of the present embodiment, it is gloomy by obtaining target in preset time period The Atmospheric CO of different height in woods layer2Concentration and Atmospheric CO2In δ13C isotope values, and Forest ecosystem is obtained according to it System breathing δ13C isotope values, obtain the CO that tree branches breathing to be measured is produced in target forest in preset time period2In δ13C The CO that soil respiration to be measured is produced in isotope value and target forest2In δ13C isotope values, and given birth to according to its determination forest The proportionate relationship of each breathing component, obtains the Carbon flux of soil to be measured, by the Carbon flux generation of the soil to be measured of acquisition in state system Enter the proportionate relationship of each breathing component in forest ecosystem, obtain forest ecosystem breathing Carbon flux, that is, use and stablize same Each component during the plain technology in position is breathed to forest ecosystem (including:Soil respiration and plant respiration) split and determine it Proportionate relationship, is then derived by determining soil respiration Carbon flux and obtains forest ecosystem breathing Carbon flux, thus, no Can only overcome case method it is short, can not to forest ecosystem breathe Carbon flux carry out Continuous Observation shortcoming, additionally it is possible to overcome Microclimate method may not apply to the restrictive condition of complicated landform, can be in the forest ecosystem of MODEL OVER COMPLEX TOPOGRAPHY more The Carbon flux accurately determined between air and plant changes, and technical support is provided for Global Scale carbon cycle process.
In order to preferably illustrate, Capital Circle ecological observation station in Haidian District, Beijing City is located at enforcement place below, Exemplified by April 1st, 2016 was measured to April 28, the forest ecosystem breathing Carbon flux assay method of the present embodiment can With including:
First, build air sampling probe 1 in forest as shown in Figure 2, be laid in respectively 0.05m in forest, 2m, 5m, 8m, 12m and 18m height, and branch air chamber and soil air chamber are installed, gas collection pipeline is connected to carbon isotope via Port Multiplier 5 Analyzer 6, determines the Atmospheric CO of different height in forest layer2Concentration and Atmospheric CO2In δ13C isotope values, and The CO that tree branches breathing to be measured is produced in forest2In δ13The CO that soil respiration to be measured is produced in C isotope values and forest2 In δ13C isotope values, and laying LI-8100 soil carbon fluxes measuring system 4 is installed in forest, determine soil to be measured Carbon flux.
Specifically, the carbon isotope analyser 6 can carry out instrument calibration with multi-point calibration method, and Port Multiplier 5 is each Passage is accessed 5 minutes, wherein, it is within first 2 minutes buffering stabilization time, is within latter 3 minutes the practical measurement time.
2nd, according to the Atmospheric CO of different height in the target forest of acquisition layer2Concentration and Atmospheric CO2In δ13The same positions of C Element value, using keeling curve methods, Atmospheric CO is drawn in excel softwares2In δ13C isotope values and Atmospheric CO2It is dense The scatter diagram reciprocal of degree;Using linear equation y=ax+b, the scatter diagram is fitted, forest ecosystem is obtained and exhales Inhale δ13C isotope values δ13Ceco;Wherein, y is δ13C isotope values, x is Atmospheric CO2Concentration inverse, b is forest ecology System breathes δ13C isotope values δ13Ceco, a is equation parameter.During numerical fitting, if Atmospheric CO2Change in concentration is smaller, Then fitting effect is poor, therefore in the present embodiment, what is gained atmosphere data carried out in once fitting, linear equation within every 7 days cuts It is forest ecosystem breathing δ away from b13C isotope values δ13Ceco;The CO produced is breathed for tree branches to be measured2In δ13C isotope values δ13CtreeThe CO produced with soil respiration to be measured2In δ13C isotope values δ13Csoil, carry out once within every 7 days Average value processing.The simulation of the present embodiment concrete numerical value is as shown in table 1 below with observed result.
Table 1
Time δ13Ctree δ13Csoil δ13Ceco
4.1-4.7 -27.65 -24.35 -26.45
4.8-4.14 -29.12 -25.35 -27.66
4.15-4.21 -28.11 -25.67 -26.98
4.22-4.28 -28.95 -26.13 -27.15
3rd, the CO produced is breathed according to the tree branches to be measured of acquisition2In δ13C isotope values and soil respiration to be measured The CO of generation2In δ13C isotope values, by isotope conservation principle and the principle of mass conservation, are obtained in forest ecosystem The proportionate relationship f of each breathing componentsoil
Wherein, the isotope conservation principle is:
δ13Ceco×Reco13Ctree×Rtree13Csoil×Rsoil (1)
RecoCarbon flux, δ are breathed for forest ecosystem13Cecoδ is breathed for forest ecosystem13C isotope values, Rtree Carbon flux, δ are breathed for tree branches to be measured13CtreeThe CO produced is breathed for tree branches to be measured2In δ13C isotope values, RsoilFor soil respiration Carbon flux to be measured, δ13CsoilThe CO produced for soil respiration to be measured2In δ13C isotope values;
The principle of mass conservation is:
Reco=Rtree+Rsoil (2)
The proportionate relationship f of each breathing component in the forest ecosystem of acquisitionsoilFor:
Wherein, the unit of each component breathing Carbon flux is μm olm in forest ecosystem-2·s-1, isotope abundance The unit of value is ‰.
Data in the table 1 determined by April 1 to April 28, this step can obtain April 1 to the date in April 28 Between, the f of every 7 dayssoilRespectively:0.467 (1-7 days April), 0.493 (8-14 days April), 0.504 (15-21 days April) and 0.461 (22-28 days April).
4th, using the LI-8100 soil carbon fluxes measuring system 4 being connected with the soil air chamber 3, treated described in acquisition The Carbon flux of soil is surveyed, during obtaining April 1 to April 28, the Carbon flux of every soil to be measured on the 7th is respectively:5.62g·m-2、5.87g·m-2、6.29g·m-2And 6.61gm-2
5th, according to formulaUnderstand forest ecosystem breathing Carbon fluxApril 1 can be obtained During day to 28 days, the forest ecosystem breathing Carbon flux R of every 7 daysecoRespectively 12.04gm-2、11.90g·m-2、 12.47g·m-2And 14.34gm-2
The present embodiment forest ecosystem breathing Carbon flux assay method, can not only overcome case method it is short, can not be right Forest ecosystem breathing Carbon flux carries out the shortcoming of Continuous Observation, additionally it is possible to overcome microclimate method to may not apply to intricately The restrictive condition of shape, can more accurately be determined between air and plant in the forest ecosystem of MODEL OVER COMPLEX TOPOGRAPHY Carbon flux changes, and technical support is provided for Global Scale carbon cycle process.
One of ordinary skill in the art will appreciate that:Realizing all or part of step of above method embodiment can lead to The related hardware of programmed instruction is crossed to complete.Foregoing program can be stored in a computer read/write memory medium.The journey Sequence upon execution, performs the step of including above-mentioned each method embodiment;And foregoing storage medium includes:ROM, RAM, magnetic disc Or CD etc. is various can be with the medium of store program codes.
Finally it should be noted that:Various embodiments above is merely illustrative of the technical solution of the present invention, rather than its limitations; Although the present invention is described in detail with reference to foregoing embodiments, it will be understood by those within the art that:Its The technical scheme described in foregoing embodiments can still be modified, or it is special to which part or whole technologies Levy carry out equivalent substitution;And these modifications or replacement, the essence of appropriate technical solution is departed from the right of the present invention Claimed scope.

Claims (9)

1. a kind of forest ecosystem breathes Carbon flux assay method, it is characterised in that including:
Obtain the Atmospheric CO of different height in target forest layer in preset time period2Concentration and Atmospheric CO2In δ13C isotopes Value;
According to the Atmospheric CO of different height in the target forest of acquisition layer2Concentration and Atmospheric CO2In δ13C isotope values, are obtained Forest ecosystem breathes δ13C isotope values;
Obtain the CO that tree branches breathing to be measured is produced in target forest in preset time period2In δ13C isotope values and target The CO that soil respiration to be measured is produced in forest2In δ13C isotope values;
The CO produced is breathed according to the tree branches to be measured of acquisition2In δ13The CO that C isotope values and soil respiration to be measured are produced2 In δ13C isotope values, determine the proportionate relationship of each breathing component in the forest ecosystem;
Obtain the Carbon flux of the soil to be measured;
The Carbon flux of the soil to be measured of acquisition is substituted into the proportionate relationship of each breathing component in the forest ecosystem, obtained gloomy Woods ecosystem respiration Carbon flux.
2. according to the method described in claim 1, it is characterised in that different in target forest layer in the acquisition preset time period The Atmospheric CO of height2Concentration and Atmospheric CO2In δ13C isotope values, including:
Coordinate Port Multiplier using carbon isotope analyser, to the Atmospheric CO of different height in target forest layer in preset time period2's Concentration and Atmospheric CO2In δ13C isotope values carry out layering measure;
Wherein, the carbon isotope analyser is connected by the Port Multiplier with multiple air sampling probes, and the air takes Sample probe is set in advance in different height in forest to be measured layer, the air sampling probe by gas collection pipeline with it is described many Road device is connected.
3. method according to claim 2, it is characterised in that coordinate Port Multiplier gloomy to target using carbon isotope analyser The time of different height progress layering measure includes in woods layer:Buffer stabilization time and practical measurement time;
The buffering stabilization time, the practical measurement time, the buffering stabilization time was less than described actual survey rear preceding Fix time.
4. according to the method described in claim 1, it is characterised in that treat assize in target forest in the acquisition preset time period The CO that the breathing of wooden branch bar is produced2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13C is same Position element value, including:
Coordinate Port Multiplier using carbon isotope analyser, in preset time period, respectively to treating assize in target forest Branch air chamber on wooden branch bar and the soil air chamber on soil to be measured in target forest are measured, and obtain preset time The CO that tree branches breathing to be measured is produced in target forest in section2In δ13Soil to be measured is exhaled in C isotope values and target forest Inhale the CO produced2In δ13C isotope values;
Wherein, the carbon isotope analyser is connected with the branch air chamber and the soil air chamber respectively by the Port Multiplier Connect, the branch air chamber and the soil air chamber are connected by gas collection pipeline with the Port Multiplier.
5. method according to claim 4, it is characterised in that coordinate Port Multiplier to branch gas using carbon isotope analyser The time that room and soil air chamber are measured includes:Buffer stabilization time and practical measurement time;
The buffering stabilization time, the practical measurement time, the buffering stabilization time was less than described actual survey rear preceding Fix time.
6. according to the method described in claim 1, it is characterised in that different height in the target forest layer according to acquisition Atmospheric CO2Concentration and Atmospheric CO2In δ13C isotope values, obtain forest ecosystem breathing δ13C isotope values, including:
According to the Atmospheric CO of different height in the target forest of acquisition layer2Concentration and Atmospheric CO2In δ13C isotope values, are utilized Keeling curve methods, draw Atmospheric CO2In δ13C isotope values and Atmospheric CO2Concentration scatter diagram reciprocal;
Using linear equation y=ax+b, the scatter diagram is fitted, forest ecosystem breathing δ is obtained13C isotope values δ13Ceco
Wherein, y is δ13C isotope values, x is Atmospheric CO2Concentration inverse, b be forest ecosystem breathe δ13C isotope values δ13Ceco, a is equation parameter.
7. according to the method described in claim 1, it is characterised in that trees branch to be measured in target forest in the preset time period The CO that bar breathing is produced2In δ13The CO that soil respiration to be measured is produced in C isotope values and target forest2In δ13C isotopes Value is by repeatedly measuring acquisition of averaging.
8. according to the method described in claim 1, it is characterised in that described to breathe what is produced according to the tree branches to be measured of acquisition CO2In δ13The CO that C isotope values and soil respiration to be measured are produced2In δ13C isotope values, determine the forest ecosystem In each breathing component proportionate relationship, including:
The CO produced is breathed according to the tree branches to be measured of acquisition2In δ13The CO that C isotope values and soil respiration to be measured are produced2 In δ13C isotope values, by isotope conservation principle and the principle of mass conservation, obtain each breathing component in forest ecosystem Proportionate relationship;
Wherein, the isotope conservation principle is:
δ13Ceco×Reco13Ctree×Rtree13Csoil×Rsoil,
RecoCarbon flux, δ are breathed for forest ecosystem13Cecoδ is breathed for forest ecosystem13C isotope values, RtreeTo be to be measured Tree branches breathe Carbon flux, δ13CtreeThe CO produced is breathed for tree branches to be measured2In δ13C isotope values, RsoilTo treat Survey soil respiration Carbon flux, δ13CsoilThe CO produced for soil respiration to be measured2In δ13C isotope values;
The principle of mass conservation is:
Reco=Rtree+Rsoil
The proportionate relationship f of each breathing component in the forest ecosystem of acquisitionsoilFor:
9. method according to claim 4, it is characterised in that the Carbon flux of the acquisition soil to be measured, including:
Using the soil carbon flux measuring system being connected with the soil air chamber, the Carbon flux of the soil to be measured is obtained.
CN201710108896.9A 2017-02-27 2017-02-27 Forest ecosystem breathes Carbon flux assay method Pending CN107014951A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710108896.9A CN107014951A (en) 2017-02-27 2017-02-27 Forest ecosystem breathes Carbon flux assay method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710108896.9A CN107014951A (en) 2017-02-27 2017-02-27 Forest ecosystem breathes Carbon flux assay method

Publications (1)

Publication Number Publication Date
CN107014951A true CN107014951A (en) 2017-08-04

Family

ID=59440558

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710108896.9A Pending CN107014951A (en) 2017-02-27 2017-02-27 Forest ecosystem breathes Carbon flux assay method

Country Status (1)

Country Link
CN (1) CN107014951A (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108776051A (en) * 2018-05-11 2018-11-09 山西师范大学 A kind of soil, plant evaporation steam isotope harvester
CN109212166A (en) * 2018-11-12 2019-01-15 福建师范大学 It is a kind of for studying the monitoring device of illumination Yu soil carbon flux relationship
CN111080173A (en) * 2019-12-31 2020-04-28 黑龙江工程学院 Estimation method of carbon flux of forest system
CN111257505A (en) * 2019-09-20 2020-06-09 浙江农林大学 Forest carbon flux metering system and method
CN111272948A (en) * 2020-02-12 2020-06-12 北京市环境保护科学研究院 Carbon dioxide flux testing device and method for representing natural attenuation process of pollutants
CN112162061A (en) * 2020-09-17 2021-01-01 中山大学 Evapotranspiration component space measuring and calculating method based on hydrogen-oxygen stable isotope observation
CN113358826A (en) * 2021-06-10 2021-09-07 北京林业大学 Device and method for measuring tree branch and leaf respiration
CN113945681A (en) * 2020-07-17 2022-01-18 中国科学院沈阳应用生态研究所 Survey forest ecosystem component respiratory delta13C device and method
CN114544911A (en) * 2022-02-14 2022-05-27 北京市农林科学院 Method and device for determining organic carbon delivery from plants to soil based on different ways
CN114814184A (en) * 2022-05-10 2022-07-29 中国科学院城市环境研究所 Method for measuring degradation rate of degradable plastic based on carbon 13 isotope method
CN116840002A (en) * 2023-06-07 2023-10-03 南方海洋科学与工程广东省实验室(广州) Analysis method of carbon dioxide emission source of mangrove seedling sediment-atmosphere interface

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102405404A (en) * 2009-02-02 2012-04-04 行星排放管理公司 System of systems for monitoring greenhouse gas fluxes
CN105842180A (en) * 2016-05-10 2016-08-10 北京林业大学 Device and method for determining soil respiration and carbon isotopes
CN105973817A (en) * 2016-05-10 2016-09-28 北京林业大学 Device and method for determining trunk respiration and 13C thereof
CN205749509U (en) * 2016-07-01 2016-11-30 中国林业科学研究院森林生态环境与保护研究所 A kind of device observing forest ecosystem soil respiration process

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102405404A (en) * 2009-02-02 2012-04-04 行星排放管理公司 System of systems for monitoring greenhouse gas fluxes
CN105842180A (en) * 2016-05-10 2016-08-10 北京林业大学 Device and method for determining soil respiration and carbon isotopes
CN105973817A (en) * 2016-05-10 2016-09-28 北京林业大学 Device and method for determining trunk respiration and 13C thereof
CN205749509U (en) * 2016-07-01 2016-11-30 中国林业科学研究院森林生态环境与保护研究所 A kind of device observing forest ecosystem soil respiration process

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孙伟 等: "稳定性同位素技术与Keeling曲线法在陆地生态系统碳/水交换研究中的应用,", 《植物生态学报》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108776051A (en) * 2018-05-11 2018-11-09 山西师范大学 A kind of soil, plant evaporation steam isotope harvester
CN108776051B (en) * 2018-05-11 2024-05-03 山西师范大学 Soil and plant evaporation water vapor isotope acquisition device
CN109212166B (en) * 2018-11-12 2024-04-30 福建师范大学 Monitoring device for researching relationship between illumination and soil carbon flux
CN109212166A (en) * 2018-11-12 2019-01-15 福建师范大学 It is a kind of for studying the monitoring device of illumination Yu soil carbon flux relationship
CN111257505A (en) * 2019-09-20 2020-06-09 浙江农林大学 Forest carbon flux metering system and method
CN111080173A (en) * 2019-12-31 2020-04-28 黑龙江工程学院 Estimation method of carbon flux of forest system
CN111080173B (en) * 2019-12-31 2020-11-03 黑龙江工程学院 Estimation method of carbon flux of forest system
CN111272948A (en) * 2020-02-12 2020-06-12 北京市环境保护科学研究院 Carbon dioxide flux testing device and method for representing natural attenuation process of pollutants
CN113945681A (en) * 2020-07-17 2022-01-18 中国科学院沈阳应用生态研究所 Survey forest ecosystem component respiratory delta13C device and method
CN112162061A (en) * 2020-09-17 2021-01-01 中山大学 Evapotranspiration component space measuring and calculating method based on hydrogen-oxygen stable isotope observation
CN112162061B (en) * 2020-09-17 2021-08-03 中山大学 Evapotranspiration component space measuring and calculating method based on hydrogen-oxygen stable isotope observation
CN113358826A (en) * 2021-06-10 2021-09-07 北京林业大学 Device and method for measuring tree branch and leaf respiration
CN114544911B (en) * 2022-02-14 2024-03-29 北京市农林科学院 Method and device for determining organic carbon delivery to soil based on plants in different ways
CN114544911A (en) * 2022-02-14 2022-05-27 北京市农林科学院 Method and device for determining organic carbon delivery from plants to soil based on different ways
CN114814184A (en) * 2022-05-10 2022-07-29 中国科学院城市环境研究所 Method for measuring degradation rate of degradable plastic based on carbon 13 isotope method
CN116840002A (en) * 2023-06-07 2023-10-03 南方海洋科学与工程广东省实验室(广州) Analysis method of carbon dioxide emission source of mangrove seedling sediment-atmosphere interface
CN116840002B (en) * 2023-06-07 2024-01-09 南方海洋科学与工程广东省实验室(广州) Analysis method of carbon dioxide emission source of mangrove seedling sediment-atmosphere interface

Similar Documents

Publication Publication Date Title
CN107014951A (en) Forest ecosystem breathes Carbon flux assay method
Rothfuss et al. Reviews and syntheses: Isotopic approaches to quantify root water uptake: a review and comparison of methods
Ma et al. An inorganic CO2 diffusion and dissolution process explains negative CO2 fluxes in saline/alkaline soils
Zhao et al. Combining LPJ-GUESS and HASM to simulate the spatial distribution of forest vegetation carbon stock in China
Winton Do climate models underestimate the sensitivity of Northern Hemisphere sea ice cover?
Théroux‐Rancourt et al. The light response of mesophyll conductance is controlled by structure across leaf profiles
Yang et al. Climatic reconstruction at the Miocene Shanwang basin, China, using leaf margin analysis, CLAMP, coexistence approach, and overlapping distribution analysis
Roden et al. The enigma of effective path length for 18O enrichment in leaf water of conifers
Meunier et al. Measuring and modeling hydraulic lift of Lolium multiflorum using stable water isotopes
Foereid et al. The sensitivity of carbon turnover in the Community Land Model to modified assumptions about soil processes
CN105547933B (en) A kind of air monitoring method and apparatus
CN105675468A (en) Evaluation method of shale organic matter porosity
Blunier et al. Planetary fertility during the past 400 ka based on the triple isotope composition of O 2 in trapped gases from the Vostok ice core
JP4928746B2 (en) Estimating the composition ratio of carbon dioxide emissions
CN106885892A (en) The method and device split to forest ecosystem evapotranspiration
Pekin et al. Plant diversity is linked to nutrient limitation of dominant species in a world biodiversity hotspot
Schrader et al. Non-stomatal exchange in ammonia dry deposition models: comparison of two state-of-the-art approaches
Hayes et al. Thorium distributions in high‐and low‐dust regions and the significance for iron supply
Wilding et al. Spatial variability: enhancing the mean estimate of organic and inorganic carbon in a sampling unit
Feig et al. Soil biogenic emissions of nitric oxide from a semi-arid savanna in South Africa
Ferretti et al. Variables related to nitrogen deposition improve defoliation models for European forests
Huang et al. Changes and net ecosystem productivity of terrestrial ecosystems and their influencing factors in China from 2000 to 2019
Chen et al. Evapotranspiration partitioning based on field‐stable oxygen isotope observations for an urban locust forest land
Yan et al. Assessing productivity and carbon sequestration capacity of subtropical coniferous plantations using the process model PnET-CN
Ishidoya et al. Secular change in atmospheric Ar∕ N 2 and its implications for ocean heat uptake and Brewer–Dobson circulation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170804