Airborne LiDAR Reveals a Vast Archaeological Landscape at the Nan Madol World Heritage Site

: An airborne LiDAR survey of the Nan Madol World Heritage Site and adjacent Temwen Island revealed a complex, irrigated cultivation system, the ﬁrst found in the Central and Eastern Caroline Islands. This informs the goals of the sustainable conservation project, funded by the U.S. Department of State Ambassadors Fund for Cultural Preservation, that inspired the survey, and expands understanding of Nan Madol and its place in the network of Paciﬁc island interaction and trade. Fieldwork veriﬁed the presence, across Temwen, of low, wet, cultivable areas, many of which are connected by water channels or separated by earthen berms. The berms themselves may also have been cultivated. In complexity, labor investment, and organization, the system is comparable to Nan Madol itself, the largest archaeological site in Micronesia, with structures on about 100 artiﬁcial islets built of stone and coral on a reef ﬂat. Constructed over a millennium, Nan Madol was the seat of the Saudeleur Dynasty, which persisted from about 1200 to 1600 CE. The cultivation system appears to have been able to provide ample food for consumption, feasting, and redistribution or trade. If the landscape alteration described here proves to date to the time of the Saudeleur Dynasty, it will o ﬀ er many avenues of research into the economic basis of Nan Madol’s regional dominance.


Introduction
An airborne LiDAR survey of the UNESCO World Heritage Site of Nan Madol and Temwen Island has revealed a vast archaeological landscape consisting of a network of cultivation features connected by a complex water management system. This landscape is adjacent to the core area of Nan Madol.
The LiDAR survey was completed to inform a sustainable conservation and management plan that is funded by the U.S. Department of State and conducted by Cultural Site Research and Management (CSRM) Foundation. Data products and models generated from the LiDAR data will guide the physical conservation and management of Nan Madol and enhance site interpretation. They also enrich understanding of the scope, organization, and influence of the Saudeleur Dynasty, which ruled from Nan Madol. The revealed features have the potential to inform the chronology and pattern of Riesenberg as having had specific or noteworthy functions. Some islets, for example, served as residences or burial places for Saudeleur elites, while others hosted rituals that enhanced the Dynasty's socio-political cohesion. Meggers and Evans's work supported Riesenberg's ethnographic research [25]. Riesenberg considered the Saudeleur Dynasty "a highly developed socio-political system that permitted the planning of the structures, organization of manpower to build them and to raise and feed the mass of labor necessary for construction" ( [24], p. 5). Subsequent research has corroborated this interpretation, with many other researchers noting the high level of social complexity implied by the construction of Nan Madol [2,[4][5][6][7][8]10,14,15,22,23,26,27]. Riesenberg's comment about the necessity of feeding laborers anticipates the findings on Temwen Island reported here.

Project Objectives
The goal of the CSRM Foundation project is to develop a Sustainable Conservation Plan that will identify steps needed to remove Nan Madol from the List of World Heritage in Danger. This includes recommending physical interventions to mitigate threats to the site. Chief among these threats are (1) destructive vegetation; (2) siltation of the canals between islets, which has encouraged an environment that nurtures destructive vegetation and makes accessing the islets difficult; (3) erosion of islets and damage to structures caused by tidal action and storm surges (both of which become more destructive as the sea level rises); and (4) human activities that damage islets and structures. To be sustainable, conservation measures must include an interpretive plan and a management plan, which are necessary to site maintenance after initial treatments have been made.
LiDAR data was collected in order to assess the condition of Nan Madol's islets and structures; to assist in the identification of vegetation that threatens the structural integrity of the site; and to model Temwen Island's hydrology, in order to understand how it might be contributing to the siltation of the site's canals. Siltation and the proliferation of vegetation have progressed to such an extent that today, visitors can see little to suggest why Nan Madol was once called "The Venice of

Project Objectives
The goal of the CSRM Foundation project is to develop a Sustainable Conservation Plan that will identify steps needed to remove Nan Madol from the List of World Heritage in Danger. This includes recommending physical interventions to mitigate threats to the site. Chief among these threats are (1) destructive vegetation; (2) siltation of the canals between islets, which has encouraged an environment that nurtures destructive vegetation and makes accessing the islets difficult; (3) erosion of islets and damage to structures caused by tidal action and storm surges (both of which become more destructive as the sea level rises); and (4) human activities that damage islets and structures. To be sustainable, conservation measures must include an interpretive plan and a management plan, which are necessary to site maintenance after initial treatments have been made.
LiDAR data was collected in order to assess the condition of Nan Madol's islets and structures; to assist in the identification of vegetation that threatens the structural integrity of the site; and to model Temwen Island's hydrology, in order to understand how it might be contributing to the siltation of the site's canals. Siltation and the proliferation of vegetation have progressed to such an extent that today, visitors can see little to suggest why Nan Madol was once called "The Venice of the Pacific." LiDAR's potential for distinguishing among vegetation, ground, and other features, and for yielding Remote Sens. 2019, 11,2152 4 of 24 precise elevation measurements, figured prominently into the CSRM Foundation's decision to collect the data that led to this paper.
CSRM Foundation planned and oversaw the airborne LiDAR survey of Nan Madol and Temwen Island, which took place on 11 and 12 March 2019. It was undertaken by Aerial Surveying Inc., a Hawaii-based company specializing in LiDAR acquisitions. CSRM Foundation oversaw the preliminary data processing, much of which was done by Aerial Surveying Inc., and either oversaw or did all the data processing and analysis that followed.
As anticipated, it proved possible to classify and process the LiDAR data to develop several kinds of useful models. Some models reveal the landscape of Temwen Island and Nan Madol beneath the vegetation that blankets them, while others display the vegetative canopy itself. This was pertinent to the conservation objectives of the LiDAR survey, but it also resulted in the detection of the cultivation features on Temwen Island that are discussed below. In addition to providing the kinds of research opportunities mentioned above, these features are likely relevant to a myriad of archaeological and conservation-related matters, such as the original selection of the construction site for Nan Madol, the socio-economic system in which the islets emerged, and the alteration of Temwen hydrology that has contributed to the siltation of Nan Madol's canals.

Terminology
As the rest of this paper will reinforce, visual inspection of the LiDAR-derived products; fieldwork; and consultation of the literature on Pacific island archaeology, anthropology, and cultivation and agriculture indicate that the features revealed by the LiDAR survey described here were created by people to enhance Temwen Island's capacity for food production through the intensive cultivation of select species. But important details, such as which crops were grown in which areas and what methods were used to cultivate them, have not been ascertained with certainty. For that reason, the terms we use to describe the features revealed by LiDAR are empirically descriptive and nontechnical. For example, it is likely that the earthen berms described in Section 3 were used to cultivate certain crops, while the low, wet areas occurring between such berms were planted with others. Both berms and wet areas could be called fields, but that term has technical connotations we will avoid here.
The word agriculture is especially connotation laden, particularly in archaeology. Though its application to the features discussed here would be defensible, we opt for the arguably more general term cultivation.

