Introducing winter rice cropping by using non-saline tidal water influx in western basins of South 24 Parganas, India

A population exceeding 3.8 million people in the western region of 24-Parganas (South) is directly or indirectly reliant on agriculture as their primary source of livelihood. The agricultural trend shows a clear lack of multi-cropping with a drop of nearly 30% in rice cultivation during the winter season. Nearly 50% of the region is directly dependent on canals. The introduction of tidal water in the canal network provides an exceptionally economical and highly effective mode of irrigation water supply. The primary aim of the study was to identify the cartographic characteristics and channel hydraulics in the summer season. It was noted that the canals have a wide discharge range of 0.03–540.03 m3/s, average evaporation loss of 9.07 mm/day with a seepage loss ranging from 0.04 to 6.36 m3/s. The tidal water ingress quantity was calculated to be 4.17 Mm3, 5.32 Mm3, 1.88 Mm3 at Diamond Harbour sluice (Sl.), Kulpi Sl. and Kholakhali Sl. respectively. It was denoted that the augmentation of tidal backwater six times monthly would suffice the winter crop water requirement for the majority of the basins. This would result in the production of 172.13 kt which was previously 17.6 kt resulting in an increase of production by 878.01%. The per capita income would also be increased by nearly 978% for the season, resulting in the macro-socioeconomic upliftment of the region.

The South 24 Parganas (S24 Pgns.) district of West Bengal, India has been found to have a record of an intense scarcity of rainfall (42 mm-11 mm-11 mm-24 mm) during the months November-December-January-February respectively 1 . This results in the need for irrigation water within the region for this extended period for effective agriculture and other relatable activities. Mono crop is the existing main crop pattern in the region due to these varied and seasonal rainfall activities throughout the year. Multi-cropping could enhance the productivity of the region. In the region, nearly 50% of the agricultural produce is directly dependent on canals. In the region, the tidal backwater is not only used for multi-cropping irrigation purposes and pisciculture but also to serve as a resource base for various other activities and opportunities like fodder production, brick manufacturing, and other microscale industries which directly and indirectly affect the socio-economic conditions of the occupants within the region. It also provides an effective way to utilize the excess freshwater which would otherwise be wasted into the Bay of Bengal if not managed properly. The global demand for water, energy, and food is increasing the stress on the resources which are utilized for supplying these basic commodities 2 . This system absolutely nullifies the requirement of any extra resources, further accentuating the supply of irrigation water and food supply proliferating the efficiency of the clean technology. Just a few studies give insight about different sorts of back-water inflow estimation done over the past decades. The backwater function was determined with the assistance of numerical integration technique 3 . The backwater surface profiles of rectangular channels at gradually varying flow conditions had also been described 4 . The estimation technique for the water level, when the stream is under steady change has also been introduced by using an expanded form of the Bernoulli's equation was applied 5 . A computational program has also been likewise discussed for composite backwater profiles in trapezoidal channels 6 . An arrangement for the augmentation of irrigation systems through an abundance of tidal back-water convergence for the region has been previously provided 7 .
The examination endeavored to depict not just the improvement of the financial return of the general population of the Magrahat basin yet additionally to make arrangements for a multi-cropping pattern of farming/ cultivation including pisciculture options in the dry season for the sustenance of their employment. They further conjectured in light of the data and field studies that the Diamond Harbor Creek khal is an estuarine canal with a regular tidal influx. Particularly amid flash floods, the catchments used to be mostly flooded. To combat this issue, in the early 1960s, a strategy was set up by the Irrigation Department of the Government of West Bengal to regulate the passage of tidal ingress through a major sluice gate (Diamond Harbour sluice) and development of similar sluice gates in each estuarine canals (Kholakhali sluice, Hara sluice, Kulpi sluice) to check the surge issue in the territory as likewise to reserve the tidal backwater for use in lean periods. A proposal for the Magrahat basin irrigation management based on mathematical estimation had been examined 8 . The work was done purely on theoretical data collected from several sources and calculations were performed henceforth. Smaller reservoirs and branch canals are not considered here for the research as it would be very time consuming and difficult to reach certain sections physically. So, there might be a small difference (almost negligible) in the calculated results related to the water availability in the study area. The farmlands in the region do not usually directly receive water from the canals. There are several entrance points cut through the canal linings to receive water from the canals by the effect of tidal influx from the canals. It helps in a controlled yet ample water supply for cultivation. The water supply from the canals is then distributed via smaller channels or "nullah(s)" to the cultivation fields. The water is often even diverted to and stored in temporary storage ponds for future use when the canals do not receive enough water. The temporary storage ponds also are used for pisciculture, adding to the productivity of the reservoirs.
On several occasions whilst conducting the extensive field survey the people were noted to be curious and at several acquaintances enquired if the canals would be dredged or would be cleaned. It was uncanny how the population even being so dependent on the canal network was actively destroying its very backbone by dumping waste, sewage, and encroaching the very network that is so important for them. The lack of awareness and knowledge about the means to protect their asset is the most likely reason for such. The possibility to irrigate the total cultivatable command area of 523.18 km 2 was examined with the sessional feeding of the tidal backwater for the period of water requirement i.e. from ground preparation to the last watering. The study would create provision for increasing the winter rice cultivation and also greatly supplement in enhancing the economic return of the population in the region.

