Estimates of potential new production in the Java-Sumatra upwelling system

The Java-Sumatra upwelling is one of the most important upwelling systems in the Indian Ocean, with maximum upwelling intensity in July through August. To estimate the nitrate supplied by upwelling, we developed a three-dimensional hydrodynamic model to calculate the mean vertical speed and determine the depth of upwelling. We used in-situ vertical nitrate profiles to assess nitrate concentration in the upwelled waters, and calculated the nitrate supply as the product of nitrate concentration and vertical transport obtained from the numerical model. The calculated result represents potential new production generated in the upwelling region. We found that on the event time scale (monthly) of Java-Sumatra upwelling, water brought to the surface originated from locations 100-m deep, giving a nitrate supply of 93.77×10 3mol/s and potential new production of 1.02×10 14gC/a.

Satellite observations of coastal upwelling in the northern Arafura Sea

Coastal upwelling is significant for marine ecosystems by lifting nutrient-rich deep waters into the euphotic zone,thereby increasing primary and secondary productivity.The satellite observations show that the northern Arafura Sea(NAS),especially in the coastal region,features high chlorophyll-a(chl-a)concentrations,implying a strong coastal upwelling.However,coastal upwelling in the NAS has not received much attention.Based on a semi-automatic image processing technology,the seasonal and interannual variability of coastal upwelling in the NAS are investigated in this study using satellite-observed sea surface temperature(SST)and wind data.The results suggest that there are seasonal coastal upwelling events in the NAS modulated by upwelling-favorable southeast monsoon(SEM).The annual mean days,mean area,and annual mean intensity of coastal upwelling events during the SEM season are 92 days,6514 km^(2),and-5.31×10^(5),respectively,while the corresponding values during the northwest monsoon(NWM)season are 32 days,5569 km^(2),and-1.41×10^(5).It is also found that the SEM coastal upwelling in the NAS displays prominent interannual variability.The strong upwelling events are found in 2010,2013,2016,and 2017 when the southeast monsoon winds were weaker.Further analysis suggests that at the interannual scale,the upwelling index(UI)averaged in the SEM season is negatively correlated with that of three upwelling indicators.This can be attributed to the limitation of onshore geostrophic flow which is evidenced by the negative correlation between the UI and the alongshore difference in sea surface height.This study highlights the important role of the southeast monsoon in the temporal variability of coastal upwelling in the NAS.

Extrapolation of upwelling irradiance just beneath the ocean surface from the in-water radiometric profile measurements

Accurate measurements of upwelling irradiance just beneath the ocean surface,E_(u)(λ,0^(-)),can be used to calculate ocean optical parameters,and further develop retrieval algorithms for remotely sensing water component concentrations.Due to the effects of sea surface waves,perturbation from instrument platform(ship),and instrument self-shading,E_(u)(λ,0^(-))is often difficult to be accurately measured.This study presents a procedure for extrapolating the E_(u)(λ,0^(-))from the in-water radiometric profile measurements.Using the optical profile data from 13 bands(ranging from 381 to 779 nm)measured by 45 casts in the Ligurian Sea during 2003–2009,the E_(u)(λ,0^(-))was extrapolated from in-water upwelling irradiance measurements between the initial shallow depth,Z_(0),and an optimal bottom depth,Z_(1),by three linear models(linear,2-degree polynomial,and exponential)and two nonlinear models(LOESS and spline).The accumulated errors of extrapolated E_(u)(λ,0^(-))at each wavelength for the five models were calculated.It was found that the optimal Z_(1) depth for the linear and exponential models was at the depth of80%of E_(u)(λ,Z_(0)),50%of E_(u)(λ,Z_(0))for the 2-degree polynomial model,40%of E_(u)(λ,Z_(0))for the LOESS model,and 15%of E_(u)(λ,Z_(0))for the spline model.The extrapolated E_(u)(λ,0^(-))derived from the five models was in good agreement with the calculated true E_(u)(λ,0^(-)).In all bands,the 2-degree polynomial model achieved the highest accuracy,followed by the LOESS model.In the short band of 381–559 nm,the linear and exponential models had the third-best performance,and the spline model performed worst within this range.For the red band of 619–779 nm,the accuracies of the exponential and spline models had the third highest performance,and the linear model produced lowest accuracy.Hence,the 2-degree polynomial model was an optimal procedure for extrapolation of E_(u)(λ,0^(-))from the in-water radiometric profile measurements.

