A method and device for analyzing the joint toxicity of prototheca reticulata and microplastics

Abstract

The application provides a kind of net primary angle algae and microplastic combined toxicity analysis method, comprising: constructing simulated water environment;Obtain the original data of various populations in the simulated water environment;Carry out 7-day ecotoxicology experiment, obtain the change data of various populations in the simulated water environment;According to the original data and change data of various populations, the combined toxicity analysis model is used to analyze the combined toxicity of the simulated water environment;The change of each type of population is predicted by using population prediction model, and population prediction trend chart is generated respectively.This analysis method has good portability, wide application range, good effect;It can also be used for the prevention and treatment of other water environment or other environmental harmful factors.The application provides a kind of combined toxicity analysis device, which integrates data acquisition, data calculation and processing, combined toxicity analysis, is easy to operate, has comprehensive functions, good prediction effect and high accuracy.

Description

Technical Field

[0001] This invention relates to a method and apparatus for combined toxicity analysis of Protoceratops reticulata and microplastics, belonging to the field of marine biology. Background Technology

[0002] How to effectively control red tides has become a goal that people are constantly exploring. At present, various control methods have been proposed at home and abroad. Exploring effective, economical and environmentally friendly methods to control the occurrence of red tides is an important task of current environmental science and an urgent need for current economic development and environmental protection. Analytical methods and models are a way to simulate the changing trends of relevant parameters in an ecosystem through mathematical modeling. Analytical methods and models for red tides provide a possible way to control them. Summary of the Invention

[0003] The purpose of this invention is to overcome the shortcomings of the prior art and provide a method for simulating the changing trends of biological populations through mathematical modeling, analyzing the harm of harmful factors to the ecosystem, and using the results of analysis and simulation to assist in decision-making and take measures for governance.

[0004] To achieve the above objectives, the present invention is implemented using the following technical solution:

[0005] On one hand, this invention provides a method for combined toxicity analysis of *Protoceratops reticulata* and microplastics, comprising: constructing a simulated aquatic environment containing microplastics and *Protoceratops reticulata* populations based on the actual water quality of the target ecosystem; obtaining raw data of various populations in the simulated aquatic environment; the raw data including: the original total dry weight of various populations, the total dry weight of all organisms, and the weight of each population; conducting ecotoxicological experiments on the simulated aquatic environment, analyzing the daily changes in the total dry weight of various populations in the simulated aquatic environment, and obtaining change data of various populations in the simulated aquatic environment; wherein, the change data of various populations includes: the average lifespan of each population in the simulated aquatic environment, the average daily decrease in dry weight of each population, and the average daily decrease rate of total dry weight of each population; and performing combined toxicity analysis on the simulated aquatic environment using a combined toxicity analysis model based on the raw data and change data of various populations in the simulated aquatic environment.

[0006] Furthermore, the analysis of the daily changes in the total dry weight of various populations in the simulated aquatic environment includes: collecting dead organisms deposited in the simulated aquatic environment daily; placing the collected dead organisms into centrifuge tubes, repeatedly centrifuging and washing with water, and then drying them; and using flow cytometry to analyze the dry weight of various populations in the dried dead organisms.

[0007] Furthermore, the various species populations in the simulated aquatic environment include: population A, population B, and population D; population A is the dominant species population; population B is the key species population; and population D is the species population that needs protection.

[0008] Furthermore, the combined toxicity analysis model is as follows:

[0009] ;

[0010] In the formula, W is the combined toxicity analysis value; The weight of population class A; The daily average decline rate of total dry weight in category A population; The weights for population type B; The daily average decline rate of total dry weight in the B-type population; The weights of the D-class population; α represents the daily average decline rate of total dry weight in population D; α, β, and λ are species differentiation inhibition factors corresponding to populations A, B, and D, respectively; and 7 represents the duration of the ecotoxicological experiment, which is 7 days.

[0011] Furthermore, the weight of class A population Weights of Class B populations and the weight of the D-class population The calculation formula is:

[0012] ;

[0013] ;

[0014] ;

