Pacific Island Ecosystems Research Center

MEPS Marine Ecology Progress Series Contact the journal Facebook Twitter RSS Mailing List Subscribe to our mailing list via Mailchimp HomeLatest VolumeAbout the JournalEditorsTheme Sections MEPS 210:223-253 (2001) - doi:10.3354/meps210223 Our evolving conceptual model of the coastal eutrophication problem James E. Cloern* US Geological Survey, MS496, 345 Middlefield Rd., Menlo Park, California 94025, USA *E-mail: jecloern@usgs.gov ABSTRACT: A primary focus of coastal science during the past 3 decades has been the question: How does anthropogenic nutrient enrichment cause change in the structure or function of nearshore coastal ecosystems? This theme of environmental science is recent, so our conceptual model of the coastal eutrophication problem continues to change rapidly. In this review, I suggest that the early (Phase I) conceptual model was strongly influenced by limnologists, who began intense study of lake eutrophication by the 1960s. The Phase I model emphasized changing nutrient input as a signal, and responses to that signal as increased phytoplankton biomass and primary production, decomposition of phytoplankton-derived organic matter, and enhanced depletion of oxygen from bottom waters. Coastal research in recent decades has identified key differences in the responses of lakes and coastal-estuarine ecosystems to nutrient enrichment. The contemporary (Phase II) conceptual model reflects those differences and includes explicit recognition of (1) system-specific attributes that act as a filter to modulate the responses to enrichment (leading to large differences among estuarine-coastal systems in their sensitivity to nutrient enrichment); and (2) a complex suite of direct and indirect responses including linked changes in: water transparency, distribution of vascular plants and biomass of macroalgae, sediment biogeochemistry and nutrient cycling, nutrient ratios and their regulation of phytoplankton community composition, frequency of toxic/harmful algal blooms, habitat quality for metazoans, reproduction/growth/survival of pelagic and benthic invertebrates, and subtle changes such as shifts in the seasonality of ecosystem functions. Each aspect of the Phase II model is illustrated here with examples from coastal ecosystems around the world. In the last section of this review I present one vision of the next (Phase III) stage in the evolution of our conceptual model, organized around 5 questions that will guide coastal science in the early 21st century: (1) How do system-specific attributes constrain or amplify the responses of coastal ecosystems to nutrient enrichment? (2) How does nutrient enrichment interact with other stressors (toxic contaminants, fishing harvest, aquaculture, nonindigenous species, habitat loss, climate change, hydrologic manipulations) to change coastal ecosystems? (3) How are responses to multiple stressors linked? (4) How does human-induced change in the coastal zone impact the Earth system as habitat for humanity and other species? (5) How can a deeper scientific understanding of the coastal eutrophication problem be applied to develop tools for building strategies at ecosystem restoration or rehabilitation? KEY WORDS: Coastal eutrophication · Nutrient enrichment · Estuaries · Nitrogen · Phosphorus · Coastal ecosystems Full text in pdf format PreviousNextExport citation RSS - Facebook - Tweet - linkedIn Cited by Published in MEPS Vol. 210. Online publication date: January 26, 2001 Print ISSN: 0171-8630; Online ISSN: 1616-1599 Copyright © 2001 Inter-Research.

Future drought is projected to occur under warmer temperature conditions as climate change progresses, referred to here as global-change-type drought, yet quantitative assessments of the triggers and potential extent of drought-induced vegetation die-off remain pivotal uncertainties in assessing climate-change impacts. Of particular concern is regional-scale mortality of overstory trees, which rapidly alters ecosystem type, associated ecosystem properties, and land surface conditions for decades. Here, we quantify regional-scale vegetation die-off across southwestern North American woodlands in 2002-2003 in response to drought and associated bark beetle infestations. At an intensively studied site within the region, we quantified that after 15 months of depleted soil water content, >90% of the dominant, overstory tree species (Pinus edulis, a piñon) died. The die-off was reflected in changes in a remotely sensed index of vegetation greenness (Normalized Difference Vegetation Index), not only at the intensively studied site but also across the region, extending over 12,000 km2 or more; aerial and field surveys confirmed the general extent of the die-off. Notably, the recent drought was warmer than the previous subcontinental drought of the 1950s. The limited, available observations suggest that die-off from the recent drought was more extensive than that from the previous drought, extending into wetter sites within the tree species' distribution. Our results quantify a trigger leading to rapid, drought-induced die-off of overstory woody plants at subcontinental scale and highlight the potential for such die-off to be more severe and extensive for future global-change-type drought under warmer conditions.