Variation of the coastal upwelling off South Java and their impact on local fishery resources

There is a vast upwelling area induced by the southeast monsoon in the waters off South Java,making the region an important fishing ground.Climate events can affect the variation of upwelling,but oceanographers have different understandings on the extent to which climate events control upwelling in this area,which leads to a lack of basis for studies on the evaluation and mechanisms of the variability of fishery resources in the region.The correlation between environmental parameters,including surface temperature(SST),chlorophyll-a(Chl-a)concentration,and climate event indices in South Java from 2003 to 2020 was analyzed.Results show that the Indian Ocean Dipole(IOD)has a greater influence on the interannual variability of upwelling intensity than ENSO.During the IOD,variations in equatorial latitudinal winds excite different types of Kelvin waves that anomalously deepen or shallow the thermocline,which is the main cause of anomalous variations in upwelling,independent of variations in the local wind field.A correlation between the interannual variability in upwelling and the annual catches was revealed,showing that climatic events indirectly affect fishery resources through upwelling effects.During positive IOD/El Niño periods,strong upwelling delivers more nutrients to the surface layer,which favors fish growth and reproduction,resulting in higher annual catches.A negative IOD/La Niña,on the other hand,leads to weaker upwelling and fewer nutrients into the surface waters.Fish tend to move in deeper waters,making traditional fishing methods less efficient and consequently lower annual catches.

Studies on the role of marine bacteria in the upwelling ecosystem——Ⅰ.Bacterial biomass in the upwelling regions

-Six cruises were carried out off the south bank of Fujian - Taiwan during the period of December 21, 1987 to November 15, 1988 to estimate the contribution of bacterial biomass carbon (BBC) to the totai particulate organic pools using epifluorescent microscopic technique. The results show that the standing crop of bacteria fluctuated from 0. 95 to 66. 60 mg /m3 (dry weight). Upwelling phenomena appeared in the region around Nanpeng Island in summer while in the region of Waixie in all seasons. The average value of BBC was 27. 60(±6. 08)mg/m3and 21. 32 (±2. 34) mg/m3 respectively. The seasonal and spatial distribution is discussed in relation to environmental factors as well as upwelling phenomena. The role of bacteria in the flow of material and energy in the upwelling ecosystem is emphasized.

Receiver function study of the crustal structure of Northeast China:Seismic evidence for a mantle upwelling beneath the eastern flank of the Songliao Basin and the Changbaishan region

Northeast China is a unique place to study intra-plate volcanism. We analyzed P-wave receiver function data recorded by 111 permanent broadband seismic stations in Northeast China. The results show that the crustal thickness varies from 27.9 km beneath the eastern flank of the Songliao Basin to 40.7 km beneath the Great Xing＇an Range region. The large depth variations of the Moho can be largely but not completely explained by surface topography. The residual Moho depth calculated based on the Airy＇s isostasy model indicates that the Moho is dynamically uplifted by 〈3 km beneath the eastern flank of the Songliao Basin and the Changbaishan region. We suggest that a mantle upwelling, which has been proposed by several recent seismic studies, might have caused the uplift.

UPWELLING AND SEDIMENTATION DYNAMICS II. A SIMPLE MODEL

It was found that the location and distributional pattern of certain sea bottom mud patches in the Yellow Sea and East China Sea are identical with the geographical location and flow pattern of the upwelling above. Following a general description by Hu, (1984) of the role of upwelling in . sedimentation, this study presents a two - dimensional numerical model to further examine, the sedimentation dynamics in upwelling areas. The model shows that advection by ocean circulation, downward settling due to gravity and turbulent mixing are the main mechanisms controlling the "movement of suspended matters, and that the upwelling :on the shelf favors deposition near the bottom. It is concluded that on the shelf (1) wherever there is upwelling, mud sediment must exist; and . (2) the closer to the center of upwelling, the greater the concentration of suspended matters in the water, the greater the depositing rate near the bottom and the finer the sediment on the sea bed.

Evolution of a coastal upwelling event during summer 2004 in the southern Taiwan Strait

