Implication of community-level ecophysiological parameterization to modelling ecosystem productivity:a case study across nine contrasting forest sites in eastern China

Parameterization is a critical step in modelling ecosystem dynamics.However,assigning parameter values can be a technical challenge for structurally complex natural plant communities;uncertainties in model simulations often arise from inappropriate model parameterization.Here we compared five methods for defining community-level specific leaf area(SLA)and leaf C:N across nine contrasting forest sites along the North-South Transect of Eastern China,including biomass-weighted average for the entire plant community(AP_BW)and four simplified selective sampling(biomass-weighted average over five dominant tree species[5DT_BW],basal area weighted average over five dominant tree species[5DT_AW],biomass-weighted average over all tree species[AT_BW]and basal area weighted average over all tree species[AT_AW]).We found that the default values for SLA and leaf C:N embedded in the Biome-BGC v4.2 were higher than the five computational methods produced across the nine sites,with deviations ranging from 28.0 to 73.3%.In addition,there were only slight deviations(<10%)between the whole plant community sampling(AP_BW)predicted NPP and the four simplified selective sampling methods,and no significant difference between the predictions of AT_BW and AP_BW except the Shennongjia site.The findings in this study highlights the critical importance of computational strategies for community-level parameterization in ecosystem process modelling,and will support the choice of parameterization methods.

Role of cancer stem cell ecosystem on breast cancer metastasis and related mouse models

Breast cancer metastasis is responsible for most breast cancer-related deaths and is influenced by many factors within the tumor ecosystem,including tumor cells and microenvironment.Breast cancer stem cells(BCSCs)constitute a small population of cancer cells with unique characteristics,including their capacity for self-renewal and differentiation.Studies have shown that BCSCs not only drive tumorigenesis but also play a crucial role in promoting metastasis in breast cancer.The tumor microenvironment(TME),composed of stromal cells,immune cells,blood vessel cells,fibroblasts,and microbes in proximity to cancer cells,is increasingly recognized for its crosstalk with BCSCs and role in BCSC survival,growth,and dissemination,thereby influencing metastatic ability.Hence,a thorough understanding of BCSCs and the TME is critical for unraveling the mechanisms underlying breast cancer metastasis.In this review,we summarize current knowledge on the roles of BCSCs and the TME in breast cancer metastasis,as well as the underlying regulatory mechanisms.Furthermore,we provide an overview of relevant mouse models used to study breast cancer metastasis,as well as treatment strategies and clinical trials addressing BCSC-TME interactions during metastasis.Overall,this study provides valuable insights for the development of effective therapeutic strategies to reduce breast cancer metastasis.

Integrating a novel irrigation approximation method with a process-based remote sensing model to estimate multi-years'winter wheat yield over the North China Plain

Accurate estimation of regional winter wheat yields is essential for understanding the food production status and ensuring national food security.However,using the existing remote sensing-based crop yield models to accurately reproduce the inter-annual and spatial variations in winter wheat yields remains challenging due to the limited ability to acquire irrigation information in water-limited regions.Thus,we proposed a new approach to approximating irrigations of winter wheat over the North China Plain(NCP),where irrigation occurs extensively during the winter wheat growing season.This approach used irrigation pattern parameters(IPPs)to define the irrigation frequency and timing.Then,they were incorporated into a newly-developed process-based and remote sensing-driven crop yield model for winter wheat(PRYM–Wheat),to improve the regional estimates of winter wheat over the NCP.The IPPs were determined using statistical yield data of reference years(2010–2015)over the NCP.Our findings showed that PRYM–Wheat with the optimal IPPs could improve the regional estimate of winter wheat yield,with an increase and decrease in the correlation coefficient(R)and root mean square error(RMSE)of 0.15(about 37%)and 0.90 t ha–1(about 41%),respectively.The data in validation years(2001–2009 and 2016–2019)were used to validate PRYM–Wheat.In addition,our findings also showed R(RMSE)of 0.80(0.62 t ha–1)on a site level,0.61(0.91 t ha–1)for Hebei Province on a county level,0.73(0.97 t ha–1)for Henan Province on a county level,and 0.55(0.75 t ha–1)for Shandong Province on a city level.Overall,PRYM–Wheat can offer a stable and robust approach to estimating regional winter wheat yield across multiple years,providing a scientific basis for ensuring regional food security.

Developing a process-based and remote sensing driven crop yield model for maize(PRYM–Maize) and its validation over the Northeast China Plain

