Method for analyzing nitrogen and phosphorus pollutants in land-sea system

Abstract

The application discloses a land-sea system nitrogen and phosphorus pollutant collaborative source analysis method, and particularly relates to the fields of environmental treatment and water pollution tracing, and comprises a comprehensive method based on a multi-isotope tracing technology, an isotope mixing model and a verification of a nitrogen and phosphorus molar ratio, obtains isotope characteristics of various forms of nitrogen sources and phosphorus sources in a land-sea system, combines hydrological and water quality monitoring data, quantitatively analyzes the multi-form contributions of river water, groundwater, coastal land sources and estuary internal sources to estuary nitrogen and phosphorus pollutants, and solves the problems of system deficiency and accuracy deficiency of the prior art in complex multi-source pollutant tracing, and provides a basis for integrated treatment of a basin-estuary-bay.

Description

Technical Field

[0001] This invention relates to the fields of environmental governance and water pollution source tracing, and more specifically, to a method for synergistic source apportionment of nitrogen and phosphorus pollutants in land-sea systems. Background Technology

[0002] Nitrogen and phosphorus pollution is a significant challenge for the ecological and environmental governance of the land-sea system. Its source apportionment involves land-based sources (rivers, groundwater), marine sources (estuarine sources), and complex migration and transformation processes. Existing research mainly focuses on the single source or single form of nitrogen and phosphorus pollutants, usually using pollution discharge statistics or model simulation methods. It lacks a holistic consideration of the land-sea coordinated system. In addition, most technologies can only estimate the flux distribution of pollutants and are difficult to accurately identify the contributions of groundwater input, estuarine release, and other hidden sources of pollutants. They are particularly limited when dealing with complex multi-source pollution.

[0003] Currently, there are no patents for river-estuary nitrogen and phosphorus pollutant source tracing technology based on integrated land and sea management. Among existing technologies, such as the watershed nitrogen and phosphorus pollutant analysis method and system based on multi-isotope joint tracing disclosed in CN114384224A, it can only calculate the pollutant flux transported by rivers flowing into the sea to the estuary, but cannot accurately identify the dynamic impact of groundwater input and coastal land sources. In particular, it lacks analysis methods for the contribution of pollution released within the estuary. For phosphorus pollutants, traditional methods usually calculate their source based on a single phosphate oxygen isotope, but due to the complexity of the morphology and incomplete data, the analysis results have great uncertainty. These problems limit the application of existing technologies in large-scale complex environmental governance. Summary of the Invention

[0004] To overcome the aforementioned deficiencies of the prior art, embodiments of the present invention provide a method for synergistic source apportionment of nitrogen and phosphorus pollutants in land-sea systems. This method utilizes a multi-isotope tracing technology system and an isotope mixing model to address the problems mentioned in the background section.

[0005] To achieve the above objectives, the present invention provides the following technical solution: a method for co-source apportionment of nitrogen and phosphorus pollutants in land-sea systems, comprising:

[0006] S1. Obtain the isotopic characteristics of nitrogen and phosphorus sources of various forms in the land-sea system. Utilize elemental and stable isotope tracing technology and isotope mixing models to quantitatively analyze the contributions of river water, groundwater, coastal terrestrial sources, and estuarine endogenous sources to different forms of nitrogen and reactive phosphate in estuarine seawater. Different forms of nitrogen include ammonia nitrogen, nitrate nitrogen, nitrite nitrogen, dissolved organic nitrogen, and particulate organic nitrogen.

[0007] S2. Based on hydrological monitoring data and water quality monitoring data, quantitatively calculate the flux of nitrogen and phosphorus pollutants of various forms in rivers into the sea and their spatiotemporal variation patterns; based on the flux of nitrogen and phosphorus pollutants in rivers and the results of source apportionment of estuary seawater, determine the flux of nitrogen and phosphorus pollutants input into groundwater.

[0008] S3. Quantitatively analyze the contribution ratio of different forms of nitrogen and phosphorus pollutants in rivers and groundwater;

