Search

Excessive food waste contributes to the greenhouse effect. Although low-carbon technology has the potential to reduce carbon emissions and generate renewable energy, the high costs may deter the adoption. To promote low-carbon food waste treatment, local authorities could develop carbon reduction policies targeted at the food waste treatment industry, which could lead to a game between local authorities and food waste treatment enterprises. Hence, a bi-level multi-objective optimization model is put forward to balance the decisions of local authorities and those of food waste treatment enterprises. Enterprises can choose which carbon reduction policies to implement, and local authorities determine the intensity of specific policies. After testing this model on Shenzhen’s food waste treatment system, the model’s validity and feasibility were confirmed. It was discovered that by identifying suitable carbon reduction policies, the model can help enterprises achieve more than 50 % carbon reduction. Moreover, different enterprises may choose different carbon reduction policies to implement. Optimal food waste disposal arrangements need flexible adjustments. The constructed model can serve as a decision-making tool for city managers to reduce carbon emissions in food waste treatment.

Staphylococcus aureus (S. aureus) is well known for its biofilm formation ability and is responsible for serious, chronic refractory infections worldwide. We previously demonstrated that advanced glycation end products (AGEs), a hallmark of chronic hyperglycaemia in diabetic tissues, enhanced biofilm formation by promoting eDNA release via sigB upregulation in S. aureus, contributing to the high morbidity and mortality of patients presenting a diabetic foot ulcer infection. However, the exact regulatory network has not been completely described. Here, we used pull-down assay and LC-MS/MS to identify the GlmS as a candidate regulator of sigB in S. aureus stimulated by AGEs. Dual-luciferase assays and electrophoretic mobility shift assays (EMSAs) revealed that GlmS directly upregulated the transcriptional activity of sigB. We constructed NCTC 8325 ∆glmS for further validation. qRT-PCR analysis revealed that AGEs promoted both glmS and sigB expression in the NCTC 8325 strain but had no effect on NCTC 8325 ∆glmS. NCTC 8325 ∆glmS showed a significant attenuation in biofilm formation and virulence factor expression, accompanied by a decrease in sigB expression, even under AGE stimulation. All of the changes, including pigment deficiency, decreased haemolysis ability, downregulation of hla and hld expression, and less and sparser biofilms, indicated that sigB and biofilm formation ability no longer responded to AGEs in NCTC 8325 ∆glmS. Our data extend the understanding of GlmS in the global regulatory network of S. aureus and demonstrate a new mechanism by which AGEs can upregulate GlmS, which directly regulates sigB and plays a significant role in mediating biofilm formation and virulence factor expression.

Accurate and reliable mangrove datasets are essential for the protection and management of mangrove ecosystems. Therefore, the evaluation of the current mangrove datasets and understanding the differences among them are critical. This study takes the Guangxi Beibu Gulf (GBG) and Guangdong-Hong Kong-Macao Greater Bay Area (GBA) as the study areas and analyzes the agreement and accuracy of eight mangrove datasets from 2000 to 2020 using area comparison, spatial agreement analysis, and absolute accuracy evaluation. The results show that; 1) significant differences exist in mangrove area and spatial distribution among the different mangrove datasets, with the percentage of high agreement areas ranging from 10% to 42%. 2) The overall accuracy of the evaluated mangrove datasets ranges from 56.2% to 95.6%, and the classification accuracy of mangrove datasets in inland areas is lower than the overall level. 3) There are regional differences in the quality of mangrove datasets, with the agreement and accuracy of mangrove datasets in the GBG being greater than those in the GBA. 4) Fine-scale mangrove mapping based on high-resolution remote sensing images, such as Sentinel, and global mangrove mapping based on the Google Earth Engine (GEE) cloud platform should be emphasized in the future. The findings of this study can provide guidance for data users to select appropriate mangrove datasets and a reference for future mangrove mapping research.