Study area
This study develops on a much broader area, covering five basins in the region the Magrahat basin, Keorapukur basin, Kholakhali basin, Hara-Hatuganj basin, and Kulpi basin (57 km 2 ) of the S24 Pgns. play a major role in imparting freshwater via a system of integrated channel networks within the catchment for effective irrigation purposes. The Diamond Harbour sluice is one of the primary sources in the Magrahat basin that provides tidal backwater to the Diamond Harbour Creek which in turn provides water for the Dasani Khal, Nazra Khal, Upper Hatuganj Khal, Sangrampur Khal, Upper Dhanpota Khal, Joynagar Khal, Kata Khal, Usti Nainan Outfall Channel, Suryapur Inner Channel, Lower Keorapukur Khal, Hotor Khal, and Lower Srichanda Khal.
The Kholakhali Khal which receives water from the Kholakhali sluice is the primary channel in the Kholakhali basin. The Hara Khal along with Mid Hatuganj Khal and Mid Dhanpota Khal receiving water from the Hara sluice plays the principal role in the tidal backwater irrigation process of Hara-Hatuganj basin. The Kulpi basin is primarily irrigated with water received from Kulpi sluice gate entering the Kulpi Khal, and a subsidiary water network formed by the Lower Hatuganj Khal and Lower Dhanpota Khal. Sustainable development is a subtle equilibrium between environmental protection, economic development, and societal upliftment which can be achieved only by efficient use of resources 9 . The study area was identified using Google earth software and the demarcations of the basins and canals in the study area were made for easy understanding of the region. The study area was subdivided into five basins that receive water directly from the Diamond Harbour sluice, Kholakhali sluice, Hara sluice, and the Kulpi sluice. The five basins in the S24 Pgns. where predominant tidal-backwater irrigation practices are being performed are mapped in Fig. 1.

