Our objective was to determine if the occurrence and abundance of Pacific loons (Gavia pacifica Lawrence) and horned grebes (Podiceps auritus L.) on 123 wetlands of Yukon Flats National Wildlife Refuge in east central Alaska were related to the limnological characteristics of those wetlands. Aquatic bird-wetland use surveys were conducted in conjunction with limnological sampling from May through September 1985–87 and May through August 1989. Results from logistic regression analysis demonstrated a significant association between the probability of wetland use by Pacific loons and shoreline length, water color, calcium and total phosphorus. Wetland use by horned grebes was related to shoreline length, pH, and chlorophyll. Aquatic bird abundance was then used as a Poisson response variable and modeled as a function of wetland limnological characteristics. Our results indicate that Pacific loon abundance was adequately modeled by linear and quadratic functions of shoreline length, color, pH, calcium and total phosphorus. Horned grebe abundance could not be modeled with this approach. The statistical techniques known collectively as generalized linear models provided a framework for the development of models for aquatic bird use of wetlands. Our results, however, indicate that while this approach shows promise, a better understanding of how to model aquatic bird abundance is needed. We then identify problems in model development and suggest avenues for future research.