Wetlands provide vital ecological services for both humans and environment, necessitating continuous, refined and up-to-date mapping of wetlands for conservation and management. In this study, we developed an automated and refined wetland mapping framework integrating training sample migration method, supervised machine learning and knowledge-driven rules using Google Earth Engine (GEE) platform and open-source geospatial tools. We applied the framework to temporally dense Sentinel-1/2 imagery to produce annual refined wetland maps of the Dongting Lake Wetland (DLW) during 2015–2021. First, the continuous change detection (CCD) algorithm was utilized to migrate stable training samples. Then, annual 10 m preliminary land cover maps with 9 classes were produced using random forest algorithm and migrated samples. Ultimately, annual 10 m refined wetland maps were generated based on preliminary land cover maps via knowledge-driven rules from geometric features and available water-related inventories, with Overall Accuracy (OA) ranging from 81.82% (2015) to 93.84% (2020) and Kappa Coefficient (KC) between 0.73 (2015) and 0.91 (2020), demonstrating satisfactory performance and substantial potential for accurate, timely and type-refined wetland mapping. Our methodological framework allows rapid and accurate monitoring of wetland dynamics and could provide valuable information and methodological support for monitoring, conservation and sustainable development of wetland ecosystem.

As the second batch of international wetland cities, Wuhan and Nanchang, both provincial capital cities in China, have abundant wetland resources. An important sign of achievement of protecting the urban wetland areas is the international wetland city designation. Understanding the growth and changes of wetlands in international wetland cities is necessary for wetland protection and management. Thus, it is crucial to conduct proper wetland mapping in these international wetland cities. By studying wetland cities, sustainable planning can be provided to promote the coordinated development of wetlands and cities. In this research, Google Earth Engine and random forest machine learning were used. According to the “Data-Information-Knowledge-Wisdom” research framework, we carried out mapping of diverse wetlands in Wuhan and Nanchang with a resolution of 10 m in 2015 and 2020. Then, using the findings of wetland mapping, we examined changes in the wetland landscapes of the two cities. Finally, this study examined changes in international wetland city indicators over this time frame. The research results are as follows. 1) Our wetland mapping results in 2015 and 2020 achieved good accuracy, with an overall accuracy of over 0.90 and a kappa coefficient of over 0.85. 2) The total wetland area in both cities increased. Nanchang grew by 91.11 km2, whereas Wuhan grew by 290.68 km2. Most restored wetland areas were far from urban construction areas. In the two cities, the fragmentation of wetlands decreased, the diversity of wetlands increased, and the growth rate of wetlands was high. 3) Wetland rates rose from 17.79% to 19.07% in Nanchang and from 19.74% to 23.12% in Wuhan, according to mapping results between 2015 and 2020. The wetland protection rate in Nanchang remained unchanged, but the wetland protection rate in Wuhan decreased. Wuhan needs to strengthen the protection of increased wetlands. In addition, the study found that most of the increased areas of wetlands were previously cultivated land. The international wetland mapping framework of this study can be easily implemented in other regions of the world.

A recombinant inbred line population including 371 lines was developed by a high kernel number per spike (KNPS) genotype T1208 and a low KNPS genotype Chuannong18 (CN18). A genetic linkage map consisting of 11,583 markers was constructed by the Wheat55K SNP Array. The quantitative trait loci (QTLs) related to KNPS were detected in three years. Eight, twenty-seven, and four QTLs were identified using the ICIM-BIP, ICIM-MET, and ICIM-EPI methods, respectively. One QTL, QKnps.sau-2D.1, which was mapped on chromosome 2D, can explain 18.10% of the phenotypic variation (PVE) on average and be considered a major and stable QTL for KNPS. This QTL was located in a 0.89 Mb interval on chromosome 2D and flanked by the markers AX-109283238 and AX-111606890. Moreover, KASP-AX-111462389, a Kompetitive Allele-Specific PCR (KASP) marker which closely linked to QKnps.sau-2D.1, was designed. The genetic effect of QKnps.sau-2D.1 on KNPS was successfully confirmed in two RIL populations. The results also showed that the significant increase of KNPS and 1000-kernel weight (TKW) was caused by QKnps.sau-2D.1 overcoming the disadvantage due to the decrease of spike number (SN) and finally lead to a significant increase of grain yield. In addition, within the interval in which QKnps.sau-2D.1 is located in Chinese Spring reference genomes, only fifteen genes were found, and two genes that might associate with KNPS were identified. QKnps.sau-2D.1 may provide a new resource for the high-yield breeding of wheat in the future.