Methodology
A detailed survey was done to perceive the previous research work in literature like journals, conference proceedings, book chapters, and several statistical data reports. It was observed unfortunately that scanty work was carried out on the topic and on the study area. Government records on soil exploration and Land Use-Land Cover (LULC) maps of the study area were collected from Survey of India, Canals Division, and Magrahat drainage division of West Bengal Irrigation and Waterways Division. Field surveys of the region were carried out to verify the positioning of the canals with the maps collected. Satellite maps were extracted from United-States-Geological-Survey (USGS) Earth Explorer and the LULC map was prepared using ArcGIS v10.4 software and the canal paths were denoted on the satellite imagery as shown in Fig. 1. For a basic and clear understanding of the basin, the LULC map was prepared using ArcGIS v10.4 software. LANDSAT 8 OLI (Operational Land Imager) dataset was utilized having a 30 m spatial resolution of the year 2019. The accuracy of the LULC map was assessed using the Kappa coefficient. The overall accuracy and the Kappa coefficient value were 0.983 and 0.950, respectively. Later it was also validated with Google Earth imagery. The tidal timings and the tidal bore schedule were collected from Kolkata Port Trust Authorities, Ministry of Shipping, Government of India (GOI). Subsequent field visits were carried out for the collection of data (from November to February) to understand the characteristics of the channels during the tide and ebb effect. The cross-sectional (C/S) area of the channels was measured using echo-sounder, model GPSMAP 585, manufactured by Garmin, USA for before tide (BT) and after tide (AT) conditions. The readings were plotted in AutoCAD 2012 and the C/S area was calculated. Float method and piezometric depth variance measurement were carried out to find out the tidal velocity of the canals. The length of each canal was measured using the length scale conversion method on the corrected satellite imagery. The Area-velocity method was applied to find out the discharge of the canals at different points of testing. The actual discharge flowing through the canals differ from the calculated discharge as there is a prevalent www.nature.com/scientificreports/ conveyance loss in the canals. The quantity of available water can only be estimated after making deductions for conveyance loss. The conveyance loss was calculated as a summation of the seepage loss and evaporation loss as the operational loss could be considered negligible in this case. Irrigation systems are designed to perform at maximum efficiency conditions. It implies that the exchange of water at the base expense and with the least water loss 10 . The irrigated regions are extensive and henceforth require the transfer of water for irrigation purposes to travel longer distances. It leads to losses amid the transfer procedure. The canal lining is usually treated to limit such losses. However, the lining gradually deteriorates with the progression of time where cracks emerge that in the end lead to the collapse in the layers of lining. These cracks and fissures of covering joints are exit routes for water so necessary measures should be taken on time to diminish these losses 11 . The Moritz equation (1967) was converted to a metric scale for this work. Seepage rates were estimated using Eq. (1) 12 .
where S 1 is the seepage losses (m 3 .day -1 .km -1 ); Q is the canal discharge (m 3 /s); V is the mean velocity of the canal (m/s); C is the coefficient directly varying on the function of the material (soil type) used to construct the canal lining. The modified A. N. Kostiakov equation (1978) was also used for re-estimating the seepage losses from earthen canals (Eq. 2) 13 .
where; S 2 is the seepage losses (m 3 /s); L 1 is the length of the canal (km); s is the water losses per km of canal length (%) was calculated using Eq. (3), where A and m are empirical constants depending on soil permeability.
Soil samples were collected from the canals for testing the soil type of the canals. The grading of the soil samples was performed by the sieve analysis method. Soil samples collected from the canals were tested using the variable head permeability test to determine the permeability of the underlying soil. Variable head permeameter was primarily used to measure the permeability of relatively less pervious soil. The coefficient of permeability, K (cm/s) was calculated using Eq. (4), where h 1 is the initial head (cm); h 2 is the final head (cm); t is the time interval (min); a is the C/S area of the standpipe (cm 2 ); A is the C/S area of the specimen (cm 2 ); L 2 is the length of specimen. (1) Scientific Reports | (2021) 11:553 | https://doi.org/10.1038/s41598-020-80797-x www.nature.com/scientificreports/ Apart from seepage loss, the canals are also subjected to evaporation loss. Evaporation loss is dependent on several climatological parameters. Considering the several climatological parameters the evaporation loss has been calculated using Meyer's equation (5).
where E L is the evaporation (mm/day); K M is the coefficient with a value of 0.36 for large deep waters and 0.50 for small, shallow waters; e w is the saturated vapor pressure at the water surface (mm of mercury); e a is the actual vapor pressure of overlying air at specified height (mm of mercury); H r is the relative humidity (Eq. 6); u 9 is the monthly mean wind velocity at 9 m above ground level (GL) (km/h) In this scenario, the wind velocity was collected from the Indian Meteorological Department (IMD), GOI reports, which were collected at a height of 1 m above ground level, so it had to be converted to find the wind velocity at 9 m above GL. The wind velocity was calculated using the 1/7 power law (Eq. 7).
where and u 1 is the wind velocity at 1 m above GL (km/h).
The total available water volumes at ebb and tide conditions were calculated from field data. The total basin area, gross cultivable area (GCA), and cultivable command area (CCA) were calculated from the LULC data. The climatic and meteorological data were extracted from CLIMWAT 2.0 software for Sagar Island recording station as it was the nearest to the climatic conditions of the region. The data were then fed into the CROPWAT 8.0 software. The CROPWAT software uses rainfall data to indicate the effective rainfall for the growth of crops in the region. The effective rainfall was calculated in the software by the United States Department of Agriculture Soil Conservation Service Method (USDA-SCS). Effective rainfall is required as not all the quantity of rainfall that the catchment area receives can be used due to surface runoff and deep percolation below the root zone. Other parameters as the max-min temperature, wind speed, relative humidity, and average sun hours were also provided for the region and were also input into the software for consideration. The crop characteristics were taken from the CROPWAT database. The two primary factors which CROPWAT soil template considers are the available water storage capacity of the soil and the daily final capacity of rainfall infiltration. The available water storage capacity is the same for both black soil and clayey loam soil i.e. 200 mm/m 14 . The daily final capacity of rainfall infiltration for black soil is 30 mm/day and for clayey loam is 30.48 mm/day, which can also be considered as identical 15 . Another secondary and notable consideration in the template is the maximum root depth. The template has an allowance for 900 cm whilst the maximum root depth for rice ranges from 55 to 65 cm, so abundant allowance has been considered on using the template 16,17 . Henceforth, the black clay soil template was used in this part of the research. Now, from the water availability and crop requirements, the deficiency was calculated to note the ideal number of times the sluice gates need to open to allow the tidal influx to occur in the region to increase the production rate of Winter Rice in the region. It is to be noted that the available water was calculated from the water entering the canal network from Hooghly River at the definitive sluice gates that are opened to allow the water to enter during the high tide and is closed after the peak period of the tide. The period of the high tide ranges from a period of 4.5-5.5 h so an average of 5 h has been considered as the mean to calculate the average tidal water influx into the basins. The tidal water entering was then reproduced by two, four, and six times to calculate the average monthly water availability for a specific month and basin to find out the water availability in the event of tidal water ingress. Propositions have been made according to the results for multicropping according to net-water availability. A concise flowchart of the methodology has been shown in Fig. 2.

