Measuring ecological preference for growth of the most influential weeds in weed community structure associated with agronomic and horticultural crops in Nile Delta, Egypt

The studying of the performance of weeds and weed assemblages associated with different crops through phytosociological researches derives its importance mainly from weeds adverse impact on crop productivity. Hopefully, the current study is bene cial for developing a sustainable long term weed control and soil management strategy. The objective of the present work is to measure three ecological preferences: crop diversi cation, soil type and crop seasonality. A sample area was selected comprising farmland of Nile Delta and its adjoining east and west territories, Egypt. Sixty-four sites in 12 governorates comprising 30 agroecosystems were monitored. 150 species were designated as most in uential weed-species out of 555 species recorded. Ecological preference for species for crop seasonality was evident through the results of Agglomerative hierarchical clustering (AHC). Three weed communities or weed assemblages’ groups (WAG) were identi ed. Their diversity evaluated and their linear correlation strength was determined. The delity measurements showed that 60-species have delity to WAG A; 46 to WAG B; 44 to WAG C. Measuring the ecological preference of species for crop diversi cation and soil type indicated that the growth activity of 36-species was more affected by crop diversi cation; 63-species preferential for growth in ne grained soil (FGS) and 87-species in coarse grained soil (CGS). By employing Redundancy analysis (RDA) with variation partitioning (VP), the variability in species frequency in WAG that can be explained by ecological preference of species for the three ecological preferences was determined and variation was partialized out to estimate their partial, linear effect. In fact, weed control is a recalcitrant issue and the current study revealed that a portion of the solution is in a successful crop diversi cation plan. The successful selection for a competitive crop that can be taken seriously as an adequate weed controlling mechanism in a crop rotation technique will


Introduction
The studying of the performance of weeds and weed assemblages associated with different crops through phytosociological researches derives its importance mainly from weeds adverse impact on crop productivity. The adverse impact of weeds on crop productivity was a matter of discussion in weed science researches Torun &Uygur, 2018). Several researchers concluded that the harmful effect on crop productivity increased by presence of more invasive species Kolar & Lodge, 2001). Such studies revealed that these invasive species appear to have speci c traits or speci c combinations of traits that allow them to outcompete not only the native species but the cultivated crops as well. For example, from the common invasive species traits: a higher ecological competence, higher phenotype plasticity, faster growth, more dispersal ability, faster renewal capability after local re regimens, …etc. The extent of persistence and resistance of these weeds against weed control plans as in the case of other weeds depend mainly on the genetically inherited characters and the mutations and gene ow which contribute to genetic variability and provide resistant alleles . The harmful effect of weeds on crop yield and the increase of food requirements encouraged ecologists to pay further attention and produce more researches. Several models describing the relationship between weed abundance and crop yield were produced, they have been useful for predicting yield loss and assisting with developing management guidelines (Zimdahl, 2004). That coincided with the production of "weed abundance and crop performance models" which have also been useful for characterising the effects of weeds on crop yields under different management systems. This approach allows researchers to evaluate the effectiveness of suites of tactics applied together . In addition, more advanced simulation models based on weed-crop competition relationships were developed and have further assisted management and research Van Ittersum et al., 2003;. For such programs to be visible, accurate information on the weed ora and the distribution, frequency, abundance and phenology of weed species and weed communities are prerequisite (Frick and Thomas, 1992;Ghersa and Holt, 1995). Such kind of data may also be valuable for understanding weed communities and for creating a higher biodiversity in arable land . Despite that weed assemblages associated with either agronomic or horticultural crops was a subject of study in several phytosociological researches Mashaly et al., 2002and Mahgoub, 2017, 2019a, 2019, however, the objective of this work was not adequately addressed. The present study aimed to measure ecological preference for growth of the most in uential weeds in weed community structure associated with agronomic and horticultural crops in Nile Delta, Egypt. Hopefully that is bene cial to help for the development of a sustainable long term weed control and soil management strategy.
A sample area was selected comprising farmland of Nile Delta and its adjoining east and west territories, Egypt. The delta is the major constituent of the farmed lands of Egypt and the governmental reports indicated that ≈ 61% of the priority reclaimable land through the Nile waters is located on the fringes of the Delta region where soil, in parts of these areas, is loamy in nature; cultivation can be relatively successful . Three ecological preferences which represent the main drivers for weed growth were selected: crop seasonality, crop diversi cation and soil type. No one denies the importance of crop seasonality as one of the main drivers affecting weed community structure. Through phytosociological researches, ecologists concluded that seasonality is an evident feature of weed plant growth, as it is equally evident in cropping rotation (Mahgoub, 2017 and a crop rotation is accompanied by a weed-ora rotation (El Hadidi & Kosinová, 1971). The cultivation of non-traditional crops throughout crop diversi cation, especially in reclaimed land, has become one of the means of improving agricultural economic policy in many countries, including Egypt. The European Conference on Crop Diversi cation in 2019 re ects the importance of this concept. Hence, it is worth to questionnaire about the preference of weeds for growth in response for this ecological preference. The effect of soil type and properties on vegetation structure was explicit through several phytosociological researches.
From different environmental factors, soil properties are of prime importance as they directly in uence plant growth and distribution . In his phytosociological study on the vegetation structure in Isthmus of Suez and adjoining farmland east Nile delta,  concluded that the variation of plant species composition and weed community structure was strongly spatially structured by soil type and properties. It is often that ecologists use multivariate analysis and ordination techniques to describe the relationship between a set of samples based on their attributes (Streibig, 1979;and Clarke & Warwick, 2001a) and in the current study, both multivariate analyses and ordination techniques employed.

The study area
Geography: The area of study is teardrop shaped (C. ≈ 27,231 Km 2 (10,514 mi 2 ), Figure 1). It comprised farmland of Nile Delta and adjoining east and west territories. The Nile Delta is one of the world's largest river deltas and is a rich agricultural region. The base of the Nile Delta represents the northern side of the study area, it covers some 240 km (150 mi) of Mediterranean coastline, from Alexandria in the west (coordinates: 31°12'27.3"N, 29°55'09.0"E;Latitude: 31.207592, Longitude: 29.919168) to Port Said in the east (31°16'09.6"N, 32°18'04.1"E;Lat, Long: 31.269344, 32.301133). From north to south the delta is approximately 160 km (99 mi) in length, from Baltim (31°33'19.7"N,31°05'45.3"E;Lat,Long: 31.555461,31.095909) to Cairo (30°02'40.9"N, 31°14'07.3"E;Lat, Long: 30.044698, 31.235365). The Suez Canal borders the study area from the east and the Western Desert from the west. The sample area comprises several types of ecological habitats: fertile lands of the Nile Delta; farmland facing the Nile including those of the Delta and its adjoining east and west territories; farmland at the fringes of salt marshes; those facing Mediterranean Sea; nearby the littoral sand dunes and reclaimed land in west-delta region and near to the western desert in west side, and at Isthmus of Suez and east-delta region in east side.
Climate: The Meteorological records of the area were obtained as a courtesy from the Egyptian Meteorological Authority "EMA", they include the records of the monthly averages of rainfall (mm), temperature (°C), relative humidity and evaporation from 1931 to 2017. The records revealed that from 100-200 mm (4-8 in) of rain falls on the delta area during an average year, and most of these falls in the winter months. The northernmost part is the wettest region of the sample area. With the onset of the winter season, at the end of October through the beginning of May, the northwestern delta area is exposed to 21 squalls. They are usually accompanied with lightning thunderstorms, low temperature, high wind speeds, and cloud bursts, which produce heavy rain showers. The Delta has relatively moderate temperatures, with highs usually not surpassing 31 °C (88 °F) in the summer. The area experiences its hottest temperatures in July and August, with a maximum average of 34 °C (93 °F). Winter temperatures are normally in the range of 9 °C (48 °F) at nights to 19 °C (66 °F) in the daytime. Minimum temperature records usually don't seem to reach freezing. With cooler temperatures and some rain, the Nile Delta region becomes quite humid during the winter months, but the records of the southern meteorological stations for its adjoining east and west territories, away from the water bodies of the Mediterranean, were an express of drier atmosphere and higher rates of evaporation.
Soil type: According to the geological map of Egypt (ESDAC, 2017) and other geological researches (Elbasiouny &Elbehiry, 2019 and, the following Quaternary formations were identi ed in the Nile Delta: recent Nile alluvium (Nile deposits); super cial deposits (sand dunes and sabkha deposits) and undivided quaternary (wadi and playa deposits, and raised beaches). The rst formation included cultivated fertile land of Nile Delta and the other two included reclaimed land in the north, east and west territories of the study area. The surface of the Delta in the south is relatively smooth if compared to its surface in the north (Abu AL-Izz, 1977). Soils of the Nile Delta are mostly heavy in the texture and rather compact at the surface and the humus status of the soils is fairly well (EL-GABALY et al., 1969). Thus, all soils with exception of the northernmost part, are manmade and are regarded as anthropogenic variants of the Gleysols and Fluvisols . The chemical and mechanical analyses for the sites and localities of the study area were quoted from those recorded by the Ministry of Agriculture and Land Reclamation (MALR), Department of soil survey, Egypt, for respect of its authority and farmers property. The depth of soil horizon pro les were 0 -30, 30-60, 60-90, 90-120 cm., in them the following soil properties were measured: 1) soil texture expressed as percentage for clay, silt, clay + silt, ne sand, coarse sand, 2) Water holding capacity (100 gm soil %), 3) Hydrolytic conductivity (cm./hour), 4) Soil reaction (PH), 5) Main salts in water saturated soil extract: for cations (Ca, Mg, Na, K), and for anions (CO3, HCO3, Cl), 6) Total soluble salts (%) and 7) Calcium carbonate content (CaCO3). The results of these analyses revealed that a soil type may dominate a sampling site (district), though another sampling site may comprise one or sometimes several soil types. Accordingly, the weight average was calculated to express the soil properties that dominate a certain area.
The surveyed sampling sites were clustered into two basic classes of soil type based on texture and a simple particle-size analysis according to USCS (The Uni ed Soil Classi cation System, ASTM International). The ecological preference for weed growth was rated according to soil texture from ne to coarse. Thus, the ecological preference of weed growth was measured in two basic soil types: Coarsegrained soils (CGS) which contain 50 percent or less of nes (including gravel soil and sandy soil), and Fine-grained soils (FGS) which contain more than 50 percent nes (silt and clay). Accordingly, out of the 64 sites surveyed a total of 40 sites were designated FGS and 24 sites CGS.