A coastal upwelling event in the southern Taiwan Strait （STWS） was investigated using intensive cruise surveys （four repeated transects in a month） and satellite data in July and early August 2004. The extensive upwelling-associated surface cold water was first observed in early July （-2.0×10a km2） along the STWS coast. Then, the cold surface water reduced in size by -50% with decreased chlorophyll concentrations after 15 days, indicating the weakening of the upwelling event. At the end of July, the cold surface water disappeared. The temporal variations of the surface cold water and the 3-D hydrography around Dongshan Island are thought to be mainly attributed to the weakened upwelling-favorable southwestern wind, the asymmetric spatial structure of the wind field and the intrusion of warm water from the northern South China Sea.

Upwelling Process of the Western Himalaya Mountains:Height and Velocity Estimation Evidenced by Formation and Evolution of the Zanda Basin,Tibet,China

Based on field geological survey, stratigraphic section measurement and indoor comprehensive investigation, the Zanda Basin＇s tectonic location in the Himalaya Plate was ascertained, and the formation and evolution of the Zanda Basin during the Pliocene to Early Pleistocene was classified as six stages： （a） primary rift-faulting stage, （b） quick rift-faulting Stage, （c） intensive rift-faulting stage, （d） stasis stage, （e） secondary rift-faulting stage, and （f） secondary quick rift-faulting stage. Based on this six-staged formation-evolution theory of the Zanda Basin, the upwelling process of the Western Himalaya Mountains from the Pliocene to Early Pleistocene was classified as the following five stages： （a） slow upwelling stage （5.4-4.4 Ma）, （b） mid-velocity upwelling stage （4.4-3.5 Ma）, （c） quick upwelling stage （3.5-3.2 Ma）, （d） upwelling-ceasing stage （3.2-2.7 Ma）, and （e） quick upwelling stage （2.7 Ma）. Research has shown that in the duration from the Early Pliocene （4.7 Ma） to the End of Pliocene （2.67 Ma）, which lasted 2.03 million years, the Himalaya Mountains had uplifted 1500 m at a velocity of 0.74 mm/a; this belongs to a mid-velocity upwening. During the 1.31 million years in the Early Stage of the Early Pleistocene, the Himalaya Mountains had risen up another 1500 m at a velocity of 1.15 mm/a; this is a rather quick upwelling. All of these data have shown that the upwelling of the Western Himalaya Mountains is along a complicated process with multi-stages, multi-velocities, and non-uniformitarian features.

The response of spring phytoplankton assemblage to diluted water and upwelling in the eutrophic Changjiang(Yangtze River) Estuary

A comprehensive study on the phytoplankton standing stocks, species composition and dominant species in the eutrophic Changjiang（Yangtze River） Estuary（CE） was conducted to reveal the response of phytoplankton assemblage to Changjiang Diluted Water（CDW） and upwelling in the spring. Phytoplankton presented peak standing stocks（13.03 μg/L of chlorophyll a, 984.5×103 cells/L of phytoplankton abundance） along the surface isohaline of 25. Sixty-six species in 41 genera of Bacillariophyta and 33 species in 19 genera of Pyrrophyta were identified, as well as 5 species in Chlorophyta and Chrysophyta. Karenia mikimotoi was the most dominant species, followed by Prorocentrum dentatum, Paralia sulcata, Pseudo-nitzschia delicatissima and Skeletonema costatum. A bloom of K. mikimotoi was observed in the stratified stations, where the water was characterized by low nitrate, low phosphate, low turbidity, and specific ranges of temperature（18–22 °C） and salinity（27–32）. K.mikimotoi and P. dentatum accumulated densely in the upper layers along the isohaline of 25. S. costatum was distributed in the west of the isohaline of 20. Benthonic P. sulcata presented high abundance near the bottom,while spread upward at upwelling stations. CDW resulted in overt gradients of salinity, turbidity and nutritional condition, determining the spatial distribution of phytoplankton species. The restricted upwelling resulted in the upward transport of P. sulcata and exclusion of S. costatum, K. mikimotoi and P. dentatum. The results suggested that CDW and upwelling were of importance in regulating the structure and distribution of phytoplankton assemblage in the CE and the East China Sea.

The cruise observation of turbulent mixing in the upwelling region east of Hainan Island in the summer of 2012

