Finding Coastal Megaclast Deposits: A Virtual Perspective

: Coastal megaclast deposits are dominated by detrital particles larger than 1 m in size. These attract signiﬁcant attention of modern researchers because of the needs of sedimentary rock nomenclature development and interpretation of storm and tsunami signatures on seashores. If so, ﬁnding localities that exhibit coastal megaclast deposits is an important task. Field studies do not o ﬀ er a quick solution, and, thus, remote sensing tools have to be addressed. The application of the Google Earth Engine has permitted to ﬁnd four new localities, namely Hondarribia in northern Spain (Biscay Bay), the Ponza Island in Italy (Tyrrhenian Sea), the Wetar Island in eastern Indonesia (Banda Sea), and the Humboldt o Coredo Bay at the Colombia / Panama border (eastern Paciﬁc). In these localities, coastal megaclast deposits consisting of blocks (1–10 m in size) and some megablocks ( > 10 m in size) are delineated and preliminary described in regard to the dominant size of particles, package density, mode of occurrence, etc. The limitations of such virtual surveys of coastal megaclast deposits are linked to an insu ﬃ ciently high resolution of satellite images, as well as ‘masking’ e ﬀ ects of vegetation cover and cli ﬀ shadows. However, these limitations do not diminish the importance of the Google Earth Engine for ﬁnding these deposits. Consideration of some tourism-related information, including photos captured by tourists and bouldering catalogues, facilitates search for promising areas for subsequent virtual surveying of megaclast distribution. It is also established that the Google Earth Engine permits quantitative analysis of composition of coastal megaclast deposits in some areas, as well as to register decade-long dynamics or stability of these deposits, which is important to interpret their origin. The current opportunities for automatic detection of coastal megaclast deposits seem to be restricted.


Introduction
Coastal sediment dynamics is always highly complex, and, thus, it provides a lot of issues for geoscience investigations. An interest in nomenclature development for large clasts that started near the end of the 20th century [1] and attention to coastal hazards facilitated by the Indian Ocean Tsunami of 2004 [2] have shaped a new international research direction, namely megaclast studies [3]. Megaclasts of ocean and sea coasts have been investigated most intensively. These studies have been conducted in so different places of the world as Baja California in Mexico [4,5], North Eleuthera of the Bahamas [6,7], and Rabat in Morocco [8]. Although the works have tended to focus on only some regions like the Mediterranean (e.g., [9][10][11][12][13]), the global evidence of large clast accumulations has become huge already, and it continues growing. This evidence requires extension and generalization for further conceptual treatment. Previous reviews of rocky shores [14] and megaclasts [3] confirm such research is promising. Undoubtedly, coastal megaclast deposits of Quaternary age are of special importance because the evidence of such deposits from the earlier geological periods remains scarce [3,15]. account that megaclasts are >1 m in size, the resolution of the available satellite images permits finding them in almost all regions and describing them preliminarily when image resolution is especially high. A significant advantage of this approach is a no-cost, quick, and geographically unrestricted search for megaclast localities, including those in difficult-to-access places. Anyway, the efficacy of the Google Earth Engine in the search for coastal megaclast deposits needs testing, and the latter was addressed in the present study. The Google Earth Pro 7.3.2.5776 version (free software) of the Google Earth Engine was employed for this purpose. Images of the maximum appropriate resolution are preferred. In some cases, resolution can be increased a bit else, but this leads to smooth contours of natural objects, i.e., to poor-visibility of megaclasts.
In order to test the application of the Google Earth Engine to coastal megaclast deposits, four localities of the latter were considered. These were found as a result of tentative visual coastline surveying with the Google Earth Engine in those areas where the occurrence of megaclasts seems to be possible, but has not been reported earlier. The localities were Hondarribia in northern Spain (the Biscay Bay coast), the Ponza Island of Italy (Tyrrhenian Sea), the Wetar Island in eastern Indonesia (Banda Sea), and the Humboldt o Coredo Bay and vicinities near the border between Colombia and Panama (eastern Pacific coast; Figure 1). These localities represent different geographical domains, namely Atlantic Europe, the Mediterranean, Southeast Asia, and Central America. For each locality, two plots are selected provisionally for more representative approach testing. When possible, these plots were selected so to demonstrate local differences of coastal megaclast deposits. Although this paper focused on finding localities, preliminary qualitative descriptions of the studied coastal megaclast deposits were provided, as this study explored the very possibility of their 'visibility' with the Google Earth Engine. The descriptions avoid details relevant to remote sensing techniques, but emphasize on the general character of the deposits, which is the primary interest of sedimentologists. In other words, the descriptions were addressed to sedimentologists, not specialists in remote sensing. Importantly, these localities were selected to demonstrate the utility of satellite images with a different resolution. it has been used efficiently in landslide susceptibility mapping [42,43]. The most elementary function of the Google Earth Engine is visual surveying of the Earth surface at an appropriate resolution. Taking into account that megaclasts are >1 m in size, the resolution of the available satellite images permits finding them in almost all regions and describing them preliminarily when image resolution is especially high. A significant advantage of this approach is a no-cost, quick, and geographically unrestricted search for megaclast localities, including those in difficult-to-access places. Anyway, the efficacy of the Google Earth Engine in the search for coastal megaclast deposits needs testing, and the latter was addressed in the present study. The Google Earth Pro 7.3.2.5776 version (free software) of the Google Earth Engine was employed for this purpose. Images of the maximum appropriate resolution are preferred. In some cases, resolution can be increased a bit else, but this leads to smooth contours of natural objects, i.e., to poor-visibility of megaclasts. In order to test the application of the Google Earth Engine to coastal megaclast deposits, four localities of the latter were considered. These were found as a result of tentative visual coastline surveying with the Google Earth Engine in those areas where the occurrence of megaclasts seems to be possible, but has not been reported earlier. The localities were Hondarribia in northern Spain (the Biscay Bay coast), the Ponza Island of Italy (Tyrrhenian Sea), the Wetar Island in eastern Indonesia (Banda Sea), and the Humboldt o Coredo Bay and vicinities near the border between Colombia and Panama (eastern Pacific coast; Figure 1). These localities represent different geographical domains, namely Atlantic Europe, the Mediterranean, Southeast Asia, and Central America. For each locality, two plots are selected provisionally for more representative approach testing. When possible, these plots were selected so to demonstrate local differences of coastal megaclast deposits. Although this paper focused on finding localities, preliminary qualitative descriptions of the studied coastal megaclast deposits were provided, as this study explored the very possibility of their 'visibility' with the Google Earth Engine. The descriptions avoid details relevant to remote sensing techniques, but emphasize on the general character of the deposits, which is the primary interest of sedimentologists. In other words, the descriptions were addressed to sedimentologists, not specialists in remote sensing. Importantly, these localities were selected to demonstrate the utility of satellite images with a different resolution.