Wetland carbon storage plays an essential role in the global carbon cycle. However, in recent decades, intensive human activities and rapid urbanization have reduced wetland C stocks in the Guangdong–Hong Kong–Macau Greater Bay Area (GBA). Long-term assessment of carbon storage in the wetland ecosystems of the GBA is needed for promoting regional sustainable development. Therefore, we proposed a framework integrating the integrated valuation of ecosystem services and tradeoffs (InVEST) model and water inundation frequency, which was retrieved via multiple index water detection rule based on Google Earth Engine platform, to estimate yearly carbon storage in wetlands across the GBA at a 30-m spatial resolution during 1995–2020. The incorporation of water inundation frequency in estimating carbon storage in wetlands presents features of long time series and higher temporal frequency. The results showed that: first, the carbon storage of wetlands estimated by the InVEST model decreased from 33.38 to 16.64 Tg·C between 1995 and 2020 in the GBA, mainly owning to the loss of paddy fields. Second, there is a significant exponentially decreasing relationship between water inundation frequency and carbon storage density within the potential wetland extent, with the R2 reaching 0.998, which denotes huge potential for estimating wetland carbon storage; third, annual wetland carbon storage estimated were divided into three main stages: a rapidly increasing period (1995–2004), a falling period (2005–2013), and a slightly fluctuating period (2015–2020). This research may provide references for supporting decision-making in implementing wetland-related sustainable development and carbon-neutrality policies in the GBA.

Understanding the variation regularity of water extent can provide insights into lake conservation and management. In this study, inter- and inner-annual variations of water extent during the period of 1987–2020 were analyzed to understand the temporal and spatial distribution characteristics of Dongting Lake. We applied the Multiple Index Water Detection Rule to extract the Dongting Lake water extent quickly and accurately based on Google Earth Engine platform, and then assessed the extraction accuracy. The water surface analysis results showed that (1) based on sustainable development goals (SDG) 6.6.1, the trend of water extent showed the downward fluctuating trend from 1987 to 2020, with the overall average water extent being 1894.48 km². (2) Among the monthly average water area, the largest extent was 2477.14km² (July) and the smallest was 848.14 km² (January). Among the seasonal mean water area, summer was the largest, with an area of 2438.06 km², and winter was the smallest at 967.34 km². (3) For the water inundation frequency, seasonal water bodies accounted for the largest proportion, with 1577.85 km²; the nonwater area was the smallest, with the area of 573.02 km²; and the permanent water area was 1086.21 km². Through the analysis of the historical water body extent of the long time series of Dongting Lake, this study reflected support for SDG, for which the research idea and design can help us understand the importance and feasibility of the SDG 6.6.1 indicator.

To verify whether a low-frequency electromagnetic field (LFE field) can help reduce structural damage during the freeze–thaw process and maintain shelf life, Mongolian cheese was frozen at −10, −20, and −30 °C, then thawed at microwave or room temperature. Results showed that LFE field-assisted frozen treatment could reduce ice crystal size and protect the protein matrix structure of cheese. Frozen–thawed cheese retained 96.5% of its hardness and showed no significant difference from the fresh one in elasticity, cohesion, and chewiness. Frozen cheese showed similar but slower ripening behavior during storage, suggesting a potential application of the LFE field in the frozen storage of high-protein foods.