[0015] In the formula, To simulate the total dry weight of all organisms in the aquatic environment; The original total stem weight of population A; The original total stem weight of the B-type population; The original total stem weight of the D-type population.

[0016] Furthermore, the daily average decline rate of total dry weight in category A populations Daily average decline rate of total dry weight in Category B populations Daily average decline rate of total dry weight in Category D population The calculation formula is:

[0017] ;

[0018] ;

[0019] ;

[0020] In the formula, The average daily dry weight decrease for Category A populations; The average daily dry weight decrease for population type B; For population type D, the daily average dry weight decrease To simulate the total dry weight of all organisms in the aquatic environment; The original total stem weight of population A; The original total stem weight of the B-type population; The original total dry weight of the D-type population; The rate of decline in life expectancy for individuals in the A-class population; The rate of decline in life expectancy for individuals in the B-type population; The lifespan decline rate of individuals in the D-class population.

[0021] Furthermore, the population lifespan decline rate = mean population lifespan in simulated aquatic environment / mean population lifespan in normal environment; wherein, the mean population lifespan in simulated aquatic environment is obtained from ecotoxicological experiments.

[0022] Furthermore, the analysis method also includes: using a population prediction model to predict the changes of populations A, B, C and D in the simulated aquatic environment over a preset n days, and generating population prediction trend charts respectively; where n=1, 2, 3...; population C is a reticulate protoceratops population.

[0023] Furthermore, the population prediction model is as follows:

[0024] Total dry weight of C-type population on day n Calculation formula: ;

[0025] Total dry weight of type A population on day n Calculation formula: ;

[0026] Total stem weight of type B population on day n Calculation formula: ;

[0027] Total dry weight of D-type population on day n Calculation formula: ;

[0028] Day n total population weight Calculation formula: ;

[0029] In the formula, The original total stem weight of population A; The original total stem weight of the B-type population; The original total dry weight of the C-type population; The original total dry weight of the D-type population; The average daily dry weight decrease for Category A populations; The average daily dry weight decrease for population type B; The average daily dry weight growth of the C-type population; The daily average dry weight decrease was for population type D.

[0030] On the other hand, the present invention provides a combined toxicity analysis device for *Procambarus reticulata* and microplastics, comprising:

[0031] Module: Used to construct a simulated aquatic environment containing microplastics and Protoceratops reticulata populations based on the actual water quality of the target ecosystem;

[0032] Data Acquisition Module: Used to collect raw data of various populations in the simulated aquatic environment; the raw data includes: the raw total dry weight of various populations, the total dry weight of all organisms, and the weight of each population.

[0033] Change Data Analysis Module: This module is used to analyze the daily changes in the total dry weight of various species in the simulated aquatic environment based on ecotoxicological experiments, and to obtain change data for various species in the simulated aquatic environment. The change data for various species include: the average lifespan of individuals in various species in the simulated aquatic environment, the average daily decrease in dry weight of various species, and the average daily decrease rate of total dry weight of various species.

[0034] Toxicity Analysis Module: Used to perform joint toxicity analysis on the simulated aquatic environment based on the original and change data of various populations in the simulated aquatic environment and using a joint toxicity analysis model.

[0035] Compared with the prior art, the beneficial effects achieved by the present invention are as follows:

[0036] (1) The method provided by the present invention has low implementation cost. Only one toxicology experiment needs to be completed to analyze and predict the hazards of reticulate protoceratops and microplastics in the marine environment, and the results of the analysis and prediction can be used to assist in decision-making and take measures for treatment.

[0037] (2) The method provided by the present invention has good effect, good portability, and wide application range. It can be used for the prevention and treatment of other water environment or other environmental harmful factors.

[0038] (3) The device provided by the present invention integrates data acquisition, data calculation and processing, and joint toxicity analysis. It is simple to operate, has comprehensive functions, good prediction effect, and high accuracy. Attached Figure Description

[0039] Figure 1 This is a flowchart of the combined toxicity analysis in an embodiment of the present invention;

[0040] Figure 2 This is the population prediction trend diagram for type C in Embodiment 1 of the present invention;

