Spatiotemporal trends of urban-induced land use and land cover change and implications on catchment surface imperviousness

Urbanization, changes in land use and land cover (LULC), and an increase in population collectively have significant impacts on urban catchments. However, a vast majority of LULC studies have been conducted using readily available satellite imagery, which often presents limitations due to its coarse spatial resolution. Such imagery fails to accurately depict the surface characteristics and diverse spectrum of LULC classifications contained within a single pixel. This study focused on the highly urbanized Dry Creek catchment in Adelaide, South Australia and aimed to determine the impact of urbanization on spatiotemporal changes in LULC and its implications for the land surface condition of the catchment. Very high spatial resolution imagery was utilized to examine changes in LULC over the past four decades. Support Vector Machine-learning-based image classification was utilized to classify and identify the changes in LULC over the study area. The classification accuracy showed strong agreement, with a kappa value greater than 0.8. The findings of this analysis showed that extensive urban development, which expanded the built-up area by 34 km2, were responsible for the decline in grass cover by 43.1 km2 over the last 40 years (1979–2019). Moreover, built-up areas, plantation, and water features, in contrast to grass cover, have demonstrated an increasing trend during the study period. The overall urban expansion over the study period was 136.6%. Urbanization intensified impervious area coverage, increasing the runoff coefficient, equivalent impervious area, and curve number by 60.6%, 60.6%, and 7.9%, respectively, while decreasing the retention capacity by 38.6%. These modifications suggest a potential variability in catchment surface runoff, prompting the need for further research to understand the surface runoff changes brought by the changes in LULC resulting from urbanization. The findings of this study can be used for land use planning and flood management.


