Comparison of ERA5 turbulent fluxes at the air-sea interface with measurements from a wave-following platform

Turbulent fluxes at the air-sea interface were estimated with data collected in 2011-2020 with a low-profile platform named OCARINA during eight experiments in five regions:2011,2015,and 2016 in the Iroise Sea;2012 in the tropical Atlantic;2014 in the Chilie-Peru upwelling;2017 and 2018 in the Mediterranean Sea,and 2018 and 2020 in Barbados.The observations were carried out with moderate winds(2-10 m s^(-1))and average wave heights of 1.5 m.In this study,the authors used the fluxes calculated by the bulk method using OCARINA-sampled data as the input.These data can validate the fluxes estimated from ERAS reanalysis data.The OCARINA and ERA5 data were taken concomitantly.To do this,the authors established an algorithm to extract the OCARINA data as closely as possible to the reanalysis data in time and position.The measurements of the OCARINA platform can conclude on the relevance of the widely used reanalysis data.

Seasonal response of nitrogen exchange fluxes to crab disturbance at sediment-water interface in coastal tidal wetlands

Coastal wetlands are hotspots for nitrogen(N)cycling,and crab burrowing is known to transform N in intertidal marsh soils.However,the underlying mechanisms remain unclear.This study conducted field experiments and used indoor control test devices to investigate the seasonal response of nitrogen to crab disturbance at the sediment-water interface in coastal tidal flat wetlands.The results showed that crab disturbance exhibited significant seasonality with large seasonal differences in cave density and depth.Due to crab disturbance,nitrogen fuxes at the sediment-water interface were much greater in the box with crabs than in the box without crabs.In summer,NH-N showed a positive flux from the sediment to the overlying water,but NO2-N and NOg-N showed positive fluxes from the sediment to the overlying water only in early stages.In winter,NH-N showed a positive flux from the sediment to the overlying water,but NO-N and NO,-N both exhibited positive and negative fluxes.These results indicated that the presence of crab burrows can cause the aerobic layer to move downward by approximately 8-15 cm in summer and directly promote nitrification at the sediment surface.

Impact of Initial Soil Conditions on Soil Hydrothermal and Surface Energy Fluxes in the Permafrost Region of the Tibetan Plateau

Accurate initial soil conditions play a crucial role in simulating soil hydrothermal and surface energy fluxes in land surface process modeling.This study emphasized the influence of the initial soil temperature(ST)and soil moisture(SM)conditions on a land surface energy and water simulation in the permafrost region in the Tibetan Plateau(TP)using the Community Land Model version 5.0(CLM5.0).The results indicate that the default initial schemes for ST and SM in CLM5.0 were simplistic,and inaccurately represented the soil characteristics of permafrost in the TP which led to underestimating ST during the freezing period while overestimating ST and underestimating SLW during the thawing period at the XDT site.Applying the long-term spin-up method to obtain initial soil conditions has only led to limited improvement in simulating soil hydrothermal and surface energy fluxes.The modified initial soil schemes proposed in this study comprehensively incorporate the characteristics of permafrost,which coexists with soil liquid water(SLW),and soil ice(SI)when the ST is below freezing temperature,effectively enhancing the accuracy of the simulated soil hydrothermal and surface energy fluxes.Consequently,the modified initial soil schemes greatly improved upon the results achieved through the long-term spin-up method.Three modified initial soil schemes experiments resulted in a 64%,88%,and 77%reduction in the average mean bias error(MBE)of ST,and a 13%,21%,and 19%reduction in the average root-mean-square error(RMSE)of SLW compared to the default simulation results.Also,the average MBE of net radiation was reduced by 7%,22%,and 21%.

A data assimilation-based forecast model of outer radiation belt electron fluxes