[0041] Figure 3 These are the population prediction trend diagrams for categories A, B, and D in Embodiment 1 of the present invention;

[0042] Figure 4 This is a trend chart of all population predictions in the simulated aquatic environment in Embodiment 1 of the present invention;

[0043] Figure 5 This is the population prediction trend diagram for type C in Embodiment 2 of the present invention;

[0044] Figure 6 These are the population prediction trend diagrams for categories A, B, and D in Embodiment 2 of the present invention;

[0045] Figure 7 This is a trend chart of all population predictions in the simulated aquatic environment in Embodiment 2 of the present invention; Detailed Implementation

[0046] The present invention will be further described below with reference to the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and should not be used to limit the scope of protection of the present invention.

[0047] This embodiment provides a method for combined toxicity analysis of *Procambarus reticulata* and microplastics, comprising: constructing a simulated aquatic environment; obtaining raw data of various populations in the simulated aquatic environment; conducting ecotoxicological experiments to obtain change data of various populations in the simulated aquatic environment; performing combined toxicity analysis on the simulated aquatic environment based on the raw data and change data of various populations using a combined toxicity analysis model; and predicting the changes of various populations using a population prediction model, and generating population prediction trend maps respectively.

[0048] Example 1

[0049] The application scenario of this embodiment is a seaside fishery. In this ecosystem (S), fish (H) are economic fish species, and their initial food is the key species (B type) *Trichogramma foldis*. The dominant species (A type) in this fishery is *Dichroa rubescens*, and *Chlorella vulgaris* is the food for rotifers and *Dichroa rubescens*. The species that needs protection (D type) is *Chlorella vulgaris*. There is a precedent for the outbreak and reproduction of *Procambarus reticulata* (C type) in this sea area.

[0050] In this embodiment, a simulated aquatic environment was prepared in the laboratory: 20 L of seawater was prepared, and the concentration of *Procambarus reticulatus* in the seawater was set to 1 × 10⁻⁶. 5 The concentration of Chlorella in seawater was set at 1 × 10⁻⁶ cells / mL. 5 Cells / mL, and simultaneously introduced 1µm diameter PS microplastics with an abundance of 0.5 µg / mL. 13,000 a-type species, Japanese tiger daphnia, and 7,000 b-type species, *Trichogramma foldis*, were placed in 20 L of seawater.

[0051] The seawater was thoroughly shaken and divided into four equal portions (20 L each). Two portions were used to record raw population data in the simulated aquatic environment. Raw data for various populations in the simulated aquatic environment were obtained by: placing each portion of seawater into centrifuge tubes, repeatedly centrifuging, washing with clean water, drying, and weighing to obtain the raw total biological dry weight. s The dry weight of populations A, B, C, and D was analyzed using flow cytometry. The average of the two sets of data was used to obtain the original total dry weight of populations A, B, C, and D.

[0052] The weights of populations A, B, and D are calculated using the following formula:

[0053] Weight of Class A population ;

[0054] Weight of Class B population ;

[0055] and the weight of the D-class population ;

[0056] The raw data of various populations in the simulated aquatic environment are recorded in Table 1:

[0057] Population raw data numerical values Original total biological dry weight 100g Original dry weight of population C 15g Original dry weight of Category A population 40g Original dry weight of population B 20g Original dry weight of population type D 25g Weight of Class A population 0.4 Weight of Class B population 0.2 Weight of Class D population 0.25

[0058] Two additional samples were used for ecotoxicology experiments. During the 7-day ecotoxicology experiment, two samples of dead organisms deposited at the bottom of the seawater were collected daily for analysis to obtain data on the changes in populations A, B, C, and D in the simulated aquatic environment. The specific method was as follows: the dead organisms collected daily were placed in centrifuge tubes, repeatedly centrifuged, washed with water, and then dried. Flow cytometry was used to analyze the dry weight of populations A, B, C, and D in the dead organisms. The average of the two sets of data was taken to obtain the daily average decrease in dry weight for populations A, B, and D. Simultaneously, the average lifespan of individuals in the simulated aquatic environment was observed under a microscope, and the lifespan decline rate of populations A, B, and D was calculated using the formula: Population lifespan decline rate = Average lifespan of individuals in the simulated aquatic environment / Average lifespan of individuals in the normal environment. Further calculations were then performed.