Introduction
Population growth and increased migration to urban areas have contributed to have increased urbanization affecting the LULC dynamics of an area (Aburas et al. 2019;Gashu and Gebre-Egziabher 2018;Shawul and Chakma 2019), consequently influencing the quantity of runoff and sediment transport (Bledsoe 2002).Urbanization is a worldwide phenomenon characterized by increasing built-up features and the conversion of natural landscapes for residential, commercial, and industrial land use activities (Li et al. 2018;Pradhan 2017).Urban sprawl transforms non-urban land into built-up areas (roads, rooftops, parking lots, and other urban surfaces), which significantly increases the proportion of impervious surfaces (Ding et al. 2022;Li et al. 2018) causing considerable flooding and waterlogging in urban catchments (Moniruzzaman et al. 2020).This results in decreased infiltration, higher surface runoff, and flooding in downstream areas (Feng et al. 2021).
The analysis of LULC changes is necessary to determine the extent and magnitude of urban expansion (Grigorescu et al. 2021).Studying changes in LULC is essential for evaluating its environmental impact in terms of both spatial and temporal dimensions (Gashu and Gebre-Egziabher 2018;Shawul and Chakma 2019), including the influence on river discharge and sediment yield.Remote sensing plays an important role in detecting LULC changes (Jamali et al. 223 Page 2 of 16 2022; Jat et al. 2008;Rojas et al. 2013;Sidhu et al. 2018) and their impact on runoff (Demeke and Andualem 2018), sediment yield, and channel morphology.IDRISI Selva (IDRISI), Sentinel Application Platform (SNAP), Earth Resources Data Analysis System (ERDAS Imagine), Environment for Visualizing Images (ENVI), ArcGIS Pro and Google Earth Engine are the most common remote sensing tools used for LULC classification (Jamali et al. 2022;Jat et al. 2008;Sidhu et al. 2018).
Several image classification methods, including the ISO cluster unsupervised classification, Maximum Likelihood (ML), random tree, Support Vector Machine (SVM), and K-nearest neighbor, can be used to produce a LULC map and assess the impact of urbanization (Chandra and Bedi 2021).Unlike unsupervised classification, supervised classification allows the user to select image pixel samples to classify data into groups, such as special subcategories, by providing training samples to cluster new inputs into predefined groups, map features according to their respective groups, and to provide high accuracy (Hu and Dong 2018;Nguyen 2020).The most frequent approach in land use classification and change detection is supervised image classification, which is a user-guided classification utilizing training samples.The most successful technique for detecting changes in land use has been identified as SVM (Kesikoglu et al. 2019).In comparison with other classification techniques, recent studies have demonstrated that supervised image classification using an SVM classifier has the highest accuracy in LULC analysis (Basheer et al. 2022;Dabija et al. 2021;Thanh Noi and Kappas 2017;Xie and Niculescu 2021).The SVM classifier uses statistical learning theory and is an advanced and powerful machine learning classification technique (Chandra and Bedi 2021;Mohammed and Sulaiman 2012).It can handle segmented raster data and common image inputs with less noise susceptibility (Ding et al. 2022;Lamine et al. 2018).In addition, this approach makes use of object-based classification, which involves image segmentation using spatial, spectral, and size information.This approach more closely mimics real-world features and yields more accurate categorization results for high-resolution imageries (Blaschke et al. 2014;Li and Shao 2014;Peña et al. 2014).
The accuracy of LULC change detection studies are dependent on spatial scale or map resolution (Manandhar et al. 2010).LULC change analysis has been conducted in several studies by employing different satellite imagery resolutions.Some studies have utilized a coarse spatial resolution of 250-m from a Moderate Resolution Imaging Spectroradiometer (MODIS) sensor (Yin et al. 2014), and 60 m (Ding et al. 2022); while, others have worked with moderate resolutions of 15 m (Mansour et al. 2020), and 30 m (Hu and Dong 2018;Hussien et al. 2022;Malede et al. 2022;Silva et al. 2018).The spatial resolution of imagery affects the accuracy of classification and representativeness of the resulting LULC map for the intended purposes (Chen et al. 2004;Fisher et al. 2018).The provision of high-resolution spatial imagery of urban areas (Chen et al. 2004;Zhen et al. 2013) provides an opportunity to clearly identify the pervious and impervious surface conditions of the study area.Unlike coarse resolution pixels (such as 30 × 30 m), finer resolution pixels (2.5 m and finer) can be utilized to distinguish road, roof, grass, and tree characteristics that can be present within a 30 × 30 m ground spatial position.Satellite sensors such as Landsat, Sentinel, and/or ASTER have high spectral resolution and multispectral bands such as Red Edge (RE), Near Infrared (NIR), Shortwave Infrared (SWIR), and Thermal Infrared (TIR), which are effective for moisture/water extraction and vegetation health (Amani et al. 2018).Compared to high spectral resolution, their spatial resolution is inadequate for distinguishing various land use aspects that depict surface imperviousness.Therefore, this study employed high-resolution SPOT5 pan-sharpened satellite and aerial imagery to quantify the spatial and temporal extent of LULC change and its implications on the surface imperviousness of a catchment.
It is crucial to manage and mitigate the ecological and geomorphological effects of urban development by comprehending the effects of urbanization on fluvial systems and sedimentation through the analysis of spatiotemporal changes in urbanization (Chin et al. 2013;Ciupa & Suligowski 2020;O'Driscoll et al. 2010).Surface imperviousness is comparatively high in established urban catchments and has been associated with higher runoff and lower levels of total suspended sediment discharge (Myers et al. 2021).Therefore, in recent years, determining how urbanization the affects hydrology of the catchment and channel stability has gained significant importance (O'Driscoll et al. 2010).
The Dry creek catchment is an urbanized catchment in the metropolitan area of Adelaide, South Australia, with significant land use changes due to urban sprawl since the 1970s (Wilkinson 2005).These changes have resulted in alterations to the surface imperviousness and hydrological system of the catchment.A recent modeling study focused on the Parafield stormwater harvesting scheme, a part of the Dry Creek catchment, has revealed that urban development within the catchment could potentially result in a 20% increase in impervious area (Clark et al. 2015).However, to date, no research has been undertaken within the catchment to explore the implications of urban-induced LULC changes and their influence on the surface imperviousness of the area.Therefore, this study aimed to examine the long-term spatiotemporal trends of urban-induced LULC change and its implication on surface imperviousness over the urbanized Dry Creek catchment between 1979 and 2019 using very high-resolution multispectral imagery.

