CN103902833B - A kind of estimation method of coastal ocean hypoxic condition - Google Patents

Abstract

The present invention discloses the estimation method of a kind of coastal ocean hypoxic condition, first obtains the data such as working region marine topography, water body spring layer, bottom dissolved oxygen and remote sensing top layer chlorophyll；Working region is carried out stress and strain model, by month and grid to historical data resampling, average statistical and minima；Obtain and be positioned at the time series data that the upper and lower buoy group of coastal ocean spring layer is measured in real time, calculate equivalence chlorophyll test value；With buoy point and the dissolved oxygen value of other mesh points of statistical correlation models coupling equivalence chlorophyll test value evaluation work region；Calculate hypoxia scope with the dissolved oxygen value of grid each point, calculate hypoxia degree with working region minimum dissolved oxygen value.The present invention utilizes hypoxia single-point buoy observed data, in conjunction with chlorophyll Ocean Color Remote Sensing view data, gives hypoxia scope and the hypoxia degree of working region.

Description

A kind of estimation method of coastal ocean hypoxic condition

Technical field

The present invention relates to the measuring method of a kind of coastal ocean hypoxic condition, particularly relate to a kind of integrated buoy and observe distant with satellite Sense data carry out the measuring method of target area bottom dissolved oxygen, the most in an economical manner, by setting up statistics warp Testing model, the top layer chlorophyll information utilizing satellite remote sensing to obtain calculates target area water bottom dissolved oxygen levels, timely and effective Ground obtains the relevant information such as hypoxia scope and degree.

Background technology

Anoxia (Hypoxia) refers in water environment that the content of oxygen is in reduced levels or oxygen and is consumed in a large number.Generally by dissolved oxygen (Dissolved oxygen, the DO) concentration water body less than 3.0mg/L is referred to as hypoxia water body.The hypoxia phenomenon of immediate offshore area, Usually occur in river mouth and offshore shelf marine site, the coastal waters eutrophication that it causes with mankind's activity, and enter Haihe River water each season The hydrographic factors such as run-off, nutrient substance, monsoon and the upper up-flow of joint are closely related, belong to seasonal hypoxia.Inshore hypoxia Phenomenon started just to have become as from the fifties in last century affects one of disaster that the marine eco-environment develops in a healthy way, occurrence frequency from Only monitor hypoxia before nineteen fifty to occur in nearly 20 sea areas, to global range at the beginning of 21 century, be no less than 400 There is hypoxia in sea area, area coverage is more than 245,000km², it is concentrated mainly on that population density is high, nutrient substance discharge capacity is big Inshore, the Northern Hemisphere (Diaz and Rosenberg, 2008).