Results
The LULC map of the region showed a difference in data from actual field data. Henceforth, ground-truthing was done as an additional reference and to prepare the adjusted LULC map, and data is depicted in Fig. 3 and Table 1 respectively. The field-collected data for canal length, the cross-sectional area before tide (C/S area BT), the cross-sectional area after tide (C/S area AT) and calculated discharge has been shown in Table 2.
It was noticed during the field survey that the canals were unlined non-prismatic channels. The results as depicted in Table 3 show that the canal substrate is primarily clayey or loamy in nature. Moritz suggested the value of the co-efficient for clay and clayey loam as 0.41. Henceforth, following Eq. (1), the seepage value of the canals was calculated and the results are depicted in Table 3.
The values of h 1 , h 2 , t are variables for every soil type and the tested results and calculations are shown in Table 4. The values of a (0.196 cm), A (78.53 cm), and L (13 cm) are constant for all the calculations, in this case. The permeability (P) of the canal linings shows that the canals have a low permeability rate. Henceforth, the value of constants 'A' is 0.7 and 'm' is 0.3. The seepage rates of each canal according to modified Kostiakov's equation have been depicted in Table 4. The seepage values retrieved using the Moritz equation and modified Kostiakov's equation although quite near but the values differ to some extent. Henceforth, for the calculation of losses, the average of the values was taken into account. This would provide a singular value, reducing confusion and ease further calculations. The mean seepage value has been calculated in Table 5. Tables 6 and 7 show the canal surface area with temperature and the relative humidity with wind velocity respectively. Now according to Eqs. (6), (7), and (5) the values e a and u 9 and E L were calculated to be 7.02 mm of Hg, 11.58 kmph, and 9.07 mm/day respectively. The mean seepage loss and the evaporation loss were deducted from the mean discharge to calculate the actual discharge for the canals, as shown in Table 8. The discharge of the  Since the research concentrates on maximizing the production of the winter rice in the region, the CCA was considered as the total area where rice is currently cultivated. The observed trend shows that whilst for monsoon rice 38.27% of the CCA is used but for winter rice only 8.249% of the CCA is used. This is a drop of nearly 30% in the growing of rice in the winter season.
As mentioned earlier, the canals receive the bulk of their water from tidal ingress. Hara-Hatuganj basin and Kulpi basin were considered co-jointly as it was noted during the field visits and also in satellite imagery that the canal network for the two basins was so intensely interwoven that differentiating the basins would provide an error as the likelihood of water distributing itself easily is very high. As the tidal velocity of the points were previously measured in the field visits. The tidal discharge was calculated as shown in Table 9. The discharge into the Keorapukur basin is from the water entering from the Magrahat basin and is thus lost to the Magrahat basin. So, the water entering the Keorapukur basin has been deducted from the water calculated to be available for the Magrahat basin. The high tide is noticed to last an average of 5 h per event, for an average of two events per month. The total monthly availability was calculated henceforth as depicted in Table 10.
This was then deducted from the average monthly water requirement of the crops which had been calculated earlier in Table 10 to find out the deficiency (if any) and the deficiency percentage (if any) for the basins for that particular month. It was also noted to take into account that the efficiency of the system remains high, i.e. the number of times the gates need to be opened be kept minimum so as to also reduce the operational cost. Furthermore, excess water entering the canal network could also result in potential flooding of the area of the basins adjoining the riverbank, as was noted earlier. The details of the calculations are shown in Table 11. The results show that the Keorapukur basin has a deficiency of more than 98% even after the tidal water influx happening more than six times, which is a major lag for the region. It is suggested that the basin opts for other means of crop cultivation. It is to be also noted that the recent advancements of the socio-economic and geo-political trend of the region are pointing towards a drastic change in the livelihood of the people. The beneficiary (i.e. inhabitants of the region) are changing to a more suburban kind of livelihood due to recent progress of the Kolkata city establishments distributing further in those regions. The shifting of several government administrative units from the central region of the Kolkata city to the regions and adjoining indicate a massive turn in the dependence of livelihood of the people in the region. The trends suggest that the ever-expanding Kolkata city is most likely to grasp the region into an urban-suburban kind of region where agriculture will not be a major key player, and the canal network might serve only as a part of an efficient drainage system. For the remaining other basins after careful analysis of the calculated data, the following plan of action is being suggested so as to get the maximum productivity in Table 12.
It is to be noted that even after allowing tidal water influx six times in November at the Diamond Harbour sluice there was a deficiency of 53.21% or 36.3 m 3 winter irrigation water. This could be easily avoided by channeling water from other basins in the network to enter the basin and also by supplementing water from the Piyali River in the east of the basin during November to fulfill that need. It was noted from the population data recorded in Census 2011 that the rural population has an average family with one earning member per family has a size of 5-12 members. The proper irrigation management of the study area would directly benefit nearly 3.8 million people. Since the current yield rate of winter rice in the region is 3290 kg/ha, the inclusion of the plan would result in the production of 172.13 kt which was previously only 17.6 kt. The increased area of cultivation could likewise provide an opportunity to employ more than 4.6 million new farmers, which will greatly supplement the agricultural sector not just within the region but also the whole nation. The average income of a farmer per harvest was nearly INR 892, which greatly dissuaded farmers in winter rice harvesting. The current plan could increase the income to INR 8724 per harvest per farmer, which is nearly ten times or a 978% uptrend.