Spatial dynamics of crop yield provide useful information for improving the production. High sensitivity of crop growth models to uncertainties in input factors and parameters and relatively coarse parameterizations in conventional remote sensing(RS) approaches limited their applications over broad regions. In this study, a process-based and remote sensing driven crop yield model for maize(PRYM–Maize) was developed to estimate regional maize yield, and it was implemented using eight data-model coupling strategies(DMCSs) over the Northeast China Plain(NECP). Simulations under eight DMCSs were validated against the prefecture-level statistics(2010–2012) reported by National Bureau of Statistics of China, and inter-compared. The 3-year averaged result could give more robust estimate than the yearly simulation for maize yield over space. A 3-year averaged validation showed that prefecture-level estimates by PRYM–Maize under DMCS8, which coupled with the development stage(DVS)-based grain-filling algorithm and RS phenology information and leaf area index(LAI), had higher correlation(R, 0.61) and smaller root mean standard error(RMSE, 1.33 t ha^(–1)) with the statistics than did PRYM–Maize under other DMCSs. The result also demonstrated that DVS-based grain-filling algorithm worked better for maize yield than did the harvest index(HI)-based method, and both RS phenology information and LAI worked for improving regional maize yield estimate. These results demonstrate that the developed PRYM–Maize under DMCS8 gives reasonable estimates for maize yield and provides scientific basis facilitating the understanding the spatial variations of maize yield over the NECP.

Exploring the Spatiotemporal Changes and Driving Forces of Ecosystem Services of Zhejiang Coasts,China,Under Sustainable Development Goals

Ecosystem services(ESs)refer to the continuous provisioning of ecosystem goods and services that benefit human beings.Over recent decades,rapid urbanization has exerted significant pressure on coastal ecosystems,resulting in biodiversity and habitat loss,environmental pollution,and the depletion of natural resources.In response to these environmental challenges,the Sustainable Development Goals(SDGs)were proposed.Given the pressing need to address these issues,understanding the changes in ESs under the SDGs is crucial for formulating specific ecological strategies.In this study,we first analyzed land use and cover change in the Zhejiang coasts of China during 2000–2020.Then,we investigated the spatiotemporal configuration of ESs by integrating carbon storage(CS),soil retention(SR),habitat quality(HQ)and water yield(WY)using the Integrated Valuation of Ecosystem Services and Tradeoffs(InVEST)model.The driving mechanisms of ESs,which varied by space and time,were also explored using the Geo-detector method.The results revealed that,over the past two decades:1)the Zhejiang coasts have experienced a significant increase of 2783.72 km^(2) in built-up land areas and a continuous decrease in farmland areas due to rapid urbanization;2)owing to higher precipitation,extensive vegetation cover,and reduced anthropogenic disturbances,forests emerge as a crucial land use type for maintaining ecosystem services such as HQ,CS,WY,and SR;3)ESs have generally declined across the entire Zhejiang coasts,with a significant decrease observed in the northern areas and an increase in the southern areas spatially;4)the expansion of built-up land areas emerged as the primary factor affecting ecosystem services,while the vegetation factor has been increasingly significant and is expected to become predominant in the near future.Our study provides insights of understanding of ecosystem service theory and emphasizing the importance of preserving biodiversity for long-term sustainable development,and valuable scientific references to support the ecological management decision-making for local governments.

Spatiotemporal variations in ecosystem services and their trade-offs and synergies against the background of the gully control and land consolidation project on the Loess Plateau,China

Studying the spatiotemporal variations in ecosystem services and their interrelationships on the Loess Plateau against the background of the gully control and land consolidation(GCLC)project has significant implications for ecological protection and quality development of the Yellow River Basin.Therefore,in this study,we took Yan'an City,Shaanxi Province of China,as the study area,selected four typical ecosystem services,including soil conservation service,water yield service,carbon storage service,and habitat quality service,and quantitatively evaluated the spatiotemporal variation characteristics and trade-offs and synergies of ecosystem services from 2010 to 2018 using the Integrated Valuation of Ecosystem Services and Trade-offs(InVEST)model.We also analysed the relationship between the GCLC project and regional ecosystem service changes in various regions(including 1 city,2 districts,and 10 counties)of Yan'an City and proposed a coordinated development strategy between the GCLC project and the ecological environment.The results showed that,from 2010 to 2018,soil conservation service decreased by 7.76%,while the other three ecosystem services changed relatively little,with water yield service increasing by 0.56% and carbon storage service and habitat quality service decreasing by 0.16% and 0.14%,respectively.The ecological environment of Yan'an City developed in a balanced way between 2010 and 2018,and the four ecosystem services showed synergistic relationships,among which the synergistic relationships between soil conservation service and water yield service and between carbon storage service and habitat quality service were significant.The GCLC project had a negative impact on the ecosystem services of Yan'an City,and the impact on carbon storage service was more significant.This study provides a theoretical basis for the scientific evaluation of the ecological benefits of the GCLC project and the realization of a win-win situation between food security and ecological security.