Numerical methods, such as finite element or finite difference, have been widely used in the past decades for modeling solid mechanics problems by solving partial differential equations (PDEs). Differently from the traditional computational paradigm employed in numerical methods, physics-informed deep learning approximates the physics domains using a neural network and embeds physics laws to regularize the network. In this work, a physics-informed neural network (PINN) is extended for application to linear elasticity problems that arise in modeling non-uniform deformation for a typical open-holed plate specimen. The main focus will be on investigating the performance of a conventional PINN approach to modeling non-uniform deformation with high stress concentration in relation to solid mechanics involving forward and inverse problems. Compared to the conventional finite element method, our results show the promise of using PINN in modeling the non-uniform deformation of materials with the occurrence of both forward and inverse problems.

Wetland cities were proposed by the International Wetland Convention Organization for the protection of urban wetlands. Few studies have performed land cover classifications for internationally recognized wetland cities or explored what contribution the classification results can make to the establishment of additional wetland cities to date. Based on Sentinel-1 and 2 data, this study used a 10-fold random forest method to classify the land cover of the first six wetland cities recognized in China. A land cover dataset, which had a resolution of 10 m and included four wetland types, was obtained and the wetland area and protected wetland areas of the six cities were calculated. The results showed that (1) the classification accuracy of six cities was good, the overall accuracy was above 90%, and the Kappa coefficient was above 0.88. (2) Cropland or forested areas were the most common non-wetland land coverage type in wetland cities and accounted for more than 20% or 40% of the land coverages, while water was the most common wetland type and accounted for more than 2% of the land coverages. From 2015 to 2020, the built area in most cities increased, while cropland and forest decreased significantly. (3) The wetland rate was 6.68–37.56% and the wetland protection rate was 49.48–73.74% in the six wetland cities. From 2015 to 2020, the wetland rate of the six cities were relatively stable, and the wetland protection rate of inland cities (Yinchuan, Changde, Harbin and Changshu) increased significantly, while those of coastal cities (Haikou and Dongying) decreased, which might be related to the change in coastline. Therefore, we found that the wetlands in these cities were well protected. Land cover classification for wetland cities can provide a reference for using remote sensing techniques used to monitor internationally wetland cities while also supporting the creation of additional wetland cities.

Doxorubicin (DOX) leads to myocardial cell damage and irreversible heart failure, which limits the clinical application of DOX. Naringin has biological functions of inhibiting inflammation, oxidative stress and apoptosis. Our aim was to investigate whether Naringin could prevent DOX-related cardiotoxicity in mice. Naringin was administered by gavage and mice were intraperitoneally injected with doxorubicin (1 mg/kg/day) for 15 days. H&E, Masson, TUNEL and others experiments were examined. NRVMs and H9C2 cells were treated with Naringin and DOX in vitro. Using IF, ELISA and Western Blot to detect the effect of Naringin and ECHS1 on cells. The results showed that Naringin could prevent DOX related cardiac injury, inhibit cardiac oxidative stress, inflammation and apoptosis in vivo and in vitro. Inhibition of ECHS1 could interfere the effect of Naringin on DOX-induced myocardial injury. Naringin may provide a new cardiac protective tool for preventing the cardiotoxicity of anthracycline drugs.

Starch biosynthesis is a complex and highly controlled process that requires coordinated activities among multiple enzymes. Starch synthase III (SSIII) is the largest protein in the starch synthase complex and its function is to lengthen long-chain amylopectin in starch synthesis. It potentially affects the activity of other key enzymes in starch synthesis through protein–protein interactions; therefore, its function and regulation play a predominant role in starch synthesis. In this review, we summarized the main research of SSIII including its biochemical characteristics, structural features, expression atlas, and regulation means. Structural features and expressional analysis indicated that SSIIIa is the main functional protein in maize endosperm rather than SSIIIb-a and SSIIIb-b, even though they are similar in the tertiary structures. The regulation investigation of SSIIIa showed that there are 13 transcription factors that control the transcription of SSIIIa. Interaction network analysis showed that SSIIIa could be involved with ten other key enzymes in starch synthesis. In conclusion, this review considerably extends our understanding of SSIII and provides the theoretical basis for improving starch synthesis by SSIII in maize.