Field sampling design and data collection
Strati ed sampling technique (Müller-Dombois & Ellenberg, 1974: pp. 177-209) was used as an ecological sampling design method. A total of 64 georeferenced sampling sites were located randomly in 12 governorates, in the Nile Delta and its adjoining east and west territories to represent 30 agroecosystems which were de ned at the level of cultivated crops (agronomic and horticultural crops) in the different ecological habitats (refer to Figure 1). In each sampling site, eld plots (relevés) for the cultivated crops were surveyed, each of which 1000-1500 m 2 . Field plots were visited and the species were recorded in the different seasons through sequential seasonal excursions during 2018. Three visits during the winter half of the year from December to May (in January, March, and April) and the other three visits in the summer half of the year from June to November (in June, August, and September). The weed assemblages associated with the cultivated crops were recorded and the presence of species was taken to indicate degree of ecological success and sociological performance. Presence estimate percentage PE% for species in each crop was calculated in the raw data set. One hundred and fty species were categorized as "the most in uential weed species" in the weed community structure associated with agricultural and horticultural crops out of the 555 species recorded. These weeds had PE% ≥ 30% in at least in one of the agroecosystems monitored. The values of PE% were calculated as an average from species records: in March & April for the eight winter crops; in August & September for the six summer crops; in August for the two melon crops and the highest average of their PE% values for the two records groups of March-April and of August-September for perennial and orchard crops. For those weed species which are not included in these records but monitored in other months, their highest PE % record was accepted in these months. The monitored crops were traditionally classi ed according to their season (period through which the crop occupies the eld) and their position in the yearly agricultural rotation. Agronomic crops included: eight winter crops (WC): Clover, Wheat, Broad bean, Flax, Potato, Artichoke, Tomato "winter" and Vegetables "winter"; nine summer crops (SC): Cotton, Maize, Rice, Tomato "summer", Peanut, Soybean or Soya bean, Watermelon, Sweet-melon and Vegetables "summer" and two perennial crops (PC): Taro and Sugarcane. Horticultural crops (Orchards (OC)) included eleven genera: Citrus, Guava, Pear, Banana, Mango, Prunus, Apple, Olive, Pomegranate, Date palm and Grape. The rst two letters were used as a crop code during analyses and if similar with others, the subsequent ones were used (for crop codes and botanical nomenclature, refer to Appendix 1). Some crops have two types of cultivations, e.g. Tomato plantings which are cultivated as a winter crop in sandy soils and as a summer crop in clay and clay loamy soils. It should be also pointed out that several genera monitored in orchards crops (e.g. Citrus and Prunus) included more than one cultivated species, subgenus, or varieties.
The ecological preference of species for soil type, crop seasonality and crop diversi cation were measured, based on their values of Fr (frequency, Cain, 1959; and Ab (abundance, Damgaard, 2009;Verberk, 2011). Given these importance values, one can link the maximum species growth activity with its ecological preferences for a certain soil type in a certain crop type.
The comparison of species importance values (Fr and Ab) measured in crops cultivated in ne grain soils sites (FGS) with those measured in coarse grain soils sites (CGS), was taken as an express of the ecological preference of species for soil type. Indicator value of species (IndVal) was calculated and an indicator value index was elaborated (Dufrêne & Legendre, 1997; see also Legendre & Legendre, 1998). It comprised species' indicator values in FGS and CGS. The ecological preference of a weed for growth in one of the soil types rather than in the other type was determined by its higher scores (refer to Synoptic, Appendix 1). An arbitrarily threshold level of 25% chosen in the indicator value index, and the best indicator species for each soil type were selected from those species that have an indicator value (IndVal) > 25% and a statistically signi cant P-value (α = 0.05) of their calculated X 2 (Chi-square statistic, , see also Ryabko et al., 2004. Through the comparison of these values the effectiveness of the best indicator species was evaluated. For crop seasonality, the combined effects of the following items were considered as an express for the ecological preference of species: frequency and abundance values in winter crops (WC), summer crops (SC) and perennial crops (PC, OC); species seasonality measurement (SS) designated as all-the-year-round weeds (A), winter weeds (W), early-appearing winter weeds (Ws), summer weeds (S), early-appearing summer weeds (Sw) and species seasonal bias percentage (Sb%); . For the ecological preference of species for crop diversi cation, the importance values for species (Fr and Ab) of pool all groups in addition to their Cv (Coe cient of variation, Hendricks & Robey, 1936, see also Krishnamoorthy & Meesook, 2013); were taken as an express of the variable. The results were deposited in the synoptic table, Appendix 1, and the maximum values between the subsets of the aforementioned items were denoted by bold red text. The methodology and mathematical formulation of the aforementioned items were given in Table 1. Species were ranked in the synoptic table in a descending order according to their frequency values. During analyses, four letters code for species was used by combining rst two letters of genus name and species name and if similar with others, the subsequent ones were used (refer to synoptic table, Appendix 1). The Botanical Nomenclature of the recorded species have been updated from that appeared in Students' Flora of Egypt,  and in checklist of  to a more recent Plant List, created by the Collaboration between the Royal Botanic Gardens (Kew), Missouri Botanical Garden (MO) and other collaborated institutions, Version 1.1, September 2013. Voucher specimens of each recorded species were collected and identi ed earlier in Cairo University Herbarium (CAI),, where they deposited as Herbarium specimens and numbered by a serial collecting number (MAHGOUB'S collecting number)..

Diversity and Biostatistics
The following software were used during Multivariate analyses (MVA):: VEGAN packages in R environment, ver. 3.6, 2019 (Oksanen et al., 2013); JUICE, ver. 7.1, 2020 (Tichý, 2001); PAST, ver. 4, 2020 and CANOCO, ver. 4.5 (ter Braak, 1988c). Agglomerative Hierarchical Clustering (AHC) was employed, as a clustering technique using Euclidean distance as a measure of dissimilarity and Ward's method (Minimum-variance clustering) as an agglomeration criterion , to classify the weed assemblages associated with the monitored crops into weed assemblages' groups (WAG) based on the variation in their oristic composition. The data were standardized prior to analyses. The identi ed WAG were named after the two most dominant species in each group (Whittaker, 1962). Their diversity was evaluated at different levels to estimate the extent of variation in vegetative structure. The following diversity indices were measured: Species richness (S) "Taxa_S" (Magurran et al., 2004, see also ; Shannon-Wiener diversity index (H) "Shannon_H" , see also ; Equitability (E) "Equitability_J" , see also  and Dominance (D) "Dominance_D" , see also Harper, 1999. The strength of linear correlation between them was determined by Person's correlation coe cient "r" , see also . The phicoe cient of association (Φ) between species and the target vegetation units or WAG was used as a statistical measure of association or delity . Quantitative delity measure was applied using sum of all covers of the species in data set and target vegetation unit. The statistical signi cance of the phi coe cient was tested by calculating the chi-square statistic (Sokal & Rohlf, 1995: 697, 736), at a level of signi cance of 0.05 (X 2 : α = 0.05, P* > 3.84). The rule of thumb introduced by Lepš & Šmilauer (2003) was applied and detrended correspondence analysis (DCA;  had been used to estimate the compositional gradient in the vegetation data of the present study, to be < 4 SD (standard deviation units); thus, redundancy analysis (RDA; Van den Wollenberg, 1977) is the appropriate ordination method to perform direct gradient analysis (ter Braak & Prentice, 1988). RDA, extracted and summarised the variation in species records (response variables) that can be explained by soil type, crop seasonality and crop diversi cation as explanatory variables. The soil physical properties and chemical properties which dominate the site were measured in FGS and CGS and were used as a comparison criterion during RDA. Whittaker's species turnover (β W ; Whittaker, 1972, see also Tuomisto, 2010) was calculated depending on the presence-absence data to quantify the amount of species turnover among the two subunits of winter and summer half of the year, and the parameters were used as a comparison criterion during RDA, as well. A Monte Carlo permutation test (499 permutations;ter Braak, 1990) was used to test for signi cance of the eigenvalue of the rst axis (RDA axis 1). The results of the analysis have been represented by an RDA correlation biplot for species vs environmental variables ( Figure 3). Variation partitioning (VP; , Zelený, 2019 was employed with RDA and the variation was partialized out to determine the partial, linear effect of the three ecological preferences. The variation explained was compared by adjusted R 2 values. The signi cance of fractions of interest was tested by applying ANOVA, permutation test for RDA under reduced model, 999 permutations. The results of VP analysis were displayed as a Venn diagram (Figure 4).

Results
The ecological preference of species for crop seasonality was evident through the results of AHC ( Figure   2) which identi ed three weed communities or weed assemblages' groups (WAG). The weed assemblages associated with winter crops were clustered in WAG A or group Sonchus oleraceus -Rumex dentatus, summer crops in WAG B or group Bassia indica -Cyperus rotundus, perennial crops and orchards in WAG C or group Cynodon dactylon -Chenopodium murale. However, in addition to the dominant species in WAG, some other common species have characterized some sites or a number of relevés and formed alliances with less common types within the main identi ed group. According to the results of tabular comparison technique, the following associations were recognized at the level of alliance. In WAG A: Polypogon monspeliensis -Chenopodium album association and Melilotus indicusMalva parvi ora associations in winter crops cultivated in CGS, Convolvulus arvensisBeta vulgaris and Vicia sativaBrassica nigra associations in those cultivated in FGS. In WAG B: Echinochloa colona -Dinebra retro exa and Aster squamatus -Amaranthus viridis associations in summer crops cultivated in FGS and Phyla nodi ora Portulaca oleracea and Panicum repensAmaranthus graecizans associations in those cultivated in CGS. In WAG C: Urtica urensSetaria viridis association in perennial crops and orchards cultivated in FGS and Erigeron bonariensis -Digitaria sanguinalis and Euphorbia peplusSisymbrium irio associations in those cultivated in CGS. The importance values of species presented in synoptic table (Appendix 1) indicated that most of them (≈ 87%) were more affected by crop seasonality. They restrict their growth activity or showed their best performance either during the winter half of the year or during the summer half of the year. However, their likely association with the monitored agroecosystems differed. The delity measurements showed that a total of sixty species have delity to WAG A; 46 to WAG B; 44 to WAG C. The less affected species by the variable of crop seasonality were the all-the-yearround weeds (A = 19 sp.).. They had convergent records in winter crops and in summer crops, or all-theyear-round in perennial crops and orchards. Their delity measurement indicated that, 4 species had a greater degree of joint delity with WAG A, 11 species with WAG B and 4 species with WAG C. The most common of them was Cynodon dactylon and Cyperus rotundus, the two dominants of WAG B and WAG C, respectively. Their Sb% and Ab values indicated that the latter ourished more than the rst in crops cultivated during the summer half of the year, while the rst was more abundant during the above denoted season in perennial crops and orchards. The remaining 131-species were more affected by the variable. The rst group of them included 40-species designated as Ws and Sw (23 sp. and 17 sp., respectively), of them 18-species had statistically signi cant delity measurements (X 2 : α = 0.05, P* > 3.841). These species showed their best performance either in winter crops or in summer crops and had some tangible growth in the early cultivated crops during the corresponding other half of the year. This group included the other four dominant species of the identi ed WAG. The second group included 94species designated as W ( = 64 sp.) and S = ( = 27 sp.).. Some of them had statistically signi cant positive delity values (Φ) for their speci ed groups (WAG), while others showed statistically signi cant negative values for their association with non-speci c groups. This group showed the highest ecological preference for crop seasonality and their records in their preferable season of growth were as twice or more than in the corresponding other half of the year. The degree of ecological preference of these species in response to the variable were recognized by comparing their values of seasonal bias (Sb%).
The values of the diversity indices re ected the variability in the weed community structure of the identi ed WAG (A-C). The weed assemblages' group for winter crops (WAG A) scored the highest total species richness (S = 136 sp.) with the lowest equitability value (E = 0.92), while the weed assemblages' group for summer crops (WAG B) scored the highest value of dominance (D = 0.02) and the lowest species richness (S = 119 sp.).. The weed assemblages' group for orchards and perennial crops (WAG C) scored the highest Equitability value and the highest Shannon-wiener diversity index value (H = 4.7, E = 0.96). Pearson's correlation coe cient value for the paired data of weed assemblages associated with WAG A (winter crops) and WAG B (summer crops), was not statistically signi cant (r AB = 0.015 & P = 0.853). However, it was statistically signi cant (P < 0.01) for the paired data of weed assemblages associated with each of both groups and those associated with perennial-crops and orchards; either at the level of WAG (r AC = 0.450; r BC = 0.392); or at the level of the four categories of agronomic and horticultural crops presented in synoptic table (r WP = 0.608; r WO = 0.428; r SP = 0.594; r SO = 0.363; r PO = 0.661).
The results presented in synoptic table (Appendix 1) indicated also that the vast majority of the very common weeds (Fr for "pool all groups" = 80-100%) and the other less common weeds were affected by crop diversi cation. 136-species out of the recorded 150 most in uential weeds had scored coe cient of variation (C v) which exceeded 50%, but the more common species with the wider ecological amplitude were the less affected by the variable. The species Fr and Ab values had showed a medium strength of negative linear relationship with C v. The more affected weeds (36-species) scored a very high coe cient of variation (C v) which exceeded 100%. Although some of these species were recorded in all crop categories while others restrict their occurrence to one category or more, all of them showed a remarkable growth activity in one crop or more in certain soil type during their favorite growth season, while they were modestly represented or not recorded in others. Such weed performance was recorded in 13 out of the 30 agroecosystems monitored. Two species (Brni, Cipu) were recorded in Clover elds planted in FGS; 2 sp. The indicator value of species (IndVal) indicated that 63-species showed ecological preference for growth in crops cultivated in ne grained soil (FGS) and 87-species in those crops cultivated in coarse grained soil (CGS). The six dominant species of the identi ed WAG (A-C) were the less affected by the variable. The other less common species were more affected with different degrees, of them 61-species had scored species' indicator value (IndVal) in their preferred soil type twice as high as in the other soil type (22-species IndVal: FGS ≥ 2 CGS and 39-species IndVal: CGS ≥ 2 FGS). The calculated chi-square statistic (X 2 ) revealed that twenty-two species were statistically signi cant with P-value less than 0.05, of them a total of 13 species (5 sp. in FGS and 8 sp. in CGS) were considered as best indicator species (indicator value > 25%).
RDA results indicated that the constrained variance is much higher than the unconstrained one. The analysis explained 79.16% of the total variance. The permutation test indicated that both the model and eigenvalue "RAD axis 1" were statistically signi cant with α = 0.01. The rst axis explained 50.6% of the variation; with the second axis they explained a total of 69.3%. The results of RDA displayed in the correlation biplot ( Figure 3) indicated that "RDA axis 1" was positively correlated with soil content of sand ( ne sand and coarse sand) while it was negatively correlated with soil content of clay, thus it can be interpreted as the "sand-clay" gradient. "RDA axis 2" was positively correlated with soil content of "CaCO3" while it was negatively correlated with soil content of "Cl", thus it can be interpreted as the "calcium carbonate-chloride" gradient. The sets of object scores and explanatory variable scores indicated that the impact of crop diversi cation on the weed community structure associated with agronomic and horticultural crops cultivated in CGS is more than that in those cultivated in FGS. The angles between the vectors of the three ecological preferences re ected their (linear) correlation.
Conceiving the objects ordination, we can realize that "the most in uential species" have shown an ecological preference of variable growth in response to the three ecological preferences. The wider ecological amplitude species at the level of the major identi ed groups (WAG) were the less affected and most of them were ordinated around the centered axis. The vast majority of those at level of alliance within each WAG (associations) and the best indicator species were more affected and they were ordinated further far away. The codes of these species were denoted with a bold blue text in RDA correlation biplot (refer to synoptic table, Appendix 1, for corresponding botanical nomenclature). The right-angled projections of the object points of the six dominant species for WAG onto the vectors of the three environmental variables indicated their variability to attain maximum growth activity in response to ecological preferences. Four of them showed ecological preference for crops grown in FGS. However, the projection of the ordinated point of Rumex dentatus onto the variable vector of soil type, indicated that the variable's value of FGS realised for that object had higher value of the variable, relative to objects of the other dominants (Cyperus rotundus, Sonchus oleraceus and Chenopodium murale, respectively). Their coordinates also showed resistance and persistence of Rumex dentatus more than other species of growth in soils that contain a high concentration of "Cl" or in those which are irrigated with water that has a high concentration of chlorides. The other two species (Cynodon dactylon and Bassia indica) showed ecological preference for crops grown in CGS characterized by lower WHC, higher PH, HC, CaCO3, HCO3, Na, Ca, Mg, K & TSS, but the Bermuda grass excelled in soils characterized by more content of ne sand.
Visualizing the RDA biplot we can realize that the above-mentioned species associations occurred as a result of more speci c environmental conditions (ecological niche) within the major groups (WAG). Their species objects were ordinated closer together and had convergent response variable scores. Six-species were identi ed as the best indicator species in WAG A, seven-species were identi ed in WAG B and ninespecies were identi ed in WAG C (refer to IndVal * > 25% in synoptic table, Appendix 1). Their effectiveness as best indicator species for certain ecological conditions could be recognized by their coordinates.
Partitioning the variation in response matrix (species records) between the three explanatory matrices (ecological preferences), had indicated that the partial linear effect for crop seasonality on the response data was greater than that of the other ecological preferences. The adjusted R 2 values of the marginal and conditional effect of crop seasonality were 76.23% & 25.71%, respectively, while those of soil type and crop diversi cation were 46.75% & 2.17%; 47.69% & 0.35%, respectively (refer to Figure 4). The results indicated a strong correlation between the three explanatory variables and the shared variation explained by the them was much higher (Adj.R 2 = 40.6%) than the shared variation explained by only two of them (Adj.R2 = 0.4% -6.34%). ANOVA results revealed that all simple (marginal) and conditional (partial) effects of the three predictors were signi cant at P < 0.001.

Discussion
The quest to bring new lands under cultivation is considered as the cornerstone for the agricultural policy in several countries. In Egypt, desert comprise about 95% of the total land surface and except for Nile Delta and the Fayium Oasis, only a narrow strip along the Nile is cultivated and the population is concentrated in these areas. It is not surprisingly that the land reclamation was not only the cornerstone of the Egyptian agricultural policy but it was also on the political agenda since the 1950s. It has been used as a remedy for a variety of problems e.g. desertification, alleviation of poverty, and reduction of graduate unemployment (Adriansen, 2009). The expansion in cultivation of non-traditional crops, especially in reclaimed land, has become one of the means of improving agricultural economic policy in the country. It usually concerns the shift of resources from low value agriculture to high value agriculture . However, it is not only a move from occasional and low value crop(s) to high value crop(s), but also a need-based situation speci c nonstop and vibrant idea. In the past few years, crop diversi cation played an essential role in the agricultural strategy and a diverse array of nontraditional crops has been cultivated. Hence, crop diversi cation was recognized as a new main driver of weed growth. However, crop cultivation is determined by several ecological factors and usually it is structured mainly by climatic conditions, edaphic qualities which characterize the soil itself (e.g. drainage, soil texture or chemical properties such as pH, soil salinity, …etc.) and crop seasonality which depends on crop genetic adaptation. Several phytosociological researches also indicated that the structure of vegetation groups is determined by the combined effects of a whole range of ecological factors Gholinejad et al., 2012). Consequently, the topic of the present study included the measurement of preference for growth of the most in uential weeds in three ecological preferences: crop diversi cation, soil type and crop seasonality.
The values of diversity indices for the identi ed WAG (A-C) re ected the effectiveness of variability in the ecological preference of species for growth in response to the combined effect of the three ecological preferences. Their values varied from: S = 119 sp. -136 sp., H = 4.4-4.7, D = 0.01-0.02, E = 0.92-0.96. WAG A was the more diverse group, the others showed lower heterogeneity in species composition and WAG B scored the highest "D" while WAG C scored the highest "E". Such increase in species diversity in WAG A (winter crops) is mainly attributed to more records of winter annuals in those crops cultivated in CGS in the northernmost part of Nile Delta and its adjoining east and west territories. Their appearance in ecological niche, matched to the presence of their speci c environmental conditions which related rst with the plentiful of water resources that are mainly available in arid habitats through high ratios of rainfall . The results deposited in the synoptic table and those of RDA with VP (Adj.R2 = 40.6%), also indicated that the three ecological preferences were strongly correlated and the dominance of a species in one or more crops rather than it was in the others was a conclusion for its ecological preference for growth in response to the combined effect of these ecological preferences.
The species with the wider ecological amplitude were the less affected. Such difference in performance and wider ecological amplitude is often caused by phenotypic plasticity and heterogeneity ). However, the results of the analyses revealed that even the six dominant species in WAG (A-C), which had the widest ecological amplitude showed a noticeable variable preference to attain their maximum growth activity. Their indicator value (IndVal) indicated that four of them showed a preference for growth in the crops grown in FGS, while the other two showed a preference for growth in those cultivated in CGS. The impact of soil texture on above-ground net primary productivity (ANPP) of plant communities had been discussed by several researchers (Lane et al., 1998;Le HouŽrou, 1984). They showed a variable performance as a response for crop seasonality and crop diversi cation, as well. Although, their frequency values indicated that they were recorded in all agroecosystems (Fr = 100%), their abundance values in the different crop categories ranged from 46.9 to 73.3% and their coe cient of variation ranged from Cv = 4.9 to 61.4%. Such vulnerability in their performance within their identi ed WAG (A-C) was the main reason for the record of eleven plant associations in these three macro-environmental classes through a tabular comparison technique. The species objects for these associations were ordinated closer together in RDA biplot. Their parameters revealed a more speci c environmental conditions (ecological niche) through which the best indicators species were identi ed. Some species thrive well at the same soil conditions (Ellenberg et al., 1992), and the increase of convergent ecological conditions between sampling areas affects their oristic composition which depends on the species phenotypic plasticity .
The results of the current study also showed that the ecological preference of the most in uential species (150 sp.) for growth in response to crop seasonality had the priority over the other ecological preferences. Most of them (≈ 87%) were more affected by the variable and VP with RDA indicated that the variation explained by crop seasonality (marginal effect) was 76.23% of the cumulative variance (79.16%). The ecological preference of species for crop seasonality was explicit through the results of the other analyses, as well. AHC separated the weed assemblages for winter crops into WAG A, summer crops in WAG B, perennial crops and orchards in WAG C, of whom 33 species had delity values (Φ) statistically signi cant for their identi ed groups. The results of Pearson's correlation coe cient indicated that the linear correlation between WAG A and WAG B was not statistically signi cant. In Egypt, 2 agronomic crops are usually grown in a seasonal sequence: a winter crop and a summer crop. However, the ecological preference of species for crop seasonality involves their ecological preference for crop diversi cation which depends on the type of crop cultivated as a conclusion of crop rotation techniques.
Usually crop rotation is preferred rather than monocropping techniques in the vast majority of the surveyed sites. The results of VP with RDA indicated that the variation explained by the variable's conditional effect was 25.71% and the shared variation explained by crop seasonality and crop diversi cation was 46.94%. The measurement of the coe cient of variation had indicated that 136 species scored C v > 50%, of them 34 species were more affected by crop diversi cation (C v ≥ 100%). The angles between the vectors in RDA correlation biplot re ect such higher linear correlation between the two variables. The aforementioned results coincided with the previously cited conclusions of Mahgoub (2017) and El Hadidi & Kosinová (1971 which indicated that seasonality is an evident feature of weed plant growth, as it is equally evident in cropping rotation and a crop rotation is accompanied by a weedora rotation. The impact of crop type was also discussed by  and Fried et al. (2008) and they declared that the type of crop has the most signi cant impact on species composition. In his textbook, Modern Weed Control, Crafts (1975) cites several crops as potential weed smothering crops, and rotating crops with ones that kill weeds by choking them out can be a very effective method of weed control. It is a way to avoid the use of herbicides, and to gain the bene ts of crop rotation (Lotz et al., 1991).
However, the cultivation of a certain crop type is a conclusion of what is available from natural resources as, soil type, quantity of farmland suitable for cultivation, plentiful of water, prevailing climatic conditions, ecological amplitude of the crop, human requirements, …etc. . The results of the RDA con rmed the adoption of the type of crops grown through the diversi cation of crops to the type of soil, which affected the owering composition of the groups of herbs associated with the cultivated crop. The higher strength of linear correlation observed in the RDA correlation biplot between crop diversi cation and CGS more than with FGS can be mainly attributed to that traditional crops were widely cultivated in the old farmlands (FGS) of the Nile Delta; while crop diversi cation through cultivation of non-traditional crops was preferred in reclaimed land (CGS) of the adjoining east and west territories. A total of 63species showed a preference for growth in the crops grown in FGS of them 22-species had IndVal: FGS ≥ 2 CGS and 87-species in those crops grown in CGS of them 39-species had IndVal: CGS ≥ 2 FGS.
Moreover, 22-species had statistically signi cant P-value of the calculated X 2 and 13-species of them were considered as best indicator species. VP with RDA had revealed that the shared variation explained by crop diversi cation and soil type was 41%.  had declared that soil structure-the spatial arrangement of individual particles, their aggregates, and of pores-plays a multifaceted key role in the factors determining crop and vegetation performance. Furthermore, Hamblin and Hamblin (1985) and Letey (1985) had indicated the effect of soil texture and properties on vegetation structure. It was also explicit that the preference for maximum growth activity of species is correlated with management practices and other soil physical and chemical properties. The coordinates of several species on RDA axis 2 indicated their ecological preference for growth in calcareous soils. The addition of CaCo3 in the form of agricultural lime (ag lime) affects soil pH and soil nutrient availability. Calcareous soil PH is > 7.0 and adding ag lime in suitable amounts has several advantages (supply Ca & Mg, increase P availability, improve N xation by legumes, enhanced N mineralization and nitri cation, neutralize soil acidity caused by high concentration of Cl, …etc.). Gregorich et al. (1994), revealed that soil organic matter (SOM) is important in maintaining several soil properties and Voroney and Angers (1995), indicated that it can be controlled by management practices, including choice of cropping, management of crop residues and methods and intensity of tillage. In addition, Zhang et al., 2006; revealed that the pattern of assemblages within communities depend on both soil variables (e.g. N, P, K and organic matter) and other related environmental conditions as the topographic variables. It is explicit that the heterogeneous field conditions are inherent in most agricultural fields (Cook & Bramley, 1998;. Such heterogeneity can be manifested in several items as soil fertility , hydrological properties  and weed communities . Despite that several researches revealed that weed community structure is affected by many factors, as farm management practices (Derksen et al., 1994;Thomas and Frick, 1993), crop type , crop seasonality (El-Demerdash et al., 1997;Mahgoub, 2017), soil characteristics (Fried et al., 2008;; however, much of the solution may be successful implementation of a crop diversi cation plan accompanied by robust weed control plans (weeding, hoeing, ploughing, Stale seed bed, farming practices, etc.).

Conclusion
Weed control is a recalcitrant issue and the current study revealed that a portion of the solution is in a successful crop diversi cation plan, especially for crops grown for organic certi cation. The successful selection for a competitive crop that can be taken seriously as an adequate weed controlling mechanism in a crop rotation technique will cause the strongest population reduction of harmful weeds on infested farmland.

References:
Abu  Climate: The Meteorological records of the area were obtained as a courtesy from the Egyptian Meteorological Authority "EMA", they include the records of the monthly averages of rainfall (mm), temperature (°C), relative humidity and evaporation from 1931 to 2017. The records revealed that from 100-200 mm (4-8 in) of rain falls on the delta area during an average year, and most of these falls in the winter months. The northernmost part is the wettest region of the sample area. With the onset of the winter season, at the end of October through the beginning of May, the northwestern delta area is exposed to 21 squalls. They are usually accompanied with lightning thunderstorms, low temperature, high wind speeds, and cloud bursts, which produce heavy rain showers. The Delta has relatively moderate temperatures, with highs usually not surpassing 31 °C (88 °F) in the summer. The area experiences its hottest temperatures in July and August, with a maximum average of 34 °C (93 °F). Winter temperatures are normally in the range of 9 °C (48 °F) at nights to 19 °C (66 °F) in the daytime. Minimum temperature records usually don't seem to reach freezing. With cooler temperatures and some rain, the Nile Delta region becomes quite humid during the winter months, but the records of the southern meteorological stations for its adjoining east and west territories, away from the water bodies of the Mediterranean, were an express of drier atmosphere and higher rates of evaporation.
Soil type: According to the geological map of Egypt (ESDAC, 2017) and other geological researches (Elbasiouny &Elbehiry, 2019 and, the following Quaternary formations were identi ed in the Nile Delta: recent Nile alluvium (Nile deposits); super cial deposits (sand dunes and sabkha deposits) and undivided quaternary (wadi and playa deposits, and raised beaches). The rst formation included cultivated fertile land of Nile Delta and the other two included reclaimed land in the north, east and west territories of the study area. The surface of the Delta in the south is relatively smooth if compared to its surface in the north (Abu AL-Izz, 1977). Soils of the Nile Delta are mostly heavy in the texture and rather compact at the surface and the humus status of the soils is fairly well (EL-GABALY et al., 1969). Thus, all soils with exception of the northernmost part, are manmade and are regarded as anthropogenic variants of the Gleysols and Fluvisols . The chemical and mechanical analyses for the sites and localities of the study area were quoted from those recorded by the Ministry of Agriculture and Land Reclamation (MALR), Department of soil survey, Egypt, for respect of its authority and farmers property.
The depth of soil horizon pro les were 0 -30, 30 -60, 60 -90, 90 -120 cm., in them the following soil properties were measured: 1) soil texture expressed as percentage for clay, silt, clay + silt, ne sand, coarse sand, 2) Water holding capacity (100 gm soil %), 3) Hydrolytic conductivity (cm./hour), 4) Soil reaction (PH), 5) Main salts in water saturated soil extract: for cations (Ca, Mg, Na, K), and for anions (CO3, HCO3, Cl), 6) Total soluble salts (%) and 7) Calcium carbonate content (CaCO3). The results of these analyses revealed that a soil type may dominate a sampling site (district), though another sampling site may comprise one or sometimes several soil types. Accordingly, the weight average was calculated to express the soil properties that dominate a certain area.
The surveyed sampling sites were clustered into two basic classes of soil type based on texture and a simple particle-size analysis according to USCS (The Uni ed Soil Classi cation System, ASTM International). The ecological preference for weed growth was rated according to soil texture from ne to coarse. Thus, the ecological preference of weed growth was measured in two basic soil types: Coarsegrained soils (CGS) which contain 50 percent or less of nes (including gravel soil and sandy soil), and Fine-grained soils (FGS) which contain more than 50 percent nes (silt and clay). Accordingly, out of the 64 sites surveyed a total of 40 sites were designated FGS and 24 sites CGS.