Because radiation belt electrons can pose a potential threat to the safety of satellites orbiting in space,it is of great importance to develop a reliable model that can predict the highly dynamic variations in outer radiation belt electron fluxes.In the present study,we develop a forecast model of radiation belt electron fluxes based on the data assimilation method,in terms of Van Allen Probe measurements combined with three-dimensional radiation belt numerical simulations.Our forecast model can cover the entire outer radiation belt with a high temporal resolution(1 hour)and a spatial resolution of 0.25 L over a wide range of both electron energy(0.1-5.0 MeV)and pitch angle(5°-90°).On the basis of this model,we forecast hourly electron fluxes for the next 1,2,and 3 days during an intense geomagnetic storm and evaluate the corresponding prediction performance.Our model can reasonably predict the stormtime evolution of radiation belt electrons with high prediction efficiency(up to~0.8-1).The best prediction performance is found for~0.3-3 MeV electrons at L=~3.25-4.5,which extends to higher L and lower energies with increasing pitch angle.Our results demonstrate that the forecast model developed can be a powerful tool to predict the spatiotemporal changes in outer radiation belt electron fluxes,and the model has both scientific significance and practical implications.

Sea spray induced air-sea heat and salt fluxes based on the wavesteepnessdependent sea spray model

Sea spray,which comprises amounts of small ocean droplets,plays a significant role in the air-sea coupling,atmospheric and oceanic dynamics,and climate.However,it remains arduous to arrive at estimates for the efficiency and accuracy of the sea spray induced air-sea heat and salt fluxes.This is because the microphysical process of sea spray evolution in the air is of extreme complexity.In this study,we iteratively calculated the sea spray induced air-sea heat and salt fluxes at various weather condition.To do so,we implemented one novel wave-steepness-dependent sea spray model into a bulk air-sea fluxes algorithm and utilized other sea spray models as comparisons.Based on the improved wave-dependent bulk turbulent algorithm,we observed that despite the negative contribution of sea spray to the sensible heat fluxes,the sea spray positively contributes to the air-sea latent heat fluxes,leading to an overall increase in the total air-sea heat fluxes.The additional heat fluxes caused by sea spray may be the missing critical process that can clarify the discrepancies observed between measured and modelled Tropical Cyclone’s development and intensification.In addition to heat fluxes,we observed that sea spray has significant impacts on the air-sea salt fluxes.As the sea salt particles are one of the main sources of the atmosphere aerosol,our results imply that sea spray could impact global and regional climate.Thus,given the significance of sea spray on the air-sea boundary layer,sea spray effects need to be considered in studies of air-sea interaction,dynamics of atmosphere and ocean.

Evaluating Parameterizations for Turbulent Fluxes over the Landfast Sea-Ice Surface in Prydz Bay, Antarctica

It is crucial to appropriately determine turbulent fluxes in numerical models.Using data collected in East Antarctica from 8 April to 26 November 2016,this study evaluates parameterization schemes for turbulent fluxes over the landfast seaice surface in five numerical models.The Community Noah Land Surface Model with Multi-Parameterizations Options(Noah_mp)best replicates the turbulent momentum flux,while the Beijing Climate System Model(BCC_CSM)produces the optimum sensible and latent heat fluxes.In particular,two critical issues of parameterization schemes,stability functions and roughness lengths,are investigated.Sensitivity tests indicate that roughness lengths play a decisive role in model performance.Based on the observed turbulent fluxes,roughness lengths over the landfast sea-ice surface are calculated.The results,which can provide a basis for setting up model parameters,reveal that the dynamic roughness length(z0m)increases with the increase of frictional velocity(u*)when u*≤0.4 m s^(−1) and fluctuates around 10^(−3 )m when u*>0.4 m s^(−1);thermal roughness length(z0t)is linearly related to the temperature gradient between air and sea-ice surface(ΔT)with a relation of lg(z0t)=−0.29ΔT−3.86;and the mean water vapor roughness length(z0q)in the specific humidity gradient(Δq)range ofΔq≤−0.6 g kg^(−1) is 10^(−6) m,3.5 times smaller than that in the range ofΔq˃−0.6 g kg^(−1).

Cross-satellite calibration of high-energy electron fluxes measured by FengYun-4A based on Arase observations