Sustainable Development Goals (SDGs) target 11.a is a good vision for the coordinated development of the economy, society and environment in urban agglomerations. However, there was an extreme lack of indicators, data or case studies for SDG target 11.a, since it is a vague “process target”, which is not conducive to the implementation of SDG target 11.a. It is important to propose a quantitative, convenient, and local policies relevant method to promote the realization or to test the implementation effects of SDG target 11.a. Combined with socio-economic data and land use data, this study uses the methods of comprehensive evaluation model, coupling and coordination degree, and comparative advantage degree methods to study the pattern evolution, coordination characteristics and advantageous areas of production–living–ecological (PLE) functions in the Guangxi Beibu Gulf Urban Agglomeration (GBG_UA) from 1995 to 2019. The results showed that, (1) considering the spatiotemporal distribution of PLE functions, the study area has a relatively stable ecological function as well as fluctuating production and living functions. Considering the coordination characteristics of PLE functions, high–high and low–low clustering effects were observed, and primary coordination maintained the highest proportion, accounting from 55.26% in 1995 to 71.05% in 2019, indicating the SDG target 11.a level in the GBG_UA was poor. Considering the advantageous areas for PLE functions, the region mostly comprises single-function advantageous areas and a few multifunction advantageous areas, including 20 single-function advantage counties (accounting for 52%), 15 dual-function advantage counties (accounting for 39%), and three multi-function advantage counties (accounting for 7.8%), which indicates the lack of diversified land use structures in this region. (2) Optimization suggestions for the coordinated development and realization of SDG target 11.a for the GBG_UA were provided. Suggestions were made based on the radiation and driving role of Nanning city to guide the coordinated development of surrounding counties (districts). Suggestions were also made to improve the design of the integrated transportation network as well as to optimize allocation according to the resource endowment of land and to realize an upgraded ecology as well as agricultural products and services. (3) The evaluation of PLE functions is a quantitative and convenient method that can optimize national and regional development planning and test the implementation effects of SDG target 11.a. This study offers foundational knowledge for the realization of SDG target 11.a in the GBG_UA and provides a reference for the research and implementation of SDG target 11.a in other regions around the world.

Atrial fibrillation (AF) is a highly prevalent cardiac arrhythmia that leads to numerous adverse outcomes including stroke, heart failure, and death. Hyperuricemia is an important risk factor that contributes to atrium injury and AF, but the underlying molecular mechanism remains to be elucidated. In this review, we discussed the scientific evidence for clarifying the role of hyperuricemia in the pathogenesis of AF. Experimental and Clinical evidence endorse hyperuricemia as an independent risk factor for the incidence of AF. Various in vivo and in vitro investigations showed that hyperuricemia might play a critical role in the pathogenesis of AF at different UA concentrations through the activation of oxidative stress, inflammation, fibrosis, apoptosis, and immunity.