Discussions
Satellite imagery models are helpful tools to comprehend the conduct of hydrological systems. Quite unfortunately, the parameterization of these hydrological models requires a wide range of information, from various sources, and from various orders (e.g., meteorological, geo-morphological). In a basin-scale hydrological demonstration, the conventional technique for model initiates with getting contour map, DEM, LULC maps, soils maps, and meteorological information. Usually, the specialist's past involvement with these datasets determines if they would be best suited for utilization in a specific scenario 9 . Additional multifaceted complexity is that different scientific networks have varying information sources and control strategies, which makes the interpretation of    www.nature.com/scientificreports/ these data by researchers unfamiliar to the basin physically, exceedingly unreliable. Whilst preparing the LULC with satellite imagery it was noticed that the imagery could be misleading on several occasions. It was noticed that the actual field-collected data for LULC is much more efficient and accurate than satellite imagery. Usually, the researcher's involvement with a region's LULC determines the suitability for utilization in a specific scenario. Frequently more updated, or progressively reasonable information items do not exist. In the current scenario familiarity with the actual field conditions helped in creating a fairly accurate recognition of the basin characteristics. It was also noted during the survey that the canals in various places were subject to encroachment, heavy silting, and wastes (plastic, domestic garbage, slaughterhouse waste) that were dumped into the canals which resulted in complete impediment of canal hydraulics in certain regions of the drainage basin.

Conclusions
A deep insight into the irrigation requirements, water availability, socio-economic statistics, geographical conditions, the agricultural scenario, and meteorological parameters provides a notion that though the water entering the canals as a result of tidal backwater influx might theoretically seem to suffice the irrigation requirements for the winter rice crop, the intensification and high-density production in the region is not prevalent mostly due   www.nature.com/scientificreports/ to the lack of good distribution. An improved canal network would provide a much more increase in the overall productivity of crops in the region. The impact of encroachment, dumping of waste along the lack of maintenance makes the canal network less efficient. Also, control over the tidal backwater returning into the river during the low tide results in a huge loss of water which could be controlled by some minor adjustments and tweaks in the system operation. The model increased the winter cropping great extent and likewise the farmers' income by nearly ten times. Nevertheless, apart from cultivation, the excess water could be a major boost to several other factors like small scale industries and aquaculture, which will greatly improve the socio-economic conditions of an entire region. The influence of tidal backwater feeding in the non-monsoon period not only allows for multi-cropping and aquaculture but also acts as an active waste and flood water drainage network for the region. The maximization of sustainable use of natural resources is the key to the economic development of any region. The study provides an extensive outlook on the water availability, drainage, and irrigation water supply characteristics of the region but the model can be easily normalized by replicating the model framework whilst changing the specific variable parameters unique to any location to effectively derive results.