Materials and Methods
In this section, we describe the methods of LiDAR data collection and analysis and the field methods used for validating the LiDAR data and ground-truthing the features visible in the LiDAR imagery.

LiDAR Data Collection
The LiDAR data were collected with a (full-waveform) ALIS560 sensor mounted to a Britten-Norman Islander, a small, two-engine, fixed-wing aircraft. The surveyed area ( Figure 2) is 8.69 square kilometers and covers all of Temwen Island and Nan Madol. Mean flight altitude was 457.2 m, mean ground speed was 46.3 m per second, the scan angle was 60 • (±30 • ), and the total swath width was 577 m. There was total coverage between scan lines, making the overlap percentage 100%. The data was output as WGS84 UTM Zone 57N, with elevations measured in meters above mean sea level (MSL), and using EGM2008 as the geoid model. The survey yielded 108.55 million recorded data points or returns. Each return consists of a 3D (x, y, z) spatial measurement, along with several other measurements, including intensity of return and red, green, and blue reflectances, that will be useful for future analyses. The mean number of returns per square meter was 11.93.

LiDAR Data Analysis
The combined point cloud for the entire region of interest was classified into ground and non-ground returns via the LP360 commercial software package, which is commonly used within the LiDAR community. The mean number of returns classified as having been reflected from the ground itself (called ground returns) was 2.7 per square meter. Once the point cloud was classified, it was used to create a Digital Terrain Model (DTM), a bare-earth representation of the ground surface, with 0.5-meter post spacing. This DTM took the form of a georeferenced raster data product (sometimes referred to as a 2.5D product, as opposed to the original point cloud itself, which is 3D; see [28]) with 0.5 × 0.5-meter cells ( Figure 3). LP360 was used to create the DTM as well.
Multiple products derived from the DTM were generated to highlight the archaeological features. We briefly describe seven of them here. The first derived product is a composite shaded relief image composed of 32,310 individual shaded relief products that were generated using a highly parallelized custom application written in C++. Azimuth angles from zero to 359 and elevation angles from 0 to 90, which effectively represent a solar illumination dome over the region of interest, were used to ensure that subtle terrain features that might be visible only under specific geometric conditions could be highlighted. The shaded relief results of each azimuth and elevation angle combination were summed and then divided by the total number of products used, creating a normalized output. This derived product can be generated using standard commercial and open-source GIS software packages through relatively straightforward Python scripting, but the process would likely be very time consuming due to the single-threaded nature of those applications, unless some batch job orchestration system, taking advantage of multiple CPU cores, multiple computers, or both, was used. That is one reason why simpler versions of this product, like the multidirectional hillshade discussed below (Figure 4), are generally created instead.

LiDAR Data Analysis
The combined point cloud for the entire region of interest was classified into ground and non-ground returns via the LP360 commercial software package, which is commonly used within the LiDAR community. The mean number of returns classified as having been reflected from the ground itself (called ground returns) was 2.7 per square meter. Once the point cloud was classified, it was used to create a Digital Terrain Model (DTM), a bare-earth representation of the ground surface, with 0.5-meter post spacing. This DTM took the form of a georeferenced raster data product (sometimes referred to as a 2.5D product, as opposed to the original point cloud itself, which is 3D; see [28]) with 0.5 × 0.5-meter cells ( Figure 3). LP360 was used to create the DTM as well.
Multiple products derived from the DTM were generated to highlight the archaeological features. We briefly describe seven of them here. The first derived product is a composite shaded relief image composed of 32,310 individual shaded relief products that were generated using a highly parallelized custom application written in C++. Azimuth angles from zero to 359 and elevation angles from 0 to 90, which effectively represent a solar illumination dome over the region of interest, were used to ensure that subtle terrain features that might be visible only under specific geometric conditions could be highlighted. The shaded relief results of each azimuth and elevation angle combination were summed and then divided by the total number of products used, creating a normalized output. This derived product can be generated using standard commercial and open-source GIS software packages through relatively straightforward Python scripting, but the process would likely be very time consuming due to the single-threaded nature of those applications, unless some batch job orchestration system, taking advantage of multiple CPU cores, multiple computers, or both, was used. That is one reason why simpler versions of this product, like the multidirectional hillshade discussed below (Figure 4), are generally created instead. The second derived product is a Sky-View Factor (SVF) image [29] that was generated using a highly parallelized custom application written in C++. Instead of using a solar illumination dome (top-down), SVF evaluates how visible each cell in the supplied terrain model would be with respect to its neighbors (bottom-up) and assigns it a score between zero and one. While both composite shaded relief and SVF use the concept of a dome around a point of interest, SVF can sometimes be more effective at highlighting subtle landscape features. The custom application mentioned above is not required to create the product; the developers of SVF provide software on their website [30].
The third derived product is a multidirectional hillshade that was generated using the Hillshade function in ArcGIS Pro. A subset of the Surface functions, the multidirectional hillshade generates grayscale 3D terrain representations that combine light from six different directions, enhancing the visibility of terrain features. This tool proved particularly important for revealing the intricate patterns of flat areas, berms, and water channels on Temwen Island, many of which are subtle, low-relief features.
The fourth derived product is an estimated canopy height raster. It was generated by calculating the differences between the digital surface model (DSM, derived from the so-called first returns) and those provided by the digital terrain model (DTM, which is often called the bare earth model, derived from the ground, or last, returns). It models the estimated height of the vegetation canopy that covers Temwen Island and Nan Madol. As discussed in subsection 4.1 Conservation Objectives, this product is particularly useful for identifying potentially harmful vegetation.
The fifth derived product is a raster of flow directions from each cell to the steepest downslope cell, generated with a traditional eight-direction (D8), pour-point flow modeling algorithm. This product can be represented as a 3D scene. It provides a visual model of the hydrology of Temwen Island and Nan Madol. The second derived product is a Sky-View Factor (SVF) image [29] that was generated using a highly parallelized custom application written in C++. Instead of using a solar illumination dome (top-down), SVF evaluates how visible each cell in the supplied terrain model would be with respect to its neighbors (bottom-up) and assigns it a score between zero and one. While both composite shaded relief and SVF use the concept of a dome around a point of interest, SVF can sometimes be more effective at highlighting subtle landscape features. The custom application mentioned above is not required to create the product; the developers of SVF provide software on their website [30].
The third derived product is a multidirectional hillshade that was generated using the Hillshade function in ArcGIS Pro. A subset of the Surface functions, the multidirectional hillshade generates grayscale 3D terrain representations that combine light from six different directions, enhancing the visibility of terrain features. This tool proved particularly important for revealing the intricate patterns of flat areas, berms, and water channels on Temwen Island, many of which are subtle, low-relief features.
The fourth derived product is an estimated canopy height raster. It was generated by calculating the differences between the digital surface model (DSM, derived from the so-called first returns) and those provided by the digital terrain model (DTM, which is often called the bare earth model, derived from the ground, or last, returns). It models the estimated height of the vegetation canopy that covers Temwen Island and Nan Madol. As discussed in Section 4.1 Conservation Objectives, this product is particularly useful for identifying potentially harmful vegetation.
The fifth derived product is a raster of flow directions from each cell to the steepest downslope cell, generated with a traditional eight-direction (D8), pour-point flow modeling algorithm. This product can be represented as a 3D scene. It provides a visual model of the hydrology of Temwen Island and Nan Madol. The sixth derived product is a stream network created by analyzing the DTM with a Flow Accumulation tool, which reveals the number of upslope cells that flow into each cell. The product included in this paper has a threshold value of 100, presenting a model that only includes cells with more than 100 cells flowing into them. This threshold is a common choice as a default threshold for exploratory stream network analysis. The resulting model highlights real and potential streams and watercourses.
The seventh derived product, based on the DTM, simulates a sea level 1.4 meters above MSL. It shows the impact that the higher sea level would have on low-lying terrain. All values below 1.4 meters MSL in the DTM were set to that value, after which a new SVF image was generated to highlight subtle terrain features. The DTM was then colorized, emphasizing variability between 1.4 and 10 meters MSL, and overlaid on the SVF image.