[0059] Daily average decline rate of total dry weight in Category A population: ;

[0060] Daily average decline rate of total dry weight in category B population: ;

[0061] Daily average decline rate of total dry weight in category D population: ;

[0062] The data on the changes in various populations in the simulated aquatic environment are recorded in Table 2:

[0063] Population change data Value / g Daily average dry weight growth of Category C population 3g The average daily dry weight of Category A population decreased 4g The average daily dry weight decrease in Category B populations 1g The average daily dry weight decrease in the D-type population 2.5g The rate of decline in lifespan of individuals in Category A populations 40% The rate of decline in lifespan of individuals in the B population 25% The rate of decline in lifespan of individuals in category D populations 30% Daily average growth rate of total dry weight of C-type population 0.2 Daily average growth rate of total dry weight of Category A population 0.117 Daily average growth rate of total dry weight of Category B population 0.073 Daily average growth rate of total dry weight of D-type population 0.113

[0064] Based on the data in Tables 1 and 2, the combined toxicity analysis values ​​of the simulated aquatic environment were calculated using a combined toxicity analysis model. The combined toxicity analysis model is as follows:

[0065] ;

[0066] In the formula, α, β, and λ are species differentiation inhibition factors corresponding to populations A, B, and D, respectively.

[0067] Different organisms have varying degrees of tolerance to combined toxicity. Populations with higher tolerance have lower inhibitory factors. The species-differentiated inhibitory factor is determined by selecting two parameters from the population life table: net reproductive rate and generation time. Under simulated aquatic conditions, the net reproductive rate of population A is calculated separately. aR0 Net reproductive rate of population type B bR0 Net reproductive rate of population type D dR0 Generation time of type A population aT Generation time of type B population bT Generation time of type D population dT We selected the one that best reflects the differences between populations of types A, B, and D as the quantitative standard to quantify the inhibitory factor of species differentiation.

[0068] In this embodiment, organism A aR0 A 20% decrease, B organism bR0 A 30% decrease, D-Bio dR0 A 40% decrease in net reproductive rate significantly reflects the differences between populations A, B, and D, and is therefore selected as the quantitative standard. α / β / λ = aR0 / bR0 / dR0 α is 1, β is 1.5, and λ is 2.

[0069] In this embodiment, the combined toxicity analysis value W=0.8764. Based on the value of W, further analytical indicators can be formulated to classify the combined toxicity hazards of Protoceratops reticulata and microplastics. Referring to the analytical indicators in Table 3 of this embodiment, the combined toxicity hazard level of Protoceratops reticulata and microplastics in this coastal fishing ground is moderate.

[0070] value range Hazard level Recommendations 0 ~ 0.6 Mild Strengthen supervision 0.6 ~ 1 moderate Active governance Greater than 1 Severe Isolate the area to prevent further damage.

[0071] Based on the data in Tables 1 and 2, a population prediction model is used to predict the changing trends of populations A, B, C, and D in the simulated aquatic environment over the next 30 days. The population prediction model is as follows:

[0072] Total dry weight of C-type population on day n Calculation formula: ;

[0073] Total dry weight of type A population on day n Calculation formula: ;

[0074] Total stem weight of type B population on day n Calculation formula: ;

[0075] Total dry weight of D-type population on day n Calculation formula: ;

[0076] Day n total population weight Calculation formula: ;

[0077] In the formula, The original total stem weight of population A; The original total stem weight of the B-type population; The original total dry weight of the C-type population; The original total dry weight of the D-type population; The average daily dry weight decrease for Category A populations; The average daily dry weight decrease for population type B; The average daily dry weight growth of the C-type population; The daily average dry weight decrease was for population type D.