Results
Four small case studies were undertaken in order to demonstrate the efficacy of the virtual examination of coastal megaclast deposits with the Google Earth Engine. Two examples deal with satellite images of exceptionally high resolution, and two other examples deal with satellite images of appropriate but somewhat limited resolution. In each case, it is also attempted to stress the relative importance of the new evidence to the regional knowledge of coastal megaclast deposits. For this purpose, the basic published information was summarized.

Example 1: Hondarribia (Spain; Biscay Bay)
The locality is situated in the vicinities of the town of Hondarribia in northern Spain (Figure 1), where the rocky shore of the Biscay Bay hosts a lot of megaclasts that form significant accumulations. Their origin is linked to the destruction of the Paleogene marine rocks exposed in cliffs and directly along the coastline [44]. Apparently, the main factors of megaclast accumulation are gravity processes and wave abrasion.
On the eastern plot, the coastal megaclast deposits were represented by the dense accumulation of angular blocks with the size of 1-5 m; the maximum size reached 7-8 m ( Figure 2). Importantly, the deposits were spatially heterogeneous-the mean size of blocks increased locally (for instance, near the shoreline and close to parent rock exposures). Although the deposits covered almost the entire area, some 'islands' of parent rock exposures and semi-detached megaclasts were also visible. On the western plot, angular megaclasts were larger in size: these include blocks and megablocks; the size of the largest particle exceeded 25 m (Figure 3). Megablocks were more numerous in the right half of the plot. Parent rock exposures and semi-detached megaclasts occurred in the left half. On both plots, clasts >1 m in size were distinguished unequivocally, and the evidence was enough to confirm the presence of true coastal megaclast deposits. The resolution of the satellite view permitted us to distinguish blocks from megablocks; the presence of boulders was evident, although these could not be registered individually. Dense package of the deposits was well visible. Size and shape of individual blocks can be registered when their size exceeds 1-2 m. The spatial position and orientation of stones did not leave an impression of significant re-working, i.e., primary, undisturbed coastal megaclast deposits could be hypothesized.

Results
Four small case studies were undertaken in order to demonstrate the efficacy of the virtual examination of coastal megaclast deposits with the Google Earth Engine. Two examples deal with satellite images of exceptionally high resolution, and two other examples deal with satellite images of appropriate but somewhat limited resolution. In each case, it is also attempted to stress the relative importance of the new evidence to the regional knowledge of coastal megaclast deposits. For this purpose, the basic published information was summarized.