[0009] S4. Based on the source apportionment results of nitrogen and phosphorus pollutants in various forms from rivers and groundwater, calculate the input flux of nitrogen and phosphorus pollutants in different forms from each source into the bay seawater through river water and groundwater.

[0010] S5. Based on the above results, preliminary source apportionment results for various forms of nitrogen and phosphorus pollutants in the estuary are formed;

[0011] S6. Based on the source apportionment results of different forms of nitrogen pollutants and the percentage composition of each form of nitrogen from each pollution source, verify the accuracy of the nitrogen pollutant source apportionment results.

[0012] S7. Based on the source apportionment results of reactive phosphate and the nitrogen-phosphorus ratio of each pollution source, verify the accuracy of the phosphate source apportionment results;

[0013] S8. Based on the source apportionment results of phosphorus content percentage and reactive sulfate in various forms from different pollution sources, determine the input of each pollution source to dissolved organic phosphorus and particulate organic phosphorus.

[0014] S9. Based on the concept of integrated land and sea management, complete the source analysis of nitrogen and phosphorus pollutants of various forms in the estuary.

[0015] In a preferred embodiment, according to step S1, the isotopic characteristics of the nitrogen and phosphorus sources are obtained and their contributions are analyzed; a plan is formulated. The isotopic contribution ratio of nitrogen or phosphorus form i in pollution source j, corresponding to region k; weighted by each sample point. Isotope abundance ratio Calculate the contribution of source j;

[0016]

[0017] in Let i be the isotopic abundance weight of morphology i at the m-th sample point, taking into account the characteristics of region k; The characteristic ratio of isotopes; The total flow contribution of the m-th sample point is calculated by combining the monitoring data of sample points in region k; N, M, and P represent the number of pollution sources, the number of sample points, and the number of pollution types, respectively.

[0018] In a preferred embodiment, the inflow flux and its spatiotemporal variation are quantitatively calculated according to step S2; T is then determined. i (t,z) represents the total sea flux of nitrogen or phosphorus form i at time t and depth z; Φ is proposed. i (x,t,z) represents the three-dimensional distribution of pollutant concentration. Combined with the dynamic change of flow rate Q(x,t), the contribution of each location to the total flux is calculated.

[0019]

[0020] Where Φ i (x,t,z) is the concentration distribution function of morphology i at location x, time t, and depth z; Q(x,t) is the flow distribution of the river at location x and time t; K z (x) is the depth correction coefficient, which represents the weight of the influence of different depth layers at location x on the flow; L is the total length of the river.

[0021] In a preferred embodiment, a groundwater input calibration model is established according to step S3; a formulation is made. The input flux of nitrogen or phosphorus form i from groundwater below depth z;

[0022]

[0023] T i (z) represents the total inflow flux into the ocean at depth z for morphology i; C i,j The proportion of pollution source j to form i; L j The total load of pollution source j; denoted as the total pollutant load of morphology i below depth layer z; J represents the total number of pollution sources.

[0024] In a preferred embodiment, the contribution ratios of different sources are quantitatively analyzed according to step S4; F is proposed. i,j The quantitative contribution ratio of pollution source j to nitrogen or phosphorus form i;

[0025]

[0026] Where Θ j Ψ is a regional characteristic factor of pollution source j; j Θ represents the environmental impact factor of pollution source j; K represents the total number of pollution sources; k Ψ is a regional characteristic factor of pollution source k; k Environmental impact factors of pollution source k;

[0027] Based on step S5, preliminary source analysis results are generated; R is proposed. i,j The preliminary analysis results of the nitrogen or phosphorus form i for pollution source j;

[0028] R i,j =F i,j ·T i

[0029] Where F i,j T represents the quantitative contribution of pollution source j to nitrogen or phosphorus form i; i Let i be the total flux into the sea for form i.