Integrating GPT into the Journalism and Communication Industries:The Technological Logic,Ecosystem Reconstruction,and Ethical Chal enges

Generative AI,represented by GPT(Generative Pre-trained Transformer),is now leading the technological revolution and is reconstructing the journalism and communication industries'ecosystems because of its powerful generative capacity and diverse range of outputs.While GPT is busy revolutionizing and innovating the production of news content,working patterns,and operation modes,it has also given rise to ethical concerns in regard to news authenticity,data security,humanistic values,and other related aspects.Therefore,it is imperative to initiate strategies and approaches,such as establishing a mechanism for verifying information authenticity,enhancing data security and privacy regulations,and instituting an ethical supervision and governance framework for AI,in order to facilitate the systematic advancement of AI-based news production while reinstating public trust.

Response of ecosystem carbon storage to land use change from 1985 to 2050 in the Ningxia Section of Yellow River Basin,China

Regional sustainable development necessitates a holistic understanding of spatiotemporal variations in ecosystem carbon storage(ECS),particularly in ecologically sensitive areas with arid and semi-arid climate.In this study,we calculated the ECS in the Ningxia Section of Yellow River Basin,China from 1985 to 2020 using the Integrated Valuation of Ecosystem Services and Tradeoffs(InVEST)model based on land use data.We further predicted the spatial distribution of ECS in 2050 under four land use scenarios:natural development scenario(NDS),ecological protection scenario(EPS),cultivated land protection scenario(CPS),and urban development scenario(UDS)using the patch-generating land use simulation(PLUS)model,and quantified the influences of natural and human factors on the spatial differentiation of ECS using the geographical detector(Geodetector).Results showed that the total ECS of the study area initially increased from 1985 until reaching a peak at 402.36×10^(6) t in 2010,followed by a decreasing trend to 2050.The spatial distribution of ECS was characterized by high values in the eastern and southern parts of the study area,and low values in the western and northern parts.Between 1985 and 2020,land use changes occurred mainly through the expansion of cultivated land,woodland,and construction land at the expense of unused land.The total ECS in 2050 under different land use scenarios(ranked as EPS>CPS>NDS>UDS)would be lower than that in 2020.Nighttime light was the largest contributor to the spatial differentiation of ECS,with soil type and annual mean temperature being the major natural driving factors.Findings of this study could provide guidance on the ecological construction and high-quality development in arid and semi-arid areas.

Multi-scenario Simulation of the Impact of Land Use Change on the Ecosystem Service Value in the Suzhou-Wuxi-Changzhou Metropolitan Area,China

As the most economically developed metropolitan area in China’s Yangtze River Delta,the rapid changing land use patterns of Suzhou-Wuxi-Changzhou(Su-Xi-Chang) metropolitan area have profound impacts on the ecosystem service value(ESV).Based on the patterns of land use change and the ESV change in Su-Xi-Chang metropolitan area from 2000 to 2020,we set up four scenarios:natural development scenario,urban development scenario,arable land protection scenario and ecological protection scenario,and simulated the impact of land use changes on the ESV in these scenarios.The results showed that:1) the area of built-up land in the Su-XiChang metropolitan area increased significantly from 2000 to 2020,and the area of other types of land decreased.Arable land underwent the highest transfer-out area,and was primarily converted into built-up land.The total ESV of Su-Xi-Chang metropolitan area increased initially then declined from 2000–2020,and the value of almost all individual ecosystem services decreased.2) Population density,GDP per area,night lighting intensity,and road network density can negatively impact the ESV.3) The total ESV loss under the natural development and urban development scenarios was high,and the expansion of the built-up land and the drastic shrinkage of the arable land contributed to the ESV decline under both scenarios.The total ESV under arable land protection and ecological protection scenarios increases,and therefore these scenarios are suitable for future land use optimization in Su-Xi-Chang.This study could provide a certain reference for land use planning and allocation,and offer guidance for the rational allocation of land resources.