The turbulent mixing in the upwelling region east of Hainan Island in the South China Sea is analyzed based on in situ microstructure observations made in July 2012. During the observation, strong upwelling appears in the coastal waters, which are 3℃ cooler than the offshore waters and have a salinity 1.0 greater than that of the offshore waters. The magnitude of the dissipation rate of turbulent kinetic energy ε in the upwelling region is O（10–9 W/kg）, which is comparable to the general oceanic dissipation. The inferred eddy diffusivity Kρ is O（10–6 m2/s）, which is one order of magnitude lower than that in the open ocean. The values are elevated to Kρ≈O（10–4 m2/s） near the boundaries. Weak mixing in the upwelling region is consistent with weak instability as a result of moderate shears versus strong stratifications by the joint influence of surface heating and upwelling of cold water.The validity of two fine-scale structure mixing parameterization models are tested by comparison with the observed dissipation rates. The results indicate that the model developed by Mac Kinnon and Gregg in 2003 provides relatively better estimates with magnitudes close to the observations. Mixing parameterization models need to be further improved in the coastal upwelling region.

Phytoplankton pigment pattern in the subsurface chlorophyll maximum in the South Java coastal upwelling system, Indonesia

Upwelling occurs on the coast of Java between June and October, forced by local alongshore winds associated with the southeasterly monsoon. This causes variations in phytoplankton community composition in the upwelling zone compared with the surrounding offshore area. Based on pigments analysis with subsequent calculations of group contributions to total chlorophyll a(Chl a) using CHEMTAX, we studied the distribution and composition of phytoplankton assemblages in the subsurface chlorophyll maximum along the south coast of Java and the influence of upwelling. Nineteen phytoplankton pigments were identified using high-performance liquid chromatography, and CHEMTAX analysis associated these to ten major phytoplankton groups. The phytoplankton community in the coastal area influenced by upwelling was characterized by high Chl a and fucoxanthin concentrations, indicating the dominance of diatoms. In contrast, in the offshore area, the Chl a and fucoxanthin concentrations declined to very low levels and the community was dominated by haptophytes represented by 19′-Hexanoyloxyfucoxanthin. Accordingly, microphytoplankton was found to be the major size class in the coastal area influenced by upwelling, while nanophytoplankton was most abundant in the offshore area. Low concentrations of other accessory pigments indicated less contribution from dinoflagellates,prasinophytes, chlorophytes and cryptophytes. Photo-pigment indices revealed that photosynthetic carotenoids(PSCs) were the largest component of the pigment pool, exceeding the proportion of Chl a, with the average PSCTP up to 0.62. These distribution trends can mainly be explained by phytoplankton adaption strategies to upwelling and subsurface conditions by changing species composition and adjusting the pigment pool.

Chemical characteristics and estimation on the vertical flux of N，P，Si in upwelling area of Taiwan Strait

Based on the data obtained from oceanographic surveys at the central and northern parts of the TaiwanStrait during the period from 1983 to 1988, it is showed that the nutrient content in the coastal and bottom waters isrelatively high with an upward and offward decrease, and that the upwelling in the Fujian coastal area controls the distribution of nutrients in summer.The upwelling water is characterized by 2. 29 μmol/dm￣3 of nitrate; 2. 83 μmol/dm￣3of silicate; 0. 20 μmol/dm￣3 of dissolved inorganic phosphate,respectively. Significant correlation between nutrientsand dissolved oxygen content and its degree of saturation, temperature, and salinity, respectively are found in the upweiling area. Average of N/P ratio in the area(15. 9) approaches to Redfield ratio. Their vertical flux is estimated tobe 23. 6 mg/(m￣2·d) for - P, 223 mg/(m￣2·d) for NO-N, 302 mg/(m￣2·d) for SiO-Si, respectively,which is the main source Of nutients in the areas in summer. The mean of the flux of P and NO-N is about86 % and 73% necessary for phytoplankton in the euphotic zone.

Degree of nutrient consumption as an aging index of upwelling or vertically mixed water in the northern Taiwan Strait

In order to determine the relative aging status of upwelled or vertically mixed water in the northern Taiwan Strait, a recently developed method of calculating the degree of nutrient consumption (DNC) was employed. Upwelling was detected at a water depth of about 75 m in summer. Bottom waters in the aphotic zone and newly-upwelled waters in the euphoric zone were found to be low in terms of their DNC. In general, a low DNC was noted alongside the other traditional upwelling indicators, such as lower temperature and degree of oxygen saturation, but higher salinity, apparent oxygen utilization, nutrient contents and chlorophyll-a. Enhanced vertical mixing, but without an apparent upwelling signal, was detected near the same location in winter.