We have produced the first 30 m resolution global land-cover maps using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data. We have classified over 6600 scenes of Landsat TM data after 2006, and over 2300 scenes of Landsat TM and ETM+ data before 2006, all selected from the green season. These images cover most of the world's land surface except Antarctica and Greenland. Most of these images came from the United States Geological Survey in level L1T (orthorectified). Four classifiers that were freely available were employed, including the conventional maximum likelihood classifier (MLC), J4.8 decision tree classifier, Random Forest (RF) classifier and support vector machine (SVM) classifier. A total of 91,433 training samples were collected by traversing each scene and finding the most representative and homogeneous samples. A total of 38,664 test samples were collected at preset, fixed locations based on a globally systematic unaligned sampling strategy. Two software tools, Global Analyst and Global Mapper developed by extending the functionality of Google Earth, were used in developing the training and test sample databases by referencing the Moderate Resolution Imaging Spectroradiometer enhanced vegetation index (MODIS EVI) time series for 2010 and high resolution images from Google Earth. A unique land-cover classification system was developed that can be crosswalked to the existing United Nations Food and Agriculture Organization (FAO) land-cover classification system as well as the International Geosphere-Biosphere Programme (IGBP) system. Using the four classification algorithms, we obtained the initial set of global land-cover maps. The SVM produced the highest overall classification accuracy (OCA) of 64.9% assessed with our test samples, with RF (59.8%), J4.8 (57.9%), and MLC (53.9%) ranked from the second to the fourth. We also estimated the OCAs using a subset of our test samples (8629) each of which represented a homogeneous area greater than 500 m × 500 m. Using this subset, we found the OCA for the SVM to be 71.5%. As a consistent source for estimating the coverage of global land-cover types in the world, estimation from the test samples shows that only 6.90% of the world is planted for agricultural production. The total area of cropland is 11.51% if unplanted croplands are included. The forests, grasslands, and shrublands cover 28.35%, 13.37%, and 11.49% of the world, respectively. The impervious surface covers only 0.66% of the world. Inland waterbodies, barren lands, and snow and ice cover 3.56%, 16.51%, and 12.81% of the world, respectively.

Multi-Resolution Land Characterization 2001 (MRLC 2001) is a second-generation Federal consortium designed to create an updated pool of nation-wide Landsat 5 and 7 imagery and derive a second-generation National Land Cover Database (NLCD 2001). The objectives of this multi-layer, multi-source database are twofold: first, to provide consistent land cover for all 50 States, and second, to provide a data framework which allows flexibility in developing and applying each independent data component to a wide variety of other applications. Components in the database include the following: (1) normalized imagery for three time periods per path/row, (2) ancillary data, including a 30 m Digital Elevation Model (DEM) derived into slope, aspect and slope position, (3) per pixel estimates of percent imperviousness and percent tree canopy, (4) 29 classes of land cover data derived from the imagery, ancillary data, and derivatives, (5) classification rules, confidence estimates, and metadata from the land cover classification. This database is now being developed using a Mapping Zone approach, with 66 Zones in the continental United States and 23 Zones in Alaska. Results from three initial mapping Zones show single-pixel land cover accuracies ranging from 73 to 77 percent, imperviousness accuracies ranging from 83 to 91 percent, tree canopy accuracies ranging from 78 to 93 percent, and an estimated 50 percent increase in mapping efficiency over previous methods. The database has now entered the production phase and is being created using extensive partnering in the Federal government with planned completion by 2006.

Rodríguez, J. P., T. D. Beard, Jr., E. M. Bennett, G. S. Cumming, S. Cork, J. Agard, A. P. Dobson, and G. D. Peterson. 2006. Trade-offs across space, time, and ecosystem services. Ecology and Society 11(1): 28. https://doi.org/10.5751/ES-01667-110128

The two largest earthquakes of the past 40 years ruptured a 1600-kilometer-long portion of the fault boundary between the Indo-Australian and southeastern Eurasian plates on 26 December 2004 [seismic moment magnitude (Mw) = 9.1 to 9.3] and 28 March 2005 (Mw = 8.6). The first event generated a tsunami that caused more than 283,000 deaths. Fault slip of up to 15 meters occurred near Banda Aceh, Sumatra, but to the north, along the Nicobar and Andaman Islands, rapid slip was much smaller. Tsunami and geodetic observations indicate that additional slow slip occurred in the north over a time scale of 50 minutes or longer.