[0078] Figure 2 The predicted results for Category C populations indicate that the control of the *Procambarus reticulata* red tide should be carried out before day 20. After day 30, the red tide has already broken out, and environmental management departments need to take measures to prevent the spread of the impact and minimize the losses.

[0079] Figure 3 (a) The subplot shows the predicted results for population A. The trend of population A indicates that protection of dominant species A should be carried out in the first 10 days. After the 10th day, the population of dominant species A will decline significantly, and the nature of the ecosystem will change. Environmental management departments need to take measures to prevent the spread of impact and minimize losses.

[0080] Figure 3 (b) Subplot shows the predicted results for the B population. The trend of the B population change indicates that protection of critical species B should be carried out in the first 20 days. After the 20th day, the population of critical species B will decline sharply, and the ecosystem will face collapse. Environmental management departments need to take measures to prevent the spread of impact and minimize losses.

[0081] Figure 3(c) The subplot shows the predicted results of the D-type population. The trend of the D-type population indicates that the protection of species D should be carried out in the first 10 days. After the 10th day, the population of species D that needs protection has declined sharply and is in danger. The environmental management department needs to take measures to transfer the species to a new sea area for protection.

[0082] Figure 4 For the prediction results of ecosystem S, Figure 4 The changing trends indicate that after day 6, the total biomass of the ecosystem exceeds the original total biomass. After day 10, the second derivative of the predicted curve is large, indicating that Protoceratops reticulata begins to reproduce rapidly.

[0083] In summary, taking measures to treat the target water quality before the 10th day can achieve better results.

[0084] Example 2

[0085] The application scenario of this embodiment is a marine protected area, in which the ecosystem (S) includes dominant fish species (Category A), key food species (Category B), other species unaffected by combined toxicity (Category E), rare and protected fish species (Category D), and Protoceratops reticulata (Category C).

[0086] In this embodiment, a simulated seawater environment was constructed by on-site sampling. The sample was taken by boat to the protected area, and 10 representative areas were selected. Two 100 L seawater samples were taken from each area. One seawater sample from each area was mixed to obtain two 1 m³ simulated seawater environments. One of the two simulated water environments was used to measure and record the original population data in the simulated water environment, and the other was used to conduct a 7-day ecotoxicological experiment.

[0087] In this embodiment, the method for obtaining the original data of various populations in the simulated aquatic environment is the same as in Embodiment 1.

[0088] The raw data of various populations in the simulated aquatic environment are recorded in Table 4:

[0089] Population raw data numerical values Original total biological dry weight 100g Original dry weight of population C 10g Original dry weight of Category A population 40g Original dry weight of population B 30g Original dry weight of population type D 20g Weight of Class A population 0.4 Weight of Class B population 0.3 Weight of Class D population 0.2

[0090] The method for obtaining change data of various populations in the simulated aquatic environment is the same as in Example 1.

[0091] The data on the changes in various populations in the simulated aquatic environment are recorded in Table 5:

[0092] Population change data Value / g Daily average dry weight growth of Category C population 1g The average daily dry weight of Category A population decreased 4g The average daily dry weight decrease in Category B populations 0.5g The average daily dry weight decrease in the D-type population 1g The rate of decline in lifespan of individuals in Category A populations 30% The rate of decline in lifespan of individuals in the B population 25% The rate of decline in lifespan of individuals in category D populations 20% Daily average growth rate of total dry weight of C-type population 0.2 Daily average growth rate of total dry weight of Category A population 0.113 Daily average growth rate of total dry weight of Category B population 0.111 Daily average growth rate of total dry weight of D-type population 0.07

[0093] Based on the data in Tables 4 and 5, the combined toxicity analysis values ​​of the simulated aquatic environment were calculated using the combined toxicity analysis model, which is the same as in Example 1.

[0094] In this embodiment, population type A aT A 20% decrease, Category B population bT A 30% decrease, Class D population dT A 40% decrease, and generation time significantly reflects the differences between populations A, B, and D, thus it was chosen as the quantitative standard, α / β / λ= aT / bT / dT α is 1, β is 1.5, and λ is 2.

