In response to the need and an intergovernmental commission for a high resolution and data-derived global ecosystem map, land surface elements of global ecological pattern were characterized in an ecophysiographic stratification of the planet. The stratification produced 3,923 terrestrial ecological land units (ELUs) at a base resolution of 250 meters. The ELUs were derived from data on land surface features in a three step approach. The first step involved acquiring or developing four global raster datalayers representing the primary components of ecosystem structure: bioclimate, landform, lithology, and land cover. These datasets generally represent the most accurate, current, globally comprehensive, and finest spatial and thematic resolution data available for each of the four inputs. The second step involved a spatial combination of the four inputs into a single, new integrated raster dataset where every cell represents a combination of values from the bioclimate, landforms, lithology, and land cover datalayers. This foundational global raster datalayer, called ecological facets (EFs), contains 47,650 unique combinations of the four inputs. The third step involved an aggregation of the EFs into the 3,923 ELUs.