Abstract Shifts in the timing of spring phenology are a central feature of global change research. Long‐term observations of plant phenology have been used to track vegetation responses to climate variability but are often limited to particular species and locations and may not represent synoptic patterns. Satellite remote sensing is instead used for continental to global monitoring. Although numerous methods exist to extract phenological timing, in particular start‐of‐spring (SOS), from time series of reflectance data, a comprehensive intercomparison and interpretation of SOS methods has not been conducted. Here, we assess 10 SOS methods for North America between 1982 and 2006. The techniques include consistent inputs from the 8 km Global Inventory Modeling and Mapping Studies Advanced Very High Resolution Radiometer NDVIg dataset, independent data for snow cover, soil thaw, lake ice dynamics, spring streamflow timing, over 16 000 individual measurements of ground‐based phenology, and two temperature‐driven models of spring phenology. Compared with an ensemble of the 10 SOS methods, we found that individual methods differed in average day‐of‐year estimates by ±60 days and in standard deviation by ±20 days. The ability of the satellite methods to retrieve SOS estimates was highest in northern latitudes and lowest in arid, tropical, and Mediterranean ecoregions. The ordinal rank of SOS methods varied geographically, as did the relationships between SOS estimates and the cryospheric/hydrologic metrics. Compared with ground observations, SOS estimates were more related to the first leaf and first flowers expanding phenological stages. We found no evidence for time trends in spring arrival from ground‐ or model‐based data; using an ensemble estimate from two methods that were more closely related to ground observations than other methods, SOS trends could be detected for only 12% of North America and were divided between trends towards both earlier and later spring.

The Thermal Emission Spectrometer (TES) investigation on Mars Global Surveyor (MGS) is aimed at determining (1) the composition of surface minerals, rocks, and ices; (2) the temperature and dynamics of the atmosphere; (3) the properties of the atmospheric aerosols and clouds; (4) the nature of the polar regions; and (5) the thermophysical properties of the surface materials. These objectives are met using an infrared (5.8‐ to 50‐μm) interferometric spectrometer, along with broadband thermal (5.1‐ to 150‐μm) and visible/near‐IR (0.3‐ to 2.9‐μm) radiometers. The MGS TES instrument weighs 14.47 kg, consumes 10.6 W when operating, and is 23.6×35.5×40.0 cm in size. The TES data are calibrated to a 1‐σ precision of 2.5 −6 ×10 −8 W cm −2 sr −1 /cm −1 , 1.6×10 −6 W cm −2 sr −1 , and ∼0.5 K in the spectrometer, visible/near‐IR bolometer, and IR bolometer, respectively. These instrument subsections are calibrated to an absolute accuracy of ∼4×10 −8 W cm −2 sr −1 /cm −1 (0.5 K at 280 K), 1–2%, and ∼1–2 K, respectively. Global mapping of surface mineralogy at a spatial resolution of 3 km has shown the following: (1) The mineralogic composition of dark regions varies from basaltic, primarily plagioclase feldspar and clinopyroxene, in the ancient, southern highlands to andesitic, dominated by plagioclase feldspar and volcanic glass, in the younger northern plains. (2) Aqueous mineralization has produced gray, crystalline hematite in limited regions under ambient or hydrothermal conditions; these deposits are interpreted to be in‐place sedimentary rock formations and indicate that liquid water was stable near the surface for a long period of time. (3) There is no evidence for large‐scale (tens of kilometers) occurrences of moderate‐grained (>50‐μm) carbonates exposed at the surface at a detection limit of ∼10%. (4) Unweathered volcanic minerals dominate the spectral properties of dark regions, and weathering products, such as clays, have not been observed anywhere above a detection limit of ∼10%; this lack of evidence for chemical weathering indicates a geologic history dominated by a cold, dry climate in which mechanical, rather than chemical, weathering was the significant form of erosion and sediment production. (5) There is no conclusive evidence for sulfate minerals at a detection limit of ∼15%. The polar region has been studied with the following major conclusions: (1) Condensed CO 2 has three distinct end‐members, from fine‐grained crystals to slab ice. (2) The growth and retreat of the polar caps observed by MGS is virtually the same as observed by Viking 12 Martian years ago. (3) Unique regions have been identified that appear to differ primarily in the grain size of CO 2 ; one south polar region appears to remain as black slab CO 2 ice throughout its sublimation. (4) Regional atmospheric dust is common in localized and regional dust storms around the margin and interior of the southern cap. Analysis of the thermophysical properties of the surface shows that (1) the spatial pattern of albedo has changed since Viking observations, (2) a unique cluster of surface materials with intermediate inertia and albedo occurs that is distinct from the previously identified low‐inertia/bright and high‐inertia/dark surfaces, and (3) localized patches of high‐inertia material have been found in topographic lows and may have been formed by a unique set of aeolian, fluvial, or erosional processes or may be exposed bedrock.