As one of the most open and dynamic regions in China, the Guangdong-Hong Kong-Macao Greater Bay Area (GBA) has been urbanizing rapidly in recent decades. The surface water in the GBA also has been suffering from urbanization and intensified human activities. The study aimed to characterize the spatiotemporal patterns and assess the losses and gains of surface water caused by urbanization in the GBA via long time-series remote sensing data, which could support the progress towards sustainable development goals (SDGs) set by the United Nations, especially for measuring SDG 6.6.1 indicator. Firstly, utilizing 4750 continuous Landsat TM/ETM+/OLI images during 1986–2020 and the Google Earth Engine cloud platform, the multiple index water detection rule (MIWDR) was performed to extract surface water extent in the GBA. Secondly, we achieved surface water dynamic type classification based on annual water inundation frequency time-series in the GBA. Finally, the spatial distribution and temporal variation of urbanization-induced water losses and gains were analyzed through a land cover transfer matrix. Results showed that (1) the average minimal and maximal surface water extents of the GBA during 1986–2020 were 2017.62 km2 and 6129.55 km2, respectively. The maximal surface water extent fell rapidly from 7897.96 km2 in 2001 to 5087.46 km2 in 2020, with a loss speed of 155.41 km2 per year (R2 = 0.86). (2) The surface water areas of permanent and dynamic types were 1529.02 km2 and 2064.99 km2 during 2000–2020, accounting for 42.54% and 57.46% of all water-related areas, respectively. (3) The surface water extent occupied by impervious land surfaces showed a significant linear downward trend (R2 = 0.98, slope = 36.41 km2 per year), while the surface water restored from impervious land surfaces denoted a slight growing trend (R2 = 0.86, slope = 0.99 km2 per year). Our study monitored the long-term changes in the surface water of the GBA, which can provide valuable information for the sustainable development of the GBA urban agglomeration. In addition, the proposed framework can easily be implemented in other similar regions worldwide.

β2-adrenergic receptor (β2AR) can be phosphorylated by G protein receptor kinases and second messenger-dependent kinases. Here, the authors demonstrate that these phosphorylation events are specific to functionally distinct and spatially segregated subpopulations of β2AR that co-exist in a single cell.

Lung cancer has the highest rate of incidence and mortality among all cancers. Most chemotherapeutic drugs used to treat lung cancer cause serious side effects and are susceptible to drug resistance. Therefore, exploring novel therapeutic targets for lung cancer is important. In this study, we evaluated the potential of TMEM16A as a drug target for lung cancer. Homoharringtonine (HHT) was identified as a novel natural product inhibitor of TMEM16A. Patch-clamp experiments showed that HHT inhibited TMEM16A activity in a concentration-dependent manner. HHT significantly inhibited the proliferation and migration of lung cancer cells with high TMEM16A expression but did not affect the growth of normal lung cells in the absence of TMEM16A expression. In vivo experiments showed that HHT inhibited the growth of lung tumors in mice and did not reduce their body weight. Finally, the molecular mechanism through which HHT inhibits lung cancer was explored by western blotting. The findings showed that HHT has the potential to regulate TMEM16A activity both in vitro and in vivo and could be a new lead compound for the development of anti-lung-cancer drugs.

Y2O3 codopedwith different Zr4+ (10% mol, 30%mol, and 49%mol) and constant Er3+ (2% mol) ions were successfully prepared via the laser annealing method. According to the XRD data, the lattice parameters decrease with the increase of Zr4+, and the grain size increases with Zr4+. The crystal space group switches from Ia-3 to Fm-3 m at a Zr4+ concentration of 49% mol. As the concentration of doped Zr4+ increases, the upconversion emission of Er3+ decreases. The nonequivalent substitution of Y3+ by Zr4+ creates lattice defects, which affect the fluorescence intensity of Er3+ ions in the material. The temperature sensing behavior of Y88%Zr10%Er2%O1.55 crystals based on the fluorescence intensity ratio between 4S3/2→4I15/2 and 4I9/2→4I15/2 transitions of Er3+ is studied. The results show that after codoping with Zr4+, the 4S3/2→4I15/2 and 4I9/2→4I15/2 transitions satisfy the Boltzmann distribution formula at temperatures above 333K. And the temperature sensing relative sensitivity can reach 2.96% K-1. [ABSTRACT FROM AUTHOR]