[0030] In a preferred embodiment, the accuracy of nitrogen pollutant source apportionment is verified according to step S6; V is formulated. i,j The error ratio of nitrogen pollutant speciation i to source j;

[0031]

[0032] Where R i,j These are preliminary analysis results; These are the results of actual measurements and analysis.

[0033] In a preferred embodiment, the accuracy of phosphorus pollutant source apportionment is verified according to step S7; Q is then determined. i,j The accuracy factor for phosphorus pollutant speciation i and source j is determined by verifying and analyzing the results through the nitrogen-phosphorus ratio and measured phosphorus load.

[0034]

[0035] in The measured load of phosphorus pollutant form i in source j; N i,j is the nitrogen-to-phosphorus ratio verification value for source j relative to form i; k is the index of the pollution source; K is the total number of pollution sources; By measuring phosphorus load from all sources and its corresponding nitrogen-to-phosphorus ratio N i,k We obtain the overall verification benchmark by weighted summation.

[0036] In a preferred embodiment, the sources of other forms of phosphorus are determined according to step S8; L is formulated. i,j Input amounts of dissolved organophosphates and particulate organophosphates for source j to form i; calibration of other phosphorus forms based on reactive phosphate results;

[0037] L i,j =Q i,j ·α i,j

[0038] Q i,j α is an accuracy factor for determining the source j of phosphorus pollutant speciation i; i,j The conversion factor is used to represent the rate at which reactive phosphate is converted to other forms of phosphorus.

[0039] Complete the final source parsing according to step 9; formulate S i Precise source apportionment results for nitrogen or phosphorus pollutants in final form i

[0040]

[0041] Where R i,j This is the preliminary source resolution result; L i,j J represents the input amount of dissolved organophosphates and particulate organophosphates from source j to form i; J represents the total number of pollution sources.

[0042] The technical effects and advantages of this invention are as follows:

[0043] 1. Achieving precise source tracing of multi-form nitrogen and phosphorus pollutants in land-sea systems: This invention establishes a comprehensive and precise source apportionment method by constructing a multi-isotope tracer system and isotope mixing model, targeting the multi-form characteristics of nitrogen and phosphorus pollutants such as ammonia nitrogen, nitrate nitrogen, and dissolved organic phosphorus, and taking into account land-based, marine-based, and internal release factors. Compared with traditional technologies, this invention can not only quantitatively analyze the contribution of rivers, groundwater, and estuaries to the multi-form nitrogen and phosphorus pollutants, but also verify the analysis results by combining the nitrogen and phosphorus molar ratio and composition of different pollution sources, thereby improving the accuracy of analysis and providing a scientific basis for the pollution control of complex watersheds and estuaries.

[0044] 2. Dynamic quantification of spatiotemporal variations in river and groundwater pollution fluxes: This invention combines hydrological and water quality monitoring data to construct a dynamic quantification model of pollutant fluxes, which can accurately reflect the spatiotemporal evolution characteristics of nitrogen and phosphorus pollutant fluxes in rivers and groundwater. By quantitatively calculating the input of different forms of pollutants in different periods and regions, this technology provides multi-scale data support for dynamic pollution control. Compared with existing static methods, the dynamic characteristics of this invention significantly improve the adaptability and targeting of pollutant control strategies.

[0045] 3. Enhance groundwater input calibration and latent pollution identification: Addressing the shortcomings of traditional technologies in identifying latent groundwater inputs, this invention, based on a residual quantity estimation algorithm and a dynamic calibration model, can efficiently identify the dynamic input flux of groundwater and estuarine pollutants and quantify their contribution to the total pollution load; it improves the sensitivity of analytical models to latent pollution sources and provides key data support for comprehensive pollution control.

[0046] 4. Collaborative analysis of nitrogen and phosphorus pollutants and establishment of cross-validation mechanism: This invention constructs a collaborative analysis system for nitrogen and phosphorus pollutants, incorporating the nitrogen-phosphorus molar ratio and content of various forms of pollutants into the verification process, thereby achieving multi-level cross-validation of the analysis results. Compared with traditional methods, this invention has a stronger systematic approach to verification, ensuring the scientific validity and credibility of source apportionment results, and is particularly suitable for research and application in complex polluted environments.