Field sampling design and data collection
Strati ed sampling technique (Müller-Dombois & Ellenberg, 1974: pp. 177-209) was used as an ecological sampling design method. A total of 64 georeferenced sampling sites were located randomly in 12 governorates, in the Nile Delta and its adjoining east and west territories to represent 30 agroecosystems which were de ned at the level of cultivated crops (agronomic and horticultural crops) in the different ecological habitats (refer to Figure 1). In each sampling site, eld plots (relevés) for the cultivated crops were surveyed, each of which 1000 -1500 m 2 . Field plots were visited and the species were recorded in the different seasons through sequential seasonal excursions during 2018. Three visits during the winter half of the year from December to May (in January, March, and April) and the other three visits in the summer half of the year from June to November (in June, August, and September). The weed assemblages associated with the cultivated crops were recorded and the presence of species was taken to indicate degree of ecological success and sociological performance. Presence estimate percentage PE% for species in each crop was calculated in the raw data set. One hundred and fty species were categorized as "the most in uential weed species" in the weed community structure associated with agricultural and horticultural crops out of the 555 species recorded. These weeds had PE% ≥ 30% in at crop in clay and clay loamy soils. It should be also pointed out that several genera monitored in orchards crops (e.g. Citrus and Prunus) included more than one cultivated species, subgenus, or varieties.
The ecological preference of species for soil type, crop seasonality and crop diversi cation were measured, based on their values of Fr (frequency, Cain, 1959; and Ab (abundance, Damgaard, 2009;Verberk, 2011). Given these importance values, one can link the maximum species growth activity with its ecological preferences for a certain soil type in a certain crop type.
The comparison of species importance values (Fr and Ab) measured in crops cultivated in ne grain soils sites (FGS) with those measured in coarse grain soils sites (CGS), was taken as an express of the ecological preference of species for soil type. Indicator value of species (IndVal) was calculated and an indicator value index was elaborated (Dufrêne & Legendre, 1997; see also Legendre & Legendre, 1998). It comprised species' indicator values in FGS and CGS. The ecological preference of a weed for growth in one of the soil types rather than in the other type was determined by its higher scores (refer to Synoptic, Appendix 1). An arbitrarily threshold level of 25% chosen in the indicator value index, and the best indicator species for each soil type were selected from those species that have an indicator value (IndVal) > 25% and a statistically signi cant P-value (α = 0.05) of their calculated X 2 (Chi-square statistic, , see also Ryabko et al., 2004. Through the comparison of these values the effectiveness of the best indicator species was evaluated. For crop seasonality, the combined effects of the following items were considered as an express for the ecological preference of species: frequency and abundance values in winter crops (WC), summer crops (SC) and perennial crops (PC, OC); species seasonality measurement (SS) designated as all-the-year-round weeds (A), winter weeds (W), early-appearing winter weeds (Ws), summer weeds (S), early-appearing summer weeds (Sw) and species seasonal bias percentage (Sb%); . For the ecological preference of species for crop diversi cation, the importance values for species (Fr and Ab) of pool all groups in addition to their Cv (Coe cient of variation, Hendricks & Robey, 1936, see also Krishnamoorthy & Meesook, 2013); were taken as an express of the variable. The results were deposited in the synoptic Botanical Nomenclature of the recorded species have been updated from that appeared in Students' Flora of Egypt,

Diversity and Biostatistics
The following software were used during Multivariate analyses (MVA): VEGAN packages in R environment, ver. 3.6, 2019 (Oksanen et al., 2013); JUICE, ver. 7.1, 2020 (Tichý, 2001); PAST, ver. 4, 2020  and CANOCO, ver. 4.5 (ter Braak, 1988c). Agglomerative Hierarchical Clustering (AHC) was employed, as a clustering technique using Euclidean distance as a measure of dissimilarity and Ward's method (Minimum-variance clustering) as an agglomeration criterion , to classify the weed assemblages associated with the monitored crops into weed assemblages' groups (WAG) based on the variation in their oristic composition. The data were standardized prior to analyses. The identi ed WAG were named after the two most dominant species in each group (Whittaker, 1962). Their diversity was evaluated at different levels to estimate the extent of variation in vegetative structure. The following diversity indices were measured: Species richness (S) "Taxa_S" (Magurran et al., 2004, see also ; Shannon-Wiener diversity index (H) "Shannon_H" , see also ; Equitability (E) "Equitability_J" , see also  and Dominance (D) "Dominance_D" , see also Harper, 1999. The strength of linear correlation between them was determined by Person's correlation coe cient "r" , see also . The phicoe cient of association (Φ) between species and the target vegetation units or WAG was used as a statistical measure of association or delity . Quantitative delity measure was applied using sum of all covers of the species in data set and target vegetation unit. The statistical signi cance of the phi coe cient was tested by calculating the chi-square statistic (Sokal & Rohlf, 1995: 697, 736), at a level of signi cance of 0.05 (X 2 : α = 0.05, P* > 3.84). The rule of thumb introduced by Lepš & Šmilauer (2003) was applied and detrended correspondence analysis (DCA;  had been used to estimate the compositional gradient in the vegetation data of the present study, to be < 4 SD (standard deviation units); thus, redundancy analysis (RDA; Van den Wollenberg, 1977) is the appropriate ordination method to perform direct gradient analysis (ter Braak & Prentice, 1988). RDA, extracted and summarised the variation in species records (response variables) that can be explained by soil type, crop seasonality and crop diversi cation as explanatory variables. The soil physical properties and chemical properties which dominate the site were measured in FGS and CGS and were used as a comparison criterion during RDA. Whittaker's species turnover (β W ; Whittaker, 1972, see also Tuomisto, 2010) was calculated depending on the presence-absence data to quantify the amount of species turnover among the two subunits of winter and summer half of the year, and the parameters were used as a comparison criterion during RDA, as well. A Monte Carlo permutation test (499 permutations; ter Braak, 1990) was used to test for signi cance of the eigenvalue of the rst axis (RDA axis 1). The results of the analysis have been represented by an RDA correlation biplot for species vs environmental variables ( Figure 3). Variation partitioning (VP; , Zelený, 2019 was employed with RDA and the variation was partialized out to determine the partial, linear effect of the three ecological preferences. The variation explained was compared by adjusted R 2 values. The signi cance of fractions of interest was tested by applying ANOVA, permutation test for RDA under reduced model, 999 permutations. The results of VP analysis were displayed as a Venn diagram (Figure 4).

Results
The ecological preference of species for crop seasonality was evident through the results of AHC ( Figure   2) which identi ed three weed communities or weed assemblages' groups (WAG). The weed assemblages associated with winter crops were clustered in WAG A or group Sonchus oleraceus -Rumex dentatus, summer crops in WAG B or group Bassia indica -Cyperus rotundus, perennial crops and orchards in WAG C or group Cynodon dactylon -Chenopodium murale. However, in addition to the dominant species in WAG, some other common species have characterized some sites or a number of relevés and formed alliances with less common types within the main identi ed group. According to the results of tabular comparison technique, the following associations were recognized at the level of alliance. In WAG A: Polypogon monspeliensis -Chenopodium album association and Melilotus indicus -Malva parvi ora associations in winter crops cultivated in CGS, Convolvulus arvensis -Beta vulgaris and Vicia sativa -Brassica nigra associations in those cultivated in FGS. In WAG B: Echinochloa colona -Dinebra retro exa and Aster squamatus -Amaranthus viridis associations in summer crops cultivated in FGS and Phyla nodi ora -Portulaca oleracea and Panicum repens -Amaranthus graecizans associations in those cultivated in CGS. In WAG C: Urtica urens -Setaria viridis association in perennial crops and orchards cultivated in FGS and Erigeron bonariensis -Digitaria sanguinalis and Euphorbia peplus -Sisymbrium irio associations in those cultivated in CGS. The importance values of species presented in synoptic table (Appendix 1) indicated that most of them (≈ 87%) were more affected by crop seasonality. They restrict their growth activity or showed their best performance either during the winter half of the year or during the summer half of the year. However, their likely association with the monitored agroecosystems differed. The delity measurements showed that a total of sixty species have delity to WAG A; 46 to WAG B; 44 to WAG C. The less affected species by the variable of crop seasonality were the all-the-year-round weeds (A = 19 sp.). They had convergent records in winter crops and in summer crops, or all-the-year-round in perennial crops and orchards. Their delity measurement indicated that, 4 species had a greater degree of joint delity with WAG A, 11 species with WAG B and 4 species with WAG C. The most common of them was Cynodon dactylon and Cyperus rotundus, the two dominants of WAG B and WAG C, respectively. Their Sb% and Ab values indicated that the latter ourished more than the rst in crops cultivated during the summer half of the year, while the rst was more abundant during the above The results presented in synoptic table (Appendix 1) indicated also that the vast majority of the very common weeds (Fr for "pool all groups" = 80 -100%) and the other less common weeds were affected by crop diversi cation. 136-species out of the recorded 150 most in uential weeds had scored coe cient of variation (C v ) which exceeded 50%, but the more common species with the wider ecological amplitude were the less affected by the variable. The species Fr and Ab values had showed a medium strength of negative linear relationship with C v . The more affected weeds (36-species) scored a very high coe cient of variation (C v ) which exceeded 100%. Although some of these species were recorded in all crop categories while others restrict their occurrence to one category or more, all of them showed a remarkable growth activity in one crop or more in certain soil type during their favorite growth season, while they were modestly represented or not recorded in others. Such weed performance was recorded in 13 out of the 30 agroecosystems monitored. Two species (Brni, Cipu) were recorded in Clover elds planted in FGS; 2 sp. The indicator value of species (IndVal) indicated that 63-species showed ecological preference for growth in crops cultivated in ne grained soil (FGS) and 87-species in those crops cultivated in coarse grained soil (CGS). The six dominant species of the identi ed WAG (A-C) were the less affected by the variable. The other less common species were more affected with different degrees, of them 61-species had scored species' indicator value (IndVal) in their preferred soil type twice as high as in the other soil type (22-species IndVal: FGS ≥ 2 CGS and 39-species IndVal: CGS ≥ 2 FGS). The calculated chi-square statistic (X 2 ) revealed that twenty-two species were statistically signi cant with P-value less than 0.05, of them a total of 13 species (5 sp. in FGS and 8 sp. in CGS) were considered as best indicator species (indicator value > 25%).
RDA results indicated that the constrained variance is much higher than the unconstrained one. The analysis explained 79.16% of the total variance. The permutation test indicated that both the model and eigenvalue "RAD axis 1" were statistically signi cant with α = 0.01. The rst axis explained 50.6% of the variation; with the second axis they explained a total of 69.3%. The results of RDA displayed in the correlation biplot (Figure 3) indicated that "RDA axis 1" was positively correlated with soil content of sand ( ne sand and coarse sand) while it was negatively correlated with soil content of clay, thus it can be Page 36/52 interpreted as the "sand -clay" gradient. "RDA axis 2" was positively correlated with soil content of "CaCO3" while it was negatively correlated with soil content of "Cl", thus it can be interpreted as the "calcium carbonate -chloride" gradient. The sets of object scores and explanatory variable scores indicated that the impact of crop diversi cation on the weed community structure associated with agronomic and horticultural crops cultivated in CGS is more than that in those cultivated in FGS. The angles between the vectors of the three ecological preferences re ected their (linear) correlation.
Conceiving the objects ordination, we can realize that "the most in uential species" have shown an ecological preference of variable growth in response to the three ecological preferences. The wider ecological amplitude species at the level of the major identi ed groups (WAG) were the less affected and most of them were ordinated around the centered axis. The vast majority of those at level of alliance within each WAG (associations) and the best indicator species were more affected and they were ordinated further far away. The codes of these species were denoted with a bold blue text in RDA correlation biplot (refer to synoptic table, Appendix 1, for corresponding botanical nomenclature). The right-angled projections of the object points of the six dominant species for WAG onto the vectors of the three environmental variables indicated their variability to attain maximum growth activity in response to ecological preferences. Four of them showed ecological preference for crops grown in FGS. However, the projection of the ordinated point of Rumex dentatus onto the variable vector of soil type, indicated that the variable's value of FGS realised for that object had higher value of the variable, relative to objects of the other dominants (Cyperus rotundus, Sonchus oleraceus and Chenopodium murale, respectively). Their coordinates also showed resistance and persistence of Rumex dentatus more than other species of growth in soils that contain a high concentration of "Cl" or in those which are irrigated with water that has a high concentration of chlorides. The other two species (Cynodon dactylon and Bassia indica) showed ecological preference for crops grown in CGS characterized by lower WHC, higher PH, HC, CaCO3, HCO3, Na, Ca, Mg, K & TSS, but the Bermuda grass excelled in soils characterized by more content of ne sand.
Visualizing the RDA biplot we can realize that the above-mentioned species associations occurred as a result of more speci c environmental conditions (ecological niche) within the major groups (WAG). Their species objects were ordinated closer together and had convergent response variable scores. Six-species were identi ed as the best indicator species in WAG A, seven-species were identi ed in WAG B and ninespecies were identi ed in WAG C (refer to IndVal * > 25% in synoptic table, Appendix 1). Their effectiveness as best indicator species for certain ecological conditions could be recognized by their coordinates.
Partitioning the variation in response matrix (species records) between the three explanatory matrices (ecological preferences), had indicated that the partial linear effect for crop seasonality on the response data was greater than that of the other ecological preferences. The adjusted R 2 values of the marginal and conditional effect of crop seasonality were 76.23% & 25.71%, respectively, while those of soil type and crop diversi cation were 46.75% & 2.17%; 47.69% & 0.35%, respectively (refer to Figure 4). The results indicated a strong correlation between the three explanatory variables and the shared variation explained by the them was much higher (Adj.R 2 = 40.6%) than the shared variation explained by only two of them (Adj.R2 = 0.4% -6.34%). ANOVA results revealed that all simple (marginal) and conditional (partial) effects of the three predictors were signi cant at P < 0.001.

Discussion
The quest to bring new lands under cultivation is considered as the cornerstone for the agricultural policy in several countries. In Egypt, desert comprise about 95% of the total land surface and except for Nile Delta and the Fayium Oasis, only a narrow strip along the Nile is cultivated and the population is concentrated in these areas. It is not surprisingly that the land reclamation was not only the cornerstone of the Egyptian agricultural policy but it was also on the political agenda since the 1950s. It has been used as a remedy for a variety of problems e.g. desertification, alleviation of poverty, and reduction of graduate unemployment (Adriansen, 2009). The expansion in cultivation of non-traditional crops, especially in reclaimed land, has become one of the means of improving agricultural economic policy in the country. It usually concerns the shift of resources from low value agriculture to high value agriculture . However, it is not only a move from occasional and low value crop(s) to high value crop(s), but also a need-based situation speci c nonstop and vibrant idea. In the past few years, crop diversi cation played an essential role in the agricultural strategy and a diverse array of nontraditional crops has been cultivated. Hence, crop diversi cation was recognized as a new main driver of weed growth. However, crop cultivation is determined by several ecological factors and usually it is structured mainly by climatic conditions, edaphic qualities which characterize the soil itself (e.g. drainage, soil texture or chemical properties such as pH, soil salinity, …etc.) and crop seasonality which depends on crop genetic adaptation. Several phytosociological researches also indicated that the structure of vegetation groups is determined by the combined effects of a whole range of ecological factors Gholinejad et al., 2012). Consequently, the topic of the present study included the measurement of preference for growth of the most in uential weeds in three ecological preferences: crop diversi cation, soil type and crop seasonality.
The values of diversity indices for the identi ed WAG (A-C) re ected the effectiveness of variability in the ecological preference of species for growth in response to the combined effect of the three ecological preferences. Their values varied from: S = 119 sp. -136 sp., H = 4.4 -4.7, D = 0.01 -0.02, E = 0.92 -0.96. WAG A was the more diverse group, the others showed lower heterogeneity in species composition and WAG B scored the highest "D" while WAG C scored the highest "E". Such increase in species diversity in WAG A (winter crops) is mainly attributed to more records of winter annuals in those crops cultivated in CGS in the northernmost part of Nile Delta and its adjoining east and west territories. Their appearance in ecological niche, matched to the presence of their speci c environmental conditions which related rst with the plentiful of water resources that are mainly available in arid habitats through high ratios of rainfall . The results deposited in the synoptic table and those of RDA with VP (Adj.R2 = 40.6%), also indicated that the three ecological preferences were strongly correlated and the dominance of a species in one or more crops rather than it was in the others was a conclusion for its ecological preference for growth in response to the combined effect of these ecological preferences.
The species with the wider ecological amplitude were the less affected. Such difference in performance and wider ecological amplitude is often caused by phenotypic plasticity and heterogeneity . However, the results of the analyses revealed that even the six dominant species in WAG (A-C), which had the widest ecological amplitude showed a noticeable variable preference to attain their maximum growth activity. Their indicator value (IndVal) indicated that four of them showed a preference for growth in the crops grown in FGS, while the other two showed a preference for growth in those cultivated in CGS. The impact of soil texture on above-ground net primary productivity (ANPP) of plant communities had been discussed by several researchers (Lane et al., 1998;Le HouŽrou, 1984). They showed a variable performance as a response for crop seasonality and crop diversi cation, as well. Although, their frequency values indicated that they were recorded in all agroecosystems (Fr = 100%), their abundance values in the different crop categories ranged from 46.9 to 73.3% and their coe cient of variation ranged from Cv = 4.9 to 61.4%. Such vulnerability in their performance within their identi ed WAG (A-C) was the main reason for the record of eleven plant associations in these three macro-environmental classes through a tabular comparison technique. The species objects for these associations were ordinated closer together in RDA biplot. Their parameters revealed a more speci c environmental conditions (ecological niche) through which the best indicators species were identi ed. Some species thrive well at the same soil conditions (Ellenberg et al., 1992), and the increase of convergent ecological conditions between sampling areas affects their oristic composition which depends on the species phenotypic plasticity .
The results of the current study also showed that the ecological preference of the most in uential species (150 sp.) for growth in response to crop seasonality had the priority over the other ecological preferences. Most of them (≈ 87%) were more affected by the variable and VP with RDA indicated that the variation explained by crop seasonality (marginal effect) was 76.23% of the cumulative variance (79.16%). The ecological preference of species for crop seasonality was explicit through the results of the other analyses, as well. AHC separated the weed assemblages for winter crops into WAG A, summer crops in WAG B, perennial crops and orchards in WAG C, of whom 33 species had delity values (Φ) statistically signi cant for their identi ed groups. The results of Pearson's correlation coe cient indicated that the linear correlation between WAG A and WAG B was not statistically signi cant. In Egypt, 2 agronomic crops are usually grown in a seasonal sequence: a winter crop and a summer crop. However, the ecological preference of species for crop seasonality involves their ecological preference for crop diversi cation which depends on the type of crop cultivated as a conclusion of crop rotation techniques. Usually crop rotation is preferred rather than monocropping techniques in the vast majority of the surveyed sites. The results of VP with RDA indicated that the variation explained by the variable's conditional effect was 25.71% and the shared variation explained by crop seasonality and crop diversi cation was 46.94%. The measurement of the coe cient of variation had indicated that 136 species scored C v > 50%, of them 34 species were more affected by crop diversi cation (C v ≥ 100%). The angles between the vectors in RDA correlation biplot re ect such higher linear correlation between the two variables. The aforementioned results coincided with the previously cited conclusions of Mahgoub (2017 and El Hadidi & Kosinová (1971) which indicated that seasonality is an evident feature of weed plant growth, as it is equally evident in cropping rotation and a crop rotation is accompanied by a weedora rotation. The impact of crop type was also discussed by  and Fried et al. (2008) and they declared that the type of crop has the most signi cant impact on species composition. In his textbook, Modern Weed Control, Crafts (1975) cites several crops as potential weed smothering crops, and rotating crops with ones that kill weeds by choking them out can be a very effective method of weed control. It is a way to avoid the use of herbicides, and to gain the bene ts of crop rotation (Lotz et al., 1991).
However, the cultivation of a certain crop type is a conclusion of what is available from natural resources as, soil type, quantity of farmland suitable for cultivation, plentiful of water, prevailing climatic conditions, ecological amplitude of the crop, human requirements, …etc. . The results of the RDA con rmed the adoption of the type of crops grown through the diversi cation of crops to the type of soil, which affected the owering composition of the groups of herbs associated with the cultivated crop. The higher strength of linear correlation observed in the RDA correlation biplot between crop diversi cation and CGS more than with FGS can be mainly attributed to that traditional crops were widely cultivated in the old farmlands (FGS) of the Nile Delta; while crop diversi cation through cultivation of non-traditional crops was preferred in reclaimed land (CGS) of the adjoining east and west territories. A total of 63species showed a preference for growth in the crops grown in FGS of them 22-species had IndVal: FGS ≥ 2 CGS and 87-species in those crops grown in CGS of them 39-species had IndVal: CGS ≥ 2 FGS.
Moreover, 22-species had statistically signi cant P-value of the calculated X 2 and 13-species of them were considered as best indicator species. VP with RDA had revealed that the shared variation explained by crop diversi cation and soil type was 41%.  had declared that soil structure-the spatial arrangement of individual particles, their aggregates, and of pores-plays a multifaceted key role in the factors determining crop and vegetation performance. Furthermore, Hamblin and Hamblin (1985) and Letey (1985) had indicated the effect of soil texture and properties on vegetation structure. It was also explicit that the preference for maximum growth activity of species is correlated with management practices and other soil physical and chemical properties. The coordinates of several species on RDA axis 2 indicated their ecological preference for growth in calcareous soils. The addition of CaCo3 in the form of agricultural lime (ag lime) affects soil pH and soil nutrient availability. Calcareous soil PH is > 7.0 and adding ag lime in suitable amounts has several advantages (supply Ca & Mg, increase P availability, improve N xation by legumes, enhanced N mineralization and nitri cation, neutralize soil acidity caused by high concentration of Cl, …etc.). Gregorich et al. (1994), revealed that soil organic matter (SOM) is important in maintaining several soil properties and Voroney and Angers (1995), indicated that it can be controlled by management practices, including choice of cropping, management of crop residues and methods and intensity of tillage. In addition, Zhang et al., 2006; revealed that the pattern of assemblages within communities depend on both soil variables (e.g. N, P, K and organic matter) and other related environmental conditions as the topographic variables. It is explicit that the heterogeneous field conditions are inherent in most agricultural fields (Cook & Bramley, 1998;. Such heterogeneity can be manifested in several items as soil fertility , hydrological properties  and weed communities . Despite that several researches revealed that weed community structure is affected by many factors, as farm management practices (Derksen et al., 1994;Thomas and Frick, 1993), crop type , crop seasonality (El-Demerdash et al., 1997;Mahgoub, 2017), soil characteristics (Fried et al., 2008;; however, much of the solution may be successful implementation of a crop diversi cation plan accompanied by robust weed control plans (weeding, hoeing, ploughing, Stale seed bed, farming practices, etc.).

Conclusions
Weed control is a recalcitrant issue and the current study revealed that a portion of the solution is in a successful crop diversi cation plan, especially for crops grown for organic certi cation. The successful selection for a competitive crop that can be taken seriously as an adequate weed controlling mechanism in a crop rotation technique will cause the strongest population reduction of harmful weeds on infested farmland. Figure 1 Location Map of the surveyed area. A total of 64 georeferenced sampling sites were surveyed in 12 governorates to represent 30 agronomic and horticultural crops cultivated in the different ecological habitats. Sampling sites characterized by "FGS" represented by green dots (40) while those characterized by "CGS" represented by yellow dots (24). "Map adapted from Bing Maps. Microsoft product screen shot(s) reprinted with permission from Microsoft Corporation. © 2019-2020 Microsoft".

Figure 2
Fan-shaped dendrogram representing AHC for weed assemblages associated with 30 agroecosystems surveyed (for crop legend, refer to synoptic table, Appendix 1).

Figure 4
Venn's diagram, showing partition of the variation explained by the three ecological preferences (79.16%). Adjusted R2 was used for comparison of the marginal and conditional effect of the explanatory variables and the individual fractions.

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download.