Numerical studies of upwelling in coastalareas of the East China Sea-Ⅰ The tide-induced upwelling

ＮｕｍｅｒｉｃａｌｓｔｕｄｉｅｓｏｆｕｐｗｅｌｌｉｎｇｉｎｃｏａｓｔａｌａｒｅａｓｏｆｔｈｅＥａｓｔＣｈｉｎａＳｅａ－ⅠＴｈｅｔｉｄｅ－ｉｎｄｕｃｅｄｕｐｗｅｌｌｉｎｇＬｕｏＹｉｙｏｎｇ，ＹｕＧｕａｎｇｙａｏ，ＨｕａｎｇＺｕｋｅ（ＲｅｃｅｉｖｅｄＯｃ...

Winter coastal upwelling off northwest Borneo in the South China Sea

Winter coastal upwelling off northwest Borneo in the South China Sea （SCS） is investigated by using satellite data, climatological temperature and salinity fields and reanalysis data. The upwelling forms in Decem- ber, matures in January, starts to decay in February and almost disappears in March. Both Ekman trans- port induced by the alongshore winter monsoon and Ekman pumping due to orographic wind stress curl are favorable for the upwelling. Transport estimates demonstrate that the month-to-month variability of Ekman transport and Ekman pumping are both consistent with that of winter coastal upwelling, but Ek- man transport is two times larger than Ekman pumping in January and February. Under the influence of E1 Nino-Southern Oscillation （ENSO）, the upwelling shows remarkable interannual variability： during winter of El Nino （La Nina） years, an anticyclonic （a cyclonic） wind anomaly is established in the SCS, which behaves a northeasterly （southwesterly） anomaly and a positive （negative） wind stress curl anomaly off the north- west Borneo coast, enhancing （reducing） the upwelling and causing anomalous surface cooling （warming） and higher （lower） chlorophyll concentration. The sea surface temperature anomaly （SSTA） associated with ENSO off the northwest Borneo coast has an opposite phase to that off southeast Vietnam, resulting in a SSTA seesaw pattern in the southern SCS in winter.

Upwelling and anthropogenic forcing on phytoplankton productivity and community structure changes in the Zhejiang coastal area over the last 100 years

Phytoplankton productivity and community structure in marginal seas have been altered significantly during the past three decades, but it is still a challenge to distinguish the forcing mechanisms between climate change and anthropogenic activities. High time-resolution biomarker records of two 210Pb-dated sediment cores（#34： 28.5°N, 122.272°E; CJ12-1269： 28.861 9°N, 122.515 3°E） from the Min-Zhe coastal mud area were compared to reveal changes of phytoplankton productivity and community structure over the past 100 years. Phytoplankton productivity started to increase gradually from the 1970 s and increased rapidly after the late 1990 s at Site #34; and it started to increase gradually from the middle 1960 s and increased rapidly after the late 1980 s at Site CJ12-1269. Productivity of Core CJ12-1269 was higher than that of Core #34. Phytoplankton community structure variations displayed opposite patterns in the two cores. The decreasing D/B（dinosterol/brassicasterol） ratio of Core #34 since the 1960 s revealed increased diatom contribution to total productivity. In contrast, the increasing D/B ratio of Core CJ12-1269 since the 1950 s indicated increased dinoflagellate contribution to total productivity. Both the productivity increase and the increased dinoflagellate contribution in Core CJ12-1269 since the 1950–1960s were mainly caused by anthropogenic activities, as the location was closer to the Changjiang River Estuary with higher nutrient concentration and decreasing Si/N ratios. However, increased diatom contribution in Core #34 is proposed to be caused by increased coastal upwelling, with higher nutrient concentration and higher Si/N ratios.

Glacial-interglacial productivity contrasts along the eastern Arabian Sea:Dominance of convective mixing over upwelling