Spatiotemporal Characteristics of Typical Ecosystem Services and Their Spatial Responses to Driving Factors in Ecologically Fragile Areas in Upper Yellow River,China

The identification of dominant driving factors for different ecosystem services(ESs)is crucial for ecological conservation and sustainable development.However,the spatial heterogeneity of the dominant driving factors affecting various ESs has not been adequately elucidated,particularly in ecologically fragile regions.This study employed the integrated valuation of ESs and trade-offs(InVEST)model to evaluate four ESs,namely,water yield(WY),soil conservation(SC),habitat quality(HQ),and carbon storage(CS),and then to identify the dominant driving factors of spatiotemporal differentiation of ES and further to characterize the spatial heterogeneity characteristics of the dominant driving factors in the eco-fragile areas of the upper Yellow River,China from 2000 to 2020.The results demonstrated that WY exhibited northeast-high and northwest-low patterns in the upper Yellow River region,while high values of SC and CS were distributed in central forested areas and a high value of HQ was distributed in vast grassland areas.The CS,WY,and SC exhibited decreasing trends over time.The most critical factors affecting WY,SC,HQ,and CS were the actual evapotranspiration,precipitation,slope,and normalized difference vegetation index,respectively.In addition,the effects of different factors on various ESs exhibited spatial heterogeneity.These results could provide spatial decision support for eco-protection and rehabilitation in ecologically fragile areas.

Spatiotemporal variations of ecosystem services and driving factors in the Tianchi Bogda Peak Natural Reserve of Xinjiang,China

Nature reserves play a significant role in providing ecosystem services and are key sites for biodiversity conservation.The Tianchi Bogda Peak Natural Reserve(TBPNR),located in Xinjiang Uygur Autonomous Region,China,is an important ecological barrier area in the temperate arid zone.The evaluation of its important ecosystem services is of great significance to improve the management level and ecological protection efficiency of the reserve.In the present study,we assessed the spatiotemporal variations of four ecosystem services(including net primary productivity(NPP),water yield,soil conservation,and habitat quality)in the TBPNR from 2000 to 2020 based on the environmental and social data using the Integrated Valuation of Ecosystem Services and Trade-offs(InVEST)model.In addition,the coldspot and hotspot areas of ecosystem services were identified by hotspot analysis,and the trade-off and synergistic relationships between ecosystem services were analyzed using factor analysis in a geographic detector.During the study period,NPP and soil conservation values in the reserve increased by 48.20%and 25.56%,respectively;conversely,water yield decreased by 16.56%,and there was no significant change in habitat quality.Spatially,both NPP and habitat quality values were higher in the northern part and lower in the southern part,whereas water yield showed an opposite trend.Correlation analysis revealed that NPP showed a synergistic relationship with habitat quality and soil conservation,and exhibited a trade-off relationship with water yield.Water yield and habitat quality also had a trade-off relationship.NPP and habitat quality were affected by annual average temperature and Normalized Difference Vegetation Index(NDVI),respectively,while water yield and soil conservation were more affected by digital elevation model(DEM).Therefore,attention should be paid to the spatial distribution and dynamics of trade-off and synergistic relationships between ecosystem services in future ecological management.The findings of the present study provide a reference that could facilitate the sustainable utilization of ecosystem services in the typical fragile areas of Northwest China.

Evaluation of Parametric Limitations in Simulating Greenhouse Gas Fluxes from Irish Arable Soils Using Three Process-Based Models

Globally a large number of process-based models have been assessed for quantification of agricultural greenhouse gas (GHG) emissions. Modelling approaches minimize the presence of spatial variability of biogeochemical processes, leading to improved estimates of GHGs as well as identifying mitigation and policy options. The comparative performance of the three dynamic models (e.g., DNDC v9.4, DailyDayCent and ECOSSE v5+) with minimum numbers of common input parameters was evaluated against measured variables. Simulations were performed on conventionally-tilled spring barley crops receiving N fertilizer at 135 - 159 kg·N·ha-1·yr-1 and crop residues at 3 t·ha-1·yr-1. For surface soil nitrate (0 - 10 cm), the ECOSSE and DNDC simulated values showed significant correlations with measured values (R2 = 0.31 - 0.55, p 0.05). Only the ECOSSE-simulated N2O fluxes showed a significant relationship (R2 = 0.33, p 0.05) with values measured from fertilized fields, but not with unfertilized ones. The DNDC and DailyDayCent models significantly underestimated seasonal/annual N2O fluxes compared to ECOSSE, with emission factors (EFs), based on an 8-year average, were 0.09%, 0.31% and 0.52%, respectively. Predictions of ecosystem respiration by both DailyDayCent and DNDC showed reasonable agreement with Eddy Covariance data (R2 = 0.34 - 0.41, p 0.05). Compared to the measured value (3624 kg·C·ha-1·yr-1), the ECOSSE underestimated annual heterotrophic respiration by 7% but this was smaller than the DNDC (50%) and DailyDayCent (24%) estimates. All models simulated CH4 uptake we