Study area description
This study focused on the Dry Creek catchment, which is located in the Northern part of Adelaide, South Australia.The catchment is located between 138°35′0″ and 138°46′0″ E longitude and 34°51′18″ and 34°45′28″ S latitude (Fig. 1).The study area covers about 115.4 km 2 .The inner norther suburbs of Adelaide, which encompass the Dry Creek catchment, were projected to reach a population of 315,367 by 2021 (Department of Planning Transport and Infrastructure 2019a).The topography of the catchment ranges between 0 and 424.81 m AHD (Australian Height Datum).The annual rainfall of the catchment ranges between 450 and 700 mm (The Goyder Institute for Water Research 2016).
In contrast to the upper watershed, which is primarily rural, the lower catchment is primarily dense residential, commercial, and industrial, and contains some salt fields scheduled for development and various constructed wetlands.According to the Department of Planning, Transport and Infrastructure (DPTI) (Department of Planning Transport and Infrastructure 2017), the 30-year plan for Greater Adelaide showed that the Dry Creek catchment area has been marked as a potential urban growth area.
The increase in urban growth rises the impervious surface area with a potential of flooding to the downstream environment.

Data source
Very high-resolution aerial imagery (1979( and 2019( ) and SPOT imagery (2006) ) were obtained from the South Australian Department for Environment and Water (DEW) and used to identify the various LULC types found throughout the catchment (Table 1).The study period was chosen based on the availability of very high-resolution imagery encompassing the entire catchment, as well as the period when the population began to grow, potentially resulting in urban expansion.The higher resolution (0.075 m and 0.25 m) imagery was resampled using the nearest neighbor resampling technique to a lower resolution (2.5 m) prior to undertaking the classification to make change detection simple.The data were pre-processed (resampling and segmentation) and used in the ArcGIS Pro 3.0 environment to create LULC maps.
Although raster imagery with smaller cell sizes results in a larger file size and takes longer to process, its importance in LULC classification cannot be overstated.This is due to its higher spatial resolution, which provides finer details of land surface features.The hydrologic soil group of the area was retrieved from the global hydrologic soil group and resampled to 2.5 m resolution using the nearest neighbor resampling technique to estimate the curve number  et al. 2018).The nearest neighbor resampling technique was selected since it is appropriate for categorical data and does not alter the value of the input cells.