Microbes exist in a range of metabolic states (for example, dormant, active and growing) and analysis of ribosomal RNA (rRNA) is frequently employed to identify the 'active' fraction of microbes in environmental samples. While rRNA analyses are no longer commonly used to quantify a population's growth rate in mixed communities, due to rRNA concentration not scaling linearly with growth rate uniformly across taxa, rRNA analyses are still frequently used toward the more conservative goal of identifying populations that are currently active in a mixed community. Yet, evidence indicates that the general use of rRNA as a reliable indicator of metabolic state in microbial assemblages has serious limitations. This report highlights the complex and often contradictory relationships between rRNA, growth and activity. Potential mechanisms for confounding rRNA patterns are discussed, including differences in life histories, life strategies and non-growth activities. Ways in which rRNA data can be used for useful characterization of microbial assemblages are presented, along with questions to be addressed in future studies.

The normalized difference water index (NDWI) has been successfully used to delineate surface water features. However, two major problems have been often encountered: (a) NDWIs calculated from different band combinations [visible, nearinfrared, or shortwave-infrared (SWIR)] can generate different results, and (b) NDWI thresholds vary depending on the proportions of subpixel water/non-water components. We need to evaluate all the NDWIs for determining the best performing index and to establish appropriate thresholds for clearly identifying water features. We used the spectral data obtained from a spectral library to simulate the satellite sensors Landsat ETM, SPOT-5, ASTER, and MODIS, and calculated the simulated NDWI in different forms. We found that the NDWI calculated from (green ‐ SWIR)/(green SWIR), where SWIR is the shorter wavelength region (1.2 to 1.8 mm), has the most stable threshold. We recommend this NDWI be employed for mapping water, but adjustment of the threshold based on actual situations is necessary.

LANDFIRE is a 5-year, multipartner project producing consistent and comprehensive maps and data describing vegetation, wildland fuel, fire regimes and ecological departure from historical conditions across the United States. It is a shared project between the wildland fire management and research and development programs of the US Department of Agriculture Forest Service and US Department of the Interior. LANDFIRE meets agency and partner needs for comprehensive, integrated data to support landscape-level fire management planning and prioritization, community and firefighter protection, effective resource allocation, and collaboration between agencies and the public. The LANDFIRE data production framework is interdisciplinary, science-based and fully repeatable, and integrates many geospatial technologies including biophysical gradient analyses, remote sensing, vegetation modelling, ecological simulation, and landscape disturbance and successional modelling. LANDFIRE data products are created as 30-m raster grids and are available over the internet at www.landfire.gov, accessed 22 April 2009. The data products are produced at scales that may be useful for prioritizing and planning individual hazardous fuel reduction and ecosystem restoration projects; however, the applicability of data products varies by location and specific use, and products may need to be adjusted by local users.

Abstract Pearly mussels (Unionacea) are widespread, abundant, and important in freshwater ecosystems around the world. Catastrophic declines in pearly mussel populations in North America and other parts of the world have led to a flurry of research on mussel biology, ecology, and conservation. Recent research on mussel feeding, life history, spatial patterning, and declines has augmented, modified, or overturned long-held ideas about the ecology of these animals. Pearly mussel research has begun to benefit from and contribute to current ideas about suspension feeding, life-history theory, metapopulations, flow refuges, spatial patterning and its effects, and management of endangered species. At the same time, significant gaps in understanding and apparent paradoxes in pearly mussel ecology have been exposed. To conserve remaining mussel populations, scientists and managers must simultaneously and aggressively pursue both rigorous research and conservation actions.