Assessing the Water-Energy-Food-Ecosystems Nexus in Smart Irrigation: A Potato Farming Case Study in Lebanon

This study examines the Water-Energy-Food-Ecosystems (WEFE) nexus in Lebanese agriculture, with a focus on the shift from conventional surface irrigation techniques to advanced smart irrigation systems in the Bekaa region, specifically targeting potato cultivation. The study quantitatively analyzes the interaction among water, energy, and agricultural outputs at the farm scale using the WEFE Nexus framework for scenario analysis. It evaluates variations in water productivity, environmental effects, and economic outcomes, offering a detailed view of existing practices and their sustainable improvement potential. The WEFE Nexus assessment demonstrates that smart irrigation integration significantly decreased resource usage: water consumption was reduced by 58%, diesel fuel use for irrigation dropped by 57%, and the demand for labor and fertilizers decreased by 47% and 49%, respectively. This change led to enhanced crop yields and increased resource efficiency, demonstrating the potential of smart irrigation as a transformative strategy for sustainable agriculture in Lebanon and other arid areas. Economic analysis showed that farmers could recover the costs of installing the smart irrigation system within 3 months. The findings highlight the need for further research on integrating smart irrigation with renewable energy, showing potential for sustainable agricultural development. .

A spatial-explicit dynamic vegetation model that couples carbon,water,and nitrogen processes for arid and semiarid ecosystems

Arid and semiarid ecosystems, or dryland, are important to global biogeochemical cycles. Dryland＇s community structure and vegetation dynamics as well as biogeochemical cycles are sensitive to changes in climate and atmospheric composition. Vegetation dynamic models has been applied in global change studies, but the com- plex interactions among the carbon （C）, water, and nitrogen （N） cycles have not been adequately addressed in the current models. In this study, a process-based vegetation dynamic model was developed to study the responses of dryland ecosystems to environmental changes, emphasizing on the interactions among the C, water, and N proc- esses. To address the interactions between the C and water processes, it not only considers the effects of annual precipitation on vegetation distribution and soil moisture on organic matter （SOM） decomposition, but also explicitly models root competition for water and the water compensation processes. To address the interactions between C and N processes, it models the soil inorganic mater processes, such as N mineralization/immobilization, denitrifica- tion/nitrification, and N leaching, as well as the root competition for soil N. The model was parameterized for major plant functional types and evaluated against field observations.

SIMPLE ECOSYSTEM MODEL OF THE CENTRAL PART OF THE EAST CHINA SEA IN SPRING

This paper’s simple ecological model to simulate the ecosystem variation and the vertical carbon flux in the central part of the East China Sea in spring, inter-reated the phytoplankton, zooplankton,autotrophic and heterotrophic bacterioplankton, nitrate, and dissolved organic carbon (DOC) in a run lasting 90 days. Except for DOC, because of poor observation precision,the major seasonal features of the vertical distribution for these components can be simulated by this model. The results show that spring bloom is just a short period of 1-2 weeks and that deposit carbon flux at the bottom interface is about 200 mg /m2 ·d in the first 20 days and then reaches its maximum of 1500mg/m2·d about 2 months later after the spring bloom.

Quantitative model of trophic interactions in Beibu Gulf ecosystem in the northern South China Sea

A mass-balanced model was constructed to determine the flow-energy in a community of fishes and invertebrates in the Beibu Gulf, northern South China Sea using Ecopath and Ecosim software. Input parameters were taken from the literature, except for the biomass of fish groups which was obtained from trawl surveys during October 1997 to May 1999 in the study area. The model consisted of 16 functional groups （boxes）, including one marine mammal and seabirds, each representing organisms with a similar role in the food web, and only covered the main trophic flow in the Beibu Gulf ecosystem. The results showed that the food web of Beibu Gulf was dominated by the detrital path and benthic invertebrates played a significant role in transferring energy from the detritus to higher trophic levels; phytoplankton was a primary producer and most utilized as a food source. Fractional trophic levels ranged from 1.0 to 4.08 with marine mammals occupying the highest trophic level. Using network analysis, the system network was mapped into a linear food chain and six discrete trophic levels were found with a mean transfer efficiency of 16.7% from the detritus, 16.2% from the primary producer within the ecosystem. The biomass density of the commercially utilized species estimated by the model is 8.46 t/km^2, only O. 48% of the net primary production.