Image classification
The number of LULC classes was determined by the current condition of the catchment by visual inspection of the ground information and the imagery, the spatial and spectral resolution of the imagery, the availability of various land use practices, and the intended purpose of the classification output or the fineness of the classification.Preliminary investigations were conducted by visiting the catchment and viewing high-resolution images on Google Earth Pro.Typical land use classes in urban settings include built-up areas, vegetation, water bodies and/or wetlands, green spaces, open land, and cropland (Gashu and Gebre-Egziabher 2018;Moniruzzaman et al. 2020).By visual inspection of high-resolution imagery, four land use land cover classes (built-up, grass, water, and plantation/trees) were identified as present within the Dry Creek catchment and were thus chosen for the classification (Table 2).
The image segmentation procedure, which combines neighboring pixels that are similar in color and have geometric properties, is commonly utilized during image classification using SVM.Three criteria (spectral detail, spatial detail, and minimum segment size in pixels) govern the separation of images into objects in image segmentation.A classification schema is necessary to establish the number and type of classes to be utilized for supervised classification.In the ArcGIS Pro 3.0 environment, a schema with four LULC classes was created and saved in an Esri classification schema file (.ecs), which employs JSON syntax (https:// pro.arcgis.com/ en/ pro-app/ latest/ help/ analy sis/ image-analy st/ the-image-class ifica tion-wizard.htm accessed on December 15, 2022).User-defined training samples were prepared at randomly distributed locations (Fig. 2

Classification accuracy
The assessment of image classification accuracy helps to evaluate the precision of classification results, ensuring the reliability of subsequent analysis and applications (Foody 2008;Zhang et al. 2020).To verify the accuracy of classification, validation polygons were selected through a stratified random sampling approach, based on clearly distinguishable features observed in both aerial imagery and SPOT imagery.Five hundred randomly generated points were derived from the sampled polygon (Table 3 and Fig. 3).These points were then compared with the classified image to assess the accuracy of the classification.After evaluating the classification accuracy performance indicators, a post-classification carried out, resulting an improved classification accuracy.To assess the accuracy of classifications, a confusion matrix was generated and used to calculate the key parameters such as overall accuracy The individual class accuracy was assessed using the two indices PA and UA, which were computed using Eqs. 1 and 2, respectively (Congalton and Green 2019).The former index, PA reflects the likelihood of a reference pixel being correctly classified, while the later index, UA indicates the likelihood that a pixel classified on the map depicts that category on the ground (Congalton 1991).
The OA that provides a comprehensive measure of how accurately the model predicts the correct class across all classes, is computed using Eq. 3 (Liu et al. 2007).The kappa coefficient evaluates how accurately the classification performed as compared to just randomly assigning values.The kappa coefficient measures the agreement between classification and truth values (Table 4) and has a value between 0 and 1, with 1 being perfect agreement and 0 representing no agreement.The kappa analysis is a discrete multivariate approach for testing classification accuracy or agreement between classified and reference data points (Rwanga and Ndambuki 2017).

Change detection
Post-classification comparison was used to detect changes in LULC at three different time periods.Change detection can be either a categorical change, pixel value change, or time series change.Categorical change computes the change between two thematic rasters and quantifies and explores the pixels that are changed or remained unchanged.The pixel value change computes the absolute or relative change between the two continuous raster images.Time series change uses LandTrendr to analyze changes in a multidimensional raster and compute the date on which each pixel changed.LandTrendr uses regression-and vertex-to-vertex-based trajectory fitting to detect abrupt events, such as disturbances, as well as longer duration processes, such as regrowth (Kennedy et al. 2010).Change detection was performed using a comparative analysis of spectral classifications produced independently at times t1 and t2 (Hussien et al. 2022).The percentages of LULC change detection (1)  and the rate of change were estimated using Eqs. 5 and 6, respectively.
where Δ is percent change, R c is rate of change, R ue is rate of urban expansion, N is the number of years between the two LULC maps, CA n-1 and CA n are the initial and final LULC class area coverage (in km 2 ), BA n-1 and BA n, are the built-up area at time t n-1 and t n (in km 2 ).The rate of urban sprawl/expansion between the three study periods was also determined using Eq. 7.

Impact of urbanization on surface imperviousness
Surface imperviousness is altered by urbanization; imperviousness can be evaluated in hydrological studies using the runoff coefficient, curve number, and retention capacity.In this study, built-up areas, such as buildings, parking lots, and roadways, were considered as impervious surfaces.The curve number represents the stormwater runoff capacity of a drainage basin.It is an empirical measure used in hydrology to anticipate direct runoff or infiltration from excess rainfall (Cronshey et al. 1985), and is estimated as a function of land use, soil type, and antecedent catchment moisture using the Soil Conservation Service (SCS) Technical Release 55 (TR-55) tables (SCS 1986).The ArcGIS Pro 3.0 environment was used to combine the LULC and hydrologic soil groups and apply the corresponding curve number value using the SCS TR-55 table.Equation 8 was used to calculate the composite curve number value for each year.The potential maximum retention is a measure of a watershed's ability to extract and absorb storm precipitation and is calculated using Eq. 9 from the curve number.The runoff coefficient, which quantifies the percentage of rainwater flowing out of a certain surface during  a storm event, was assigned for each year and LULC class.Equation 10 was used to generate the weighted runoff coefficient and Eq.11 for estimating the equivalent impervious area that are used to compute runoff utilizing the rational technique for individual LULC class. (8) where CN composite is the composite curve number used for runoff volume computations; C w is the weighted runoff coefficient; i is an index for the catchment subdivisions of uniform land use and soil type; CN i is the curve number for subdivision i; C i is runoff coefficient value for each subdivision; EIA is the equivalent impervious area (km 2 ), and A i is the drainage area of subdivision i (km 2 ).

Accuracy of LULC classification
The overall accuracy and kappa coefficient obtained from the confusion matrix agreed consistently demonstrated strong

Spatiotemporal land use and land cover maps
The classification results revealed that in 1979, the catchment was dominated by grass cover occupying an area of 67.65 km 2 followed by built-up and plantation with 24.9 km 2 and 22.52 km 2 areal coverage, respectively ( From the spatial map of the LULC maps, it was observed that grass was dominantly found in the eastern and central parts of the catchment in 1979.Most plantation/tree LULC classes have been found in the eastern (upper) parts of the catchment over the last four decades (Fig. 4).

Spatiotemporal changes in LULC classes
The LULC maps for 1979, 2006, and 2019 indicated a substantial trend of LULC change over the last four decades.Notably, the built-up areas expanded by 33.98 km 2 (29.46%) while the grass cover diminished by 43.09 km 2 (33.73%) between 1979 and 2019.During this period, the rate of change of grass cover changes closely corresponded to the rate of the built-up area change, with a decrement of  9 and Fig. 5).
In comparison with the plantation and water classes, the built-up area increased extensively.Regardless, grass cover has consistently declined over all time periods.The water class increased due to the construction of Mawson Lakes and other reservoirs in some parts of the catchment and pond water in the upstream parts of the catchment in mining holes.Between 2006 and 2019, the rate of change in built-up area was 1.03 km 2 /yr; whereas, the rate of change in grass cover was −1.29 km 2 /yr.
Plantation coverage increased due to the implementation of Urban Forest Biodiversity Program (UFBP, 2002) and urban greening, which includes private greening, streetscapes and transportation corridors, riparian vegetation or green and blue corridors, and urban parks.Riparian landscapes (lands near rivers and streams) help to stabilize stream channel morphology, protect streams from upland pollution sources, and divert sediment-producing activities away from the stream (McKergow et al. 2003).
Post-classification categorical LULC change detection was carried out in ArcGIS Pro 3.0.A significant transition from grass to built-up areas was noted due to increased population and urbanization, with strong demand for residential, commercial, and industrial buildings with car parks and roadways.There was also a visible transition from built-up area to grass cover.This was because places under development were entirely regarded as built-up areas and later developed with backyards and some trees, particularly in residential districts.
The spatial changes in LULC and urbanization showed that the central, northeastern, and southern parts of the catchment significantly changed from grass cover to builtup areas in the three consecutive study periods (Fig. 6).This resulted in an overall increase in impervious area throughout the catchment.The plantation/trees showed a slight increase due to the conservation practices of river buffers and parks as well as an overall urban greening.The development of Mawson Lakes Pond in downstream areas and other ponds in upstream mining areas resulted in a minor increase in the spatial distribution of water features within the catchment.During the period between 1979 and 2006, approximately 55.7% of the catchment experienced distinct categorical changes in LULC distribution (Table 10).Furthermore, between 2006 and 2019, approximately 44.08% of the catchment experienced categorical shift in LULC (Table 11).Notably, over the entire study period , approximately 62.98% of the catchments showed categorical changes between the LULC classes (Table 12).A research conducted by Manandhar et al. (2010) in the Lower Hunter district of New South Wales also reported that approximately 28% of the study area underwent changes in land use and land cover classes.
Approximately 72.68 km 2 of the catchment area experienced categorical LULC changes between 1979 and 2019, with a maximum extent of 32.26 km 2 of grass cover   converted to built-up.During the study period, only 37.02% of the catchment area remained unchanged (Table 12).In addition to categorical change identification, works such as Manandhar et al. (2010) employ post-classification change detection based on a transition matrix constructed with twodimensional tables, with the rows of the matrix illustrating the initial map class and the columns representing the subsequent map categories.Jamei et al. (2022) reported similar findings in Melbourne, where population growth and urban expansion caused a notable change in LULC.
The study found that urban expansion across the entire study period was 34 km 2 , with an average annual change rate of 3.41% specifically within the built-up area.This observation indicated that urban development in the study area has been increasing (Table 13).In a similar study conducted by Ma and Xu (2010) in Guangzhou City, China, revealed a notably higher annual rate of change reaching 9.7% during the study period between 1995 and 2002.

Population growth and contributions to LULC
According to the Australian Bureau of Statistics (ABS), Australia's population has increased by over 25% from 2001 to 2016, reaching twenty-four million people.Currently, 90% of Australians live in cities, which was 60% in 1911.Notably, overseas migration currently has accounted for slightly over 55% of Australia's population growth since 2001, making it as the predominant contributor to the country's population growth (ABS 2019).The Australian Bureau Statistics data for South Australia showed that the state's population increased by 251,420 residents between 1979 and 2006, followed by an increase 226,284 individuals between 2006 and 2021 (Table 14).
The population statistics of inner north suburbs of Adelaide, in which Dry Creek is located, has shown consistent upward trajectory from the baseline census population in 2016 to the 2036 projection (Table 15).Between 2016 and 2021, there was a noticeable increase of 11,882 people.This rate of growth is anticipated to be maintained, with the population expected to increase by 34,128 people between 2021 and 2036 (Department of Planning Transport and Infrastructure 2019b).This projected increase in population will increase demand for the number of new residential houses and other institutional, commercial, and recreational built-up areas.The increase in population, has directly contributed to the increase in built-up area coverage in the Dry Creek catchment.This demonstrates how the need for residential buildings, parking lots, and other infrastructure has been raised by an expanding population.There was a strong correlation between the changes in LULC and population changes.While most LULC classes showed a strong positive correlation, grass cover demonstrated a strong negative correlation (Fig. 7).Although, the amount of data available is limited in performing an in-depth quantitative analysis of the historical trends and correlation between changes in LULC and population, the trends presented using the available data reveal that ongoing population changes will continue to influence the changes in LULC across the Dry Creek catchment.Ding et al. (2022) confirmed a similar result, with a growth in population having a substantial positive association with an increase in urban land.Similarly, Ouedraogo et al. (2010) investigated the significant relationship between the population and areas of cropland and open woodland in Burkina Faso.From the projected population, it was observed that the population will continue to increase impacting the current LULC, which further impacts the hydrology of the catchment.

Urbanization impact on catchment surface imperviousness
The surface imperviousness of a catchment is highly interrelated to the runoff coefficient and runoff generation capacity (Feng et al. 2021).In the Dry Creek catchment, significant urbanization has increased the built-up areas.Surface imperviousness is significantly correlated with built-up areas in the catchment.Built-up areas were classified as impermeable surfaces; whereas, water, grass, and plantation cover were classified as pervious.The analysis of the Dry Creek catchment revealed a marked 136.55% rise in imperviousness, with corresponding 37.56% decrease in pervious areas (Table 16).This suggests that the conversion of other LULC classes to built-up classes as a result of urbanization contributes to surface imperviousness (Shuster et al. 2005).An increasing trend of urbanization across the catchment makes a direct influence on the runoff coefficient and equivalent impervious area (Fig. 8 and Table 16).A noticeable spatial change in the runoff coefficient and equivalent impervious area are particularly significant in the central, northern, and western parts of the catchment.The increase in the runoff coefficient and equivalent impervious area increases the volume of runoff generated from the catchment by reducing infiltration of water into the underlying soil layer thus impacting groundwater recharge.Runoff also occurs more rapidly from impervious areas than runoff from pervious areas-hence flooding downstream.The weighted  runoff coefficient of the Dry Creek catchment increased from 0.33 to 0.53 while EIA from 38.34 to 61.58 km 2 , increasing the volume of runoff estimated using the rational method.Greater runoff volume increases the risk of flooding in downstream sections of the catchment.A study conducted on the effects of urbanization on storm runoff in Kathmandu city by Bajracharya et al. (2015) also showed that urbanization increases the surface imperviousness and runoff volume.Similarly, a study by Zhang et al. (2015) in Madagascar found an increase in urbanized area by 52%, runoff depth by 58% and runoff coefficient by 5.8%.Further, the study by McMahon et al. (2003) emphasized that the assessment of extreme low and high-flow episodes' frequency, magnitude, and duration should consider the impervious surface area's contribution to streamflow generation processes.
Irrespective of the rainfall condition, the catchment's weighted curve number exhibited an increase from 80.94 to 87.37 (Table 16).Ansari et al. (2016) also observed the increase in built-up areas in Nagpur urban watershed in India raised the magnitude of curve number from 76 to 80 between 2000 and 2012.The increase in the curve number resulted in a 38.56% decline in retention capacity and increasing the risk of flooding to the surrounding environment.The hydrological cycle is fundamentally affected by the type of soil and land cover because it regulates infiltration and affects surface and groundwater fluxes (Jaafar et al. 2019;Malede et al. 2023).Overall, changes in LULC due to urbanization in Dry Creek alter the surface imperviousness, runoff coefficient, equivalent impervious area, curve number, and retention capacity of the soil, which could potentially have an impact on the hydrological conditions of the catchment, including runoff and groundwater recharge.

Conclusion
The study of LULC change caused by urbanization is crucial for analyzing the influence of urban sprawl on catchment surface imperviousness and its consequences for the hydrological processes of urban catchments.High-quality data and a reliable classification method using remote sensing are helpful in the study of changes in land use land cover.LULC classification utilizing high-quality images provide more accurate classification results, which can subsequently be used for environmental and hydrological studies.In the present study, very high-resolution aerial imagery with a resolution of 0.075 m and SPOT satellite imagery with a resolution of 2.5 m were utilized to assess the long-term spatiotemporal LULC changes caused by urban sprawl over the Dry Creek catchment.This study developed LULC maps of the area and assessed categorical changes in the catchment using an advanced and user-friendly ArcGIS Pro 3.0 environment employing the SVM method of image classification.Following an examination of the accuracy assessment performance indicators of automatic classification, image post-classification was demonstrated to improve the classification accuracy.The overall classification accuracy was found to be high, with kappa coefficient and overall accuracy results greater than 0.83 and 0.88, respectively.The SVM image classification approaches performed well in classifying high-resolution imageries, with good performance metrics.The built-up area, plantation and water expanded by 34.0 km 2 , 8.7 km 2 and 0.4 km 2 , respectively from 1979 to 2019, whereas grass cover declined by 43.09 km 2 , revealing a significant transformation in the LULC.Between 1979 and 2019, the study area has undergone substantial changes in land use and land cover.The LULC in the study area changed categorically by 62.98% during the study period.High rates of change were observed in the grass and built-up areas, with respective values of −1.08 and 0.85 km 2 /yr.Urbanization raised the impervious surface from 24.90 to 58.90 km 2 with 136.55% urban expansion.During the study period, an increase in surface imperviousness increased the runoff coefficient and curve number by 60.63% and 7.95%, respectively, while decreasing soil retention capacity by 38.56%.This undoubtedly increases the quantity of catchment surface runoff by increasing the risk of flooding to the downstream environment.The findings of this study have significant implications for the development and implementation of flood and stormwater control strategies.Therefore, future research should focus on how LULC changes caused by urbanization affect runoff and other hydrological and morphological catchment parameters.

Fig. 1
Fig. 1 Location map of dry creek catchment

Fig. 3
Fig. 3 Spatial distribution of training samples (left) and validation samples (right)

Fig. 5 Fig. 6
Fig. 5 Percentage change in each LULC class over the study period

Fig. 7
Fig. 7 Variation of LULC changes against population changes over the study period

Fig. 8
Fig. 8 Spatiotemporal distribution of curve number and EIA within Dry Creek catchment and Table 3) for each LULC class.Object-based SVM image classification was carried out by generating 500 randomly stratified training samples per class.