[0047] 5. Dynamic calibration of multiple sources of phosphorus pollutants: In view of the complexity of phosphorus pollutant source tracing, this invention establishes a dynamic calibration model for the conversion relationship between reactive phosphate and other phosphorus forms, such as dissolved organic phosphorus and particulate phosphorus. Combined with source apportionment and verification mechanisms, this method improves the comprehensiveness and accuracy of phosphorus pollutant source tracing, and provides strong technical support for refined governance.

[0048] 6. Provides a scientific basis for integrated land-sea pollution control; by integrating the input, migration, and transformation of multi-source pollution in the land-sea system, this invention proposes a source apportionment technology framework based on integrated land-sea management; this method can effectively identify the key contribution ratios of pollution from rivers, groundwater, estuaries, and marine sources, and dynamically reflect their role in governance, providing a reliable basis for the formulation of regional and large-scale integrated watershed-estuary-bay governance strategies; compared with traditional single-form or single-source apportionment technologies, this invention has significant advantages. Attached Figure Description

[0049] Figure 1 This is a flowchart of the present invention. Detailed Implementation

[0050] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0051] Refer to the instruction manual appendix Figure 1 An embodiment of the present invention provides a method for co-source apportionment of nitrogen and phosphorus pollutants in a land-sea system, comprising:

[0052] S1. Obtain the isotopic characteristics of nitrogen and phosphorus sources of various forms in the land-sea system. Utilize elemental and stable isotope tracing technology and isotope mixing models to quantitatively analyze the contributions of river water, groundwater, coastal terrestrial sources, and estuarine endogenous sources to different forms of nitrogen and reactive phosphate in estuarine seawater. Different forms of nitrogen include ammonia nitrogen, nitrate nitrogen, nitrite nitrogen, dissolved organic nitrogen, and particulate organic nitrogen.

[0053] S2. Based on hydrological monitoring data and water quality monitoring data, quantitatively calculate the flux of nitrogen and phosphorus pollutants of various forms in rivers into the sea and their spatiotemporal variation patterns; based on the flux of nitrogen and phosphorus pollutants in rivers and the results of source apportionment of estuary seawater, calculate the flux of nitrogen and phosphorus pollutants input into groundwater by percentage.

[0054] S3. Quantitatively analyze the contribution ratio of different forms of nitrogen and phosphorus pollutants in rivers and groundwater;

[0055] S4. Based on the source apportionment results of nitrogen and phosphorus pollutants in various forms from rivers and groundwater, calculate the input flux of nitrogen and phosphorus pollutants in different forms from each source into the bay seawater through river water and groundwater.

[0056] S5. Based on the above results, preliminary source apportionment results for various forms of nitrogen and phosphorus pollutants in the estuary are formed;

[0057] S6. Based on the source apportionment results of different forms of nitrogen pollutants and the percentage composition of each form of nitrogen from each pollution source, verify the accuracy of the nitrogen pollutant source apportionment results.

[0058] S7. Based on the source apportionment results of reactive phosphate and the nitrogen-phosphorus ratio of each pollution source, verify the accuracy of the phosphate source apportionment results;

[0059] S8. Based on the source apportionment results of phosphorus content percentage and reactive sulfate in various forms from different pollution sources, determine the input of each pollution source to dissolved organic phosphorus and particulate organic phosphorus.

[0060] S9. Based on the concept of integrated land and sea management, complete the source analysis of nitrogen and phosphorus pollutants of various forms in the estuary.

[0061] According to step S1, obtain the isotopic characteristics of nitrogen and phosphorus sources and analyze their contributions; formulate... The isotopic contribution ratio of nitrogen or phosphorus form i in pollution source j, corresponding to region k; weighted by each sample point. Isotope abundance ratio Calculate the contribution of source j;