We use the High-energy Electron Experiments(HEP)instrument onboard Arase(ERG)to conduct an energy-dependent cross-satellite calibration of electron fluxes measured by the High Energy Particle Detector(HEPD)onboard FengYun-4A(FY-4A)spanning from April 1,2017,to September 30,2019.By tracing the two-dimensional magnetic positions(L,magnetic local time[MLT])of FY-4A at each time,we compare the datasets of the conjugate electron fluxes over the range of 245–894 keV in 6 energy channels for the satellite pair within different sets of L×MLT.The variations in the electron fluxes observed by FY-4A generally agree with the Arase measurements,and the percentages of the ratios of electron flux conjunctions within a factor of 2 are larger than 50%.Compared with Arase,FY-4A systematically overestimates electron fluxes at all 6 energy channels,with the corresponding calibration factors ranging from 0.67 to 0.81.After the cross-satellite calibration,the electron flux conjunctions between FY-4A and Arase show better agreement,with much smaller normalized root mean square errors.Our results provide a valuable reference for the application of FY-4A high-energy electron datasets to in-depth investigations of the Earth’s radiation belt electron dynamics.

Regularity of Fluxes in Nonlinear Hyperbolic Balance Laws

This paper addresses the issue of the formulation of weak solutions to systems of nonlinear hyperbolic conservation laws as integral balance laws.The basic idea is that the“meaningful objects”are the fluxes,evaluated across domain boundaries over time intervals.The fundamental result in this treatment is the regularity of the flux trace in the multi-dimensional setting.It implies that a weak solution indeed satisfies the balance law.In fact,it is shown that the flux is Lipschitz continuous with respect to suitable perturbations of the boundary.It should be emphasized that the weak solutions considered here need not be entropy solutions.Furthermore,the assumption imposed on the flux f(u)is quite minimal-just that it is locally bounded.

基于FLUXNET站点数据的中纬度森林对地表通量的影响

利用最新FLUXNET数据集,采用森林站和非森林站(开阔地)的对比分析用以表征毁林的作用,探讨中纬度地区毁林对地表通量和波文比的季节和日内影响。结果表明,毁林使得地表通量和波文比降低,其中以感热通量的减少为主。通过对比不同森林类型毁林发现,在日尺度上,地表通量变化在正午达到最强值,混合林地表通量变化相比于常绿针叶林较强,而弱于落叶阔叶林。在季节尺度上,感热通量在常绿针叶林和混合林中减少,潜热通量则在常绿针叶林和落叶阔叶林中减少。通过对不同气候区下地表通量变化发现,毁林使得常湿冷温气候地表通量以潜热通量减少为主,而常湿温暖气候以感热通量减少为主。夏干温暖气候下冬季感热通量减少,潜热通量增加,夏季则感热通量增加,潜热通量减少。整体而言,不同森林类型毁林影响地表通量变化的符号存在差异,而背景气候也决定了地表通量变化是以感热还是潜热通量变化为主。研究揭示了中纬度毁林对地表通量的影响,可以为中纬度地区植被恢复/造林政策的制定和实施提供观测依据,也可以为模式检验提供一定参考价值。

Diffusion and reaction mechanism of limestone and quartz in fluxed iron ore pellet roasting process

The increase to the proportion of fluxed pellets in the blast furnace burden is a useful way to reduce the carbon emissions in the ironmaking process.In this study,the interaction between calcium carbonate and iron ore powder and the mineralization mechanism of fluxed iron ore pellet in the roasting process were investigated through diffusion couple experiments.Scanning electron microscopy with energy dispersive spectroscopy was used to study the elements’diffusion and phase transformation during the roasting process.The results indicated that limestone decomposed into calcium oxide,and magnetite was oxidized to hematite at the early stage of preheating.With the increase in roasting temperature,the diffusion rate of Fe and Ca was obviously accelerated,while the diffusion rate of Si was relatively slow.The order of magnitude of interdiffusion coefficient of Fe_(2)O_(3)-CaO diffusion couple was 10^(−10) m^(2)·s^(−1) at a roasting temperature of 1200℃for 9 h.Ca_(2)Fe_(2)O_(5) was the initial product in the Fe_(2)O_(3)-CaO-SiO_(2) diffusion interface,and then Ca_(2)Fe_(2)O_(5) continued to react with Fe_(2)O_(3) to form CaFe_(2)O_(4).With the expansion of the diffusion region,the sillico-ferrite of calcium liquid phase was produced due to the melting of SiO_(2) into CaFe_(2)O_(4),which can strengthen the consolidation of fluxed pellets.Furthermore,andradite would be formed around a small part of quartz particles,which is also conducive to the consolidation of fluxed pellets.In addition,the principle diagram of limestone and quartz diffusion reaction in the process of fluxed pellet roasting was discussed.