Identification of Areas of Interest for Field Verification
Fieldwork to verify the presence and ostensible cultivation function of the features revealed by the LiDAR, and to validate and calibrate the LiDAR elevation measurements, was scheduled for 16-29 June 2019. In preparation, six priority areas were selected for inspection ( Figure 5). These areas were selected chiefly on the basis of two criteria: Including particularly evident and representative examples of the features visible in the LiDAR data, and being relatively lightly covered by vegetation. The latter criterion was important to the performance of the Global Navigation Satellite System (GNSS) equipment used for verification and validation, which must have a relatively unobstructed view of the sky to achieve the most precise measurements possible (GNSS systems use global positioning system satellites as well as other constellations of positioning satellites). Temwen itself is almost entirely covered by thick vegetation, but we placed one of our areas of interest, Area 6, within Nan Madol, where some of the islets are relatively clear of vegetation. This allowed us to The sixth derived product is a stream network created by analyzing the DTM with a Flow Accumulation tool, which reveals the number of upslope cells that flow into each cell. The product included in this paper has a threshold value of 100, presenting a model that only includes cells with more than 100 cells flowing into them. This threshold is a common choice as a default threshold for exploratory stream network analysis. The resulting model highlights real and potential streams and watercourses.
The seventh derived product, based on the DTM, simulates a sea level 1.4 m above MSL. It shows the impact that the higher sea level would have on low-lying terrain. All values below 1.4 m MSL in the DTM were set to that value, after which a new SVF image was generated to highlight subtle terrain features. The DTM was then colorized, emphasizing variability between 1.4 and 10 m MSL, and overlaid on the SVF image.

Identification of Areas of Interest for Field Verification
Fieldwork to verify the presence and ostensible cultivation function of the features revealed by the LiDAR, and to validate and calibrate the LiDAR elevation measurements, was scheduled for 16-29 June 2019. In preparation, six priority areas were selected for inspection ( Figure 5). These areas were selected chiefly on the basis of two criteria: Including particularly evident and representative examples of the features visible in the LiDAR data, and being relatively lightly covered by vegetation. The latter criterion was important to the performance of the Global Navigation Satellite System (GNSS) equipment used for verification and validation, which must have a relatively unobstructed view of the sky to achieve the most precise measurements possible (GNSS systems use global positioning system satellites as well as other constellations of positioning satellites). Temwen itself is almost entirely covered by thick vegetation, but we placed one of our areas of interest, Area 6, within Nan Madol, where some of the islets are relatively clear of vegetation. This allowed us to fully exploit the GNSS system's precision and gather highly accurate geospatial measurements (see Section 3.1 for details). For later analyses, many of which will be relevant to the Sustainable Conservation Plan, these data will be quantitatively compared with the LiDAR data, to verify and validate the elevations measured for Nan Madol's islets.
Remote Sens. 2019, 11, x FOR PEER REVIEW 8 of 24 fully exploit the GNSS system's precision and gather highly accurate geospatial measurements (see subsection 3.1 for details). For later analyses, many of which will be relevant to the Sustainable Conservation Plan, these data will be quantitatively compared with the LiDAR data, to verify and validate the elevations measured for Nan Madol's islets. Another area of interest, Area 1, was chosen because it appeared, in the LiDAR data products, to play critical roles in the current hydrology on Temwen Island and downstream in Nan Madol. Areas 2-5 were chosen on the basis of the first criterion mentioned above.

Field Verification
Georeferenced LiDAR-derived images, such as the shaded relief product described in subsection 2.2, were loaded on iPads and used to navigate to the six areas of interest through the ArcGIS Collector application. Preparing these maps for use in the application was done with ArcGIS Online. ArcGIS Collector also allowed the fieldwork team to add point, line, and polygon features to the maps loaded onto the iPads while in the field. Another application, the GPS Camera 55 Field Survey, automatically recorded the location, elevation, time acquired, and direction of each photo taken with the iPads, and it allowed the team to add notes to the photos.
Pohnpei State Historic Preservation Office (HPO) personnel and Temwen Island landowners were informally surveyed for information about archaeological feature locations. HPO personnel used machetes to clear vegetation from parts of the areas of interest. Without clearing vegetation, discerning and identifying the cultivation features under Temwen's profuse plant growth would have been extraordinarily difficult. HPO and Pohnpei State Department of Land personnel identified the owners of the parcels of land where the areas of interest were located. The HPO was instrumental in getting landowner permission to access and clear vegetation from the parcels. The GNSS system was used to collect data not only for validating and calibrating the LiDAR data (in future analyses), but also to confirm and record the locations of features visible in the LiDAR-derived products. Because the system records precise geospatial data in three dimensions, these data were used, and can be used further in the future, to verify and compare not only the Another area of interest, Area 1, was chosen because it appeared, in the LiDAR data products, to play critical roles in the current hydrology on Temwen Island and downstream in Nan Madol. Areas 2-5 were chosen on the basis of the first criterion mentioned above.

Field Verification
Georeferenced LiDAR-derived images, such as the shaded relief product described in Section 2.2, were loaded on iPads and used to navigate to the six areas of interest through the ArcGIS Collector application. Preparing these maps for use in the application was done with ArcGIS Online. ArcGIS Collector also allowed the fieldwork team to add point, line, and polygon features to the maps loaded onto the iPads while in the field. Another application, the GPS Camera 55 Field Survey, automatically recorded the location, elevation, time acquired, and direction of each photo taken with the iPads, and it allowed the team to add notes to the photos.
Pohnpei State Historic Preservation Office (HPO) personnel and Temwen Island landowners were informally surveyed for information about archaeological feature locations. HPO personnel used machetes to clear vegetation from parts of the areas of interest. Without clearing vegetation, discerning and identifying the cultivation features under Temwen's profuse plant growth would have been extraordinarily difficult. HPO and Pohnpei State Department of Land personnel identified the owners of the parcels of land where the areas of interest were located. The HPO was instrumental in getting landowner permission to access and clear vegetation from the parcels. The GNSS system was used to collect data not only for validating and calibrating the LiDAR data (in future analyses), but also to confirm and record the locations of features visible in the LiDAR-derived products. Because the system records precise geospatial data in three dimensions, these data were used, and can be used further in the future, to verify and compare not only the shapes and locations of the Temwen features in horizontal, plan view, but also in terms of their elevations and profiles.