[0062]

[0063] in Let i be the isotopic abundance weight of morphology i at the m-th sample point, taking into account the characteristics of region k; The isotopic characteristic ratio represents the stable isotopic ratio of source j in the m-th sample of region k, such as δ 15 N; The total flow contribution of the m-th sample point is calculated by combining the monitoring data of sample points in region k; N, M, and P represent the number of pollution sources, the number of sample points, and the number of morphologies, respectively; the denominator in the above formula is normalized for all samples of all pollution sources to ensure that the total contribution ratio is 1.

[0064] Quantitatively calculate the inflow flux into the sea and its spatiotemporal variations based on step S2; propose T i (t,z) represents the total sea flux of nitrogen or phosphorus form i at time t and depth z; Φ is proposed. i(x,t,z) represents the three-dimensional distribution of pollutant concentration. Combined with the dynamic change of flow rate Q(x,t), the contribution of each location to the total flux is calculated.

[0065]

[0066] Where Φ i (x,t,z) is the concentration distribution function of morphology i at location x, time t, and depth z; Q(x,t) is the flow distribution of the river at location x and time t; K z (x) is the depth correction coefficient, representing the weight of the influence of different depth layers on the flow at location x; L is the total river length; dx is the integral variable over the river length, ensuring that the formula can be extended from local, segment-by-segment calculations to the entire river length; correction coefficient K z (x) is used to account for the uneven distribution of river depth.

[0067] Establish a groundwater input calibration model based on step S3; formulate... The input flux of nitrogen or phosphorus form i from groundwater below depth z;

[0068]

[0069] T i (z) represents the total inflow flux of morphology i at depth layer z, which originates from step S2; C i,j The proportion of pollution source j to form i, which comes from step S1; L j The total load of pollution source j is determined based on monitoring data; Let be the total pollutant load of morphology i below depth z; J represents the total number of pollution sources; the above formula estimates the implicit input of groundwater by subtracting the known contribution of pollution sources.

[0070] Quantitative analysis of the contribution ratios from different sources is performed based on step S4; F is proposed. i,j The quantitative contribution ratio of pollution source j to nitrogen or phosphorus form i;

[0071]

[0072] Where Θ j Ψ represents the regional characteristic factor of pollution source j, indicating the pollution distribution characteristics of that source; j Θ represents the environmental impact factor of pollution source j, determined by combining climate and hydrodynamic influences; K represents the total number of pollution sources; k Ψ is a regional characteristic factor of pollution source k, used to represent the geographical distribution and emission characteristics of pollution source k, and to quantify the impact of the environmental characteristics of the area where the pollution source is located on the diffusion and migration of pollutants; kThe environmental impact factor of pollution source k is used to adjust the relative influence of the pollution source; the above formula adjusts the contribution ratio of groundwater input through characteristic factors and impact factors to determine the specific source.

[0073] Based on step S5, preliminary source analysis results are generated; R is proposed. i,j The preliminary analysis results of the nitrogen or phosphorus form i for pollution source j;

[0074] R i,j =F i,j ·T i

[0075] Where F i,j The quantitative contribution ratio of pollution source j to nitrogen or phosphorus form i is determined based on the results of step S4; T i Let be the total flux into the sea for morphology i, based on the results of step S2; the above formula provides a preliminary analysis by combining the proportions of each source with the total flux.

[0076] Verify the accuracy of nitrogen pollutant source apportionment according to step S6; formulate V i,j The error ratio of nitrogen pollutant speciation i to source j;

[0077]

[0078] Where R i,j These are preliminary analysis results; The above formula is used to compare the relative error between the model's analytical results and the measured values.

[0079] Verify the accuracy of phosphorus pollutant source apportionment according to step S7; formulate Q. i,j The accuracy factor for phosphorus pollutant speciation i and source j is determined by verifying and analyzing the results through the nitrogen-phosphorus ratio and measured phosphorus load.