基于Ameriflux通量观测数据的Hi-GLASS潜热通量产品验证

潜热通量产品的验证与分析对于研究气候变化及能量循环具有重要意义。全球陆表高分辨率蒸散产品(high resolution global lAnd surface evapotranspiration product,Hi-GLASS ET)融合了5种传统蒸散算法,能够生产出较高精度的陆表潜热通量产品,但目前没有针对此产品的验证研究。利用Ameriflux通量观测站点的潜热通量观测值与相应的Hi-GLASS陆表潜热通量产品估算值进行对比,获取多组有效验证数据。验证结果显示,所选站点实际观测值与产品估算值的决定系数(R 2)为0.6,均方根误差(RMSE)为34.4 W/m^(2),平均偏差(Bias)为-13.4 W/m^(2),克林-古普塔效率(Kling-Gupta efficiency,KGE)为0.49,Hi-GLASS潜热通量产品具有较高的精度,算法的拟合结果较好;此外,空间分布也表明Hi-GLASS陆表潜热通量产品符合正常的自然规律。由于数据获取的局限性,仅采用了美国地区18个站点数据对产品进行验证,在其他地区仍需进一步验证。

FPGA-based position reconstruction method for neutron beam flux spatial distribution measurement in BNCT

A new measurement method for the spatial distribution of neutron beam flux in boron neutron capture therapy(BNCT)is being developed based on the two-dimensional Micromegas detector.To address the issue of long processing times in traditional offline position reconstruction methods,this paper proposes a field programmable gate array based online position reconstruction method utilizing the micro-time projection chamber principle.This method encapsulates key technical aspects:a self-adaptive serial link technique built upon the dynamical adjustment of the delay chain length,fast sorting,a coordinate-matching technique based on the mapping between signal timestamps and random access memory(RAM)addresses,and a precise start point-merging technique utilizing a circular combined RAM.The performance test of the selfadaptive serial link shows that the bit error rate of the link is better than 10-12 at a confidence level of 99%,ensuring reliable data transmission.The experiment utilizing the readout electronics and Micromegas detector shows a spatial resolution of approximately 1.4 mm,surpassing the current method’s resolution level of 5 mm.The beam experiment confirms that the readout electronics system can obtain the flux spatial distribution of neutron beams online,thus validating the feasibility of the position reconstruction method.The online position reconstruction method avoids traditional methods,such as bubble sorting and traversal searching,simplifies the design of the logic firmware,and reduces the time complexity from O(n2)to O(n).This study contributes to the advancement in measuring neutron beam flux for BNCT.

Hydrographic Characteristics and Oceanic Heat Flux in the Upper Arctic Ocean over the Alpha Ridge Observed by the DTOP Platform in 2018 and 2021

In 2018 and 2021,the Drift-Towing Ocean Profilers(DTOP)provided extensive temperature and salinity data on the upper 120m ocean through their drifts over the Alpha Ridge north of the Canada Basin.The thickness and temperature maximum of Alaska Coastal Water(ACW)ranged from 20m to 40m and-1.5℃to-0.8℃,respectively,and the salinity generally maintained from 30.2 to 32.5.Comparison with World Ocean Atlas 2018’s climatology manifested a 40m-thick and warm ACW roughly ex-ceeding the temperature maximum by 0.4–0.5℃in June–August 2021.This anomalously warm ACW was highly related to the ex-pansion of the Beaufort Gyre in the negative Arctic Oscillation phase.During summer,the under-ice oceanic heat flux F_(w)^(OHF)was elevated,with a maximum value of above 25Wm^(-2).F_(w)^(OHF)was typically low in the freezing season,with an average value of 1.2Wm^(-2).The estimates of upward heat flux contributed by ACW to the sea ice bottom F_(w)^(OHF)were in the range of 3–4Wm^(-2)in June–August 2021,when ACW contained a heat content of more than 80MJm^(-2).The heat loss over this period was driven by a weak stratification upon the ACW layer associated with a surface mixed layer(SML)approaching the ACW core.After autumn,F_(w)^(OHF)was reduced(<2 Wm^(-2))except during rare events when it elevated F_(w)^(OHF)slightly.In addition,the intensive and widespread Ekman suction,which created a violent upwelling north of the Canada Basin,was largely responsible for the substantial cooling and thinning of the ACW layer in the summer of 2021.

