A finite difference version of the Levenberg-Marquardt method for nonlinear least squares problems has been extended to include inverse problems in distributed estuarial hydraulic systems. The objective in solving the inverse problems was to establish a numerical simulation procedure for estimating bestdistributed Manning's roughness coefficients from sets of observed tide heights. As an illustration, spatially varying Manning's roughness coefficients for the Upper Delaware River Estuary system were determined for several representative sets of tide height data for the period October 1973 to June 1974. The roughness coefficients were modeled as polynomial functions of distance. Manning's n was thus found generally to vary inversely with distance from the head of tide at Trenton to Wilmington. The spatially distributed tidal-averaged Reynolds number Re was used to correlate Manning's n and Darcy-Weisbach's f. The resultant n-Re relationships displayed threi: distinct hydrodynamic flow regimes characterized as having turbulence. Both n and fwere found to be independent of Re for Re > 1.52 X 10 • but inversely related to Re for Re < 1.2 X 10L Among the numerical techniques used to simulate tidal hydraulic transients it was found that a 'hopscotch' finite difference method yielded the best compromise between computational economy and overall accuracy. PROLOGUE Inverse problems in water resources research have been gaining the fancy of many investigators in recent years. The literature is replete with claims and counterclaims of generalized methodologies and unique results. In the sense of Matalas and Maddock [1976] there is certainly an 'aura of mystique' that surrounds these debates. The bottom lines on both sides in these controversies, however, are often tinted with certain degrees of subjective value judgments and risks. Perhaps after all we are all saved a little in the end by Einstein's universal axiom which states that '... the most incomprehensible thing is that anything is comprehensible.' To be sure, inverse problem modelers do not claim complete and final understanding of what they have concocted. What they concoct, if the concocting is done with a sense of practicality, usually results in tangible measures of progress. The mere process of formulating and then solving a real world inverse problem usually constitutes a learning process in itself. The process begets new ideas and underscores deficiencies. In the present• work an attempt is made to tackle a real world inverse problem involving the spatially distributed Manning's roughness coefficients for open estuarial river systems. The work entails consideration of real world input data and hydrodynamic state equations, including a model for spatially distributed frictional energy dissipation, development of an appropriate inverse problem algorithm, and implementation of a fast, stable, and accurate numerical integration routine. The combination of these elements should constitute a useful numerical procedure for workers in the field. The reader will also ...