Dissolved oxygen is ambient parameter important in marine ecosystems, is the material conditions that in ocean, most of biologies are depended on for existence, To maintaining, marine ecosystems are extremely important.The distribution of Dissolved Oxygen in Seawater concentration, change and temperature, salinity, biological activity and Water sports etc. are closely related, significant to the state of ecological environment understanding sea area.The appearance of hypoxic disorders often results in Stratobios amount reduces, local species miniaturization, and causes Fish prevalence to increase, and bio-diversity reduces, and ultimately causes life The change of state system structure and destruction.In China, hypoxia phenomenon is the seasonal environmental problem in entrance of Changjiang River and surrounding waters thereof.90 Age Mo, Li Daoji etc. (2002) is at Changjiang Estuary domain discovery Dissolved Oxygen in Water as little as 1.0mg/l, and Dissolved Oxygen in Water 2.0 The anoxic zone area of mg/l is up to 13,700km2.Chen Zhendong etc. (2007) think, the hypoxia of entrance of Changjiang River and surrounding waters thereof is existing Relevant as being deposited to bottom water body with the Organic substance that this region source of students Organic substance and river carry.This region seasonality hypoxia is directed to (Lu Yuan is included to Changjiang Diluted Water diffusion, the invasion of stratification intensity, Taiwan warm current's water, upper up-flow, nutritive salt supply, organic matter Organic) level transports and the factor relevant (opening through et al., 2010) such as continental shelf width, landform.

The formation of hypoxia phenomenon is a complicated process, by the joint effect of many factors.It is generally believed that water body stratification is scarce The external physical background that oxygen is formed, bottom organic matter decomposition oxygen consumption is the biogeochemistry endogenous cause of ill of anoxia formation and development (Rabalais, et al., 2010;Zhu et al., 2011).The large amount of organic that deposit or water bottom contain, There is Degradation under certain condition, consume a large amount of oxygen, cause dissolved oxygen to reduce.If but water body exchange is very well, biological fall Solving the oxygen consumed can be supplemented timely, ocean also can not form hypoxia phenomenon.So biological factor is not to cause ocean low The essential condition of oxygen.In ocean, the formation of low-oxygen area also needs to the existence of physical factor, and such as, water body stratification makes bottom water body difficult Swap with the water body higher with upper strata dissolved oxygen, can not get after consuming effectively supplementing so that bottom dissolves, cause low-oxygen area Formed.This stratification includes salinity stratification, temperature stratification etc., typically in the season that waters, river mouth temperature is higher, a large amount of high temperature, Less salt, low-density fresh water cover on top layer, easily make bottom high salt, high density sea water form independent water body, cause stratification to be made With, make top layer oxygen be difficult to bottom and exchange.Therefore in ocean the formation of low-oxygen area also with the topography and geomorphology in waters, flow field, temperature The physical factors such as degree are relevant (Wang, 2009).Thus, water column (body) stratification refers to be formed halocline (or pycnocline) interface, The vertical exchange of water body is prevented from, and is the essential condition of anoxia formation.Typically in the spring that wind-force is more weak, a large amount of diluted water from River mouth enters bay, is easily formed halocline (or pycnocline), stops the vertical mixing of water body and exchange, along with the temperature rises, layer It is turned into be used in and peaks summer.The river mouth bottom water seasonal hypoxias such as entrance of Changjiang River only occur in the water body of below pycnocline, and one As be formed at the end of spring and the beginning of summer, peak in midsummer, at the beginning of autumn late summer terminate, there is obvious seasonal characteristic.

River mouth seasonal hypoxia phenomenon has year and Seasonal, and the scope of region hypoxia, thickness and hypoxia degree are by many Plant the control of envirment factor.Observe hypoxia phenomenon: the feature of Satellite Remote Sensing is that spatial coverage is big, and synchronousness is good, The resolution of space-time is higher；Boats and ships monitoring has monitoring parameter and at most and obtains the advantage that the quality of data is the highest, but its space covers Scope is few, and synchronousness is poor, and spatial and temporal resolution is relatively low；Anchor system buoy is to realize anoxic zone water body section ecology to join with environment Unique effective means of number long term monitoring, but it is restricted by platform and the big technical factor of sensor two.Thus, scientifically and rationally The spatial distribution of cognitive hypoxia and hypoxia degree, be the most effectively to utilize various monitoring measures, analyzes hypoxia development and drills Law and variation tendency thereof.

Summary of the invention

It is an object of the invention to for the deficiencies in the prior art, it is provided that the measuring method of a kind of coastal ocean hypoxic condition.

It is an object of the invention to be achieved through the following technical solutions: the measuring method of a kind of coastal ocean hypoxic condition, including Following steps:

(1) the history survey data for many years such as working region marine topography, water body spring layer, bottom dissolved oxygen and remote sensing top layer chlorophyll is obtained；

(2) working region is carried out stress and strain model, by month and grid to historical data resampling, average statistical and minima；

(3) time sequences such as being positioned at temperature that the upper and lower buoy group of coastal ocean spring layer measures in real time, salinity, dissolved oxygen, chlorophyll is obtained Column data；

(4) with the degree of depth and the temperature measured in real time, the salinity annual average calculating spring layer intensity of buoy；

(5) obtain the chlorophyll remote sensing average view data with the buoy time series same period in the range of working region, calculate equivalence chlorophyll Value；

(6) the statistical correlation model of any two points in grid is set up with Dissolved Oxygen in Water balance；

(7) with buoy point and the dissolved oxygen value of other mesh points of statistical correlation models coupling equivalence chlorophyll test value evaluation work region；

(8) dissolved oxygen value with grid each point calculates the area of hypoxia, calculates hypoxia degree with working region minimum dissolved oxygen value.

Further, in described step (6), based on the dissolved oxygen equilibrium equation of the following water body of coastal ocean spring layer, it is considered to Water body initial dissolution oxygen degree, the oxygen exchange of vertical and horizontal directions and the Organic substance oxygen consumption of water body and substrate, adjust in conjunction with history Consult reference materials, set up dissolved oxygen statistical model, calculate the oxygen consumption between diverse location with the difference of equivalence chlorophyll concentration on different mesh points Difference.

The invention has the beneficial effects as follows, The present invention gives the observation of integrated single-point buoy and satellite remote sensing date is carried out at the bottom of target area The measuring method of layer dissolved oxygen, in an economical manner, by setting up statistics empirical model, utilizes single-point buoy monitor in real time and defend The top layer chlorophyll data that star remote sensing obtains calculate target area bottom dissolved oxygen levels, obtain relevant with hypoxia timely and effectively The relevant information such as area and degree.

Accompanying drawing explanation

Fig. 1 is the framework map of the estimation method of coastal ocean hypoxic condition of the present invention；

Fig. 2 is the FB(flow block) of the measuring method of coastal ocean hypoxic condition；

Fig. 3 is the schematic diagram of Application Gridization measuring and calculating target area coastal ocean hypoxic condition.

Detailed description of the invention

Dissolved oxygen is ambient parameter important in marine ecosystems, is the material conditions that in ocean, most of biologies are depended on for existence, To maintaining, marine ecosystems are extremely important.The water body that dissolved oxygen concentration is generally less than 3.0mg/L is referred to as hypoxia water body.Hypoxia The formation of phenomenon is a complicated process, by the joint effect of many factors.It is generally believed that the formation of water hypoxia and water body Stratification is relevant with bottom organic matter decomposition oxygen consumption.The former is the external physical background that anoxia is formed, and the latter is anoxia formation and development Biogeochemistry endogenous cause of ill.

During seasonal hypoxia, there is stratification phenomena in target area water body, makes top layer oxygen be difficult to bottom dissolved oxygen and exchanges.Profit With single-point monitering buoy, combine satellite remote sensing water colour data and can calculate the dissolved oxygen levels of other positions of bottom water body, target area. Coastal ocean seasonality hypoxic condition is speculated, including hypoxia scope and hypoxia degree.As it is shown in figure 1, system includes low Oxygen oceanographic buoy, ground satellite station, the buoy information processing terminal, human-computer interaction terminal, the hypoxia information processing terminal and communication network Network is constituted.Hypoxia oceanographic buoy utilizes data link and mail server that Monitoring Data sends to the buoy information processing terminal, water Ocean Color Remote Sensing data are sent to ground satellite station by color remote sensing satellite, and ground satellite station, mail server and the information processing terminal are logical Cross the Internet to be connected with each other.

The estimation method of the present invention a kind of coastal ocean hypoxic condition, comprises the following steps:

1, the history survey data for many years such as working region marine topography, water body spring layer, bottom dissolved oxygen and remote sensing top layer chlorophyll is obtained；

Determine target operation regional extent, obtain the marine topography of working region, different the water body spring layer in month, bottom dissolved oxygen And the history survey data for many years such as remote sensing top layer chlorophyll, water body spring layer information include hypoxia occur the spring layer intensity during month and Depth of pycnocline.

2, working region is carried out stress and strain model, by month and grid to historical data resampling, average statistical and minima；

Target operation regional extent is carried out spatial gridding subdivision by longitude and latitude direction.During hypoxia occurs month, by month And mesh point carries out resampling to the depth of water, spring layer intensity, depth of pycnocline, bottom dissolved oxygen, top layer chlorophyll, calculate grid model Enclose average and the minima of interior each mesh point above-mentioned parameter；Calculate the mean depth at dissolved oxygen minima place, it is determined as with reference to water Deeply；Calculate the bottom dissolved oxygen spring layer strength mean value less than 3.0mg/L position, it is determined as with reference to spring layer intensity, obtain working region Bottom dissolved oxygen is less than hypoxia month and the persistent period of 3.0mg/L.

3, time serieses such as being positioned at temperature that the upper and lower buoy group of coastal ocean spring layer measures in real time, salinity, dissolved oxygen, chlorophyll is obtained Data；

Buoy group is positioned at target operation regional extent, it is to avoid be positioned at dissolved oxygen minimum and vicinity region in historical data； Buoy group at least includes 2 water body subsurface buoys, arranges pre-test water body depth of pycnocline in subsurface buoy, and water body subsurface buoy lays respectively at buoy institute In the both sides up and down of a depth of pycnocline position, record the degree of depth of upper and lower buoy；The measured value of subsurface buoy at least include temperature, salinity, Dissolved oxygen, chlorophyll；Measuring frequency one day more than 2 times, measurement interval is uniform, and the time series of measured value was more than 3 days, extremely Include that hypoxia speculates before and after the date each 1 day less.

4, with the degree of depth and the temperature measured in real time, the salinity annual average calculating spring layer intensity of buoy；

Calculate buoy group and be positioned at the depth difference of water body spring layer both sides subsurface buoy, calculate temperature, the annual average of salinity measured by 2 subsurface buoys, At least include that hypoxia speculates before and after the date each 1 day, it is judged that whether spring layer intensity is more than with reference to spring layer intensity.

5, obtain the chlorophyll remote sensing average view data with the buoy time series same period in the range of working region, calculate equivalence chlorophyll test value；

Obtaining the chlorophyll satellite remote sensing date with the buoy time series same period in target operation regional extent, mean data can be Monthly average value, first quarter moon meansigma methods etc., the time of the moon or first quarter moon must not be later than hypoxia first 1 day of the date of supposition, when chlorophyll satellite When remote sensing spatial resolution is inconsistent with the locus of gridding, chlorophyll data are carried out gridding resampling, calculate equivalence Chlorophyll.

Dissolved oxygen depletion amount R that Organic substance causes_oRelevant with bottom organic concentration, can be expressed as:

R₀=r_ogk_gg

Wherein, r_ogFor the carbon-to-oxygen ratio of Organic substance oxygen consumption, k_gFor phosphorus content coefficient in Organic substance, g is bottom Organic substance amount, again with top layer Chlorophyllous function is expressed:

g=h×C^chl+g'

Wherein, C^chlBeing the top layer chlorophyll data of remote sensing observations, h is the top layer chlorophyll estimation organic coefficient of bottom, and g' is one Little correction, negligible.

The zmount of oxygen consumption that identical Organic substance causes in the water body of different depth is different, and Organic substance oxygen consumption needs to carry out adjusting for depth. Thus, employing depth value calculating chlorophyll test value correction bottom Organic substance oxygen consumption of equal value:

k^chl=h×|d-d₀|/d_o×(r_ogk_g)

Wherein, k^chlIt is revised bottom organic coefficient of oxygen consuming of equal value with top layer chlorophyll, d_oFor with reference to the depth of water, d is for waiting to count Calculate the mesh point depth of water, other coefficients (hr in availability coefficient k Unified Expression formula_ogk_g)。

6, the statistical correlation model of any two points in grid is set up with Dissolved Oxygen in Water balance；

Setting up spring layer following Dissolved Oxygen in Water equilibrium equation, dissolved oxygen equilibrium equation includes the horizontal O of initial dissolution oxygen_pre(mg/L), Dissolved oxygen vertical proliferation effect O_ver(mg/L), dissolved oxygen levels spreading effect O_lat(mg/L), dissolved oxygen water body oxygen consumption O_wc (mg/L), dissolved oxygen Oxygen Consumption By Sediments O_sed(mg/L):

O_obs=O_pre+O_ver+O_lat-O_wc-O_sed

Wherein, O_obsIt it is Dissolved Oxygen in Water measured value.

SetBottom dissolved oxygen average statistical is investigated, for the arbitrary mess point a in the range of working region, with target for history Dissolved oxygen relationship expression between the buoy point b place mesh point arranged in the range of working region is as follows:

$O_{\mathrm{pre}}^a-O_{\mathrm{pre}}^b\≈{\stackrel{\‾}{O}}_{\mathrm{pre}}^a-{\stackrel{\‾}{O}}_{\mathrm{pre}}^b$

$O_{\mathrm{sed}}^a-O_{\mathrm{sed}}^b\≈{\stackrel{\‾}{O}}_{\mathrm{sed}}^a-{\stackrel{\‾}{O}}_{\mathrm{sed}}^b$

$O_{\mathrm{ver}}^a\≈O_{\mathrm{ver}}^b$

Wherein,It is the water body initial dissolution oxygen value of mesh point a,It is the water body initial dissolution oxygen value of buoy point b place mesh point,It is the water body initial dissolution oxygen average of history same period of mesh point a,It it is the history same period of buoy point b place mesh point Water body initial dissolution oxygen average,It is the Oxygen Consumption By Sediments value of mesh point a,It it is the Oxygen Consumption By Sediments of buoy point b place mesh point Value,It is the Oxygen Consumption By Sediments average of history same period of mesh point a,It it is the end of history same period of buoy point b place mesh point Matter oxygen consumption average,It is the dissolved oxygen vertical proliferation amount of mesh point a,The dissolved oxygen being buoy point b place mesh point is vertical Straight diffusing capacity.

Dissolved oxygen water body oxygen consumption and horizontal proliferation effect in following for spring layer Dissolved Oxygen in Water equilibrium equation are expressed as history investigation bottom Dissolved oxygen average statistical and its deviation delta O, then dissolved oxygen equilibrium equation can be expressed as:

$O_{\mathrm{obs}}=O_{\mathrm{pre}}+O_{\mathrm{ver}}-O_{\mathrm{sed}}+{\stackrel{\‾}{O}}_{\mathrm{lat}}+\mathrm{\Δ}{O}_{\mathrm{lat}}-{\stackrel{\‾}{O}}_{\mathrm{wc}}-\mathrm{\Δ}{O}_{\mathrm{wc}}$

Wherein,It is that history investigates Dissolved Oxygen in Water horizontal proliferation statistics of variables average, Δ O_latWith horizontal proliferation amount average statistical Departure,It is that history investigates water body oxygen consumption average statistical, Δ O_wcIt is and the departure of water body oxygen consumption average statistical.

The buoy point b place mesh point then arranged in the range of arbitrary mess the point a, with working region in target operation regional extent Between dissolved oxygen relation can be rewritten as history investigation bottom dissolved oxygen average statistical and its deviation, such as following formula:

$O_{\mathrm{obs}}^a-O_{\mathrm{obs}}^b\≈{\stackrel{\‾}{O}}_{\mathrm{obs}}^a-{\stackrel{\‾}{O}}_{\mathrm{obs}}^b+\mathrm{\Δ}{O}_{\mathrm{lat}}^a-\mathrm{\Δ}{O}_{\mathrm{lat}}^b+\mathrm{\Δ}{O}_{\mathrm{wc}}^b-\mathrm{\Δ}{O}_{\mathrm{wc}}^a$

Wherein,It is the Dissolved Oxygen in Water level of mesh point a,It is the Dissolved Oxygen in Water level of buoy point b place mesh point, It is the average statistical of the Dissolved Oxygen in Water level of mesh point a,It it is the Dissolved Oxygen in Water level of buoy point b place mesh point Average statistical,It is the departure of mesh point a and horizontal proliferation amount average statistical,It it is buoy point b place mesh point With the departure of horizontal proliferation amount average statistical,It is inclined with water body oxygen consumption average statistical of buoy point b place mesh point Residual quantity,It it is the departure of mesh point a and water body oxygen consumption average statistical.

Consumption that the oxygen process of water column is mainly caused by Organic substance and the horizontal proliferation effect relevant to Dissolved Oxygen in Water position Magnitude of deviation is relevant.At the water body oxygen consumption variable quantity of other positions and buoy, water body oxygen consumption variable quantity is expressed as:

$\mathrm{\Δ}{O}_{\mathrm{wc}}^b-\mathrm{\Δ}{O}_{\mathrm{wc}}^a=k_b^{\mathrm{Chl}}\×C_b^{\mathrm{Chl}}-k_a^{\mathrm{Chl}}\×C_a^{\mathrm{Chl}}$

Wherein,K after the mesh point adjusting for depth of buoy point b place^chl,It is the kchl after mesh point a adjusting for depth,It is The top layer chlorophyll data of buoy point b place mesh point,It is the top layer chlorophyll data of mesh point a.

Thus, at the Dissolved Oxygen in Water level of other positions and buoy, dissolved oxygen levels has a statistic correlation:

$O_{\mathrm{est}}^a=\mathrm{\Δ}{O}_{\mathrm{obs}}^b+\mathrm{\Δ}{\stackrel{\‾}{O}}_{a-b}+k_b^{\mathrm{Chl}}\×C_b^{\mathrm{Chl}}-k_a^{\mathrm{Chl}}\×C_a^{\mathrm{Chl}}+\mathrm{\Δ}{O}_{\mathrm{lat}}^a-\mathrm{\Δ}{O}_{\mathrm{lat}}^b$

Wherein,It is the Dissolved Oxygen in Water guess value of mesh point a,It is that buoy point b place mesh point observed quantity is investigated with history The difference of dissolved oxygen average statistical,It it is the history investigation bottom dissolved oxygen system between mesh point a and buoy point b place mesh point The difference of meter average.

In conjunction with remote sensing top layer chlorophyll and the dissolved oxygen statistical result of each mesh point history same period, ignore in the formula in step 5 O_latWith Δ O_lat, the coefficient k of each mesh point is calculated with above-mentioned Formula Solution, and equal with the calculated coefficient k of each mesh point It is worth as the unified coefficient k used of each point in speculating.

A kind of coastal ocean hypoxia RS statistics model observed based on single-point buoy is we show at this.

7. with buoy point and the dissolved oxygen value of other mesh points of statistical correlation models coupling equivalence chlorophyll test value evaluation work region；

Integrating step 5 and step 6, by grid cycle calculations each point equivalence chlorophyll test value and water body oxygen consumption, and carry out dissolved oxygen water Flat spreading effect operates.

Owing to lacking the direct measurement of dissolved oxygen levels spreading effect, the Δ O between mesh point a and buoy point b place mesh point_lat Difference use 0.1 mg/L be unit each mesh point is carried out mathematics smooth so that the Dissolved Oxygen in Water gradient after Ping Hua meet this The dissolved oxygen gradient condition of mesh point history same period.

8. the dissolved oxygen value with grid each point calculates the area of hypoxia, calculates hypoxia degree with working region minimum dissolved oxygen value.

The meridian L total length about 40008km of the earth, in the area of dimension A, public at a length of often degree 111.13 of warp-wise In, dimension to distance l is:

l=L×cos(A)/360

The area (square kilometre) of grid cell is:

s=12350.614×x×y×cos(A)

X is warp-wise unit grid length (spending), and y is warp-wise unit grid length (spending).

Use buoy observed data, by hypoxia degree and the length of hypoxia persistent period, set up offshore sea waters and count in real time based on buoy According to hypoxia harm instruction parameter, this parameter is with 1 to 5 extent of injury expressing hypoxia, and 1 be the most slightly, and 5 is the most serious.

The syntagmatic of the hypoxia extent of injury and persistent period see table:

Table 1

The present invention utilizes single-point buoy observed data, in conjunction with remote sensing image data, gives area and the degree of working region hypoxia.

Below, what we were explained with reference to the present invention is embodied as form.

Fig. 1 is the frame diagram of the estimation method representing a kind of coastal ocean hypoxic condition relevant with embodiment of the present invention.Framework Figure includes human-computer interaction terminal 1, Ocean Color Remote Sensing satellite 2, ground satellite station 3, hypoxia oceanographic buoy 4, buoy information processing eventually End 5, communication network 6 and the hypoxia information processing terminal 7.User, after determining working region, utilizes at human-computer interaction terminal 1 Reason data, collect and process working region history survey data.Water colour data are transferred to ground satellite station by Ocean Color Remote Sensing satellite 2 3, ground satellite station 3 preserves chlorophyll mean data by agreement filename and file path.Hypoxia oceanographic buoy 4 passes through wireless network Real-time Monitoring Data is transferred to the buoy information processing terminal 5 by network.The hypoxia information processing terminal 7 is gone through by Network Capture process History survey data, Ocean Color Remote Sensing data and buoy real-time observed data.Ground satellite station, the buoy information processing terminal, man-machine friendship Terminal and the hypoxia information processing terminal are connected with each other by the Internet mutually.

Fig. 2 is the FB(flow block) of the measuring method of coastal ocean hypoxic condition.User obtains working region history by step 1 and adjusts Looking into data, data include marine topography, water body spring layer, bottom dissolved oxygen and top layer chlorophyll etc.；Man-machine friendship is utilized by step 2 Terminal network is formatted and is calculated water body spring layer and the average of bottom dissolved oxygen and minimum mutually, determines that working region dissolved oxygen is minimum The locus of value.Utilize hypoxia oceanographic buoy 4 and the buoy information processing terminal 5 by step 3 obtain buoy measure temperature, The time series datas such as salinity, dissolved oxygen, chlorophyll, are calculated the average of above-mentioned measurement data by agreement, preserve average daily by agreement Value Data.Read buoy by step 4 and measure the annual average data of temperature, calculate spring layer intensity, and with the spring layer of floating-point point Intensity infers the spring layer intensity of working region；Read buoy by agreement and measure the annual average data of bottom dissolved oxygen, it is judged that floating-point institute Whether hypoxia occurs at point；User obtains the chlorophyll data of the working region that ground satellite station processes by month by step 5 Calculate and preserve；The equivalent chlorophyll on each nexus is calculated by step 6；The hypoxia information processing terminal calculates non-buoy by step 5 The dissolved oxygen value of place mesh point, is calculated the dissolved oxygen value of whole grid by step 7, calculates dissolved oxygen scope and low by step 8 Oxygen degree.

As it is shown on figure 3, in an embodiment, working region scope is east longitude 122-124, north latitude 29.5-32.5 degree, empty at this In the range of between, by longitude 0.25 degree, 0.3 degree of grid of dimension turns to 50 × 100 mesh points.At longitude and latitude, buoy is set, floating Mark includes 3 subsurface buoys, and it is upper and lower that subsurface buoy lays respectively at spring layer.

In an embodiment, Dissolved Oxygen in Water horizontal proliferation effect processes and uses 0.1mg/L is that unit carries out mathematics to each mesh point Smooth so that the Dissolved Oxygen in Water gradient after Ping Hua meets the dissolved oxygen gradient condition of this mesh point history same period: minimum dissolved oxygen The dissolved oxygen gradient condition in region, onshore, position for often to spend less than 4.0mg/L longitudinal, minimum dissolved oxygen position to The dissolved oxygen gradient condition in land region territory is for often to spend less than 8.0mg/L longitudinal.

The invention is not restricted to above embodiment, in the invention scope recorded in detail in the claims, all changes can be carried out More, these changes are also contained in the scope of the present invention certainly, and this is self-evident.

Claims (5)

1. the estimation method of a coastal ocean hypoxic condition, it is characterised in that comprise the steps:

(1) working region marine topography, water body spring layer, bottom dissolved oxygen and remote sensing top layer chlorophyll history survey data for many years are obtained；

(2) working region is carried out stress and strain model, by month and grid to historical data resampling, average statistical and minima；

(3) obtain and be positioned at temperature that the upper and lower buoy group of coastal ocean spring layer measures in real time, salinity, dissolved oxygen, chlorophyll time series Data；

(4) with the degree of depth and the temperature measured in real time, the salinity annual average calculating spring layer intensity of buoy；

(5) obtain the chlorophyll remote sensing average view data with the buoy time series same period in the range of working region, calculate green with top layer leaf The element organic coefficient of oxygen consuming of revised bottom of equal value；

(6) the statistical correlation model of any two points in grid is set up with Dissolved Oxygen in Water balance；

(7) with buoy point and the dissolved oxygen value of other mesh points of statistical correlation models coupling equivalence chlorophyll test value evaluation work region；

(8) calculate hypoxia scope with the dissolved oxygen value of grid each point, calculate hypoxia degree with working region minimum dissolved oxygen value.

The estimation method of a kind of coastal ocean hypoxic condition the most according to claim 1, it is characterised in that described step (1) In, described water body spring layer includes that hypoxia occurs the spring layer intensity during month and depth of pycnocline.

The estimation method of a kind of coastal ocean hypoxic condition the most according to claim 1, it is characterised in that described step (2) Particularly as follows: target operation regional extent is carried out spatial gridding subdivision by longitude and latitude direction；During hypoxia occurs month, press Month and mesh point carry out resampling to the depth of water, spring layer intensity, depth of pycnocline, bottom dissolved oxygen, top layer chlorophyll, calculate net The average of each mesh point above-mentioned parameter and minima in the range of lattice；Calculate the mean depth at dissolved oxygen minima place, it is determined as ginseng Examine the depth of water；Calculate the bottom dissolved oxygen spring layer strength mean value less than 3.0mg/L position, it is determined as with reference to spring layer intensity, obtain work Region bottom dissolved oxygen is less than hypoxia month and the persistent period of 3.0mg/L.

The estimation method of a kind of coastal ocean hypoxic condition the most according to claim 1, it is characterised in that described step (5) In, described revised bottom organic coefficient of oxygen consuming of equal value with top layer chlorophyll is obtained by following formula:

k^chl=h × | d-d₀|/d_o×(r_og k_g)

Wherein, k^chlIt is revised bottom organic coefficient of oxygen consuming of equal value with top layer chlorophyll, d_oFor with reference to the depth of water, d is for waiting to count Calculating the mesh point depth of water, h is the top layer chlorophyll estimation organic coefficient of bottom, r_ogFor the carbon-to-oxygen ratio of Organic substance oxygen consumption, k_gFor having Phosphorus content coefficient in machine thing.

The estimation method of a kind of coastal ocean hypoxic condition the most according to claim 1, it is characterised in that described step (6) In the grid set up, the statistical correlation model of any two points is:

$O_{est}^a={\mathrm{\ΔO}}_{obs}^b+\mathrm{\Δ}{\stackrel{\‾}{O}}_{a-b}+k_b^{Chl}\×C_b^{Chl}-k_a^{Chl}\×C_a^{Chl}+{\mathrm{\ΔO}}_{lat}^a-{\mathrm{\ΔO}}_{lat}^b$

Wherein,It is the Dissolved Oxygen in Water guess value of mesh point a,It is that buoy point b place mesh point observed quantity is investigated with history The difference of dissolved oxygen average statistical,It it is the history investigation bottom dissolved oxygen system between mesh point a and buoy point b place mesh point The difference of meter average,K after the mesh point adjusting for depth of buoy point b place^chl,It is the k after mesh point a adjusting for depth^chl,It is the top layer chlorophyll data of buoy point b place mesh point,It is the top layer chlorophyll data of mesh point a,It is Mesh point a and the departure of horizontal proliferation amount average statistical,It is that buoy point b place mesh point is equal with horizontal proliferation amount statistics The departure of value.