Background: Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of both hospital and community infections globally. It's important to illuminate the differences between community-acquired MRSA (CA-MRSA) and hospital-acquired MRSA (HA-MRSA), but there have been confusions on the definition, especially for the MRSA isolates identified within 48 h of admission. This study aimed to determine the molecular characteristics and virulence genes profile of CA and HA-MRSA isolates identified less than 48 h after hospital admission in our region.Methods: A total 62 MRSA isolates identified within 48 h after admission and the clinical data were collected. Antimicrobial susceptibility test (AST) of collected isolates were performed according to the guidelines of Clinical and Laboratory Standards Institute (CLSI) 2015, and staphylococcal cassette chromosome mec (SCCmec) typing, multilocus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE) and virulence gene profiling were performed to explore the molecular diversity.Results: SCCmec III and sequence type (ST) 239 were the most prevalent SCCmec type and ST in both CA and HA-MRSA groups. HA-MRSA group had higher prevalence of SCCmec III (87.2 %) and ST239 (79.5 %) compared with CA-MRSA (60.9 and 43.4 %, both P < 0.001), while the frequency of SCCmec IV (26.0 %) and ST59 (21.7 %) were higher in CA-MRSA than its counterpart (P < 0.001 and P = 0.003). MRSA-ST239-III was the predominant type in this study (61.3 %, 38/62), especially in HA-MRSA group (76.9 %, 30/39). However, CA-MRSA strains exhibited more diversity in genotypes in this study. Meanwhile, CA-MRSA tended to have lower resistant percentage to non-β-lactams antibiotics but more virulence genes carriage, especially the staphylococcal enterotoxins (SE) genes. Notably, seb gene was only detected in CA-MRSA isolates (52.2 %), likely a significant marker for CA-MRSA isolates. Panton-Valentine leukocidin gene (PVL) was highly detected in both groups, while appeared no significantly different between CA-MRSA (47.8 %) and HA-MRSA (43.6 %).Conclusions: Our findings support a difference in the molecular epidemiology and virulence genes profile of CA-MRSA and HA-MRSA. Furthermore, this study indicates a possible transmission from HA-MRSA to CA-MRSA, which may cause the overlap of the definition. [ABSTRACT FROM AUTHOR]

[Objective] As an innovative way of regulating territorial space, the planning of flexible space is of great significance for improving the scientific quality and effectiveness of spatial planning, and reducing the disturbance and impact of uncertain factors. This study aimed to explore the delineation of land use layout flexible space under uncertain conditions. [Methods] The study took Xuzhou City, a typical resource-based city, as an example. First, an interval optimization model was used to obtain the flexible quantity interval of various land use types. Then, the upper and lower limits of the interval served as the quantity demand to import the PLUS model to simulate the spatial layout of various land use types, and the difference between the two simulated layouts was designated as flexible space. [Results] The results show that: (1) There are differences in the size of the flexible quantity interval of each land use type, among which land use types that contribute more to addressing future uncertainty are rural residential land, farmland, water, urban industrial and mining land, while the smaller ones are grassland and bare land. (2) The flexible space distribution of various land use types shows a certain pattern, and the flexible space of farmland is concentrated in the hotspots of land use change around urban and rural areas, the flexible space of garden plot and woodland is distributed in suitable growth areas such as river banks and low hills, and the flexible space of urban industrial and mining land and rural residential land is mainly distributed in the outer edge areas of their original scope. (3) Based on the flexible space of each land use type, the land types were merged according to the zoning function of the planning, and the conflicting parts of the flexible space were processed according to the suitability to obtain the flexible space of different spatial planning zones. Among them, the flexible space areas of agricultural production zones, ecological protection zones, urban construction zones, and rural development zones account for 1.74%, 0.20%, 1.09%, and 3.31% of the total area, respectively. [Conclusion] The probability distribution of possible conflicts in land use space is the prerequisite for the delineation of flexible space. It is necessary to delineate flexible space of land use layout from an uncertain perspective. Based on interval optimization and PLUS model, the delineation method of flexible space proposed in this study seems feasible, which can provide support for planning decisions.