Integrating ecosystem services evaluation and landscape pattern analysis into urban planning based on scenario prediction and regression model

Urban planning has become a widely concern for minimizing the negative effects of urban expansion on terrestrial ecosystems. We developed an interdisciplinary modeling framework to evaluate the effectiveness and shortcomings of urban expansion management strategies. A three-step method was applied to Yinchuan Plain in the northwestern of China, including(1)analyzing the relationship between landscape pattern and ecosystem service values through mathematical statistics;(2) predicting landscape pattern and ecosystem services change under different scenarios based on cellular automaton model(SLEUTH-3r model); and(3) designing and validating optimized scenario through integrating historical analysis experiments and future multi-comparison suggestions. Results have suggested that landscape composition and configuration can significantly affect regional ecosystem service values, especially the connectivity and shape of landscape. Compact urban growth policy and medium environment protection policy are the appropriate setting for urban expansion plan. Optimization validation of the combined designed scenario implied the reliability of this method. Our results highlighted the significance of integrating application of landscape pattern analysis, ecosystem service value evaluation,model simulation and multi-scenario prediction in urban planning.

Analysis of Ecosystem Degradation Factors in Yuanmou Arid-Hot Valleys Based on Interpretative Structural Model

For ecological restoration and reconstruction of the degraded area, it is an important premise to correctly understand the degradation factors of the ecosystem in the arid-hot valleys. The factors including vegetation degradation, land degradation, arid climate, policy failure, forest fire, rapid population growth, excessive deforestation, overgrazing, steep slope reclamation, economic poverty, engineering construction, lithology, slope, low cultural level, geological hazards, biological disaster, soil properties etc, were selected to study the Yuanmou arid-hot valleys. Based on the interpretative structural model （ISM）, it has found out that the degradation factors of the Yuanmou arid-hot valleys were not at the same level but in a multilevel hierarchical system with internal relations, which pointed out that the degradation mode of the arid-hot valleys was ＂straight （appearance）-penetrating-background＂. Such researches have important directive significance for the restoration and reconstruction of the arid-hot valleys ecosystem.

Adjoint Assimilation in Marine Ecosystem Models and an Example of Application

This paper aims at a review of the work carried out to date on the adjoint assimilation of data in marine ecosys-tem models since 1995. The structure and feature of the adjoint assimilation in marine ecosystem models are also introduced. To illustrate the application of the adjoint technique and its merits, a 4-variable ecosystem model coupled with a 3-D physical model is established for the Bohai Sea and the Yellow Sea. The chlorophyll concentration data derived from the SeaWiFS o-cean colour data are assimilated in the model with the technique. Some results are briefly presented.

Modeling Study of Seasonal Variation of the Pelagic-Benthic Ecosystem Characteristics of the Bohai Sea

Based on experiment data of the Sino-German comprehensive investigations in the Bohai Sea in 1998 and 1999, a simple coupled pelagic-benthic ecosystem multi-box model is used to simulate the ecosystem seasonal variation. The pelagic sub-model consists of seven state variables： phytoplankton, zooplankton, TIN, TIP, DOC, POC and dissolved oxygen （DO）. The benthic sub-model includes macro-benthos, meiobenthos, bacteria, detritus, TIN and TIP in the sediment. Besides the effects of solar radiation, water temperature and the nutrient from sea bottom exudation, land-based inputs are considered. The impact of the advection terms between the boxes is also considered. Meanwhile, the effects of the micro- bial-loop are introduced with a simple parameterization. The seasonal variations and the horizontal distributions of the ecosystem state variables of the Bohai Sea are simulated. Compared with the observations, the results of the multi-box model are reasonable. The modeled results show that about 13% of the photosynthesis primary production goes to the main food loop, 20% transfers to the benthic domain, 44% is consumed by the respiration of phytoplankton, and the rest goes to DOC. Model results also show the importance of the microbial food loop in the ecosystem of the Bohai Sea, and its contribution to the annual zooplankton production can be 60%-64%.