The western continental margin of India is one of the highly productive regions in the global ocean.Primary productivity is induced by upwelling and convective mixing during the southwest and northeast monsoons respectively.Realizing the importance of high primary productivity,a sediment core was collected below the current oxygen minimum zone（OMZ） from the southwestern continental margin of India.This was dated by AMS radiocarbon and as many as 60 paleoclimate/paceoceanographic proxies,such as particle size,biogenic components,major,trace and rare earth elements（REEs） which were measured for the first time to determine sources of sediment,biogeochemical processes operating in the water column and their variations since the last glacial cycle.R-mode factor analysis of comprehensive data indicates that the dominant regulator of paleoproductivity is the southwest monsoon wind induced upwelling.Other paleoproductivity related factors identified are the marine biogenic component and biogenic detritus（as an exported component from the water column added to the bottom sediment）.All paleoproductivity components increased significantly during the marine isotope stage-1（MIS-1）compared to those accumulated from MIS-4 to MIS-2.The second group of factors identified are the terrigenous sediments with heavy minerals like zircon and ilmenite.The terrigenous sediment,in particular,increased during MIS-2 when the sea-level was lower;however,the heavy mineral component fluctuated over time implying pulsed inputs of sediment.The diagenetic fraction and reducing component are the third group of factors identified which varied with time with increased accumulation during the MIS transitions.The primary productivity along the southwestern continental margin of India seems to have been controlled principally by the upwelling during the southwest monsoon season that was weaker from MIS-4 to MIS-2,as relative to that during the MIS-1.In contrast,increased glacial productivity noticed in sediments deposited below the current oxygen minimum zone（OMZ） along the north of the study area that can be linked to entrainment of nutrients through the intensified convective mixing of surface water during the northeast monsoon.The sequestration of greenhouse gases by the western continental margin of India was higher during glacial than interglacial cycles.

Role of Ekman Transport Versus Ekman Pumping in Driving Summer Upwelling in the South China Sea

Relative roles of Ekman transport and Ekman pumping in driving summer upwelling in the South China Sea(SCS) are examined using QuikSCAT scatterometer wind data.The major upwelling regions in the SCS are the coastal regions east and southeast of Vietnam(UESEV),east and southeast of Hainan Island(UESEH),and southeast of Guangdong province(USEG).It is shown that the Ekman transport due to alongshore winds and Ekman pumping due to offshore wind stress curl play different roles in the three upwelling systems.In UESEV,Ekman pumping and Ekman transport are equally important in generating upwelling.The Ekman transport increases linearly from 0.49 Sv in May to 1.23 Sv in August,while the Ekman pumping increases from 0.36 to 1.22 Sv during the same period.In UESEH,the mean estimates of Ekman transport and Ekman pumping are 0.14 and 0.07 Sv,respectively,indicating that 33% of the total wind-driven upwelling is due to Ekman pumping.In USEG,the mean Ekman transport is 0.041 Sv with the peak occurring in July,while Ekman pumping is much smaller(0.003 on average),indicating that the upwelling in this area is primarily driven by Ekman transport.In the summers of 2003 and 2007 following El Ni o-Southern Oscillation(ENSO) events,both Ekman transport and Ekman pumping decrease in UESEV due to the abnormally weak southwest monsoon.During the same events,however,Ekman transport is slightly enhanced and Ekman pumping is weakened in UESEH and USEG.

Nutrient Structure of the Taiwan Warm Current and Estimation of Vertical Nutrient Fluxes in Upwelling Areas in the East China Sea in Summer

Based on field data for nutrients collected on the continental shelf of the East China Sea(ECS) during summer 2006, the structure and variations of nutrients in every water mass related to the Taiwan Warm Current(TWC) were analyzed. The supplementary effect of nutrient of upwelling on harmful algal blooms(HABs) in the ECS was also estimated, based on upwelling data. Then the maintenance contribution of nutrient of upwelling to HABs was assessed. The results showed that N/P ratio is fairly low in both surface and deep layers of the TWC, which possibly controls nutrient structure of the HABs-frequently-occuring areas. In upwelling areas, the rate of phosphate(PO4-P) uptake exceeds that of nitrate(NO3-N) of the TWC. The TWC may relieve PO4-P limitation during the process of HABs. Furthermore, upwelling plays an important role in providing nutrients to HABs. After estimating nutrient fluxes(NO3-N, PO4-P, Si O3-Si) in the upwelling areas along a typical section(S07), the results showed that the nutrient uptake rate is the greatest at 10-20 m below euphotic zone, sustaining the ongoing presence of HABs. The uptake rate of PO4-P is the highest among dissolved inorganic nutrients. Therefore, upwelling is most likely the main source of PO4-P supply to HABs.