[0095] In this embodiment, the combined toxicity analysis value W=0.859. Based on the value of W, further analytical indicators can be formulated to classify the combined toxicity hazards of Protoceratops reticulata and microplastics. Referring to the analytical indicators in this embodiment, the combined toxicity hazard level of Protoceratops reticulata and microplastics in this coastal fishing ground is moderate.

[0096] The data from Tables 4 and 5 were input into the population prediction model to simulate and predict the changing trends of populations A, B, C, and D in the simulated aquatic environment over the next 60 days. The population prediction model was the same as in Example 1.

[0097] Figure 5 The simulation results for Category C populations indicate that the control of the red tide of Protoceratops reticulata should be carried out between 30 and 40 days. After 40 days, the red tide has already broken out, and environmental management departments need to take measures to prevent the spread of the impact and minimize the losses.

[0098] Figure 6 (a) The subplot shows the simulation prediction results of population A. The trend of population A indicates that protection of dominant species A should be carried out in the first 20 days. After the 20th day, the population of dominant species A will decline significantly, the nature of the ecosystem will change, and environmental management departments need to take measures to prevent the spread of impact and minimize losses.

[0099] Figure 6 (b) The subplot shows the simulation prediction results of the B population. The trend of the B population change indicates that the protection of the key food species B should be carried out in the first 18 days. After the 18th day, the population of key species B will decline severely, and the ecosystem will face collapse. Environmental management departments need to take measures to prevent the impact from spreading and minimize the losses.

[0100] Figure 6 (c) The subplot shows the simulation prediction results of the D-type population. The trend of the D-type population change indicates that the protection of the rare and protected species D should be carried out before day 30. After day 30, the population of the protected species DD will decline severely and face danger. The environmental management department needs to take measures to transfer the species to a new sea area for protection.

[0101] Figure 7 The simulation prediction results for ecosystem S Figure 7The changing trends indicate that after day 16, the total biomass of the simulated aquatic environment exceeds the original total biomass, and after day 18, the second derivative of the predicted curve is large, indicating that Protoceratops reticulata has just begun to reproduce rapidly.

[0102] In summary, taking measures to treat the target water quality before the 18th day can achieve better results.

[0103] Example 3

[0104] This embodiment provides a device for combined toxicity analysis of *Procambarus reticulata* and microplastics, comprising: a construction module for constructing a simulated aquatic environment as described in Example 1 or Example 2; a data acquisition module for acquiring raw data of various populations in Example 1 or Example 2; a change data analysis module for obtaining change data of various populations in Example 1 or Example 2; a toxicity analysis module for performing combined toxicity analysis on the simulated aquatic environment using the combined toxicity analysis model described in Example 1 or Example 2; and a prediction module for predicting the future development trend of populations in the simulated aquatic environment using the population prediction model described in Example 1 or Example 2, and generating a corresponding population prediction trend map.