The Kobe earthquake struck at the edge of the densely populated Osaka‐Kyoto corridor in southwest Japan. We investigate how the earthquake transferred stress to nearby faults, altering their proximity to failure and thus changing earthquake probabilities. We find that relative to the pre‐Kobe seismicity, Kobe aftershocks were concentrated in regions of calculated Coulomb stress increase and less common in regions of stress decrease. We quantify this relationship by forming the spatial correlation between the seismicity rate change and the Coulomb stress change. The correlation is significant for stress changes greater than 0.2–1.0 bars (0.02–0.1 MPa), and the nonlinear dependence of seismicity rate change on stress change is compatible with a state‐ and rate‐dependent formulation for earthquake occurrence. We extend this analysis to future mainshocks by resolving the stress changes on major faults within 100 km of Kobe and calculating the change in probability caused by these stress changes. Transient effects of the stress changes are incorporated by the state‐dependent constitutive relation, which amplifies the permanent stress changes during the aftershock period. Earthquake probability framed in this manner is highly time‐dependent, much more so than is assumed in current practice. Because the probabilities depend on several poorly known parameters of the major faults, we estimate uncertainties of the probabilities by Monte Carlo simulation. This enables us to include uncertainties on the elapsed time since the last earthquake, the repeat time and its variability, and the period of aftershock decay. We estimate that a calculated 3‐bar (0.3‐MPa) stress increase on the eastern section of the Arima‐Takatsuki Tectonic Line (ATTL) near Kyoto causes fivefold increase in the 30‐year probability of a subsequent large earthquake near Kyoto; a 2‐bar (0.2‐MPa) stress decrease on the western section of the ATTL results in a reduction in probability by a factor of 140 to 2000. The probability of a M w = 6.9 earthquake within 50 km of Osaka during 1997–2007 is estimated to have risen from 5–6% before the Kobe earthquake to 7–11% afterward; during 1997–2027, it is estimated to have risen from 14–16% before Kobe to 16–22%.

Understanding groundwater recharge is essential for successful management of water resources and modeling fluid and contaminant transport within the subsurface. This book provides a critical evaluation of the theory and assumptions that underlie methods for estimating rates of groundwater recharge. Detailed explanations of the methods are provided - allowing readers to apply many of the techniques themselves without needing to consult additional references. Numerous practical examples highlight benefits and limitations of each method. Approximately 900 references allow advanced practitioners to pursue additional information on any method. For the first time, theoretical and practical considerations for selecting and applying methods for estimating groundwater recharge are covered in a single volume with uniform presentation. Hydrogeologists, water-resource specialists, civil and agricultural engineers, earth and environmental scientists and agronomists will benefit from this informative and practical book. It can serve as the primary text for a graduate-level course on groundwater recharge or as an adjunct text for courses on groundwater hydrology or hydrogeology. For the benefit of students and instructors, problem sets of varying difficulty are available at http://wwwbrr.cr.usgs.gov/projects/GW_Unsat/Recharge_Book/

HypoDD is a Fortran computer program package for relocating earthquakes with the double-difference algorithm of Waldhauser and Ellsworth (2000). This document provides a brief introduction into how to run and use the programs ph2dt and hypoDD to compute double-difference (DD) hypocenter locations. It gives a short overview of the DD technique, discusses the data preprocessing using ph2dt, and leads through the earthquake relocation process using hypoDD. The appendices include the reference manuals for the two programs and a short description of auxiliary programs and example data. Some minor subroutines are presently in the c language, and future releases will be in c.

Coral recruitment is a key process in the maintenance and recovery of coral reef ecosystems. While intense competition between coral and algae is often assumed on reefs that have undergone phase shifts from coral to algal dominance, data examining the competitive interactions involved, particularly during the larval and immediate post-settlement stage, are scarce. Using a series of field and outdoor seawater table experiments, we tested the hypothesis that common species of macroalgae and cyanobacteria inhibit coral recruitment. We examined the effects of Lyngbya spp., Dictyota spp., Lobophora variegata (J. V. Lamouroux) Womersley, and Chondrophycus poiteaui (J. V. Lamouroux) Nam (formerly Laurencia poiteaui) on the recruitment success of Porites astreoides larvae. All species but C. poiteaui caused either recruitment inhibition or avoidance behavior in P. astreoides larvae, while L. confervoides and D. menstrualis significantly increased mortality rates of P. astreoides recruits. We also tested the effect of some of these macrophytes on larvae of the gorgonian octocoral Briareum asbestinum. Exposure to Lyngbya majuscula reduced survival and recruitment in the octocoral larvae. Our results provide evidence that algae and cyanobacteria use tactics beyond space occupation to inhibit coral recruitment. On reefs experiencing phase shifts or temporary algal blooms, the restocking of adult coral populations may be slowed due to recruitment inhibition, thereby perpetuating reduced coral cover and limiting coral community recovery.