High-chromium vanadium-titanium magnetite all-pellet integrated burden optimization and softening-melting behavior based on flux pellets

High-chromium vanadium-titanium magnetite(HVTM)is a crucial polymetallic-associated resource to be developed.The allpellet operation is a blast furnace trend that aims to reduce carbon dioxide emissions in the future.By referencing the production data of vanadium-titanium magnetite blast furnaces,this study explored the softening-melting behavior of high-chromium vanadium-titanium magnetite and obtained the optimal integrated burden based on flux pellets.The results show that the burden with a composition of 70wt%flux pellets and 30wt%acid pellets exhibits the best softening-melting properties.In comparison to that of the single burden,the softening-melting characteristic temperature of this burden composition was higher.The melting interval first increased from 307 to 362℃and then decreased to 282℃.The maximum pressure drop(ΔPmax)decreased from 26.76 to 19.01 kPa.The permeability index(S)dropped from 4643.5 to 2446.8 kPa·℃.The softening-melting properties of the integrated burden were apparently improved.The acid pellets played a role in withstanding load during the softening process.The flux pellets in the integrated burden exhibited a higher slag melting point,which increased the melting temperature during the melting process.The slag homogeneity and the TiC produced by over-reduction led to the gas permeability deterioration of the single burden.The segregation of the flux and acid pellets in the HVTM proportion and basicity mainly led to the better softening-melting properties of the integrated burden.

THE RIEMANN PROBLEM FOR ISENTROPIC COMPRESSIBLE EULER EQUATIONS WITH DISCONTINUOUS FLUX

We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.

Electromagnetic Performance Analysis of Variable Flux Memory Machines with Series-magnetic-circuit and Different Rotor Topologies

In this paper,the electromagnetic performance of variable flux memory(VFM)machines with series-magnetic-circuit is investigated and compared for different rotor topologies.Based on a V-type VFM machine,five topologies with different interior permanent magnet(IPM)arrangements are evolved and optimized under same constrains.Based on two-dimensional(2-D)finite element(FE)method,their electromagnetic performance at magnetization and demagnetization states is evaluated.It reveals that the iron bridge and rotor lamination region between constant PM(CPM)and variable PM(VPM)play an important role in torque density and flux regulation(FR)capabilities.Besides,the global efficiency can be improved in VFM machines by adjusting magnetization state(MS)under different operating conditions.

Optimum Profiles of Endwall Contouring for Enhanced Net Heat Flux Reduction and Aerodynamic Performance

Successfully utilized non-axisymmetric endwalls to enhance turbine efficiencies(aerodynamic and turbine inlet temperatures)by controlling the characteristics of the secondary flow in a blade passage.This is accomplished by steady-state numerical hydrodynamics and deep knowledge of the field of flow.Because of the interaction between mainstream and purge flow contributing supplementary losses in the stage,non-axisymmetric endwalls are highly susceptible to the inception of purge flow exit compared to the flat and any advantage rapidly vanishes.The conclusions reveal that the supreme endwall pattern could yield a lowering of the gross pressure loss at the design stage and is related to the size of the top-loss location being productively lowered.This has led to diminished global thermal exchange lowered in the passage of the vane alone.The reverse flow adjacent to the suction side corner of the endwall is migrated farther from the vane surface,as the deviated pressure spread on the endwall accelerates the flow and progresses the reverse flow core still downstream.The depleted association between the tornado-like vortex and the corner vortex adjacent to the suction side corner of the endwall is the dominant mechanism of control in the contoured end wall.In this publication,we show that the non-axisymmetric endwall contouring by selective numerical shape change method at most prominent locations is advantageous in lowering the thermal load in turbines to augment the net heat flux reduction as well as the aerodynamic performance using multi-objective optimization.