[0105] The above description is only a preferred embodiment of the present invention. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the technical principles of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for combined toxicity analysis of *Procambarus reticulata* and microplastics, characterized in that, Includes the following steps: Based on the actual water quality of the target ecosystem, a simulated aquatic environment containing microplastics and reticulate protoceratops populations was constructed. Obtain raw data for various populations in the simulated aquatic environment; the raw data includes: the raw total dry weight of each population, the total dry weight of all organisms, and the weights of each population. Ecotoxicological experiments were conducted on a simulated aquatic environment to analyze the daily changes in the total dry weight of various species populations in the simulated aquatic environment, and to obtain change data for various species populations in the simulated aquatic environment. Among them, the change data for various species populations include: the rate of decline in individual lifespan of various species populations, the average daily decrease in dry weight of various species populations, and the average daily decline in total dry weight of various species populations. Based on the original and change data of various species in the simulated aquatic environment, a joint toxicity analysis model was used to conduct a joint toxicity analysis on the simulated aquatic environment. The simulated aquatic environment includes three population groups: Group A, Group B, and Group D. Group A consists of dominant species; Group B consists of key species; and Group D consists of species requiring protection. The combined toxicity analysis model is as follows: ； In the formula, W is the combined toxicity analysis value; The weight of population class A; The daily average decline rate of total dry weight in category A population; The weights for population type B; The daily average decline rate of total dry weight in the B-type population; The weights of the D-class population; α represents the daily average decline rate of total dry weight in population D; α, β, and λ are species differentiation inhibition factors corresponding to populations A, B, and D, respectively; where 7 represents the duration of the ecotoxicological experiment, which is 7 days. The weight of the A-class population Weights of Class B populations and the weight of the D-class population The calculation formula is: ； ； ； In the formula, To simulate the total dry weight of all organisms in the aquatic environment; The original total stem weight of population A; The original total stem weight of the B-type population; The original total dry weight of the D-type population; The average daily decline rate of total dry weight of the Class A population Daily average decline rate of total dry weight in Category B populations Daily average decline rate of total dry weight in Category D population The calculation formula is: ; ； ； In the formula, The average daily dry weight decrease for Category A populations; The average daily dry weight decrease for population type B; For population type D, the daily average dry weight decrease To simulate the total dry weight of all organisms in the aquatic environment; The original total stem weight of population A; The original total stem weight of the B-type population; The original total dry weight of the D-type population; The rate of decline in life expectancy for individuals in the A-class population; The rate of decline in life expectancy for individuals in the B-type population; The lifespan decline rate of individuals in the D-class population.

2. The method as described in claim 1, characterized in that, The analysis simulates the daily changes in the total dry weight of various populations in the aquatic environment, including: Collect dead organisms deposited in a simulated aquatic environment every day; The collected dead organisms were placed in centrifuge tubes, repeatedly centrifuged, washed with clean water, and then dried. The dry weight of various populations in dried dead organisms was analyzed using flow cytometry.

3. The method as described in claim 1, characterized in that, The population lifespan decline rate = mean lifespan of the population under simulated aquatic conditions / average lifespan of the population under normal conditions; wherein, the mean lifespan of the population under simulated aquatic conditions is obtained from ecotoxicological experiments.

4. The method as described in claim 1, characterized in that, Also includes: The population prediction model was used to predict the changes of populations A, B, C and D in a simulated aquatic environment over a predetermined period of n days, and population prediction trend charts were generated for each population. Among them, n=1, 2, 3...; population C is the Protoceratops reticulata population.

5. The method as described in claim 4, characterized in that, The population prediction model is as follows: Total dry weight of C-type population on day n Calculation formula: ; Total dry weight of type A population on day n Calculation formula: ; Total stem weight of type B population on day n Calculation formula: ; Total dry weight of D-type population on day n Calculation formula: ; Day n total population weight Calculation formula: ; In the formula, The original total stem weight of population A; The original total stem weight of the B-type population; The original total dry weight of the C-type population; The original total dry weight of the D-type population; The average daily dry weight decrease for Category A populations; The average daily dry weight decrease for population type B; The average daily dry weight growth of the C-type population; The daily average dry weight decrease was for population type D.

6. A combined toxicity analysis device for *Procambarus reticulata* and microplastics, characterized in that, The apparatus for performing the method according to any one of claims 1 to 5 comprises: Module: Used to construct a simulated aquatic environment containing microplastics and Protoceratops reticulata populations based on the actual water quality of the target ecosystem; Data Acquisition Module: Used to collect raw data of various populations in the simulated aquatic environment; the raw data includes: the raw total dry weight of various populations, the total dry weight of all organisms, and the weight of each population. Change Data Analysis Module: This module is used to analyze the daily changes in the total dry weight of various species in the simulated aquatic environment based on ecotoxicological experiments, and to obtain change data for various species in the simulated aquatic environment. The change data for various species include: the average lifespan of individuals in various species in the simulated aquatic environment, the average daily decrease in dry weight of various species, and the average daily decrease rate of total dry weight of various species. Toxicity Analysis Module: Used to perform joint toxicity analysis on the simulated aquatic environment based on the original and change data of various populations in the simulated aquatic environment and using a joint toxicity analysis model.

Images