Table 1
Details of the data used for this study DEW = Department of Environment and Water, R = Red, G = Green, B = Blue, CIR = Color infrared, ORNL DAAC = OAK RIDGE National laboratory distributed active archive center for biogeochemical dynamics, Australian bureau of statistics (Adelisardou et al. 2022;Shawul and Chakma 2019;Zhang et al. 2020' accuracy (PA) parameters(Adelisardou et al. 2022;Shawul and Chakma 2019;Zhang et al. 2020).

Table 3
Training and validation samples used for classification and accuracy assessment

Table 4
Error Congalton and Green (2019)ing the pixel count in the mapped land cover category i and the reference land cover category j adapted fromCongalton and Green (2019)n is the total pixel count

Table 5
Classification accuracy of 1979 aerial imagery * Overall accuracy

Table 8
Fig. 4 Spatiotemporal distribution of LULC classes in Dry Creek catchment −1.08 km 2 /yr for grass cover and increment of 0.85 km 2 /yr for built-up areas (Table

Table 9
Rate of LULC change over the study period

Table 10
Categorical LULC changes between 1979 and 2006

Table 11
Categorical LULC change between 2006 and 2019

Table 12
Categorical LULC change between 1979 and 2019

Table 16
Surface imperviousness, runoff coefficient and retention capacity of dry creek catchment * AMC II indicates the antecedent moisture condition II, (