Results
Products derived from the LiDAR data set not only prove themselves useful for preparing the Nan Madol Sustainable Conservation Plan, but also revealed, to all appearances and in accordance with Pohnpei's larger archaeological context, a sophisticated and extensive landscape of cultivation features hidden under Temwen Island's vegetation. Ground verification supported this interpretation. Along with notes and georeferenced photos, GNSS measurements confirmed and recorded the shapes, layouts, and locations of the features.

GNSS Accuracy
Precise cartographic data was gathered using Trimble GNSS receivers supplied by the Spatial Observation Laboratory for Archaeological Research (SOLAR) at Johns Hopkins University. A Trimble R10 receiver was used as a base station and a Trimble R2 as a rover. Differential correction of the rover data was performed using Trimble GPS Pathfinder Office software. All points collected in Area 6, which proved to offer good access to the sky, had accuracy estimates of zero to five centimeters after post-processing. On Temwen Island, where vegetation is far denser, precision was often compromised (Table 1). However, elevations obtained in areas with less tree cover have accuracies comparable to those of the measurements obtained during the LiDAR survey.

LiDAR Analyses and Models
Quantitative analyses of the LiDAR data, which were mostly restricted to statistical characterizations of the distributions of elevation measurements, and visual analysis of the products derived from it yielded several conclusions. Visualizations of the LiDAR data products display the structural integrity of Nan Madol's islets and the structures on them. The islets' and structures' outlines are clear. Many islets have been visibly, profoundly affected by periodic inundations associated with normal and extreme tidal cycles and storm surges. Many structures, even on relatively intact islets, have evidently been damaged not only by inundation, but also by vegetation on the islets and in the canals between them. The siltation of the canals seems to be related to this process ( Figure 6). This phenomenon has been noted previously by several observers on the ground [31,32] (but see also [33]). Certain products make visible the extent and density of mangrove growth in and near Nan Madol.
Remote Sens. 2019, 11, x FOR PEER REVIEW 10 of 24 Figure 6. A multidirectional hillshade overlaying a colorized DTM. Some islets and structures have visibly better structural integrity than others do. This product also reveals the extent and density of mangrove growth (the polygon encompassing an area of mangrove encroachment along Nan Madol's southwestern edge is meant to be expository, not comprehensive). The blue line denotes one of Nan Madol's principal canal systems, which runs past several of the largest and best-preserved islets.
Some of the LiDAR-derived data products, and in particular the inundation-simulation product described in subsection 2.2, make it possible to quantitatively and visually model phenomena that threaten Nan Madol, such as normal and extreme tides and storm surges (Figures 7 and 8). By continuing to work with these models and to use LiDAR elevations calibrated with the precise elevation measures collected with the GNSS system, these models can be made more accurate. With more accurate models, better assessments of conservation threats associated with wind and water can be made.
Data products that include canopy height calculations allow for the mapping and easy identification of high vegetation, which poses special threats to Nan Madol, such as falling branches and trees (Figure 9). These products allow us to identify potentially harmful species, and even especially threatening individual plants. Figure 6. A multidirectional hillshade overlaying a colorized DTM. Some islets and structures have visibly better structural integrity than others do. This product also reveals the extent and density of mangrove growth (the polygon encompassing an area of mangrove encroachment along Nan Madol's southwestern edge is meant to be expository, not comprehensive). The blue line denotes one of Nan Madol's principal canal systems, which runs past several of the largest and best-preserved islets.
Some of the LiDAR-derived data products, and in particular the inundation-simulation product described in Section 2.2, make it possible to quantitatively and visually model phenomena that threaten Nan Madol, such as normal and extreme tides and storm surges (Figures 7 and 8). By continuing to work with these models and to use LiDAR elevations calibrated with the precise elevation measures collected with the GNSS system, these models can be made more accurate. With more accurate models, better assessments of conservation threats associated with wind and water can be made.
Data products that include canopy height calculations allow for the mapping and easy identification of high vegetation, which poses special threats to Nan Madol, such as falling branches and trees ( Figure 9). These products allow us to identify potentially harmful species, and even especially threatening individual plants.     Estimated canopy height and density at Nan Dowas. Note that the model shows three trees with heights between 17.8 and 28.1 meters. As height is an important selection criterion for vegetation removal, models like this one will be important tools for vegetation management.

Field Verification of LiDAR-Derived Models
Other data products make the complex system of cultivable areas and water channels across Temwen Island clearly visible (Figures 10 and 11). In addition to simply revealing these features, which are normally hidden from satellite, aerial, and ground-level view by thick vegetation, these products make immediate suggestions about how the cultivation and hydrological systems on Temwen Island functioned and how they might relate to Nan Madol.
Inspection and verification on the ground confirmed the presence of low areas, separated by earthen berms and connected to one another by channels, that allow water to collect in them. The channels that connect these areas allow water to flow from one area to the next. Because they collect water, these low areas are conducive to the growth of certain plants: In several of the areas inspected on the ground, vegetation in the low, wet areas and water channels was distinct, in terms of density and species present, from that on the earthen berms that separate them. This suggests their potential use in the past as wet fields or pondfields for crops such as taro. Water continues to flow, at least after heavy rain (which is common on Pohnpei), through some of these low areas and the channels that connect them. The berms separating these low areas and crossed by the channels may have been cultivated as well, perhaps with cultigens that prefer drier conditions. Such cultivation on raised mounds or platforms has been noted on Pohnpei and elsewhere in Oceania [34,35].
Low areas, water channels, and berms separating them are found on both flat and sloped areas. Systems of low areas, channels, and berms on the sloped areas manifest in the LiDAR-derived products as terrace-like. By clearing vegetation-for example, along a 70-meter transect in western Temwen ( Figure 10)-these terrace-like systems, proceeding downslope and transitioning into systems of parallel, alternating low areas and berms on flat ground, were revealed. Figure 9. Estimated canopy height and density at Nan Dowas. Note that the model shows three trees with heights between 17.8 and 28.1 m. As height is an important selection criterion for vegetation removal, models like this one will be important tools for vegetation management.

Field Verification of LiDAR-Derived Models
Other data products make the complex system of cultivable areas and water channels across Temwen Island clearly visible (Figures 10 and 11). In addition to simply revealing these features, which are normally hidden from satellite, aerial, and ground-level view by thick vegetation, these products make immediate suggestions about how the cultivation and hydrological systems on Temwen Island functioned and how they might relate to Nan Madol.
Inspection and verification on the ground confirmed the presence of low areas, separated by earthen berms and connected to one another by channels, that allow water to collect in them. The channels that connect these areas allow water to flow from one area to the next. Because they collect water, these low areas are conducive to the growth of certain plants: In several of the areas inspected on the ground, vegetation in the low, wet areas and water channels was distinct, in terms of density and species present, from that on the earthen berms that separate them. This suggests their potential use in the past as wet fields or pondfields for crops such as taro. Water continues to flow, at least after heavy rain (which is common on Pohnpei), through some of these low areas and the channels that connect them. The berms separating these low areas and crossed by the channels may have been cultivated as well, perhaps with cultigens that prefer drier conditions. Such cultivation on raised mounds or platforms has been noted on Pohnpei and elsewhere in Oceania [34,35].
Low areas, water channels, and berms separating them are found on both flat and sloped areas. Systems of low areas, channels, and berms on the sloped areas manifest in the LiDAR-derived products as terrace-like. By clearing vegetation-for example, along a 70-meter transect in western Temwen ( Figure 10)-these terrace-like systems, proceeding downslope and transitioning into systems of parallel, alternating low areas and berms on flat ground, were revealed. Remote Sens. 2019, 11, x FOR PEER REVIEW 13 of 24 Figure 10. Large-scale landscape modifications in western Temwen Island. In (a), a north-facing profile of the terrace-like features, sloped areas, low areas, and berm that are depicted in the plan view in (b). In (c), the LiDAR-derived multidirectional hillshade, overlaid on the DTM, showing the locations of the profile and plan view. A day after heavy rain on Pohnpei, water was observed flowing through roughly west-east water channels, down the slope of the ridge in western Temwen. As shown in (b), the channels bisect terrace-like features. The channels run into low areas that alternate with berms on flatter ground below the ridge (the alternating pattern continues eastward). The low areas slope gently northward, bringing the water past a series of berms and crossing other west-east channels. If the terrace-like features were originally designed to draw water along their long axes, such water-drawing mechanisms are now defunct.
Fieldwork confirmed that the earthen berms, terrace-like features, and many of the lower areas below and around them are roughly rectangular and elongated in plan. At their highest points, the earthen berms that were observed on the ground are approximately 0.5-1.25 meters high. The observed terrace-like features are approximately 0.25-1.25 meters above the next terrace-like feature or low area. With one possible exception, described below, observed berms and terrace-like features have no retaining walls or other non-earthen support. The margins between them and the lower areas around them, and between one terrace-like feature and the next, are sloped. Aside from these slopes, many of the berms and low areas are level and spacious: One especially large berm in north-central Temwen, not surveyed on the ground, can be seen in the LiDAR to cover approximately 2800 square meters, or about 0.7 acres. In (c), the LiDAR-derived multidirectional hillshade, overlaid on the DTM, showing the locations of the profile and plan view. A day after heavy rain on Pohnpei, water was observed flowing through roughly west-east water channels, down the slope of the ridge in western Temwen. As shown in (b), the channels bisect terrace-like features. The channels run into low areas that alternate with berms on flatter ground below the ridge (the alternating pattern continues eastward). The low areas slope gently northward, bringing the water past a series of berms and crossing other west-east channels. If the terrace-like features were originally designed to draw water along their long axes, such water-drawing mechanisms are now defunct.
Fieldwork confirmed that the earthen berms, terrace-like features, and many of the lower areas below and around them are roughly rectangular and elongated in plan. At their highest points, the earthen berms that were observed on the ground are approximately 0.5-1.25 m high. The observed terrace-like features are approximately 0.25-1.25 m above the next terrace-like feature or low area. With one possible exception, described below, observed berms and terrace-like features have no retaining walls or other non-earthen support. The margins between them and the lower areas around them, and between one terrace-like feature and the next, are sloped. Aside from these slopes, many of the berms and low areas are level and spacious: One especially large berm in north-central Temwen, not surveyed on the ground, can be seen in the LiDAR to cover approximately 2800 square meters, or about 0.7 acres. LiDAR-derived flow direction ( Figure 12) and flow accumulation models ( Figure 13) of Temwen Island show how the cultivation features influence the landscape's hydrology. They suggest how the Temwen water management system might have functioned while it was in use, how Temwen's overall topography and hydrology formed a basis for the cultivation system, and how the system interfaced with Nan Madol. Field observations verified that the earthen berms between low areas act as retaining walls for water, altering the natural flow directions that would be expected from the underlying topography. Some channels were observed bringing water from low areas at higher elevations to low areas at lower elevations ( Figure 10). Some channels flow past several low areas along their course, suggesting a spatial and water management scheme integrating many areas of cultivation.
In the southwestern part of the island (in Area 4), a wall of basalt boulders, with a similar but smaller wall built on the retained earth behind (roughly east of) it, was found. These were the only potentially constructed stone features found during field verification. As detectable in the LiDAR data, the larger of the two is about 29 meters long. Downslope (roughly west) of these walls were more earthen terraces and channels. These walls could be naturally formed, by volcanic activity, but they appear anthropogenic, and monumental structures are not uncommon on Pohnpei, even outside Nan Madol. Throughout Pohnpei, retaining walls for terraces are typically rock faced (for example, in Awak and Wene; see [36,37]).

Discussion
This section discusses the relevance of the LiDAR-derived products to the ongoing conservation planning and activity at Nan Madol. It proceeds to address some of the many implications of the discovery of the Temwen Island cultivation system for understanding and research at Nan Madol, on Pohnpei, and in the region. Current knowledge and ongoing inquiry about the site and its local and regional contexts will be reflected in an interpretive plan, which will constitute part of the larger Sustainable Conservation Plan. LiDAR-derived flow direction ( Figure 12) and flow accumulation models ( Figure 13) of Temwen Island show how the cultivation features influence the landscape's hydrology. They suggest how the Temwen water management system might have functioned while it was in use, how Temwen's overall topography and hydrology formed a basis for the cultivation system, and how the system interfaced with Nan Madol. Field observations verified that the earthen berms between low areas act as retaining walls for water, altering the natural flow directions that would be expected from the underlying topography. Some channels were observed bringing water from low areas at higher elevations to low areas at lower elevations ( Figure 10). Some channels flow past several low areas along their course, suggesting a spatial and water management scheme integrating many areas of cultivation.
In the southwestern part of the island (in Area 4), a wall of basalt boulders, with a similar but smaller wall built on the retained earth behind (roughly east of) it, was found. These were the only potentially constructed stone features found during field verification. As detectable in the LiDAR data, the larger of the two is about 29 m long. Downslope (roughly west) of these walls were more earthen terraces and channels. These walls could be naturally formed, by volcanic activity, but they appear anthropogenic, and monumental structures are not uncommon on Pohnpei, even outside Nan Madol. Throughout Pohnpei, retaining walls for terraces are typically rock faced (for example, in Awak and Wene; see [36,37]).

Conservation Objectives
As mentioned in subsection 1.2, salient conservation threats to Nan Madol include disturbance to islets and structures by wind and water, the siltation of the site's canals, and vegetation overgrowth. The LiDAR data set presented here will help design measures to counter each of these.

Stabilization of Islets and Structures
Imagery and other data products that were generated from the LiDAR data set will be used to prioritize and guide physical interventions for stabilizing Nan Madol's islets and structures. Models derived from the data set provide convenient and quantifiable means for assessing the physical integrity of archaeological features. The use of the LiDAR's precise elevation measurements to assess and forecast vulnerability to sea level change and acute threats such as storm surges is an example (Figures 7 and 8).
Prior archaeological fieldwork suggests reasons for the differential preservation of Nan Madol islets that is apparent in the LiDAR models. The findings from this fieldwork indicate that the human occupation of areas offshore of Temwen Island predated the artificial islets and monumental structures that comprise Nan Madol [5,6,8,23]. Nan Madol's most intact, and many of its most notable, islets and structures lie along a southwest-northeast axis. These are the islets most visible in Figure 6. Athens, drawing upon his excavations and surveys at the islets [3,22,23], posits the presence of natural islets along that axis that predated Nan Madol's islets. He suggests that the augmentation of these natural islets with basalt and coral rubble inspired the construction of other, later islets both nearer to shore and farther seaward ( [5], p. 21).
Later survey and excavation by Ayres of islets closer to the Temwen Island shore show-on the basis of pottery content and architectural form-that those islets are older than ones closer to the outer reef edge. In the inland islets, very little columnar rock is used in construction. Ayres and geologist G. Goles confirmed that most of the construction boulders used in the inland islets were from the adjacent Temwen lava flows (see [38]). Also, at the islet called Dauahpeidak, for example, a concentration of pottery, shell ornaments and food remains was found in a reef sand deposit buried beneath the islet construction. Many of the sherds were water-rolled, indicating either marine reworking of a sand reef islet or materials left as a result of a reef-flat occupation in houses on pilings. This suggested that dwellings on the reef flat predated the development of major artificial islets. This would be consistent with earlier Lapita and post-Lapita residential locations known from Melanesia. Between the idea of marine reworking of a sand reef islet and that of reef-flat occupation in houses on pilings, Ayres favored the latter as a working hypothesis. In any case, he did not think that the use of basalt boulders to construct the artificial islet foundations was necessarily linked to sea level change [6,7]. Columnar basalt was used as facing on some of the earlier islets, but might have been placed later to enhance the islets' prestige. Athens [23] and others [6,7,13,21,39,40] tell us that Nan Madol's islets differ greatly in construction technique, features, pattern and shape, and the oral accounts associated with them. In Figure 6, the islets along the southwest-northeast axis display the greatest structural integrity and are at higher elevations above MSL than are those to the landward and seaward. These islets have been identified by numerous ethnographers, ethnohistorians, and archaeologists as the most significant with regard to the roles they played in the Saudeleur socio-political system (see e.g., [16,20,41,42]). That these islets are relatively well preserved is consistent with the construction scenarios offered by Athens and Ayres. Such assessments assist in supporting decisions about prioritizing islets and structures for conservation treatments. These treatments could include vegetation removal and physical stabilization of structures. Erosion has caused deflation of both inland and seaward islets, and most of the landward islets have been damaged by siltation and tidal flooding [39,40]. These damaging influences make an exact count of islets difficult; only one new islet has been recorded since the 1970s.

Discussion
This section discusses the relevance of the LiDAR-derived products to the ongoing conservation planning and activity at Nan Madol. It proceeds to address some of the many implications of the discovery of the Temwen Island cultivation system for understanding and research at Nan Madol, on Pohnpei, and in the region. Current knowledge and ongoing inquiry about the site and its local and regional contexts will be reflected in an interpretive plan, which will constitute part of the larger Sustainable Conservation Plan.

Conservation Objectives
As mentioned in Section 1.2, salient conservation threats to Nan Madol include disturbance to islets and structures by wind and water, the siltation of the site's canals, and vegetation overgrowth. The LiDAR data set presented here will help design measures to counter each of these.

Stabilization of Islets and Structures
Imagery and other data products that were generated from the LiDAR data set will be used to prioritize and guide physical interventions for stabilizing Nan Madol's islets and structures. Models derived from the data set provide convenient and quantifiable means for assessing the physical integrity of archaeological features. The use of the LiDAR's precise elevation measurements to assess and forecast vulnerability to sea level change and acute threats such as storm surges is an example (Figures 7 and 8).
Prior archaeological fieldwork suggests reasons for the differential preservation of Nan Madol islets that is apparent in the LiDAR models. The findings from this fieldwork indicate that the human occupation of areas offshore of Temwen Island predated the artificial islets and monumental structures that comprise Nan Madol [5,6,8,23]. Nan Madol's most intact, and many of its most notable, islets and structures lie along a southwest-northeast axis. These are the islets most visible in Figure 6. Athens, drawing upon his excavations and surveys at the islets [3,22,23], posits the presence of natural islets along that axis that predated Nan Madol's islets. He suggests that the augmentation of these natural islets with basalt and coral rubble inspired the construction of other, later islets both nearer to shore and farther seaward ( [5], p. 21).
Later survey and excavation by Ayres of islets closer to the Temwen Island shore show-on the basis of pottery content and architectural form-that those islets are older than ones closer to the outer reef edge. In the inland islets, very little columnar rock is used in construction. Ayres and geologist G. Goles confirmed that most of the construction boulders used in the inland islets were from the adjacent Temwen lava flows (see [38]). Also, at the islet called Dauahpeidak, for example, a concentration of pottery, shell ornaments and food remains was found in a reef sand deposit buried beneath the islet construction. Many of the sherds were water-rolled, indicating either marine reworking of a sand reef islet or materials left as a result of a reef-flat occupation in houses on pilings. This suggested that dwellings on the reef flat predated the development of major artificial islets. This would be consistent with earlier Lapita and post-Lapita residential locations known from Melanesia. Between the idea of marine reworking of a sand reef islet and that of reef-flat occupation in houses on pilings, Ayres favored the latter as a working hypothesis. In any case, he did not think that the use of basalt boulders to construct the artificial islet foundations was necessarily linked to sea level change [6,7]. Columnar basalt was used as facing on some of the earlier islets, but might have been placed later to enhance the islets' prestige.
Athens [23] and others [6,7,13,21,39,40] tell us that Nan Madol's islets differ greatly in construction technique, features, pattern and shape, and the oral accounts associated with them. In Figure 6, the islets along the southwest-northeast axis display the greatest structural integrity and are at higher elevations above MSL than are those to the landward and seaward. These islets have been identified by numerous ethnographers, ethnohistorians, and archaeologists as the most significant with regard to the roles they played in the Saudeleur socio-political system (see e.g., [16,20,41,42]). That these islets are relatively well preserved is consistent with the construction scenarios offered by Athens and Ayres. Such assessments assist in supporting decisions about prioritizing islets and structures for conservation treatments. These treatments could include vegetation removal and physical stabilization of structures. Erosion has caused deflation of both inland and seaward islets, and most of the landward islets have been damaged by siltation and tidal flooding [39,40]. These damaging influences make an exact count of islets difficult; only one new islet has been recorded since the 1970s.

Hydrology and Siltation of Nan Madol's Canals
The development of the landscape of cultivation features on Temwen Island altered the island's hydrology and the flow of water into and through the canals at Nan Madol. As a result, rainfall runoff into Nan Madol may be lower than what would be expected on a natural Temwen Island landscape. This may have been the case even while the cultivation system was being maintained. At present, there can be little doubt that water flow through the network of water channels is blocked because of the degradation of the network, due to causes such as heavy vegetation on Temwen Island, historic landscape modifications, and ground disturbance by pigs.
With physical treatments such as vegetation clearance and targeted maintenance of the water management system, it may be possible to increase the flow of runoff into and through Nan Madol's canals. This would offer a minimally invasive way of clearing silt, which is linked with problematic vegetation overgrowth, from the canals. The LiDAR data set will continue to be used to model the hydrology and vegetation of Temwen Island and Nan Madol to assess the feasibility of testing such treatments.

Vegetation Management
LiDAR models show the location and density of mangroves ( Figure 6). These models will be used to guide the clearance of mangroves from selected areas. Priority areas include the canals that skirt the most notable of Nan Madol's islets. Ideally, mangroves will not be completely removed from these canals but allowed to remain in place along the canals' edges, much as they often do along the banks of natural rivers. Removing mangroves completely could destabilize the walls of the islets along which the canals flow. Mangroves also provide some protection from tidal cycles and storm surges and serve important ecological roles. Allowing the mangroves to grow unchecked, however, impedes water flow through the canals, encouraging siltation and further overgrowth, as well as making the canals impassable.
Other LiDAR models have been created to highlight trees likely to fall or from which limbs might drop ( Figure 9). Falling trees and limbs have damaged structures and the islets themselves at Nan Madol. These models will be used during the initial removal of threatening vegetation.
The LiDAR imagery shows a large area in the ravine on the south side of Temwen Island that has been used in recent decades as a field for growing swamp taro. This large field is now almost completely abandoned; local residents say that many members of the family that once grew swamp taro there have left Pohnpei. The remaining few members now grow taro only in an area adjacent to the house in which they live. Clearing this area of vegetation might assist efforts to further document the Temwen cultivation system as well as to increase water flow into Nan Madol.

Implications and Interpretation Objectives
LiDAR has revealed a massive and complex cultivation system on Temwen Island. It is this advantage of LiDAR, that it can reveal whole archaeological landscapes hidden under heavy vegetation, that has led to its comparison with radiocarbon dating as a groundbreaking technological advance in archaeology. By virtue of its abilities to provide data sets that can reveal an enormous range of archaeologically relevant features and are, geospatially, not merely representative but nearly complete, the importance of LiDAR is difficult to overstate (see [43,44]). Further LiDAR surveys, not only of other parts of Pohnpei but of other Pacific islands, could contextualize the Temwen system and offer new synoptic perspectives.
In plan, the low, water-collecting areas; water channels; and berms on Temwen are similar to the systems of terraces described by Sand in New Caledonia, which he writes are "considered as the most elaborate of the Western Pacific" ( [45], p. 167). Within that system, associated with the Traditional Kanak Cultural Complex from the second millennium CE, he documents terraces, water channels, and pondfield sites. In describing the construction and use of these systems, Sand cites Captain James Cook's observation that "not an Inch [sic] of ground is lost" ( [46], p. 538). The system has been described as "certainly the most technically complex pondfield irrigation within Melanesia" ( [47], p. 80; cited in [45], p. 167). Such complexity is in accord with what is observed in the LiDAR images of Temwen Island. The Temwen system also bears a resemblance to some Polynesian terracing systems, including the Kohala field system on Hawai'i Island [48] and recently described slope terraces on Tutuila in American Samoa [49].
As mentioned in Section 3, retaining walls for terraces on Pohnpei are typically rock faced [37]. That Temwen's cultivation features are almost entirely earthen is therefore unusual. This is perhaps explained in part by the paucity of surface rock on Temwen Island. There is a notable exception about half a kilometer to the east of the abovementioned walls, where a profusion of basalt boulders lies on and near an approximately south-facing slope down to the lagoon. (These boulders could have been rolled down the slope into Nan Madol.) The basalt walls are at the head of a ravine that runs down to the northwestern edge of Nan Madol, the same ravine mentioned in Section 4.1.3. This ravine might have brought fresh water to a point at which it could have been relatively easily collected by the site's inhabitants. The flow of water through the ravine and into the site might also have encouraged water flow through the site itself, potentially clearing sediment from the canals ( Figure 13).
The LiDAR models presented here provide the first view of an intensive cultivation system, including water control features, of such scale in the Central and Eastern Caroline Islands. This finding is relevant to ongoing discussions of the development of social complexity throughout the Pacific. Previous localized archaeological surveys and excavations on Pohnpei documented terraces and alignments that were interpreted as agricultural features [34,[50][51][52][53][54], some with water control mechanisms (see [34], p. 151). Yet a large, complex system with integrated water channels had not been posited.
Without the view provided by LiDAR, archaeologists who have done research on Pohnpei, and in particular at Nan Madol, have generally emphasized the importance of breadfruit (see e.g., [55], p. 85), yams, and reef fish, with taro contributing little to the Pohnpeian diet. Though the role of bananas is acknowledged, it is understudied, likely because of their lack of prestige in the traditional food system or their limited archaeological visibility [56]. None of these foods can be effectively preserved once harvested, with the exception of breadfruit, and less commonly banana, which can be fermented [34,[55][56][57]. This makes year-round cultivation and harvest of crops necessary. Food security has long been seen as having shaped Pohnpeian culture (and Micronesian culture more broadly) in several defining ways, including by encouraging the improvement of inter-island connections [53,58] and social systems that emphasize the production of large amounts of food for feasting [18,54,59]. But food production on a scale made possible by the cultivation system described here has not been previously suggested. The Temwen cultivation system can be considered landesque capital in the sense used by Håkansson and Widgren: A choice of agricultural techniques that serve important economic, social, and ritual purposes [60]. The term was originally used by Sen to refer simply to capital goods that replace land (as opposed to laboresque capital, which replaces labor) [61]. In this strictly economic sense, too, the Temwen system fits the concept: Spriggs maintains that irrigated systems for growing taro can generate yields at least four times higher than what is possible with swidden agriculture ( [62], cited in [63]).
The concept of landesque capital has been used by geographers and anthropologists to explore food storage as it relates to exchange in the service of food security and the prestige economy. Under certain conditions, for example, concentrations of landesque capital could have generated tradable surpluses of foods such as taro [64,65]. Such tradable surpluses would have constituted wealth, or, as Bayliss-Smith and Hviding put it, "social and nutritional value to be saved, accumulated, and later mobilized" ( [64], p. 425).
Without extensive archaeological fieldwork or the application of technologies like LiDAR, it can be easy to overlook the physical evidence even of large investments in landesque capital. In 1568, the first Europeans to arrive at the island of Guadalcanal, in the Solomon Islands, saw a planned agricultural landscape. In 1610, a Captain de Queirós, on his way to Vanuatu, wrote to the King of Spain that he had observed gardens "artificially severed (divided) into beds, bordered, and paled (fenced)" ( [66], p. 221; quoted in [63], p. 75). Yet when Europeans next recorded what they saw on Guadalcanal in the nineteenth century, there was little trace of terracing or irrigation [63]. It was only in 2000 that an archaeologist, David Roe, discovered remnants of stone walls, terraces and irrigated fields [67]. On Vanuatu, Bayliss-Smith and Hviding have examined historical documentation by de Queirós that engineered landscapes for taro production existed in the Tavai, Apuna, and Raovi valleys and their inland catchments [63].
The existence of the Temwen system supports the hypothesis of Ayres et al. that Temwen Island served as an agricultural area to support Nan Madol, and that evidence of agricultural intensification would be found there ( [40], p. 9). To investigate this, Ayres and colleagues conducted fieldwork in 2011 that documented several features, such as stone walls that were suspected to have been an element of terraces. They report finding 13 boulder alignments thought to be related to terracing, seven stone platforms or enclosures, and four stone walls, as well as other agricultural features, such as yam enclosures and breadfruit fermentation pits ( [40], p. 42; [52], p. 60).
These findings were preceded by others that are consistent with the interpretation of the Temwen Island features as relating to cultivation and possibly coeval with Nan Madol's fluorescence. Ayres and Haun found water control features, terraces, enclosures, pits, and walls in Pohnpei's Awak Valley [36] (see also [37]). Many of these features are agricultural. Among them are terraces along stream courses that the researchers suggested might be wet terraces for taro cultivation.

Conclusions
The LiDAR survey discussed here was done to advance the goals of the U.S. Department of State-funded Sustainable Conservation Plan for Nan Madol. Chief among these goals is to ameliorate threats to the Nan Madol World Heritage Site so that the site can be removed from the List of World Heritage in Danger. The discovery of an extensive and sophisticated cultivation system influences decisions that must be made about conservation treatments and a site management plan. This management plan will provide for site monitoring, cyclical maintenance, interpretation, facilities, visitor services, and site and visitor protection.
It is likely that the Temwen system revealed by LiDAR is contemporaneous with the islets and structures of the Nan Madol World Heritage Site. The system is comparable to the islets and structures of Nan Madol in terms of the complexity and magnitude of the effort that would have been required to build and maintain it. The investment of labor, for example, and the level of social organization necessary to construct and maintain the Temwen system would appear to be similar to what would have been required to construct and maintain Nan Madol. If the Temwen system is a creation of the Saudeleur Dynasty, it emphasizes the high level of social organization associated with the dynasty; regardless, it suggests many avenues for future research.
The LiDAR models presented here comprise the most convincing evidence thus far for a formal, complex cultivation system incorporating an intensive network of growing areas and water management features in the Central and Eastern Carolines. This is the region of Micronesia thought to have been initially settled by Lapita-descendant communities. Further investigations on Temwen Island, and comparative studies, will be needed to determine how this development might inform the chronology and pattern of human interaction among the islands in the region.
Dateable samples from the field system might prove informative not only regarding how the development of the field system and that of Nan Madol's islets and structures compare chronologically, but also where a model for the field system might have originated. Athens states that the pottery found at Nan Madol is best described as belonging to the Late Lapita Plain Ware pottery tradition, which he considered to be unlike pottery found elsewhere in Micronesia, based on shape and decorative technique: A simple globular form for pots and bowls, "with perhaps some rim notching or punctuation" ( [5], p. 29). He sees this form and decoration persisting over the centuries, despite changes in temper (p. 26). Athens further states that Late Lapita Plain Ware is found only in the southeast Solomon and New Hebrides islands (the latter now Vanuatu; see also [24]). Lebot's observation that the kava sub-species found on Pohnpei is most similar to that on Vanuatu and the Bismarcks constitutes another link between Pohnpei and that region [68].
Comparative studies in the future could also address the spread of stone architecture through the Pacific. As Carson notes ( [69], p. 2), "the stone city of Nan Madol in Pohnpei has long been regarded as the 'Venice of the Pacific,' with its canals leading through a labyrinthine network of long-abandoned stonework architecture, but rarely is this discussed in relation to the other monumental stonework traditions that developed throughout Micronesia and the Pacific all about the same time." McCoy and colleagues have presented evidence that Nan Madol's stone architecture is among the Pacific's oldest ( [2], p. 295). It remains to be determined how Nan Madol might be related to the spread of Pacific stone architecture, and how the discovered landscape of cultivation features might further this line of inquiry.
The cultivation system reported on here could be an important reference for future research on the development of such systems in the Pacific. Though oral histories indicate that Nan Madol received food as tribute from across Pohnpei, the resident population could not have been adequately fed by food grown or obtained on site. The site itself could not have produced enough food for food-focused ritual and feasting. To the degree that it provided food, and food security, to Nan Madol, the Temwen cultivation system carries enormous implications for the site.
The Temwen system would appear from our inspection on the ground to have formed wet fields suitable for growing taro. Soil testing, including analyses of pollen samples, phytoliths, and plant macroremains, is likely needed to confirm this. Devoting a substantial proportion of the Temwen system to taro would have resulted in an enormous and relatively secure store of food. Such a supply of taro could have been tapped not only for redistribution or trade [64], but also for feasting that enhanced the prestige of the rulers. If a suitable method can be devised to preserve or store it, food constitutes a form of wealth: Food is essential and, contingent on the preservation or storage method, it can be reapplied or reallocated. Under certain circumstances, it can be traded for, among other desirable materials and objects, prestige goods (some of which are themselves foods).
A deeper understanding of the Temwen cultivation system and its relation to Nan Madol requires research on how the system functioned: How the features were built or maintained, and how they harmonized with cultivation methods and environmental dynamics. This calls for studies such as sampling and analyzing the construction of the features, investigating the chronology of the Temwen system's development and how it relates to Nan Madol's growth, and comparing Temwen with early field systems in other parts of the Pacific. Comparison or analogy with traditional Pohnpeian cultivation methods, following Carson's method in Samoa [70], could also be informative.
Detailed knowledge of the Temwen system will allow future researchers to draw conclusions about Nan Madol's food security and wealth; how these advantages contributed to the site's development, construction, and maintenance; and the economic and political strategies of Nan Madol's rulers. With its socio-political, ritual, and economic importance, Nan Madol may have served as a symbol of affluence and stability to other islands and atolls in the region. As such, in times of need, such as those triggered by typhoons or droughts, Nan Madol may have been looked to by regional neighbors or clan members from elsewhere in Micronesia for aid or trade. This idea is supported by Petersen's argument that the presence of members of clans on numerous islands and atolls across Micronesia created a network of support and trade for voyagers and clan members in need [58]. If such a model holds, then the food security and wealth that supported Nan Madol would have made it an important node at the center of regional social and economic networks.