CN115758801A - High-resolution weather-driven planning carbon sink numerical evaluation method, system and terminal - Google Patents
High-resolution weather-driven planning carbon sink numerical evaluation method, system and terminal Download PDFInfo
- Publication number
- CN115758801A CN115758801A CN202211584270.2A CN202211584270A CN115758801A CN 115758801 A CN115758801 A CN 115758801A CN 202211584270 A CN202211584270 A CN 202211584270A CN 115758801 A CN115758801 A CN 115758801A
- Authority
- CN
- China
- Prior art keywords
- meteorological
- land
- data
- model
- vegetation
- 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
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a high-resolution weather-driven planning carbon sink numerical evaluation method, a system and a terminal, belonging to the technical field of numerical simulation, wherein the method comprises the following steps: generating high-resolution meteorological driving data with the spatial resolution of 1km-5 km; and simulating the first net productivity of the vegetation in the research area range by utilizing the land model based on the meteorological driving data and the underlying surface data. The method fully considers the influence of meteorological conditions and underlying surface data on the carbon sink effect of the vegetation, carries out 30m-100m resolution high-precision vegetation NPP simulation suitable for the planning field by combining the generated high-resolution meteorological driving data with the modifiable high-resolution land utilization type data driving a land model, quantitatively evaluates the carbon sink effect difference generated by different planning land scenes, further carries out an optimized space planning scheme under the carbon neutralization target, reasonably arranges construction indexes and implements a time sequence plan to provide decision reference, and has higher practicability.
Description
Technical Field
The invention relates to the technical field of numerical simulation, in particular to a high-resolution meteorological-driven planning carbon sink numerical evaluation method, a high-resolution meteorological-driven planning carbon sink numerical evaluation system and a high-resolution meteorological-driven planning carbon sink numerical evaluation terminal.
Background
"carbon peak-to-peak" and "carbon neutralization" are important subjects of socioeconomic development for a long period of time in the future, and the carbon sink effect of characterizing vegetation by Net Primary Productivity (NPP) is one of the mainstream methods for carrying out balance on carbon sinks in greenbelts such as forests, grasslands and the like at home and abroad at present.
At present, carbon sink capacity is estimated by a remote sensing process model method and a GIS space analysis module, wherein the remote sensing process model method comprises the steps of establishing a soil carbon flux ecological excitation remote sensing model to estimate regional carbon reserve through a light energy utilization model and a soil foundation respiration model, establishing a soil carbon flux remote sensing inversion model to estimate regional carbon reserve through establishing a spectral index, estimating vegetation first net productivity through a remote sensing image, and establishing a soil carbon flux ecological mechanism remote sensing model to estimate soil carbon sink capacity through the soil foundation respiration model; the GIS method obtains the change of the input and output of the ecosystem carbon in different land utilization modes by establishing the quantitative relation among the land utilization property, the carbon emission coefficient and the energy carbon sequestration capacity.
In the mainstream vegetation first net productivity estimation method, both methods have certain limitations, on one hand, the remote sensing model estimation method and the GIS space analysis module estimation method are highly dependent on remote sensing image data and are influenced by data precision, and the high-resolution estimation of the land scale for planning is difficult to achieve, on the other hand, the vegetation growth is related to the land utilization type, the meteorological condition, the soil type and other information, and the vegetation net productivity is affected by the geographical position of the area and the meteorological condition.
Disclosure of Invention
The invention aims to overcome the problems in the prior art and provides a high-resolution meteorological-driven planning carbon sink numerical evaluation method, a high-resolution meteorological-driven planning carbon sink numerical evaluation system and a high-resolution meteorological-driven planning carbon sink numerical evaluation terminal.
The purpose of the invention is realized by the following technical scheme: the method for evaluating the change of the land carbon sink for planning based on meteorological driving comprises the following steps:
generating high-resolution meteorological driving data with the spatial resolution of 1km-5 km;
and simulating the first net productivity of the vegetation in the research area range by utilizing the land model based on the meteorological driving data and the underlying surface data.
In one example, high resolution meteorological drive data of 1km to 5km is obtained from meteorological modeling using analytic or reanalyzed or forecasted field data.
In one example, the simulating a first net productivity of vegetation over the area of interest further comprises:
carrying out format conversion and interpolation processing on the meteorological drive data to generate an initial field and a boundary field;
and rewriting the driving data processing program of the land model, and compiling the driving data processing program into the driving data processing program of the meteorological model for outputting meteorological driving data, thereby generating an initial field file and a boundary field file which are required by the land model for simulating the vegetation first net productivity.
In one example, the method further comprises an underlying surface data updating step of:
mapping the land utilization type in the existing land utilization type data with the land utilization type in the meteorological model and the land model;
acquiring the data of the existing underlying surface by utilizing the latitude and longitude range of the grid, and calculating the area proportion of different land utilization types in each grid to further obtain a land utilization type matrix;
and updating static data files of the meteorological model and the land model according to the land utilization type matrix, and further completing the updating of the underlying surface data.
In one example, the deriving the land use type matrix includes:
and cutting the land utilization type data according to the longitude and latitude of the boundary of the meteorological model simulation grid points, calculating land utilization type ratios in the current grid points according to the quantity of certain specific land utilization type data, and further finishing different land utilization type ratios of all the grid points in the meteorological model simulation grid to obtain a land utilization type matrix.
In an example, the method further comprises a planning scenario underlying data updating step of:
and replacing the increment of the corresponding land utilization type in the simulation grid corresponding to the planning scene area by the updated underlying surface data, reducing other land utilization types, recalculating the occupation ratios of different land utilization types in each grid, further updating static data files of the meteorological model and the land model, and realizing the conversion from the planning scene to the model simulation scene.
In one example, the method further comprises a simulation result analysis processing step of:
extracting vegetation first net productivity information and vegetation type information according to a simulation result of the vegetation first net productivity, and storing the vegetation first net productivity information and the vegetation type information as a simplified result file;
generating projection information according to the longitude and latitude information of the simulation grid, the actual longitude and latitude information and the standard longitude STAND _ LON information, and generating a spatial distribution map by combining a simulation result of the first net productivity of the vegetation.
It should be further noted that the technical features corresponding to the above-mentioned method examples can be combined with each other or substituted to form a new technical solution.
The invention also comprises a meteorological-driven system for evaluating the change of carbon sink for planning land, which has the same technical concept as the meteorological-driven method for evaluating the change of carbon sink for planning land, which is formed by combining any one or more of the above examples, and comprises the following steps:
the meteorological model is used for generating high-resolution meteorological driving data with the spatial resolution of 1km-5 km;
and the land model is used for simulating the first net productivity of the vegetation in the research area range according to the meteorological driving data and the underlying surface data.
The system further comprises a meteorological element extraction unit for performing the steps of:
carrying out format conversion and interpolation processing on the meteorological drive data to generate an initial field and a boundary field;
and rewriting the driving data processing program of the land model, and compiling the driving data processing program into the driving data processing program of the meteorological model for outputting meteorological driving data, thereby generating an initial field file and a boundary field file which are required by the land model for simulating the vegetation first net productivity.
It should be further noted that the technical features corresponding to the above-mentioned system examples can be combined with each other or replaced to form a new technical solution.
The invention also includes a storage medium having stored thereon computer instructions operable to perform the steps of the method for estimating changes in geo-carbon sequestration for meteorological-driven planning based formed by any one or more of the above-described example compositions.
The invention also includes a terminal comprising a memory and a processor, wherein the memory stores computer instructions executable on the processor, and the processor executes the computer instructions to execute the steps of the method for estimating carbon sink change for planning land based on meteorological actuation formed by any one or more of the above examples.
Compared with the prior art, the invention has the beneficial effects that:
1. in one example, the method fully considers the influence of meteorological conditions and underlying surface data on the carbon sink effect of the vegetation, the generated high-resolution meteorological driving data drives a land model to develop 30m-100 m-resolution high-precision vegetation NPP simulation suitable for the planning field, the carbon sink effect difference generated by different planning land scenes is quantitatively evaluated, and then a carbon neutral under-target optimized space planning scheme, reasonable arrangement of construction indexes and implementation of a time sequence plan are developed to provide decision reference, so that the method has high practicability.
2. In one example, an initial field file and a boundary field file required by vegetation first net productivity simulation are generated by generating a land model, so that high-resolution meteorological driving data can drive the land model to carry out hectometer-level high-resolution NPP simulation, the technical blank of utilizing the high-resolution meteorological driving data to drive the land model to carry out carbon sink simulation research and application is filled, and a new thought is provided for high-precision NPP numerical simulation.
3. In one example, the influence of different planning scenes on meteorological conditions can be reflected through the underlying surface data updating step, the generation of meteorological driving data is reversely corrected, and the accuracy of the meteorological driving data is improved; and the land model is further driven by using the changed meteorological conditions, and the vegetation carbon sink change condition can be more accurately quantified and evaluated by combining the updated underlying surface data.
4. In one example, through the step of updating the data of the underlying surface of the planning scene, different underlying surface scenes can be designed by self-combining the current situation of the research area to carry out the first net productivity simulation of the multi-scene vegetation, and decision references are provided for optimizing a space planning scheme, reasonably arranging construction indexes and implementing a time sequence plan under the carbon neutralization target.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the invention without limiting the invention.
FIG. 1 is a flow chart of a method in an example of the invention;
FIG. 2 is a schematic diagram of high resolution weather-driven data generation in accordance with an example of the present invention;
FIG. 3 is a schematic diagram of a planned situational vegetation first net productivity assessment in accordance with an example of the present invention;
FIG. 4 is a spatial distribution diagram of NPP annual average values of an area of interest according to an example of the present invention;
FIG. 5 is a plot of a current site carbon sink map for a typical area at 50m resolution;
FIG. 6 is a plot of the carbon sink change evaluation for a 50m resolution plan for the same exemplary region of FIG. 5.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it should be noted that directions or positional relationships indicated by "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", and the like are directions or positional relationships described based on the drawings, and are only for convenience of description and simplification of description, and do not indicate or imply that the device or element referred to must have a specific orientation, be configured and operated in a specific orientation, and thus, should not be construed as limiting the present invention. Furthermore, ordinal words (e.g., "first and second," "first through fourth," etc.) are used to distinguish between objects, and are not limited to the order, but rather are to be construed to indicate or imply relative importance.
In the description of the present invention, it should be noted that, unless otherwise explicitly specified or limited, the terms "mounted," "connected" and "connected" are to be construed broadly and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in a specific case to those of ordinary skill in the art.
In addition, the technical features involved in the different embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
In an example, as shown in fig. 1, a method for estimating a change of a terrestrial carbon sink for planning based on meteorological driving specifically includes the following steps:
s1: generating spatial high-resolution meteorological driving data, wherein the high-resolution meteorological driving data represent high-resolution meteorological driving data with the resolution of 1km-5km, and comprise vegetation growth related elements such as air temperature, rainfall, soil temperature, soil humidity and the like;
s2: and simulating the first net productivity of the vegetation in the research area range by utilizing the land model based on the meteorological driving data and the underlying surface data.
Specifically, a land model HRLDAS (high resolution land data assimilation system) is widely used for land process simulation at present, and concerns soil-related information, soil problems, and the like. At present, the HRLDAS public driving data available in China is GLDAS, the spatial resolution of the GLDAS is 0.25 degrees and about 25Km, the resolution of related meteorological elements is low, and the requirement of high-resolution vegetation first net productivity simulation within 1 kilometer is difficult to meet. In order to solve the problem, high-resolution meteorological driving data in a research area is generated through a meteorological model such as a weather research and forecast model (WRF) to obtain a meteorological element simulation result with hourly resolution.
Furthermore, the underlying surface data comprise land utilization types, soil types, leaf area indexes and the like, when the land model simulates the first net productivity of the vegetation in the research area range, the influence of meteorological conditions and underlying surface data on the carbon sink effect of the vegetation is fully considered, the generated high-resolution meteorological driving data drives the land model to carry out high-precision vegetation NPP (first net productivity of the vegetation) simulation, the carbon sink effect difference generated by different planning land scenes is quantitatively evaluated, further, a space planning optimization scheme under the carbon neutralization target, a reasonable arrangement construction index and an implementation time sequence plan are developed to provide decision reference, and the method has high practicability.
As an option, the method can be used for developing vegetation first net productivity simulation based on historical meteorological conditions, can be used for developing carbon sink change analysis about half a year in the future based on meteorological driving data provided by a forecasting field, is flexible in simulation analysis, and can give consideration to both historical review simulation and forward-looking prediction simulation.
In one example, high resolution weather driven data of 1km to 5km is obtained via meteorological model simulation using either analytic site data (including but not limited to Final Operational Global Analysis, FNL, as published by the national environmental forecast center, NECP) or reanalysis site data (including but not limited to the fifth generation weather reanalysis data set, ERA5, for Global climate as published by the European mid-term weather forecast center, ECMWF), or forecaster site data (including but not limited to the Global forecast System, GFS, data set, as published by the national weather service). Specifically, a WRF model is used for carrying out high-resolution meteorological driving data simulation in a research area, different meteorological driving data are selected for carrying out simulation according to different applications, historical review simulation is taken as an example, 1-degree-of-resolution FNL analysis field data are selected for carrying out meteorological driving in the example, and dimension reduction is completed through step-by-step nesting, so that a meteorological driving data simulation result with the spatial resolution of 1.5km is obtained. Compared with the traditional WRF meteorological simulation, the method is different from the Lambert orthomorphic projection used by the traditional simulation, the projection type needs to be set as lat-lon (longitude and latitude projection) when the meteorological field simulation is carried out, the grid resolution is unit, and the grid resolution is not meter used by the traditional simulation; in addition, the weather analog output (wrfout) file time is set to hours, i.e., one wrfout file is output per hour.
In one example, specifically, since the HRLDAS model does not provide a processing interface for WRF data and does not support the use of the WRF model as driving data, the present invention develops a set of meteorological element extraction procedures, as shown in fig. 2, to perform the following steps before simulating the first net productivity of vegetation over a study area:
s00: carrying out format conversion and interpolation processing on the meteorological drive data to generate an initial field and a boundary field;
s01: and rewriting the driving data processing program of the land model, and compiling the driving data processing program into the driving data processing program of the meteorological model for outputting meteorological driving data, thereby generating an initial field file and a boundary field file which are required by the land model for simulating the vegetation first net productivity.
Wherein the initial field comprises 2m air temperature (T2), canopy water Content (CANWAT), surface air Temperature (TSK), soil Temperature (TSLB), soil humidity (SMOIS), snowfall (SNOW) and elevation (HGT) output by the WRF model; the boundary field comprises 2m air temperature (T2), 10m wind field (U10, V10), ground air Pressure (PSFC), 2m specific humidity (Q2), subsurface short wave radiation (SWDOWN), subsurface long wave radiation (LWDNB) and precipitation (RC, RN) which are output by the WRF model. The meteorological element extraction program extracts meteorological data output by the WRF model, reads the required meteorological elements hour by hour according to time, converts the meteorological data according to data units required by the HRLDAS model, creates folders according to variables respectively and stores the folders in a NetCDF format; and then, rewriting the HRLDAS driving data processing program create _ forcing.exe, compiling the result into wrf _ forcing.exe, enabling the result to be compatible with the initial field file and the boundary field file created by the extraction program, and generating the HRLDAS _ setup initial field file and the LDASIN boundary field file required by the HRLDAS simulation, thereby driving the HRLDAS model, reducing the model simulation result error caused by driving data, and meeting the requirement of the HRLDAS model for developing hundred-meter-level high-resolution simulation.
According to the invention, the high-resolution initial field file and the boundary field file which are required by the vegetation first net productivity simulation are generated through the land model, so that the high-resolution meteorological driving data can drive the land model to carry out the hundred-meter-level high-resolution NPP simulation, the technical blank of utilizing the high-resolution meteorological driving data to drive the land model to carry out the carbon sink simulation research and application is filled, and a new thought is provided for the high-precision NPP numerical simulation.
In one example, the HRLDAS model needs to generate static data using the geogrid module of the WRF model, and the underlying surface data is an important component of the WRF static data. Although the WRF model provides basic data, as for land utilization types, the urbanization speed of China is high, the change of the underlying surface is obvious, an underlying surface data set USGS of the WRF model is data in 1992, an MODIS data set is data in 2002, the timeliness is poor, in order to better reflect the recent urbanization development process, the influence of the underlying surface data on vegetation growth is fully considered, the NPP simulation precision is improved, and the benchmark high-resolution underlying surface data needs to be updated, and the method comprises the following substeps:
s11: mapping the land utilization type in the existing land utilization type data with the land utilization type in the meteorological model and the land model;
s12: acquiring the data of the existing underlying surface by utilizing the latitude and longitude range of the grid, and calculating the proportion of different land utilization types in each grid to further obtain a land utilization type matrix;
s13: and updating static data files of the meteorological model and the land model according to the land utilization type matrix, and further completing the updating of the underlying surface data.
Specifically, the existing land use type data may be acquired through official published information or the like disclosed by the public network. Taking a global 30m land cover data set FROM-GLC data disclosed by Qinghua university as an example, establishing a relational mapping table of land utilization types in the data set and land utilization types of a meteorological model WRF, and mapping in the following way:
table 1 land use type mapping table
For planning the situational land utilization, the land utilization mapping is established by referring to the upper table, and after the classification mapping is finished, the land utilization type matrix C is used as an option to obtain a land utilization type matrix C Xi,Yi,k The method specifically comprises the following steps:
cutting the land use type data according to the longitude and latitude of four boundaries of the simulation grid lattice point of the WRF model, and according to a certain specific land use type data Lu i The land utilization type ratio in the grid points of the current grid is calculated, and then different land utilization type ratios of all the grid points in the meteorological model simulation grid are completed, and a land utilization type matrix is obtained. More specifically, the proportion calculation formula of the kth land utilization type of the WRF model simulation grid with coordinates (Xi, yi) is as follows:
wherein n represents the number of land use types; lu denotes all land use types.
Further, updating the static data files of the meteorological model and the land model specifically includes:
formatted land use type matrix C Xi,Yi,k The method can be used for updating LANDSUSEF (land utilization type) variables in the WRF model geo _ em static data file, calculating the type with the largest proportion in different land utilization types as a dominant land utilization type, namely, the method can be used for updating LU _ INDEX (dominant land utilization type INDEX) variables, then selecting grids with the dominant land utilization type in the LU _ INDEX variables as a water body, setting the grids corresponding to the land mask LANMASK variables to be 0, and finishing the updating related to the land utilization type data.
The underlying data update in this example will be applied to both the geo _ em static data file required for the WRF model and the geo file required for the HRLDAS, which typically has a higher spatial resolution, which may reach hundreds of meters or even higher, relative to the resolution of the WRF model up to 3 km.
In the example, through the step of updating the underlying surface data, the influence of the underlying surface data on meteorological conditions in different planning situations can be reflected, the generation of meteorological driving data is reversely corrected, and the accuracy of the meteorological driving data is improved, namely, the updated land utilization type data has a certain improvement effect on the meteorological element numerical simulation; and the land model is further driven by using the changed meteorological conditions, and the vegetation carbon sink change condition can be more accurately quantified and evaluated by combining the updated underlying surface data.
In an example, the planning scenario is a land utilization planning scenario in a specific area, different planning scenarios may affect meteorological conditions, and in order to ensure simulation accuracy, the method further includes updating the underlying surface data of the planning scenario in the geo file of the HRLDAS model. Specifically, the planning scenario is stored in a GeoTIFF file form, and the new land utilization type is recorded in a grid form; reading a scenario land utilization type file, calculating in the same way as the updating step of the datum high-resolution underlying surface data, replacing the increment of the corresponding land utilization type in the simulation grid, reducing other land utilization types, recalculating the proportion of different land utilization types, calculating a dominant land utilization type and a land and water mask, and writing the result into a geo file of the HRLDAS model to realize the conversion from a planning scenario to a model simulation scenario.
As an option, when fine research is performed, the method for updating underlying surface data of the planning scene is also applied to a geo _ em file required by a WRF model, and meteorological field condition changes caused by the planning scene are simulated again, such as changes of elements such as quantized soil temperature and humidity and precipitation, so that the accuracy of evaluation work is improved.
In this example, based on the model basic data updating program, the multi-scenario underlying surface types are established according to the space planning, so that the replacement of the designated area, the designated range and the designated underlying surface types is realized, and the carbon sink influence of different underlying surface types on different space positions is quantitatively evaluated. Meanwhile, through the step of updating the data of the underlying surface under the planning situation, different underlying surface situations can be designed by self-combining the current situation of the research area to carry out the first net productivity simulation of the multi-situation vegetation, and decision reference is provided for optimizing a space planning scheme, reasonably arranging construction indexes and implementing a time sequence plan under the carbon neutralization target.
Combining the above steps results in a preferred example of developing NPP estimates for a multi-project scenario, as shown in fig. 3, when the high resolution vegetation first net productivity simulation comprises the steps of:
s1': preparing a high-resolution land utilization type covering the simulation area and the periphery, and finishing data updating of the WRF reference high-resolution underlying surface;
s2': carrying out high-resolution meteorological field simulation by using a WRF model to obtain a high-resolution meteorological field covering a research area by 1km-5 km;
and S3': reading high-resolution meteorological field data, completing HRLDAS driving data conversion, and generating HRLDAS high-resolution meteorological driving data;
s4': carrying out NPP simulation by using an HRLDAS model;
and S5': and finishing data extraction and processing to generate the GeoTIFF format NPP raster data.
Further, the NPP data difference obtained by simulation of different planning scenarios is the first net productivity difference of vegetation caused by the planning scenario, and the actual simulation process needs to be initialized and simulated for at least one month to obtain soil temperature and humidity data with higher precision.
To sum up, the invention couples the meteorological model with the land model, fully reflects the vegetation first net productivity change of the planned land vegetation type change under the influence of meteorological conditions based on the numerical simulation technology, utilizes the vegetation growth related factors such as air temperature, precipitation, soil temperature, soil humidity and the like provided by the high-resolution meteorological data, combines the vegetation data and the soil data in the high-resolution land utilization type data, simulates the plant first net productivity, and quantitatively evaluates the plant first net productivity value under different national and soil space planning situations by adjusting the vegetation coverage condition in the land utilization type data.
In one example, in order to facilitate the analysis and utilization of the simulation result in the planning field, the method further comprises the steps of analyzing and processing the simulation result, namely extracting vegetation first net productivity information and vegetation type information according to the simulation result of the vegetation first net productivity, and storing the vegetation first net productivity information and the vegetation type information as a simplified result file; and then generating projection information according to the longitude and latitude information of the simulation grid, the actual longitude and latitude information and the standard longitude STAND _ LON information, and generating a spatial distribution map by combining a simulation result of the first net productivity of the vegetation. Taking the NPP simulation result of the metropolis as an example, the specific implementation process for generating the spatial distribution map is as follows:
s61: reading an NPP simulation result file of HRLDAS, extracting NPP and IVGTYP variables in a simulation result by utilizing ncks program of the nco program package, respectively corresponding to the first net productivity of vegetation and the type of vegetation, storing the NPP and IVGTYP variables as an HRLDAS simplification result file, reducing the size of the file and improving the subsequent processing speed;
s62: reading an HRLDAS simplification result file, and generating projection information by using file information such as LAT1 (corresponding to the central latitude of the simulation grid), LON1 (corresponding to the central longitude of the simulation grid), TRUELAT1 (for Ying Lan burt projection standard latitude 1), TRULAT2 (for Ying Lan burt projection standard latitude 2), STAND _ LON (consistent with the central longitude of the grid) and the like;
s63: the NPP simulation results are read, statistics is carried out according to time periods, and a GeoTIFF file is generated by combining projection information and is used for analyzing the change condition of the subsequent carbon sink as shown in FIG. 4.
To further illustrate the technical effects of the scheme of the invention, on the basis of the current land carbon sink condition of 50m resolution of a certain typical area as shown in fig. 5, the carbon sink simulation is performed on the 50m resolution planning land of the typical area by adopting the preferred example method of the invention to obtain the carbon sink change evaluation diagram of the planning land as shown in fig. 6, and it can be seen that the high-resolution (50 m in the example) high-precision vegetation NPP simulation can be developed by driving the land model through the generated high-resolution meteorological driving data, and the method is suitable for carbon sink change evaluation under the planning situation, provides decision reference for developing a space planning scheme, reasonably arranging construction indexes and implementing a time sequence plan under the carbon neutralization target, and has high practicability.
In an example, the present invention further provides a storage medium, which has the same inventive concept as the meteorological drive-based planned land carbon sink change assessment method formed by any one or more of the above examples in combination, and on which computer instructions are stored, and the computer instructions are executed when executed to perform the steps of the meteorological drive-based planned land carbon sink change assessment method formed by any one or more of the above examples in combination.
Based on such understanding, the technical solution of the present embodiment or parts of the technical solution may be essentially implemented in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, an optical disk, or other various media capable of storing program codes.
In an example, the present invention further provides a terminal, wherein any example or any combination of examples corresponding to the above meteorological-driven planned land carbon sink change assessment method has the same inventive concept, and the terminal comprises a memory and a processor, wherein the memory stores computer instructions executable on the processor, and the processor executes the steps of the above meteorological-driven planned land carbon sink change assessment method when executing the computer instructions. The processor may be a single or multi-core central processing unit or a specific integrated circuit, or one or more integrated circuits configured to implement the present invention.
In one example, a terminal, i.e., an electronic device, is represented in the form of a general purpose computing device, and components of the electronic device may include, but are not limited to: the at least one processing unit (processor), the at least one memory unit, and a bus connecting the various system components including the memory unit and the processing unit.
Wherein the storage unit stores program code executable by the processing unit to cause the processing unit to perform steps according to various exemplary embodiments of the present invention as described in the above section "exemplary methods" of the present specification. For example, the processing unit may execute the above-described method for estimating changes in carbon sink for planning based on meteorological actuation.
The memory unit may include a readable medium in the form of a volatile memory unit, such as a random access memory unit (RAM) 3201 and/or a cache memory unit, and may further include a read only memory unit (ROM).
The storage unit may also include a program/utility having a set (at least one) of program modules including, but not limited to: an operating system, one or more application programs, other program modules, and program data, each of which, or some combination thereof, may comprise an implementation of a network environment.
The bus may be any of several types of bus structures including a memory unit bus or memory unit controller, a peripheral bus, an accelerated graphics port, a processing unit, or a local bus using any of a variety of bus architectures.
The electronic device may also communicate with one or more external devices (e.g., keyboard, pointing device, bluetooth device, etc.), with one or more devices that enable a user to interact with the electronic device, and/or with any devices (e.g., router, modem, etc.) that enable the electronic device to communicate with one or more other computing devices. Such communication may be through an input/output (I/O) interface. Also, the electronic device may communicate with one or more networks (e.g., a Local Area Network (LAN), a Wide Area Network (WAN), and/or a public network, such as the internet) via a network adapter. The network adapter communicates with other modules of the electronic device over the bus. It should be appreciated that other hardware and/or software modules may be used in conjunction with the electronic device, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems, to name a few.
Through the above description, those skilled in the art will readily understand that the exemplary embodiments described herein may be implemented by software, or by software in combination with necessary hardware. Therefore, the technical solution according to the exemplary embodiment may be embodied in the form of a software product, which may be stored in a non-volatile storage medium (which may be a CD-ROM, a usb disk, a removable hard disk, etc.) or on a network, and includes several instructions to enable a computing device (which may be a personal computer, a server, a terminal device, or a network device, etc.) to execute the method of the exemplary embodiment of the present application.
The invention also comprises a carbon sink change evaluation system for planning land based on meteorological driving, which has the same inventive concept as the carbon sink change evaluation method for planning land based on meteorological driving formed by any one or more examples, and the system comprises:
the meteorological model is used for generating high-resolution meteorological driving data with the spatial resolution of 1km-5 km; wherein the meteorological model is preferably a WRF model, for driving the land model based on the generated meteorological drive data.
A land model, preferably an HRLDAS model, for simulating a first net productivity of vegetation over the area of interest from the meteorological drive data, underlay surface data.
In an example, the system further comprises a meteorological element extraction unit for performing the steps of:
carrying out format conversion and interpolation processing on the meteorological driving data to generate an initial field and a boundary field;
and rewriting the driving data processing program of the land model, and compiling the driving data processing program into the driving data processing program of the meteorological model for outputting meteorological driving data, thereby generating an initial field file and a boundary field file which are required by the land model for simulating the vegetation first net productivity.
In one example, the system further comprises an underlying data updating unit for performing the steps of:
mapping the land utilization type in the existing land utilization type data with the land utilization type in the meteorological model and the land model;
acquiring the data of the existing underlying surface by utilizing the latitude and longitude range of the grid, and calculating the proportion of different land utilization types in each grid to further obtain a land utilization type matrix;
and updating static data files of the meteorological model and the land model according to the land utilization type matrix, and further completing updating of the underlying surface data.
In an example, the system further comprises a planning scenario underlying data updating unit for performing the steps of: and replacing the increment of the corresponding land utilization type in the simulation grid corresponding to the planning scene area by the updated underlying surface data, reducing other land utilization types, recalculating the occupation ratios of different land utilization types in each grid, further updating static data files of the meteorological model and the land model, and realizing the conversion from the planning scene to the model simulation scene.
In one example, the system further comprises a simulation result analysis processing unit for performing the steps of:
extracting vegetation first net productivity information and vegetation type information according to a simulation result of the vegetation first net productivity, and storing the vegetation first net productivity information and the vegetation type information as a simplified result file;
and generating projection information according to the longitude and latitude information of the simulation grid, the actual longitude and latitude information and the STAND _ LON information, and generating a spatial distribution map by combining a simulation result of the first net productivity of the vegetation.
The above detailed description is for the purpose of describing the invention in detail, and it should not be construed that the detailed description is limited to the description, and it will be apparent to those skilled in the art that various modifications and substitutions can be made without departing from the spirit of the invention.
Claims (10)
1. A planning land carbon sink change evaluation method based on meteorological driving is characterized by comprising the following steps: which comprises the following steps:
generating high-resolution meteorological driving data with the spatial resolution of 1km-5 km;
and simulating the first net productivity of the vegetation in the research area range by utilizing the land model based on the meteorological driving data and the underlying surface data.
2. The method for estimating carbon sink change for planning based on meteorological drive according to claim 1, wherein: and (3) obtaining high-resolution meteorological driving data of 1km-5km by adopting analytical field data or reanalysis field data or forecast field data through meteorological model simulation.
3. The method for estimating carbon sink change for planning based on meteorological drive according to claim 1, wherein: the method further comprises, prior to simulating a first net productivity of vegetation over the area of interest:
carrying out format conversion and interpolation processing on the meteorological drive data to generate an initial field and a boundary field;
and rewriting the driving data processing program of the land model, and compiling the driving data processing program into the driving data processing program of the meteorological model for outputting meteorological driving data, thereby generating an initial field file and a boundary field file which are required by the land model for simulating the vegetation first net productivity.
4. The method for estimating carbon sink change for planning based on meteorological drive according to claim 1, wherein: the method also comprises the following step of updating the underlying surface data:
mapping the land utilization type in the existing land utilization type data with the land utilization type in the meteorological model and the land model;
acquiring the data of the existing underlying surface by utilizing the latitude and longitude range of the grid, and calculating the proportion of different land utilization types in each grid to further obtain a land utilization type matrix;
and updating static data files of the meteorological model and the land model according to the land utilization type matrix, and further completing the updating of the underlying surface data.
5. The method for estimating carbon sink change for planning based on meteorological drive according to claim 4, wherein: the obtaining of the land use type matrix comprises:
and cutting the land utilization type data according to the longitude and latitude of the boundary of the meteorological model simulation grid point, calculating the land utilization type ratio in the current grid point according to the area of certain specific land utilization type data, and further finishing different land utilization type ratios of all grid points in the meteorological model simulation grid to obtain a land utilization type matrix.
6. The method for estimating carbon sink change for planning based on meteorological drive according to claim 4, wherein: the method also comprises a step of updating the underlying surface data of the planning scene:
and replacing the increment of the corresponding land utilization type in the simulation grid corresponding to the planning scene area by the updated underlying surface data, reducing other land utilization types, recalculating the occupation ratios of different land utilization types in each grid, further updating static data files of the meteorological model and the land model, and realizing the conversion from the planning scene to the model simulation scene.
7. The method for estimating carbon sink change for planning based on meteorological drive according to claim 1, wherein: the method also comprises a simulation result analysis processing step:
extracting vegetation first net productivity information and vegetation type information according to a simulation result of the vegetation first net productivity, and storing the vegetation first net productivity information and the vegetation type information as a simplified result file;
and generating projection information according to the longitude and latitude information of the simulation grid, the actual longitude and latitude information and the standard longitude information, and generating a spatial distribution map by combining a simulation result of the first net productivity of the vegetation.
8. The utility model provides a land carbon sink change evaluation system for planning based on meteorological drive which characterized in that: it comprises the following steps:
the meteorological model is used for generating high-resolution meteorological driving data with the spatial resolution of 1km-5 km;
and the land model is used for simulating the first net productivity of the vegetation in the research area range according to the meteorological driving data and the underlying surface data.
9. The system according to claim 8, wherein the system comprises: the system further comprises a meteorological element extraction unit for performing the steps of:
carrying out format conversion and interpolation processing on the meteorological drive data to generate an initial field and a boundary field;
and rewriting the driving data processing program of the land model, and compiling the driving data processing program into the driving data processing program of the meteorological model for outputting meteorological driving data, thereby generating an initial field file and a boundary field file which are required by the land model for simulating the vegetation first net productivity.
10. A terminal comprising a memory and a processor, the memory having stored thereon computer instructions executable on the processor, characterized in that: the processor when executing the computer instructions performs the steps of the method for estimating changes in terrestrial carbon interchange for planning based on meteorological driving as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211584270.2A CN115758801A (en) | 2022-12-09 | 2022-12-09 | High-resolution weather-driven planning carbon sink numerical evaluation method, system and terminal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211584270.2A CN115758801A (en) | 2022-12-09 | 2022-12-09 | High-resolution weather-driven planning carbon sink numerical evaluation method, system and terminal |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115758801A true CN115758801A (en) | 2023-03-07 |
Family
ID=85345222
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211584270.2A Pending CN115758801A (en) | 2022-12-09 | 2022-12-09 | High-resolution weather-driven planning carbon sink numerical evaluation method, system and terminal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115758801A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116842351A (en) * | 2023-09-01 | 2023-10-03 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | Coastal wetland carbon sink assessment model construction method, assessment method and electronic equipment |
| CN117574329A (en) * | 2024-01-15 | 2024-02-20 | 南京信息工程大学 | Nitrogen dioxide refined space distribution method based on ensemble learning |
| CN117576339A (en) * | 2023-09-20 | 2024-02-20 | 中国气象局上海台风研究所(上海市气象科学研究所) | Unstructured grid generation method with highest resolution of 1km |
| CN117745977A (en) * | 2024-02-20 | 2024-03-22 | 四川开澜科技有限公司 | Mountain natural wind field analysis method, system and terminal based on 3DGIS |
| CN118227989A (en) * | 2024-03-19 | 2024-06-21 | 广东省气象数据中心 | Intelligent analysis and comparison method and system for meteorological data under different pad surfaces |
-
2022
- 2022-12-09 CN CN202211584270.2A patent/CN115758801A/en active Pending
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116842351A (en) * | 2023-09-01 | 2023-10-03 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | Coastal wetland carbon sink assessment model construction method, assessment method and electronic equipment |
| CN116842351B (en) * | 2023-09-01 | 2023-11-10 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | Coastal wetland carbon sink assessment model construction method, assessment method and electronic equipment |
| CN117576339A (en) * | 2023-09-20 | 2024-02-20 | 中国气象局上海台风研究所(上海市气象科学研究所) | Unstructured grid generation method with highest resolution of 1km |
| CN117576339B (en) * | 2023-09-20 | 2024-05-24 | 中国气象局上海台风研究所(上海市气象科学研究所) | Unstructured grid generation method with highest resolution of 1km |
| CN117574329A (en) * | 2024-01-15 | 2024-02-20 | 南京信息工程大学 | Nitrogen dioxide refined space distribution method based on ensemble learning |
| CN117574329B (en) * | 2024-01-15 | 2024-04-30 | 南京信息工程大学 | Nitrogen dioxide refined space distribution method based on ensemble learning |
| CN117745977A (en) * | 2024-02-20 | 2024-03-22 | 四川开澜科技有限公司 | Mountain natural wind field analysis method, system and terminal based on 3DGIS |
| CN117745977B (en) * | 2024-02-20 | 2024-05-07 | 四川开澜科技有限公司 | Mountain natural wind field analysis method, system and terminal based on 3DGIS |
| CN118227989A (en) * | 2024-03-19 | 2024-06-21 | 广东省气象数据中心 | Intelligent analysis and comparison method and system for meteorological data under different pad surfaces |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zamani et al. | A comparison of CMIP6 and CMIP5 projections for precipitation to observational data: the case of Northeastern Iran | |
| CN115758801A (en) | High-resolution weather-driven planning carbon sink numerical evaluation method, system and terminal | |
| Lindberg et al. | Urban Multi-scale Environmental Predictor (UMEP): An integrated tool for city-based climate services | |
| Xu et al. | Urban morphology detection and computation for urban climate research | |
| Wan | MODIS land surface temperature products users’ guide | |
| Zawadzka et al. | Assessment of heat mitigation capacity of urban greenspaces with the use of InVEST urban cooling model, verified with day-time land surface temperature data | |
| Klouček et al. | How does data accuracy influence the reliability of digital viewshed models? A case study with wind turbines | |
| CN104408258A (en) | Large-scale soil organic carbon spatial distribution simulation method involving environmental factors | |
| KR20140021179A (en) | Method and apparatus for predicting daily solar radiation level | |
| Rocha et al. | Mapping evapotranspirative and radiative cooling services in an urban environment | |
| CN114357879A (en) | Future land utilization simulation method and device, storage medium and terminal equipment | |
| Berckmans et al. | Reinitialised versus continuous regional climate simulations using ALARO-0 coupled to the land surface model SURFEXv5 | |
| Mortensen | Wind resource assessment using the WAsP software (DTU Wind Energy E-0135) | |
| Ma et al. | Individual structure mapping over six million trees for New York City USA | |
| CN117421871A (en) | Offshore wind power potential evaluation method and device and computer equipment | |
| Barde et al. | Performance of the CORDEX-SA regional climate models in simulating summer monsoon rainfall and future projections over East India | |
| Shiri et al. | Estimation of CO 2 flux components over northern hemisphere forest ecosystems by using random forest method through temporal and spatial data scanning procedures | |
| Mortensen | 46200 Planning and Development of Wind Farms: Wind resource assessment using the WAsP software | |
| Sperati et al. | A new Wind Atlas to support the expansion of the Italian wind power fleet | |
| Brandi et al. | Influence of projected climate change, urban development and heat adaptation strategies on end of twenty-first century urban boundary layers across the Conterminous US | |
| Jabal et al. | Integrated Approach for Land Surface Temperature Assessment in Different Topography of Iraq | |
| CN113935622A (en) | Ecological geological environment quality evaluation method and system and electronic equipment | |
| Coelho et al. | Current and future geographical distribution of the indoor conditions for high thermal inertia historic buildings across Portugal via hygrothermal simulation | |
| Reuter et al. | Modelling climate analogue regions for a central European city | |
| Wijesekara et al. | Impact of land-use changes on the hydrological processes in the Elbow river watershed in southern Alberta |
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 |