Weather and climate extremes have been varying and changing on many different time scales. In recent decades, heat waves have generally become more frequent across the United States, while cold waves have been decreasing. While this is in keeping with expectations in a warming climate, it turns out that decadal variations in the number of U.S. heat and cold waves do not correlate well with the observed U.S. warming during the last century. Annual peak flow data reveal that river flooding trends on the century scale do not show uniform changes across the country. While flood magnitudes in the Southwest have been decreasing, flood magnitudes in the Northeast and north-central United States have been increasing. Confounding the analysis of trends in river flooding is multiyear and even multidecadal variability likely caused by both large-scale atmospheric circulation changes and basin-scale “memory” in the form of soil moisture. Droughts also have long-term trends as well as multiyear and decadal variability. Instrumental data indicate that the Dust Bowl of the 1930s and the drought in the 1950s were the most significant twentieth-century droughts in the United States, while tree ring data indicate that the megadroughts over the twelfth century exceeded anything in the twentieth century in both spatial extent and duration. The state of knowledge of the factors that cause heat waves, cold waves, floods, and drought to change is fairly good with heat waves being the best understood.

ADVERTISEMENT RETURN TO ISSUEPREVArticleNEXTWater solubility enhancements of DDT and trichlorobenzene by some surfactants below and above the critical micelle concentrationDaniel E. Kile and Cary T. ChiouCite this: Environ. Sci. Technol. 1989, 23, 7, 832–838Publication Date (Print):July 1, 1989Publication History Published online1 May 2002Published inissue 1 July 1989https://doi.org/10.1021/es00065a012RIGHTS & PERMISSIONSArticle Views2917Altmetric-Citations363LEARN ABOUT THESE METRICSArticle Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated. Share Add toView InAdd Full Text with ReferenceAdd Description ExportRISCitationCitation and abstractCitation and referencesMore Options Share onFacebookTwitterWechatLinked InReddit PDF (895 KB) Get e-Alerts Get e-Alerts

Oil-shale deposits are in many parts of the world. They range in age from Cambrian to Tertiary and were formed in a variety of marine, continental, and lacustrine depositional environments. The largest known deposit is in the Green River Formation in the western United States; it contains an estimated 213 billion tons of in-situ shale oil (about 1.5 trillion U.S. barrels). Total resources of a selected group of oil shale deposits in 33 countries are estimated at 409 billion tons of in-situ shale oil, which is equivalent to 2.8 trillion U.S. barrels of shale oil. These amounts are very conservative because (1) several deposits mentioned herein have not been explored sufficiently to make accurate estimates, and (2) some deposits were not included in this survey.

It has been proposed that Jupiter's satellite Europa currently possesses a global subsurface ocean of liquid water. Galileo gravity data verify that the satellite is differentiated into an outer H 2 O layer about 100 km thick but cannot determine the current physical state of this layer (liquid or solid). Here we summarize the geological evidence regarding an extant subsurface ocean, concentrating on Galileo imaging data. We describe and assess nine pertinent lines of geological evidence: impact morphologies, lenticulae, cryovolcanic features, pull‐apart bands, chaos, ridges, surface frosts, topography, and global tectonics. An internal ocean would be a simple and comprehensive explanation for a broad range of observations; however, we cannot rule out the possibility that all of the surface morphologies could be due to processes in warm, soft ice with only localized or partial melting. Two different models of impact flux imply very different surface ages for Europa; the model favored here indicates an average age of ∼50 Myr. Searches for evidence of current geological activity on Europa, such as plumes or surface changes, have yielded negative results to date. The current existence of a global subsurface ocean, while attractive in explaining the observations, remains inconclusive. Future geophysical measurements are essential to determine conclusively whether or not there is a liquid water ocean within Europa today.