[0080]

[0081] in The measured load of phosphorus pollutant form i in source j; N i,j is the nitrogen-to-phosphorus ratio verification value of source j for form i; k is the index of the pollution source, used to represent different pollution sources; K is the total number of pollution sources, representing the sum of all pollution sources involved in the formula; By measuring phosphorus load from all sources and its corresponding nitrogen-to-phosphorus ratio N i,k We obtain the overall verification benchmark by weighted summation.

[0082] Based on step S8, determine the sources of other forms of phosphorus; formulate L i,j Input amounts of dissolved organophosphates and particulate organophosphates for source j to form i; calibration of other phosphorus forms based on reactive phosphate results;

[0083] L i,j =Q i,j ·α i,j

[0084] Q i,j α is an accuracy factor for determining the source j of phosphorus pollutant speciation i; i,j The conversion factor is used to represent the rate at which reactive phosphate is converted to other forms of phosphorus.

[0085] Complete the final precise source analysis according to step 9; formulate S i Precise source apportionment results for nitrogen or phosphorus pollutants in final form i

[0086]

[0087] Where R i,j As preliminary source resolution results, according to step S5; L i,j The input amount of dissolved organophosphates and particulate organophosphates of source j to form i; J is the total number of pollution sources; the final result is obtained by combining preliminary analysis and supplementary amounts.

[0088] The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for co-source apportionment of nitrogen and phosphorus pollutants in land-sea systems, characterized in that, include: S1. Obtain the isotopic characteristics of nitrogen and phosphorus sources of various forms in the land-sea system. Utilize elemental and stable isotope tracing technology and isotope mixing models to quantitatively analyze the contributions of river water, groundwater, coastal terrestrial sources, and estuarine endogenous sources to different forms of nitrogen and reactive phosphate in estuarine seawater. Different forms of nitrogen include ammonia nitrogen, nitrate nitrogen, nitrite nitrogen, dissolved organic nitrogen, and particulate organic nitrogen. S2. Based on hydrological monitoring data and water quality monitoring data, quantitatively calculate the flux of nitrogen and phosphorus pollutants of various forms in rivers into the sea and their spatiotemporal variation patterns; based on the flux of nitrogen and phosphorus pollutants in rivers and the results of source apportionment of estuary seawater, determine the flux of nitrogen and phosphorus pollutants input into groundwater. S3. Quantitatively analyze the contribution ratio of different forms of nitrogen and phosphorus pollutants in rivers and groundwater; S4. Based on the source apportionment results of nitrogen and phosphorus pollutants in various forms from rivers and groundwater, calculate the input flux of nitrogen and phosphorus pollutants in different forms from each source into the bay seawater through river water and groundwater. S5. Based on the above results, preliminary source apportionment results for various forms of nitrogen and phosphorus pollutants in the estuary are formed; S6. Based on the source apportionment results of different forms of nitrogen pollutants and the percentage composition of each form of nitrogen from each pollution source, verify the accuracy of the nitrogen pollutant source apportionment results. S7. Based on the source apportionment results of reactive phosphate and the nitrogen-phosphorus ratio of each pollution source, verify the accuracy of the phosphate source apportionment results; S8. Based on the source apportionment results of phosphorus content percentage and reactive sulfate in various forms from different pollution sources, determine the input of each pollution source to dissolved organic phosphorus and particulate organic phosphorus. S9. Based on the concept of integrated land and sea management, complete the source analysis of nitrogen and phosphorus pollutants of various forms in the estuary; According to step S1, obtain the isotopic characteristics of nitrogen and phosphorus sources and analyze their contributions; formulate... The isotopic contribution ratio of nitrogen or phosphorus form i in pollution source j, corresponding to region k; weighted by each sample point. Isotope abundance ratio Calculate the contribution of source j; in Let i be the isotopic abundance weight of morphology i at the m-th sample point, taking into account the characteristics of region k; The characteristic ratio of isotopes; The total flow contribution of the m-th sample point is calculated by combining the monitoring data of sample points in region k; N, M, and P represent the number of pollution sources, the number of sample points, and the number of pollution types, respectively. Quantitatively calculate the inflow flux into the sea and its spatiotemporal variations based on step S2; propose T i (t,z) represents the total sea flux of nitrogen or phosphorus form i at time t and depth z; Φ is proposed. i (x,t,z) represents the three-dimensional distribution of pollutant concentration. Combined with the dynamic change of flow rate Q(x,t), the contribution of each location to the total flux is calculated. Where Φ i (x,t,z) is the concentration distribution function of morphology i at location x, time t, and depth z; Q(x,t) is the flow distribution of the river at location x and time t; K z (x) is the depth correction coefficient, representing the weight of the influence of different depth layers at location x on the flow; L is the total length of the river; Establish a groundwater input calibration model based on step S3; formulate... The input flux of nitrogen or phosphorus form i from groundwater below depth z; T i (z) represents the total inflow flux into the ocean at depth z for morphology i; C i,j The proportion of pollution source j to form i; L j The total load of pollution source j; denoted as the total pollutant load of morphology i below depth layer z; J represents the total number of pollution sources. Quantitative analysis of the contribution ratios from different sources is performed based on step S4; F is proposed. i,j The quantitative contribution ratio of pollution source j to nitrogen or phosphorus form i; Where Θ j Ψ is a regional characteristic factor of pollution source j; j Θ represents the environmental impact factor of pollution source j; K represents the total number of pollution sources; k Ψ is a regional characteristic factor of pollution source k; k Environmental impact factors of pollution source k; Based on step S5, preliminary source analysis results are generated; R is proposed. i,j The preliminary analysis results of the nitrogen or phosphorus form i for pollution source j; R i,j =F i,j ·T i Where F i,j T represents the quantitative contribution of pollution source j to nitrogen or phosphorus form i; i Let i be the total flux into the sea for form i.