Milk fat globule epidermal growth factor 8 alleviates liver injury in severe acute pancreatitis by restoring autophagy flux and inhibiting ferroptosis in hepatocytes

BACKGROUND Liver injury is common in severe acute pancreatitis(SAP).Excessive autophagy often leads to an imbalance of homeostasis in hepatocytes,which induces lipid peroxidation and mitochondrial iron deposition and ultimately leads to ferroptosis.Our previous study found that milk fat globule epidermal growth factor 8(MFG-E8)alleviates acinar cell damage during SAP via binding toαvβ3/5 integrins.MFG-E8 also seems to mitigate pancreatic fibrosis via inhibiting chaperone-mediated autophagy.AIM To speculate whether MFG-E8 could also alleviate SAP induced liver injury by restoring the abnormal autophagy flux.METHODS SAP was induced in mice by 2 hly intraperitoneal injections of 4.0 g/kg L-arginine or 7 hly injections of 50μg/kg cerulein plus lipopolysaccharide.mfge8-knockout mice were used to study the effect of MFG-E8 deficiency on SAPinduced liver injury.Cilengitide,a specificαvβ3/5 integrin inhibitor,was used to investigate the possible mechanism of MFG-E8.RESULTS The results showed that MFG-E8 deficiency aggravated SAP-induced liver injury in mice,enhanced autophagy flux in hepatocyte,and worsened the degree of ferroptosis.Exogenous MFG-E8 reduced SAP-induced liver injury in a dose-dependent manner.Mechanistically,MFG-E8 mitigated excessive autophagy and inhibited ferroptosis in liver cells.Cilengitide abolished MFG-E8’s beneficial effects in SAP-induced liver injury.CONCLUSION MFG-E8 acts as an endogenous protective mediator in SAP-induced liver injury.MFG-E8 alleviates the excessive autophagy and inhibits ferroptosis in hepatocytes by binding to integrinαVβ3/5.

Impact of Cattaneo-Christov Heat Flux in the Nanofluid Flow over an Inclined Permeable Surface with Irreversibility Analysis

This study discusses the magnetohydrodynamic nanofluid flow over an inclined permeable surface influenced by mixed convection, and Cattaeo-Christov heat flux. The heat transfer analysis is performed in the presence of a heat source/sink and thermal stratification. To gauge the energy loss during the process, an irreversibility analysis is also performed. A numerical solution to the envisaged problem is obtained using the bvp4c package of MATLAB. Graphs are drawn to assess the consequences of the arising parameters against the associated profiles. The results show that an augmentation in the magnetic field and nanomaterial volume fraction results in an enhancement in the temperature profile. A strong magnetic field can significantly reduce the fluid velocity. The behavior of the Skin friction coefficient against the different estimates of emerging parameters is discussed. .

Effects of hydrodynamics processes on phosphorus fluxes from sediment in large,shallow Taihu Lake

The turnover of phosphorus （P） in lake sediments, a major cause of eutrophication and subsequent deterioration of water quality, is in need of deep understanding. In this study, effects of resuspension on P release were studied in cylindrical microcosms with Yshape apparatus. The results indicated that there was a positive correlation between flux of suspended substance across sediment-water interface （Fss） and the wind speed, and an increasing Fss during each wind process followed by a steady state. The maximal Fss under fight, moderate, and strong wind conditions were 299.9±41.1,573.4±61.7, and 2093.8±215.7 g/m^2, respectively. However, flux of P across sediment-water interface （Fp） did not follow a similar pattern as Fss responding to wind intensity, which increased and reached the maximum in initial 120 rain for fight wind, then decreased gradually, with maximal flux of 9.4±1.9 mg/m^2. A rapid increase of Fp at the first 30 rain was observed under moderate wind, with maximal flux of 11.2±0.6 mg/m^2. Surprisingly, strong wind caused less Fp than under light and moderate wind conditions with maximal flux of 3.5±0.9 mg/m^2. Fss in water column declined obviously during the sedimentation process after winds, but Fp varied with wind regime. No obvious difference was detected on Fp after 8 h sedimentation process, compared with the initial value, which means little redundant P left in the water column after winds.