Example 1: Hondarribia (Spain; Biscay Bay)
The locality is situated in the vicinities of the town of Hondarribia in northern Spain (Figure 1), where the rocky shore of the Biscay Bay hosts a lot of megaclasts that form significant accumulations. Their origin is linked to the destruction of the Paleogene marine rocks exposed in cliffs and directly along the coastline [44]. Apparently, the main factors of megaclast accumulation are gravity processes and wave abrasion.
On the eastern plot, the coastal megaclast deposits were represented by the dense accumulation of angular blocks with the size of 1-5 m; the maximum size reached 7-8 m ( Figure 2). Importantly, the deposits were spatially heterogeneous-the mean size of blocks increased locally (for instance, near the shoreline and close to parent rock exposures). Although the deposits covered almost the entire area, some 'islands' of parent rock exposures and semi-detached megaclasts were also visible. On the western plot, angular megaclasts were larger in size: these include blocks and megablocks; the size of the largest particle exceeded 25 m ( Figure 3). Megablocks were more numerous in the right half of the plot. Parent rock exposures and semi-detached megaclasts occurred in the left half. On both plots, clasts >1 m in size were distinguished unequivocally, and the evidence was enough to confirm the presence of true coastal megaclast deposits. The resolution of the satellite view permitted us to distinguish blocks from megablocks; the presence of boulders was evident, although these could not be registered individually. Dense package of the deposits was well visible. Size and shape of individual blocks can be registered when their size exceeds 1-2 m. The spatial position and orientation of stones did not leave an impression of significant re-working, i.e., primary, undisturbed coastal megaclast deposits could be hypothesized.   The information from this locality fills an important gap in the knowledge of the geographical distribution of coastal megaclast deposits along the Atlantic coast of Europe. These deposits were reported previously from the British Isles [21,[27][28][29][30][31][32][33][34][35][45][46][47][48][49][50][51], northern France (first of all, Brittany) [52][53][54][55][56], and the western and southwestern coasts of the Iberian Peninsula [57][58][59][60][61]. The Hondarribia locality represents the other domain, i.e., the Biscay Bay, the coasts of which are prone to severe storms creating and transporting megaclasts, similarly to the above-mentioned regions [62]. This may be the first megaclast locality reported from the Biscay Bay coast.

Example 2: Ponza Island (Italy; Tyrrhenian Sea)
The locality corresponds to the Ponza Island (Pontine Archipelago) that is situated in the Tyrrhenian Sea, near the western coast of central Italy ( Figure 1). The Pliocene volcanic and volcaniclastic rocks are exposed in tall cliffs [63][64][65], the collapse of which leads to the formation and accumulation of numerous megaclasts. As rocky shores with cliffs dominate the island, the distribution of coastal megaclast deposits is significant there. Gravity processes and abrasion contribute to local megaclast formation.
On the southern plot, the coastal megaclast deposits were represented by the accumulation of angular blocks with the common size of 1-5 m (Figure 4). Blocks of a bigger size and rare small megablocks (up to 12-13 m in size) occurred along the very shoreline, as well as in the northern part of the plot. Density of the deposit differed: it was significant near the shoreline, but megaclasts occurred individually or in small groups in the other places where they were either mixed with boulders or lay directly on the rock surface. On the northwestern plot, a small group of fully-and partially-detached megaclasts was observed ( Figure 5). These include several megablocks with the size up to 40 m and numerous blocks in between. These megaclasts were angular, and some smoothed surfaces result from wave and wind 'polishing'. On both plots, clasts >1 m in size were distinguished easily, and the delineation of true coastal megaclast deposits was easy. The resolution of satellite images permitted us to distinguish blocks from megablocks, the presence of boulders was evident, although these could not be registered individually. Different package density was visible. Size and shape of individual particles can be registered when their size exceeded 1-2 m. However, a shadow from the tall cliff on the northwestern plot 'masked' some megaclasts and did not permit us to characterize them on an individual basis ( Figure 5). Some redeposition of megaclasts could not be excluded for the southern plot, but the view of the deposits on the northwestern plot leaves impression of a 'fresh' and then undisturbed rockfall. The information from this locality fills an important gap in the knowledge of the geographical distribution of coastal megaclast deposits along the Atlantic coast of Europe. These deposits were reported previously from the British Isles [21,[27][28][29][30][31][32][33][34][35][45][46][47][48][49][50][51], northern France (first of all, Brittany) [52][53][54][55][56], and the western and southwestern coasts of the Iberian Peninsula [57][58][59][60][61]. The Hondarribia locality represents the other domain, i.e., the Biscay Bay, the coasts of which are prone to severe storms creating and transporting megaclasts, similarly to the above-mentioned regions [62]. This may be the first megaclast locality reported from the Biscay Bay coast.

Example 2: Ponza Island (Italy; Tyrrhenian Sea)
The locality corresponds to the Ponza Island (Pontine Archipelago) that is situated in the Tyrrhenian Sea, near the western coast of central Italy ( Figure 1). The Pliocene volcanic and volcaniclastic rocks are exposed in tall cliffs [63][64][65], the collapse of which leads to the formation and accumulation of numerous megaclasts. As rocky shores with cliffs dominate the island, the distribution of coastal megaclast deposits is significant there. Gravity processes and abrasion contribute to local megaclast formation.
On the southern plot, the coastal megaclast deposits were represented by the accumulation of angular blocks with the common size of 1-5 m ( Figure 4). Blocks of a bigger size and rare small megablocks (up to 12-13 m in size) occurred along the very shoreline, as well as in the northern part of the plot. Density of the deposit differed: it was significant near the shoreline, but megaclasts occurred individually or in small groups in the other places where they were either mixed with boulders or lay directly on the rock surface. On the northwestern plot, a small group of fully-and partially-detached megaclasts was observed ( Figure 5). These include several megablocks with the size up to 40 m and numerous blocks in between. These megaclasts were angular, and some smoothed surfaces result from wave and wind 'polishing'. On both plots, clasts >1 m in size were distinguished easily, and the delineation of true coastal megaclast deposits was easy. The resolution of satellite images permitted us to distinguish blocks from megablocks, the presence of boulders was evident, although these could not be registered individually. Different package density was visible. Size and shape of individual particles can be registered when their size exceeded 1-2 m. However, a shadow from the tall cliff on the northwestern plot 'masked' some megaclasts and did not permit us to characterize them on an individual basis ( Figure 5). Some redeposition of megaclasts could not be excluded for the southern plot, but the view of the deposits on the northwestern plot leaves impression of a 'fresh' and then undisturbed rockfall.  The Mediterranean is a 'classical' region for studies of coastal megaclasts and boulder deposits [9][10][11][12][13]. Particularly, megaclast accumulations have been reported from the Western Mediterranean continental and island coasts [66][67][68]. However, Ponza Island is of special importance because of two reasons. First, it represents lenticular, almost round-island distribution of the coastal megaclast deposits. Second, the latter originates from volcanic and volcaniclastic rocks. Therefore, finding this locality with the Google Earth Engine extends the regional knowledge of coastal megaclast deposits.

Example 3: Wetar Island (Indonesia; Banda Sea)
The locality corresponds to the Banda Sea coast of the Wetar Island in eastern Indonesia ( Figure  1). There, Cenozoic volcanic, volcaniclastic, and reefal carbonate deposits crop out [69]. Presumably, their destruction by wave abrasion and weathering leads to megaclast creation; it cannot be excluded that some megaclasts were formed as a result of slope collapse and subsequent downslope   The Mediterranean is a 'classical' region for studies of coastal megaclasts and boulder deposits [9][10][11][12][13]. Particularly, megaclast accumulations have been reported from the Western Mediterranean continental and island coasts [66][67][68]. However, Ponza Island is of special importance because of two reasons. First, it represents lenticular, almost round-island distribution of the coastal megaclast deposits. Second, the latter originates from volcanic and volcaniclastic rocks. Therefore, finding this locality with the Google Earth Engine extends the regional knowledge of coastal megaclast deposits.

Example 3: Wetar Island (Indonesia; Banda Sea)
The locality corresponds to the Banda Sea coast of the Wetar Island in eastern Indonesia ( Figure  1). There, Cenozoic volcanic, volcaniclastic, and reefal carbonate deposits crop out [69]. Presumably, their destruction by wave abrasion and weathering leads to megaclast creation; it cannot be excluded that some megaclasts were formed as a result of slope collapse and subsequent downslope The Mediterranean is a 'classical' region for studies of coastal megaclasts and boulder deposits [9][10][11][12][13]. Particularly, megaclast accumulations have been reported from the Western Mediterranean continental and island coasts [66][67][68]. However, Ponza Island is of special importance because of two reasons. First, it represents lenticular, almost round-island distribution of the coastal megaclast deposits. Second, the latter originates from volcanic and volcaniclastic rocks. Therefore, finding this locality with the Google Earth Engine extends the regional knowledge of coastal megaclast deposits.

Example 3: Wetar Island (Indonesia; Banda Sea)
The locality corresponds to the Banda Sea coast of the Wetar Island in eastern Indonesia ( Figure 1). There, Cenozoic volcanic, volcaniclastic, and reefal carbonate deposits crop out [69]. Presumably, their destruction by wave abrasion and weathering leads to megaclast creation; it cannot be excluded that some megaclasts were formed as a result of slope collapse and subsequent downslope transport. This means that the origin of this new locality can be highly complex, but the role of abrasion seems to be leading.
On the eastern plot, blocks with the maximum size of up to 10 m (the prevailing size was 3-5 m) occurred sporadically ( Figure 6). These lay either individually or in small groups. The biggest stones were angular. The deposit package density was very low with regard to the distance between megaclasts (even when these occur in groups). These coastal megaclast deposits were restricted to the very shoreline. On the northeastern plot, the number of megaclasts was bigger (Figure 7). These were chiefly blocks of different size (commonly < 5 m), but there were also megablocks, the biggest of which reached 20 m in size. The particle angularity was well visible. The package density was moderate (even high in some places), but the deposits also occurred like a narrow ribbon along the shoreline. On the both plots, clasts >3 m in size were distinguished easily; the satellite image evidence is enough to confirm the presence of true coastal megaclast deposits. The resolution of the images permitted us to distinguish blocks from megablocks; boulders were not seen at all (it was unclear whether these exist). A different package density was visible without any difficulty. Size and shape of individual particles can be registered when their size exceeds 5 m. Better to say, these can be examined chiefly for large blocks and all megablocks. Dense vegetation cover of the island and, particularly, along the shoreline 'masks' the coastal megaclast deposits, but this does not preclude for megaclast accumulation tracing over the studied plots. The position and the orientation of megaclasts do not suggest against their resedimentation, and, if so, it cannot be excluded that high-energy events contributed to motion and destruction of these stones. transport. This means that the origin of this new locality can be highly complex, but the role of abrasion seems to be leading. On the eastern plot, blocks with the maximum size of up to 10 m (the prevailing size was 3-5 m) occurred sporadically ( Figure 6). These lay either individually or in small groups. The biggest stones were angular. The deposit package density was very low with regard to the distance between megaclasts (even when these occur in groups). These coastal megaclast deposits were restricted to the very shoreline. On the northeastern plot, the number of megaclasts was bigger (Figure 7). These were chiefly blocks of different size (commonly < 5 m), but there were also megablocks, the biggest of which reached 20 m in size. The particle angularity was well visible. The package density was moderate (even high in some places), but the deposits also occurred like a narrow ribbon along the shoreline. On the both plots, clasts >3 m in size were distinguished easily; the satellite image evidence is enough to confirm the presence of true coastal megaclast deposits. The resolution of the images permitted us to distinguish blocks from megablocks; boulders were not seen at all (it was unclear whether these exist). A different package density was visible without any difficulty. Size and shape of individual particles can be registered when their size exceeds 5 m. Better to say, these can be examined chiefly for large blocks and all megablocks. Dense vegetation cover of the island and, particularly, along the shoreline 'masks' the coastal megaclast deposits, but this does not preclude for megaclast accumulation tracing over the studied plots. The position and the orientation of megaclasts do not suggest against their resedimentation, and, if so, it cannot be excluded that high-energy events contributed to motion and destruction of these stones.   Coastal megaclasts have been reported earlier from many localities of Southeast Asia, including those of China, Taiwan, Thailand, Malaysia, Indonesia, and the Philippines . Strong influences of tsunamis [95,96] and super-typhoons [97] contribute to the formation of such deposits. However, all this evidence has been obtained outside eastern Indonesia and chiefly on coasts open to either the Indian or Pacific Oceans. The Wetar Island locality sheds light on the coastal megaclast deposits formed in a distinct geographical setting, namely on the coast of a relatively small, intra-island sea, which makes addition to the regional knowledge of such deposits. Therefore, the Google Earth Engine permits extension of the available knowledge of megaclast occurrence in Southeast Asia.

Example 4: Humboldt o Coredo Bay (Colombia/Panama; Eastern Pacific)
The locality corresponds to the Humdoldt o Coredo Bay and its vicinities. It is situated on the Pacific coast of Central America, exactly at the border between Colombia and Panama, north of Jurado ( Figure 1) where the Cretaceous-Eocene basement and volcanic arc complex dominate the local geological setting [98]. Wave destruction and weathering of hard parent rocks contributes to formation of multiple megaclasts. Apparently, abrasion is chiefly responsible for the origin of these deposits.
On the southern plot, small and medium blocks (typical size is <5 m) formed lenticular deposits around a small bay (Figure 8). The biggest particles were angular, but the shape of the majority of particles could not be recognized at the available resolution. The width of the ribbon of deposits differed, as well as their package density. A few huge blocks (but these did not reach the size of megablocks) lay individually in the central part of the plot. Generally, it is possible to record the deposits heterogeneity. On the northern plot, the coastal megaclast deposits were very similar to those described on the southern plot, with two exceptions: the mean size of blocks was bigger, and there were a few megablocks (up to 20 m in size or even more; Figure 9). On both plots, clasts >3 m in size were distinguished easily, and the presence of true coastal megaclast deposits was evident. The resolution of satellite images permitted us to distinguish blocks from megablocks; boulders were not seen at all, although these, presumably, exist. A different package density was visible without any difficulty. The size and shape of individual particles can be registered when their size exceeds 5 m (large blocks and megablocks). The difficulties were linked, first, to the 'masking' effect of the vegetation cover and, second, uncertainty with some features on the northern plot that could be either megaclasts or parent rock exposures (or semi-detached megaclasts). Presumably, the both difficulties will remain even in the case of a much higher resolution of satellite images. The spatial position and orientation of stones permitted us to hypothesize the absence of significant reworking, i.e., these seemed to be primary, undisturbed coastal megaclast deposits.
Previous megaclast studies on the Pacific coast of the Americas are scarce, and, particularly, focused on two geographical domains, namely the Baja California Peninsula in Mexico [4,5] and northern and central Chile [41,99,100]. If so, the information from the Humboldt o Coredo Bay locality fills significant gap in the regional knowledge of coastal megaclast deposits characterizing those of Central America. Finding this locality with the Google Earth Engine implies the existence of coastal megaclast deposits on the Pacific coast of Central America (probably, for the first time).
(large blocks and megablocks). The difficulties were linked, first, to the 'masking' effect of the vegetation cover and, second, uncertainty with some features on the northern plot that could be either megaclasts or parent rock exposures (or semi-detached megaclasts). Presumably, the both difficulties will remain even in the case of a much higher resolution of satellite images. The spatial position and orientation of stones permitted us to hypothesize the absence of significant reworking, i.e., these seemed to be primary, undisturbed coastal megaclast deposits.  Previous megaclast studies on the Pacific coast of the Americas are scarce, and, particularly, focused on two geographical domains, namely the Baja California Peninsula in Mexico [4,5] and northern and central Chile [41,99,100]. If so, the information from the Humboldt o Coredo Bay locality fills significant gap in the regional knowledge of coastal megaclast deposits characterizing those of Central America. Finding this locality with the Google Earth Engine implies the existence of coastal megaclast deposits on the Pacific coast of Central America (probably, for the first time).

Synthesis: Virtual Approach Opportinities and Limitations
The four examples given above illustrate the application of the Google Earth Engine to the identification and even preliminary description of coastal megaclast deposits. The approach seems to be promising. In all four cases, the deposits were identified and delineated without any difficulty and irrespective to the different resolution of satellite images. All study areas were new to megaclast research, and these were found on the basis of the only virtual exploration of the coastline, without consideration of any preliminary evidence. This approach can be used for visual surveying of the coasts and identification of localities suitable for subsequent field or virtual investigations (or both). Moreover, some basic characteristics of coastal megaclast deposits can be examined with the Google Earth Engine taken alone, although such a study cannot 'replace' in-depth field investigations to be undertaken further. These characteristics include the dominant particle size, package density, and spatial homogeneity versus heterogeneity of deposits. Individual, group, or lenticular occurrence of megaclasts can be also registered. As for the particle size, in all cases, it was not only possible to

Synthesis: Virtual Approach Opportinities and Limitations
The four examples given above illustrate the application of the Google Earth Engine to the identification and even preliminary description of coastal megaclast deposits. The approach seems to be promising. In all four cases, the deposits were identified and delineated without any difficulty and irrespective to the different resolution of satellite images. All study areas were new to megaclast research, and these were found on the basis of the only virtual exploration of the coastline, without consideration of any preliminary evidence. This approach can be used for visual surveying of the coasts and identification of localities suitable for subsequent field or virtual investigations (or both). Moreover, some basic characteristics of coastal megaclast deposits can be examined with the Google Earth Engine taken alone, although such a study cannot 'replace' in-depth field investigations to be undertaken further. These characteristics include the dominant particle size, package density, and spatial homogeneity versus heterogeneity of deposits. Individual, group, or lenticular occurrence of megaclasts can be also registered. As for the particle size, in all cases, it was not only possible to make a distinction of blocks from megablocks, but also to note the presence of smaller and larger blocks, i.e., to deal with subcategories of the main grain-size divisions. Some notions on the possible origin of the deposits are also possible. Importantly, the use of the Google Earth Engine permitted us a very quick survey of coasts. Apparently, one researcher can check up to a thousand kilometers of coastal zone for the presence of megaclast accumulations per 5-7 working hours. No field investigation can allow such a speed of work and so massive information gathering. Moreover, many areas will remain 'untouched' by megaclast research due to the relatively circle of the involved researchers and logistic restrictions. If so, it is better to know the position of the localities and to have the only preliminary characteristics of the deposits than to concentrate on only well-known and well-accessible localities studied for years.
As implied by the examples, virtual studies of coastal megaclast deposits with the Google Earth Engine face limitations of three kinds. First, the resolution of satellite images did not permit us to deal with boulders and other finer components of the deposits, and the presence of these components could not be proven in some cases. Second, the resolution of the satellite images available for the Hondarribia and Ponza Island localities allowed the size and shape documentation of all megaclasts, whereas the resolution of the satellite images of the Wetar Island and Humboldt o Coredo Bay localities permitted doing this for only huge blocks and megablocks, i.e., for the only part of the deposit. In the two latter cases, there were also difficulties with distinction between true megaclasts, semi-detached megaclasts, and exposures of parent rocks. This means that detailed investigations of coastal megaclast deposits and their quantitative descriptions with the Google Earth Engine are possible in only some cases. Third, there are local limitations of the approach like tall cliff shadows (the Ponza Island) or dense vegetation cover (Wetar Island and Humboldt o Coredo Bay). Anyway, all these limitations do not preclude from the documentation of the spatial occurrence of megaclast deposits along coasts of seas and oceans. When necessary, field studies can be undertaken later, and these studies can be coordinated effectively with the information obtained with the Google Earth Engine. In other words, the limitations are less important than the opportunities. Moreover, it is expected that the resolution of satellite images employed by the noted software will continue to increase in the nearest future.

Tourism Information as a Source for Geographical Justification of Virtual Surveys
The total world coastline is too lengthy taking into account continental and all island coasts, and even the Google Earth Engine does not make efficient its visual surveying for the purposes of the identification and description of megaclast deposits. The latter occurs locally, and it is sensible to focus on those areas, which seem to be promising for finding megaclasts. Some previous descriptions ('occasional' descriptions of megaclast deposits can be found in the literature) and geological, geomorphological, and geographical knowledge (for instance, the selection of coasts with tall cliffs and prone to severe storms, super-typhoons, and tsunamis) is helpful in many cases for the geographical justification of virtual surveys. However, some other, non-scientific evidence can be also considered.
The first opportunity is linked to photos provided by the users (first of all, tourists) of the Google Earth Engine, the geographical attachment of which is displayed directly on satellite images. Megaclasts and their accumulations are notable features, and these can be recognized easily on photos. Moreover, the latter can help in satellite image interpretation, i.e., for the correct description of the particle shape and package density, as well as for the distinction between true megaclasts, semi-detached megaclasts, and parent rock exposures; the content of boulders and finer components of deposits can be also registered. In the present study, the photos provided by the Google Earth Engine were helpful in the cases of the Hondarribia and Ponza Island localities. However, the number of the proper photos was limited for some remote, infrequently visited places, or these were not available at all.
The second opportunity that can be recommended for the geographical inventory of coastal boulder deposits is also relevant to tourism-linked information available online. Among various sport and tourism activities, bouldering has gained importance [101][102][103][104][105]. A lot of information on individual 'boulders' (these may be either isolated rocks, but also true megaclasts) has been accumulated. This information is easily accessible with the available on-line catalogues [106][107][108], and it permits finding coastal megaclasts and their groups, as well as their preliminary visual examination with the available images. For instance, there are typical coastal boulder deposits on the Ao Sane Beach in Phuket, Thailand [106]. These deposits consist of mixed angular blocks and boulders, and there are also angular blocks lying separately on a sandy substrate. In regard to the position of the later, resedimentation cannot be excluded. Such information can be of the utmost importance for finding promising areas to focus virtual surveys of coastal megaclast deposits.

Do Satellite Images Provide Data for the Quantitative Analysis of Coastal Megaclast Deposits?
Finding new localities of coastal megaclast deposits with the Google Earth Engine poses a question of more detailed analysis providing sedimentological information. The examples given above imply that some descriptions of such deposits are possible, but these are preliminary and too qualitative. Actually, in some cases, quantitative interpretations seem to be possible. For instance, grain-size distribution of megaclasts can be analyzed with statistical tools. There are four conditions that allow such an investigation. First, the resolution of satellite images of the study area should allow unequivocal recognition of particles with the size of 1-2 m (minimal size of megaclasts). Second, there should be not any 'masking' effects of shadows, vegetation, etc. Third, the majority of megaclasts should be oriented subhorizontally because vertical orientation does not permit the measurement of the maximum size. Fourth, megaclasts should form a relatively thin layer' alternatively, the lowermost particles covered by the other blocks cannot be measured.
From the four examples considered in this study, the only Hondarribia locality satisfied all four conditions outlined above. To demonstrate the efficacy of satellite image analysis, the western plot of the locality (Figure 3) was considered. The Google Earth Engine provides the option of the measurement of objects directly on a satellite image. The maximum size of 100 megaclasts was measured with this tool on the chosen plot. Apparently, the measured particles constituted about 70% of all megaclasts visible on this plot. The size of megaclasts varied between 1.0 and 26.2 m. The mean size was 7.8 m, and the median size was 6.3 m. These values corresponded to coarse blocks, according to the classification of Bruno and Ruban [16]. The grain-size distribution of megaclasts is shown on a histogram ( Figure 10). The majority of particles had the size of 2.5-5 m, i.e., these were medium blocks. However, coarse blocks (5-10 m) and fine megablocks (10-25 m) were the most common. Generally, these coastal boulder deposits included 37% of coarse blocks, 26% of medium blocks, 24% of fine megablocks, 12% of fine blocks, and 1% of medium megablocks, i.e., it was dominated by the particles with the size of 2.5-25 m, which indicates on a restricted sorting. Undoubtedly, this deposit also bore some finer components like boulders, gravels, and even sand, but these were invisible on the satellite image. Anyway, the dominance of the noted megaclasts in these deposits was indisputable.  [16]. The grain-size distribution of megaclasts is shown on a histogram ( Figure 10). The majority of particles had the size of 2.5-5 m, i.e., these were medium blocks. However, coarse blocks (5-10 m) and fine megablocks (10-25 m) were the most common. Generally, these coastal boulder deposits included 37% of coarse blocks, 26% of medium blocks, 24% of fine megablocks, 12% of fine blocks, and 1% of medium megablocks, i.e., it was dominated by the particles with the size of 2.5-25 m, which indicates on a restricted sorting. Undoubtedly, this deposit also bore some finer components like boulders, gravels, and even sand, but these were invisible on the satellite image. Anyway, the dominance of the noted megaclasts in these deposits was indisputable.

Satellite Evidence of a Decade-Long Stability of Coastal Megaclast Deposits
The modern studies of megaclast focus much (even over-emphasize) on their origin and, particularly, their relevance to extreme events such as storms and tsunamis [4][5][6][7][8][9][10][11][12][13][19][20][21][22][23][24][25][26][27]. Evidently, these studies deal with the dynamics of the coastline and megaclasts themselves. The Google Earth Engine gives exceptional opportunity to contribute to such investigations, as a time series of satellite images is available. This means that direct comparison of the views of some areas at different time slices is possible. Three limitations are as follows. First, 'old' satellite images are not available for all areas and the date of these images differs for different areas. Second, the resolution of 'old' images can be lower than that of the current images. Third, due to a different time of image making, shadows, vegetation, and, importantly, water level may look differently. However, the third problem can be overcome in many cases.
In order to understand the importance of the analysis of the image time series for the understanding of the dynamics of coastal boulder deposits, the new localities identified in the present work were taken as examples. 'Old' satellite images (captured 10-20 years ago) of appropriate resolution are available for three of them (Hondarribia, Ponza Island, and Wetar Island-one plot was selected for analysis in each case). However, in the only case of the Ponza Island, the maximum appropriate resolution of the 'old' image was the same as that of the 'new' image. In two other cases (Hondarribia and Wetar Island), the resolution of the 'old' image was lower, but it allowed recognition of the principal features, i.e., megaclasts, via comparison to the higher-resolution 'new' image. Megaclasts from the western plot of the Hondarribia locality occupied the same position in 2001 ( Figure 11) as in 2018 when the currently available image was captured ( Figure 3). Apparently, the number of blocks and megablocks did not change, although some of them had worse visibility on the 'old' image because of the higher sea level (apparently, the time of tide). No changes were also found between 2007 ( Figure 12) and 2017 ( Figure 5) on the northwestern plot of the Ponza Island. Moreover, the 'old' image was of better quality in this case. The shadow from a tall cliff did not 'mask' a part of the deposits, and the better seawater transparency permitted us to document that up to a quarter of the deposits stretch to the submarine environment. Finally, nothing related to megaclasts changed on the eastern plot of the Wetar Island between 2009 ( Figure 13) and 2019 ( Figure 6).      The evidence presented above implies a decade-long stability of the coastal megaclast deposits in the new localities. Megaclasts were not destroyed, replaced, or over-turned despite that the regions are susceptible to severe weather conditions [62], and tsunamis cannot be excluded there [109,110]. Even if the periodicity of the high-energy events that can be responsible for changes of coastal megaclast deposits seems to be longer than a decade, the registered stability is notable. This is especially the case of the Hondarribia locality located on a high-energy coast of the Biscay Bay, for which a two-decade long comparison (Figures 3 and 11) is possible. Apparently, severe storms of the 2000s and the 2010s [62] did not affect megaclasts of this locality.
It is sensible to add that the Google Earth Engine gives opportunity to get images with a high frequency (Sentinel 2 or SPOT 7 images). If so, some localities can be specially monitored in the case of any future severe storms or tsunamis to document the dynamics of coastal megaclast deposits.

Automatic Detection of Coastal Megaclast Deposits: the Current State of the Problem
The development of the Google Earth Engine has permitted to pose the question of the invention of advanced digital tools for automatic detection of particular features on satellite images. For instance, such tools are helpful in urban studies [111] and archaeology [112]. In geomorphological studies, two successfully tested approaches are notable. The first has been proposed by Luijendijk et al. [113] for the global-scale mapping of sandy beaches and establishing their dynamics. The second approach has been employed by Vos et al. [114] to document shoreline dynamics (some examples also deal with sandy beaches). The validity and the importance of the both above-mentioned studies are indisputable, and the availability of any similar approach would help in the global mapping of coastal megaclast deposits. It is not the purpose of the present work to The evidence presented above implies a decade-long stability of the coastal megaclast deposits in the new localities. Megaclasts were not destroyed, replaced, or over-turned despite that the regions are susceptible to severe weather conditions [62], and tsunamis cannot be excluded there [109,110]. Even if the periodicity of the high-energy events that can be responsible for changes of coastal megaclast deposits seems to be longer than a decade, the registered stability is notable. This is especially the case of the Hondarribia locality located on a high-energy coast of the Biscay Bay, for which a two-decade long comparison (Figures 3 and 11) is possible. Apparently, severe storms of the 2000s and the 2010s [62] did not affect megaclasts of this locality.
It is sensible to add that the Google Earth Engine gives opportunity to get images with a high frequency (Sentinel 2 or SPOT 7 images). If so, some localities can be specially monitored in the case of any future severe storms or tsunamis to document the dynamics of coastal megaclast deposits.

Automatic Detection of Coastal Megaclast Deposits: the Current State of the Problem
The development of the Google Earth Engine has permitted to pose the question of the invention of advanced digital tools for automatic detection of particular features on satellite images. For instance, such tools are helpful in urban studies [111] and archaeology [112]. In geomorphological studies, two successfully tested approaches are notable. The first has been proposed by Luijendijk et al. [113] for the global-scale mapping of sandy beaches and establishing their dynamics. The second approach has been employed by Vos et al. [114] to document shoreline dynamics (some examples also deal with sandy beaches). The validity and the importance of the both above-mentioned studies are indisputable, and the availability of any similar approach would help in the global mapping of coastal megaclast deposits. It is not the purpose of the present work to develop such an advanced interpretative technique, but some limitations of the latter can be discussed in the light of the present findings.
First of all, the main differences of coastal megaclast deposits from sandy beaches should be noted. These deposits are more localized and often narrow; these rarely form lengthy 'belts' like sandy beaches. The latter are often distinguished on satellite images by a white or yellow color, whereas coastal megaclast deposits may look very differently in each case. Importantly, these deposits are often not continuous particle packages, but individual stones or stone groups dispersed along the coastline (this is especially well visible at the Weater Island locality- Figures 6 and 7). To make a clear distinction between true megaclasts and semi-detached megaclasts or rock exposures is not always easy, especially when the surface of a given area is not flat. When sandy beaches are detected, the objects of study are beaches themselves, and the researchers do not to pay attention to sand particles that cannot be determined from the space. In contract, studies of coastal megaclast deposits require attention to giant particles. That is why such studies depend stronger on image resolution. The experience with the Ponza Island ( Figures 5 and 12) implies water transparency influences on the detection of megaclasts that are fully or partly drowned, whereas the evidence from the Hondarribia locality (Figures 3 and 11) indicates on the importance of the sea level, i.e., tides. As a result, visibility of coastal megaclast deposits decreases under certain conditions ( Figure 14). influences on the detection of megaclasts that are fully or partly drowned, whereas the evidence from the Hondarribia locality (Figures 3 and 11) indicates on the importance of the sea level, i.e., tides. As a result, visibility of coastal megaclast deposits decreases under certain conditions ( Figure  14). The above-said does not suggest against the necessity of techniques for automated detection of coastal megaclast deposits, but reveals serious barriers for the invention of such techniques. Moreover, it should be stressed that the studies of Luijendijk et al. [113] are aimed at making globalor regional-scale conclusions. In contrast, what do the megaclast researchers actually need are particular localities, especially providing representative examples of coastal megaclast deposits from different parts of the world. In regard to this, simple visual surveys of shorelines with the Google Earth Engine provide reasonable information.

Conclusions
The use of the Google Earth Engine has permitted to us find four new localities of coastal megaclast deposits in Atlantic Europe, the Mediterranean, Southeast Asia, and Central America. The dominance of blocks and subordinate number of megablocks, as well as their chiefly angular shape are visible on satellite images. More detailed, quantitative investigations are possible on the basis of the both field and virtual studies of these localities. Despite some limitations, the Google Earth The above-said does not suggest against the necessity of techniques for automated detection of coastal megaclast deposits, but reveals serious barriers for the invention of such techniques. Moreover, it should be stressed that the studies of Luijendijk et al. [113] are aimed at making global-or regional-scale conclusions. In contrast, what do the megaclast researchers actually need are particular localities, especially providing representative examples of coastal megaclast deposits from different parts of the world. In regard to this, simple visual surveys of shorelines with the Google Earth Engine provide reasonable information.

Conclusions
The use of the Google Earth Engine has permitted to us find four new localities of coastal megaclast deposits in Atlantic Europe, the Mediterranean, Southeast Asia, and Central America. The dominance of blocks and subordinate number of megablocks, as well as their chiefly angular shape are visible on satellite images. More detailed, quantitative investigations are possible on the basis of the both field and virtual studies of these localities. Despite some limitations, the Google Earth Engine seems to be an almost ideal instrument for quick tracing the geographical distribution of megaclast deposits along the coasts of seas and oceans, which is important for the coordination of further research. The growth of tourism and voluminous online representation of the relevant photos facilitate finding areas promising for identification of megaclast localities. The novelty of this study is triple, i.e., the potential of virtual search for coastal megaclast deposits was proven, four new localities found, and the knowledge of megaclast spatial occurrence in several large regions was summarized (it was also explained how geographical gaps in this knowledge could be filled with virtual finding of new localities). More generally, this work did not argue for 'replacing' field studies by virtual surveys. In contrast, it demonstrated how to find efficiently new localities for further field studies. In the other words, finding localities as a research task is separated from in-depth sedimentological investigations.
This paper dealt with four, almost randomly selected localities. Further virtual surveys of the coastal zones will permit to extend the relevant knowledge, as well as to improve the methods of the satellite image analysis. One of the central problems is the precise megaclast size and shape description, as satellite images show the only 2D projection of these features, which can be oriented differently. It is important problem to think how some tools proposed earlier [115][116][117] can be coupled with the use of the Google Earth Engine in coastal megaclast studies. It cannot be excluded that satellite images can be used for automatic detection of megaclast-promising areas in the future.
Author Contributions: D.A.R. is the only author of this work, and he is fully responsible for its content. The Author has read and agree to the published version of the manuscript.
Funding: This research received no external funding.