2. The method for co-source apportionment of nitrogen and phosphorus pollutants in land-sea systems according to claim 1, characterized in that: Verify the accuracy of nitrogen pollutant source apportionment according to step S6; formulate V i,j The error ratio of nitrogen pollutant speciation i to source j; Where R i,j These are preliminary analysis results; These are the results of actual measurements and analysis.

3. The method for co-source apportionment of nitrogen and phosphorus pollutants in land-sea systems according to claim 2, characterized in that: Verify the accuracy of phosphorus pollutant source apportionment according to step S7; formulate Q. i,j The accuracy factor for phosphorus pollutant speciation i and source j is determined by verifying and analyzing the results through the nitrogen-phosphorus ratio and measured phosphorus load. in The measured load of phosphorus pollutant form i in source j; N i,j is the nitrogen-to-phosphorus ratio verification value for source j relative to form i; k is the index of the pollution source; K is the total number of pollution sources; By measuring phosphorus load from all sources and its corresponding nitrogen-to-phosphorus ratio N i,k We obtain the overall verification benchmark by weighted summation.

4. The method for co-source apportionment of nitrogen and phosphorus pollutants in land-sea systems according to claim 3, characterized in that: Based on step S8, determine the sources of other forms of phosphorus; formulate L i,j For the input amounts of source j to form i of dissolved organophosphates and particulate organophosphates; based on The results of reactive phosphate were used to calibrate the sources of other phosphorus forms. 50 i,j =Q i,j ·α i,j Q i,j α is an accuracy factor for determining the source j of phosphorus pollutant speciation i; i,j The conversion factor is used to represent the rate at which reactive phosphate is converted to other forms of phosphorus. Complete the final source parsing according to step S9; formulate S... i Precise source apportionment results for nitrogen or phosphorus pollutants in final form i Where R i,j This is the preliminary source resolution result; L i,j J represents the input amount of dissolved organophosphates and particulate organophosphates from source j to form i; J represents the total number of pollution sources.

Images