scaling gridded river networks for macroscale hydrology development analysis and control of error a simple and robust river network scaling algorithm nsa is presented to rescale fine resolution networks to any coarser resolution the algorithm was tested over the danube river basin and the european continent coarse resolution networks at 2 5 5 10 and 30 min resolutions were derived from higher resolution gridded networks using nsa and geomorphometric attributes such as river order shape index and width function these parameters were calculated and compared at each resolution simple scaling relationships were found to predict decreasing river lengths with coarser resolution data this relationship can be used to correct river length as a function of grid resolution the length corrected width functions of the major river basins in europe were compared at different resolutions to assess river network performance the discretization error in representing basin area and river lengths at coarser resolutions were analyzed and simple relationships were found to calculate the minimum number of grid cells needed to maintain the catchment area and length within a desired level of accuracy this relationship among geomorphological characteristics such as shape index and width function derived from gridded networks at different resolutions suggests that a minimum of 200 300 grid cells is necessary to maintain the geomorphological characteristics of the river networks with sufficient accuracy

holocene lake level trends in the rocky mountains u s a abstract the availability of water shapes life in the western united states and much of the water in the region originates in the rocky mountains few studies however have explicitly examined the history of water levels in the rocky mountains during the holocene here we examine the past levels of three lakes near the continental divide in montana and colorado to reconstruct holocene moisture trends using transects of sediment cores and sub surface geophysical profiles from each lake we find that mid holocene shorelines in the small lakes 4 110 ha were as much as 10 m below the modern lake surfaces our results are consistent with existing evidence from other lakes and show that a wide range of settings in the region were much drier than today before 3000 2000 years ago we also discuss evidence for millennial scale moisture variation including an abruptly initiated and terminated wet period in colorado from 4400 to 3700 cal yr bp and find only limited evidence for low lake stands during the past millennium the extent of low water levels during the mid holocene which were most severe and widespread ca 7000 4500 cal yr bp is consistent with the extent of insolation induced aridity in previously published regional climate model simulations like the simulations the lake data provide no evidence for enhanced zonal flow during the mid holocene which has been invoked to explain enhanced mid continent aridity at the time the data including widespread evidence for large changes on orbital time scales and for more limited changes during the last millennium confirm the ability of large boundary condition changes to push western water supplies beyond the range of recent natural variability

improvements in large scale drainage networks derived from digital elevation models 1 this paper presents an improved algorithm for deriving drainage networks for coarse gridded distributed hydrological models based on relatively fine resolution digital elevation models the proposed algorithm aims to reproduce actual drainage networks more closely when applied in regions with meandering rivers running in parallel and for cell sizes of the order of 10 kilometers to achieve this goal the cotat algorithm of reed 2003 was improved introducing another parameter related to the minimum upstream flow path into the cell besides the area threshold the proposed algorithm was tested by applying it to the catchments of the rivers tapajos 530 000 km 2 and grande 145 000 km 2 in south america results show that drainage networks are generally improved by the proposed algorithm needing only minor manual corrections

rise and fall of late pleistocene pluvial lakes in response to reduced evaporation and precipitation evidence from lake surprise california widespread late pleistocene lake systems of the basin and range province indicate substantially greater moisture availability during glacial periods relative to modern times but the climatic factors that drive changes in lake levels are poorly constrained to better constrain these climatic forcing factors we present a new lacustrine paleoclimate record and precipitation estimates for lake surprise a closed basin lake in northeastern california we combine a detailed analysis of lake hydrography and constitutive relationships describing the water balance to determine the influence of precipitation evaporation temperature and seasonal insolation on past lake levels at its maximum extent during the last deglaciation lake surprise covered 1366 km 2 36 of the terminally draining surprise valley watershed using paired radiocarbon and 230 th u analyses we dated shoreline tufa deposits from wave cut lake terraces in surprise valley california to determine the hydrography of the most recent lake cycle this new lake hydrograph places the highest lake level 176 m above the present day playa at 15 19 0 18 calibrated ka 14 c age this significantly postdates the last glacial maximum lgm when lake surprise stood at only moderate levels 65 99 m above modern playa similar to nearby lake lahontan to evaluate the climatic factors associated with lake level changes we use an oxygen isotope mass balance model combined with an analysis of predictions from the paleoclimate model intercomparison project 3 pmip3 climate model ensemble our isotope mass balance model predicts minimal precipitation increases of only 2 18 during the lgm relative to modern compared to an 75 increase in precipitation during the 15 19 ka highstand lgm pmip3 climate model simulations corroborate these findings simulating an average precipitation increase of only 6 5 relative to modern accompanied by a 28 decrease in total evaporation driven by a 7 c decrease in mean annual temperature lgm pmip3 climate model simulations also suggest a seasonal decoupling of runoff and precipitation with peak runoff shifting to the late spring early summer from the late winter early spring our coupled analyses suggest that moderate lake levels during the lgm were a result of reduced evaporation driven by reduced summer insolation and temperatures not by increased precipitation reduced evaporation primed basin and range lake systems particularly smaller isolated basins such as surprise valley to respond rapidly to increased precipitation during late heinrich stadial 1 hs1 post lgm highstands were potentially driven by increased rainfall during hs1 brought by latitudinally extensive and strengthened midlatitude westerly storm tracks the effects of which are recorded in the region9s lacustrine and glacial records these results suggest that seasonal insolation and reduced temperatures have been underinvestigated as long term drivers of moisture availability in the western united states

rates timing and cyclicity of holocene eolian activity in north central united states evidence from varved lake sediments most of the sediment components that accumulated in elk lake northwestern minnesota during the holocene are autochthonous or biogenic delivered to the sediment water interface on a seasonal schedule preserved in distinct annual laminae varves the main allochthonous component is detrital clastic material as measured by bulk sediment concentrations of aluminum sodium potassium titanium and quartz that enters the lake mostly as eolian dust the eolian clastic influx to elk lake was considerably greater during the mid holocene 8 4 ka than it has been for the past 4000 yr when periods of increased eolian activity correspond to the time of the little ice age and the dust bowl geochemical records of eolian activity exhibit distinct cyclicities with dominant periodicities of 400 and 84 yr

automated recognition of valley lines and drainage networks from grid digital elevation models a review and a new method abstract the literature on valley and drainage line recognition from grid digital elevation models dems is reviewed two problems preventing the successful delineation of fully connected drainage networks of single line width are identified a the positioning of the ends of drainage networks and b the assignment of drainage directions to a dem in particular the assignment of drainage directions across flat areas and at closed depressions in the dem pits a new method for drainage network delineation is suggested which incorporates an alternative method for drainage direction assignment over a dem this method is found to produce realistic directions and performs well on two complex dems

middle holocene dry period in the northern midwestern united states lake levels and pollen stratigraphy abstract four cores from lake mendota in the yahara river valley of south central wisconsin provide pollen charcoal and sediment stratigraphic evidence of a middle holocene dry period in the northern midwest the lake level was lower between about 6500 and 3500 yr b p indicating that runoff from the upper yahara river basin may have ceased during the drier seasons of the year derived estimates of precipitation between 6500 and 3500 yr b p indicate a decrease of about 10 down to about 700 mm from the present precipitation value of 800 mm a warming between 6500 and 3500 yr b p is also interpreted based on a change of vegetation from rich mesophytic forest before 6500 yr b p to quercus savanna with increased charcoal abundances after that time after 3500 yr b p a closed quercus forest and decreased charcoal suggest a cooler and wetter climate the qualitative paleoecological interpretations quantitative precipitation reconstructions based on pollen from the region including the lake mendota data and estimates of hydrologic budget based on lake level changes all show indications of a dry and warm middle holocene evidence from a large regional array of sites also supports these interpretations

a methodology for preserving channel flow networks and connectivity patterns in large scale distributed hydrological models physiographic data are often used to parameterize hydrological models and in the past physiographic parameters have often been derived manually however this can be a lengthy and unreliable process particularly for application to a gridded hydrological or atmospheric model applied to large or continental scale basins an important attribute of gridded models is drainage direction current methods that determine drainage directions for large or continental scale basins by general circulation models gcms route flow using lowest neighbour algorithms these methods however do not reflect the hydrology of the basin subunit this paper proposes a method of parameterizing hydrological models with physiographic data using the arclnfo macro language to create an interface between the topographic parameterization topaz software and the watflood hydrological model the interface uses output raster data created by topaz i e drainage identification to supply physiographic parameters required by watflood the interface watpaz is an expert system based on a manual method of deriving parameters for the watflood distributed model the watpaz interface uses grouped response units to subdivide the watershed this allows large drainage basins to be subdivided at a scale that allows computational efficiency while preserving the hydrological variability of the watershed to test whether the watpaz method improves the current gcm methodology for determining drainage directions watpaz is applied on a local basin wolf creek a regional scale basin athabasca and a continental scale basin mackenzie an examination of flow directions derived from this new method with current gcm methods is carried out the results indicate that a substantial improvement is made to flow routing within the basin using the channel network to determine drainage directions for each segment

hydraulic causes for basin hydrograph skewness it has been suggested that hydrograph skewness depends on the relative dominance between hillslope and channel network transport processes where the former ones make the hydrograph positively skewed while the latter tends to make it negatively skewed more recently however the role of river hydraulics in shaping the hydrograph has been highlighted we present a set of numerical modelling experiments using a hydrodynamic model of river networks in which we investigate how channel hydraulics influence the shape of the hydrograph particularly its skewness we further investigate the influence of baseflow rainfall intensity cross section geometry and basin scale on unit hydrograph response we show that river hydraulics has a decisive role because positively skewed hydrographs may occur even when water inputs to the river network are negatively skewed and in catchments whose width function is also negatively skewed we show additional results related to the effect of the degree of non linearity in the relationship between celerity and discharge the effects of baseflow and rainfall intensity these further confirm that hydraulic factors may be decisive in determining hydrograph shape

digital terrain modelling a review of hydrological geomorphological and biological applications the topography of a catchment has a major impact on the hydrological geomorphological and biological processes active in the landscape the spatial distribution of topographic attributes can often be used as an indirect measure of the spatial variability of these processes and allows them to be mapped using relatively simple techniques many geographic information systems are being developed that store topographic information as the primary data for analysing water resource and biological problems furthermore topography can be used to develop more physically realistic structures for hydrologic and water quality models that directly account for the impact of topography on the hydrology digital elevation models are the primary data used in the analysis of catchment topography we describe elevation data sources digital elevation model structures and the analysis of digital elevation data for hydrological geomorphological and biological applications some hydrologic models that make use of digital representations of topography are also considered

application of the promethee technique to determine depression outlet location and flow direction in dem with the fast growing progress of computer technologies spatial information on watersheds such as flow direction watershed boundaries and the drainage network can be automatically calculated or extracted from a digital elevation model dem the stubborn problem that depressions exist in dems has been frequently encountered while extracting the spatial information of terrain several filling up methods have been proposed for solving depressions however their suitability for large scale flat areas is inadequate this study proposes a depression watershed method coupled with the preference ranking organization method for enrichment evaluations promethees theory to determine the optimal outlet and calculate the flow direction in depressions three processing procedures are used to derive the depressionless flow direction 1 calculating the incipient flow direction 2 establishing the depression watershed by tracing the upstream drainage area and determining the depression outlet using promethee theory 3 calculating the depressionless flow direction the developed method was used to delineate the shihmen reservoir watershed located in northern taiwan the results show that the depression watershed method can effectively solve the shortcomings such as depression outlet differentiating and looped flow direction between depressions the suitability of the proposed approach was verified

holocene climatic change in the northern midwest pollen derived estimates mapping of holocene pollen data in the midwestern united states has revealed several broadscale vegetational changes that can be interpreted in climatic terms these changes include 1 the early holocene northward movement of the spruce dominated forest and its later southward movement after 3000 yr b p and 2 the eastward movement of the prairie forest border into southwestern wisconsin by 8000 yr b p and its subsequent westward retreat after 6000 yr b p when certain basic assumptions are met multiple regression models can be derived from modern pollen and climate data and used to transform the pollen record of these vegetational changes into quantitative estimates of temperature or precipitation to maximize the reliability of the regression equations we followed a sequence of procedures that minimize violations of the assumptions that underlie regression analysis reconstructions of precipitation during the holocene indicated that from 9000 to 6000 yr b p precipitation decreased by 10 to 25 over much of the midwest while mean july temperature increased by 0 5 to 2 0 c at 6000 yr b p precipitation was less than 80 of its modern values over parts of wisconsin and minnesota after 6000 yr b p precipitation generally increased while mean july temperature decreased in the north and increased in the south the time of the maximum temperature varies within the midwest and is earlier in the north and later in the south

determination of the drainage structure of a watershed using a digital elevation model and a digital river and lake network abstract distributed hydrological models require a detailed definition of a watershed s internal drainage structure the conventional approach to obtain this drainage structure is to use an eight flow direction matrix d8 which is derived from a raster digital elevation model dem however this approach leads to a rather coarse drainage structure when monitoring or gauging stations need to be accurately located within a watershed this is largely due to limitations of the d8 approach and the lack of information over flat areas and pits the d8 approach alone is also unable to differentiate lakes from plain areas to avoid these problems a new approach using a digital river and lake network drln as input in addition to the dem has been developed this new approach allows for an accurate fit between the drln and the modelled drainage structure which is represented by a flow direction matrix and a modelled watercourse network more importantly the identification of lakes within the modelled network is now possible the proposed approach which is largely rooted in the d8 approach uses the drln to correct modelled flow directions and network calculations for dem cells overlapped by the drln flow directions are determined using drln connections only the flow directions of the other dem cells are evaluated with the d8 approach which uses a dem that has been modified as a function of distance to the drln the proposed approach has been tested on the chaudiere river watershed in southern quebec canada the modelled watershed drainage structure showed a high level of coherence with the drln a comparison between the results obtained with the d8 approach and those obtained by the proposed approach clearly demonstrated an improvement over the conventionally modelled drainage structure the proposed approach will benefit hydrological models which require data such as a flow direction matrix a river and lake network and sub watersheds for drainage structure information

validation of a new global 30 min drainage direction map abstract digital drainage direction maps are a prerequisite for analyzing the flow of water on the land surface of the earth for continental or global studies the most appropriate and most frequently used resolution for such data sets appears to be 30 longitude by latitude in this paper we present the new global drainage direction map ddm30 a 30 raster map of surface drainage directions which organizes the earth s land area into drainage basins and provides the river network topology ddm30 was generated by first upscaling two drainage direction maps ddms at higher resolutions the resulting map was then extensively corrected in an iterative manner by comparison against vectorized high resolution river maps and other geographic information finally it was co referenced to the locations of 935 gauging stations provided by the global runoff data centre grdc which again involved manual corrections ddm30 was validated against drainage basin areas from the literature against the given upstream areas of the grdc stations and most importantly against information from hydro1k a data set based on a hydrologically corrected 1 km digital elevation model which is thought to afford the best information on surface drainage currently available at the global scale in the course of the validation the quality of ddm30 was compared to three other 30 ddms the validation results show that ddm30 provides a more accurate representation of drainage directions and river network topology than the other 30 ddms

a morphology independent methodology for quantifying planview river change and characteristics from remotely sensed imagery abstract remotely sensed imagery of rivers has long served as a means for characterizing channel properties and detection of planview change in the last decade the dramatic increase in the availability of satellite imagery and processing tools has created the potential to greatly expand the spatial and temporal scale of our understanding of river morphology and dynamics to date the majority of gis and automated analyses of planview changes in rivers from remotely sensed data has been developed for single threaded meandering river systems these methods have limited applicability to many of the earth s rivers with complex multi channel planforms here we present the methodologies of a set of analysis algorithms collectively called spatially continuous riverbank erosion and accretion measurements scream scream analyzes planview river metrics regardless of river morphology these algorithms quantify both the erosion and accretion rates of riverbanks from binary masks of channels generated from imagery acquired at two time periods additionally the program quantifies the area of change between river channels and the surrounding floodplain and area of islands lost or formed between these two time periods to examine variations in erosion rates in relation to local channel attributes and make rate comparisons between river systems of varying sizes the program determines channel widths and bank curvature at every bank pixel scream was developed and tested on rivers with diverse and complex planform morphologies in imagery acquired from a range of observational platforms with varying spatial resolutions validation and verification of scream generated metrics against manual measurements show no significant measurement errors in determination of channel width erosion and bank aspects scream has the potential to provide data for both the quantitative examination of the controls on erosion rates and for the comparison of these rates across river systems ranging broadly in size and planform morphology

where is the north china south china block boundary in eastern china the qinling dabieshan belt results from the collision of the north china and south china blocks the eastern extension of the belt the sulu area consists of several stacked units from top to bottom these are 1 a weakly metamorphosed slate sandstone unit of neoproterozoic paleozoic age 2 a gneiss quartzite unit metamorphosed under high pressure p conditions 3 a marble amphibolite unit 4 an ultrahigh p unit that was very deeply subducted and then exhumed and 5 a migmatitic dome resulting from the melting of ultrahigh p rocks all these units underwent the same deformation which was characterized by a top to the northwest extensional ductile shearing during their exhumation granulite facies restites in the migmatite are not significantly different from granulitized eclogites derived from the ultrahigh p rocks the migmatite exhibits the same structural and petrologic evolution as that determined for the other units the resemblance of petrologic and structural features and the lack of ocean basin rock shows that the boundary between the north china block and south china block must be placed north of the sulu area

the treatment of flat areas and depressions in automated drainage analysis of raster digital elevation models methods developed to process raster digital elevation models dem automatically in order to delineate and measure the properties of drainage networks and drainage basins are being recognized as potentially valuable tools for the topographic parameterization of hydrological models all of these methods ultimately rely on some form of overland flow simulation to define drainage courses and catchment areas and therefore have difficulty dealing with closed depressions and flat areas on digital land surface models some fundamental assumptions about the nature of these problem topographic features in dem are implicit in the various techniques developed to deal with them in automated drainage analysis the principal assumptions are 1 that closed depressions and flat areas are spurious features that arise from data errors and limitations of dem resolution 2 that flow directions across flat areas are determined solely by adjacent cells of lower elevation and 3 that closed depressions are caused exclusively by the underestimation of dem elevations it is argued that while the first of these assumptions is reasonable given the quality of dems generally available for hydrological analysis the others are not rather it seems more likely that depressions are caused by both under and overestimation errors and that flow directions across flat areas are determined by the distribution of both higher and lower elevations surrounding flat areas two new algorithms are introduced that are based on more reasonable assumptions about the nature of flat areas and depressions and produce more realistic results in application these algorithms allow breaching of depression outlets and consider the distribution of both higher and lower elevations in assigning flow directions on flat areas the results of applying these algorithms to some real and hypothetical landscapes are presented 1998 john wiley sons ltd

temporal correspondence between pluvial lake highstands in the southwestern us and heinrich event 1 pluvial lakes were abundant in the southwestern united states during pleistocene glaciations particu larly in the great basin many of these lakes occupied closed basins therefore fluctuations of their water surface elevations are valuable sources of paleoclimate information histories of the largest lakes are well constrained whereas dozens of smaller lakes that were present in this region have received relatively little scientific attention given their dimensions these smaller lakes were climatically sensitive and can offer important information about quaternary climate variability here we present new ages for the highstands of three previously undated small lakes based on radiocarbon dating of gastropod shells recovered from beach ridges these results are combined with other published and unpublished 14 c ages to yield an extensive compilation of highstand shoreline ages for lakes of all sizes throughout the southwestern us the results indicate that although some lakes reached highstands during the last glacial maximum the strongest temporal correspondence is between highstands and heinrich event h1 these results are consistent with speleothem based reconstructions of effective moisture in the southwestern us which show increased precipitation during stadials of the last glacial cycle copyright 2012 john wiley sons ltd

how did the hydrologic cycle respond to the two phase mystery interval abstract lake estancia s transition from a big dry episode during the first half of the mystery interval to a big wet episode during the second half has equivalents in records from across the planet at the time of this transition chinese monsoons experienced pronounced weakening closed basin lakes in both the great basin of the western united states and in the southern altiplano of south america underwent a major expansion mountain glaciers in southern hemisphere middle latitudes had retreated and the rates of increase of co2 and of 18o in antarctic ice underwent a decrease finally the precipitous drop in dust rain over antarctica and the southern ocean terminated as did a similar drop in the 13c to 12c ratio in atmospheric co2 these changes are consistent with a southward shift of the thermal equator the cause of such a shift is thought to be an expansion of sea ice caused by a shutdown in deep water production in the northern atlantic this creates a dilemma because a similar southward shift is an expected consequence of the heinrich event 1 which initiated the mystery interval

paleoclimate simulations for north america over the past 21 000 years features of the simulated climate and comparisons with paleoenvironmental data maps of upper level and surface winds and of surface temperature and precipitation illustrate the results of a sequence of global paleoclimatic simulations spanning the past 21 000 years for north america we review a the large scale features of circulation temperature and precipitation that appear in the simulations from the ncar community climate model version 1 ccm 1 b the implications of the simulated climate for the past continental scale distributions of three plant taxa picea spp pseudotsuga menziesii and artemisia tridentata which are broadly representative of the vegetation across the continent and c the potential explanations in terms of atmospheric circulation or surface energy and water balance processes for mismatches between the simulations and observations most of the broad scale features of previous paleoclimatic simulations with the ncar ccm 0 for north america are present in the current simulations many of the elements of a conceptual model based on previous climate simulations that describes the controls of paleoclimatic variations across north america during the past 21 000 years are found in simulations reviewed here these include 1 displacement of the jet stream by the laurentide ice sheet to the south of its present position in both winter and summer 2 generation of a glacial anticyclone over the ice sheet at the lgm and the consequent induction of large scale sinking motions induced over eastern north america 3 changes in the strength of surface atmospheric circulation features through time including weakening of the aleutian low in winter and strengthening of the eastern pacific and bermuda high pressure systems in summer as the ice sheet decreased in size 4 development of a heat low at the surface and a strengthened ridge in the upper atmosphere over the continent at the time of the maximum summer insolation anomaly 5 increases in summer temperature earlier in regions remote from the ice sheet these increases appear earlier in the present ccm 1 simulations than in the previous ccm 0 ones however and 6 continuation of negative winter temperature anomalies into the middle holocene in general simulated surface conditions that are discordant with paleoenvironmental observations can be attributed to the simulation of particular atmospheric circulation patterns e g those that suppress precipitation or advect warm air into a region with these mismatches amplified in beringia and the southeastern united states by surface energy and water balance processes

deriving flow directions for coarse resolution 1 4 km gridded hydrologic modeling 1 the national weather service hydrology laboratory nws hl is currently testing a grid based distributed hydrologic model at a resolution 4 km commensurate with operational radar based precipitation products to implement distributed routing algorithms in this framework a flow direction must be assigned to each model cell a new algorithm referred to as cell outlet tracing with an area threshold cotat has been developed to automatically accurately and efficiently assign flow directions to any coarse resolution grid cells using information from any higher resolution digital elevation model although similar to previously published algorithms this approach offers some advantages use of an area threshold allows more control over the tendency for producing diagonal flow directions analyses of results at different output resolutions ranging from 300 m to 4000 m indicate that it is possible to choose an area threshold that will produce minimal differences in average network flow lengths across this range of scales flow direction grids at a 4 km resolution have been produced for the conterminous united states

improving hydrological information acquisition from dem processing in floodplains extraction of hydrological information from digital elevation models dems is a required step when conducting any spatially distributed hydrological modelling in particular automated methods are proposed to extract the drainage structure from the dem however a realistic river network is not always derived from conventional dem processing methods indeed inaccuracy occurs in flat areas corresponding to floodplains in these areas additional sources of information are required to extract the correct drainage direction from the dem in this study it is demonstrated that traditional approaches of dem preprocessing such as the commonly known stream burning fail to provide correct maps of drainage directions and catchment areas when the extension of flat areas is large a new method is proposed to take advantage of available imagery data this method is based on a double dem burning process dem is first burned in the main rivers using the channel network and then in the floodplain using the spatial distribution of floodplains provided by classified satellite images in this sense the method has been referred to as the floodplain burning approach or simply fb approach spatial distribution of floodplains is derived from a multitemporal sar image classification a system of equations is used to vary the elevation offset required to be burnt in each cell representing the floodplain according to the minimal distance from the channel network which is calculated by a distance transformation the fb approach was applied to a sub basin located within the larger amazon river basin the region is characterized by large floodplain extensions basin delineation maps derived from the new method were compared with those obtained from the traditional stream burning method and highlighted more realistic results copyright 2008 john wiley sons ltd

river reach length and slope estimates for large scale hydrological models based on a relatively high resolution digital elevation model application of regular grid based distributed hydrological models requires information related to river drainage networks such as flow directions flow accumulated areas basin delineation and length and slope of river reaches for every grid cell most of those data may be extracted automatically from digital elevation models dem however when dealing with large basins available dems are usually in higher spatial resolution than the model grid upscaling procedures have therefore been developed to extract low resolution flow directions and flow accumulated areas from relatively high resolution dems this paper presents a new methodology which extends the upscaling method to automatically extract length and slope of river reaches for large scale grid based hydrological models the proposed method is new and is believed to be the first attempt to produce such information automatically from dems the methodology was applied to parts of the uruguay river basin in south america using the globally wide available dem produced by the shuttle radar topography mission srtm 90m quality of results was assessed by comparing calculated river lengths with distances measured over vectorized river networks which were assumed to be correct it was shown that the proposed methodology adequately assigns river lengths to every model cell while ensuring that the whole length of the streams is considered the effect of dem resolution on calculated river length errors was analyzed by resampling the srtm dem to three different cell sizes with best results obtained for the higher resolution dem the modification of the original high resolution dem through the process known as stream burning was also tested largely improving the quality of the results finally optimal distance transforms were used for the calculation of distance increments instead of euclidean local distances giving results that were generally better the quality of results for calculated slope values could not be assessed because no reliable data were available for comparison

holocene record of precipitation seasonality from lake calcite δ18o in the central rocky mountains united states a context for recent hydroclimatic extremes and variability is provided by a 10 k y sediment carbonate oxygen isotope 18 o record at 5 100 yr resolution from bison lake 3255 m above sea level in northwestern colorado united states winter precipitation is the primary water source for the alpine headwater lake in the upper colorado river basin and lake water 18 o measurements reflect seasonal variations in precipitation 18 o holocene lake water 18 o variations are inferred from endogenic sedimentary calcite 18 o based on comparisons with historic watershed discharge records and tree ring reconstructions drought periods i e drier winters and or a more rain dominated seasonal precipitation balance generally correspond with higher calcite 18 o values and vice versa early to middle holocene 18 o values are higher implying a rain dominated seasonal precipitation balance lower more variable 18 o values after ca 3500 yr ago indicate a snow dominated but more seasonally variable precipitation balance the middle to late holocene 18 o record corresponds with records of el nino southern oscillation intensification that supports a teleconnection between rocky mountain climate and north pacific sea surface temperatures at decade to century time scales

mid holocene climates of the americas a dynamical response to changed seasonality simulations of the climatic response to mid holocene 6 ka bp orbital forcing with two coupled ocean atmosphere models foam and csm show enhancement of monsoonal precipitation in parts of the american southwest central america and northernmost south america during northern hemisphere summer the enhanced onshore flow that brings precipitation into central america is caused by a northward displacement of the inter tropical convergence zone driven by cooling of the equatorial and warming of the northern subtropical and mid latitude ocean ocean feedbacks also enhance precipitation over the american southwest although the increase in monsoon precipitation there is largely driven by increases in land surface temperature the northward shift in the equatorial precipitation band that causes enhanced precipitation in central america and the american southwest has a negative feedback effect on monsoonal precipitation in northern south america the simulations demonstrate that mid holocene aridity in the mid continent of north america is dynamically linked to the orbitally induced enhancement of the summer monsoon in the american southwest with a spatial structure wet in the southwest and dry in the mid continent similar to that found in strong monsoon years today changes in winter precipitation along the west coast of north america in central america and along the gulf coast caused by southward displacement of the westerly storm tracks indicate that changes in the northern hemisphere winter monsoon also play a role in regional climate changes during the mid holocene although the simulations with foam and csm differ in detail the general mechanisms and patterns are common to both the model results thus provide a coherent dynamical explanation for regional patterns of increased or decreased aridity shown by vegetation lake status and aeolian data from the americas

stability of the astronomical frequencies over the earth s history for paleoclimate studies the expected changes over the past 500 million years in the principal astronomical frequencies influencing the earth s climate may be strong enough to be detectable in the geological records and such effects have been inferred in several cases calculations suggest that the shortening of the earth moon distance and of the length of the day back in time induced a shortening of the fundamental periods for the obliquity and climatic precession from 54 to 35 41 to 29 23 to 19 and 19 to 16 thousand years over the last half billion years at the same time the precessional constant increased from 50 to 61 arc seconds per year the changes in the frequencies of the planetary system due to its chaotic motion are much smaller their influence on the changes of the periods of climatic precession obliquity and eccentricity of the earth s orbit around the sun can be neglected eccentricity periods used for quaternary climate studies may therefore be considered to have been more or less constant for pre quaternary times

abrupt reorganization of north pacific and western north american climate during the last deglaciation dramatic hydroclimate shifts occurred in western north america during the last deglaciation but the timing and mechanisms driving these changes are uncertain and debated and previous modeling has largely relied on linear interpolations between equilibrium snapshot simulations using a published transient climate simulation and a range of proxy records we analyze the region s climate evolution in order to identify the mechanisms governing hydroclimate shifts a rapid loss of ice around 14 000 years ago causes an abrupt reorganization of the circulation which precipitates drying and moistening of southwestern and northwestern north america respectively the atmospheric circulation transitions between two states on a timescale of decades to centuries during which time the westerly jet shifts north by about seven degrees in contrast to previous studies we find that changes in the water budget of western north america prior to this event are not attributable to variations in the position of the jet but rather to the intensity of moisture transport into the continent

omphacite fabric variation in the cabo ortegal eclogite nw spain relationships with strain symmetry during high pressure deformation petrostructural study of cabo ortegal eclogite tectonites provides an insight into the relationships between omphacite lattice preferred orientation and the fundamental deformation parameters in particular strain symmetry microstructural evidence is consistent with a tectonic origin of eclogite fabrics by coeval pyroxene plastic deformation and recrystallization under eclogite facies metamorphic conditions omphacite shape fabrics and garnet spatial distributions are anisotropic and representative of the fields of either apparent flattening or constriction omphacite lattice preferred orientations are dominated by 010 axis maxima normal to the foliation s type or flattening fabrics by 001 axis maxima parallel to the lineation l type or constriction or correspond to other fabric patterns described in the literature omphacite crystallographic fabrics are consistently asymmetric with respect to the structural framework defined by foliation and lineation this suggests that non coaxial deformation components either simple or general shear accompanied eclogite fabric development in addition to any flattening or constriction strain paths and permit shear sense determination these petrostructural features were acquired in the course of eo hercynian ductile deformation probably during deep subduction under 1 5 1 7 gpa confining pressure temperatures above 600 700 c low differential stress 15 20 mpa and slow strain rate 10 12 10 15s 1

mid holocene hydrologic model of the shingobee watershed minnesota abstract a hydrologic model of the shingobee watershed in north central minnesota was developed to reconstruct mid holocene paleo lake levels for williams lake a surface water body located in the southern portion of the watershed hydrologic parameters for the model were first estimated in a calibration exercise using a 9 yr historical record 1990 1998 of climatic and hydrologic stresses the model reproduced observed temporal and spatial trends in surface groundwater levels across the watershed mid holocene aquifer and lake levels were then reconstructed using two paleoclimatic data sets ccm1 atmospheric general circulation model output and pollen transfer functions using sediment core data from williams lake calculated paleo lake levels based on pollen derived paleoclimatic reconstructions indicated a 3 5 m drop in simulated lake levels and were in good agreement with the position of mid holocene beach sands observed in a williams lake sediment core transect however calculated paleolake levels based on ccm1 climate forcing produced only a 0 05 m drop in lake levels we found that decreases in winter precipitation rather than temperature increases had the largest effect on simulated mid holocene lake levels the study illustrates how watershed models can be used to critically evaluate paleoclimatic reconstructions by integrating geologic climatic limnologic and hydrogeologic data sets

pluvial lakes in the great basin of the western united states a view from the outcrop abstract paleo lakes in the western united states provide geomorphic and hydrologic records of climate and drainage basin change at multiple time scales extending back to the miocene recent reviews and studies of paleo lake records have focused on interpretations of proxies in lake sediment cores from the northern and central parts of the great basin in this review emphasis is placed on equally important studies of lake history during the past 30 years that were derived from outcrop exposures and geomorphology in some cases combined with cores outcrop and core records have different strengths and weaknesses that must be recognized and exploited in the interpretation of paleohydrology and paleoclimate outcrops and landforms can yield direct evidence of lake level facies changes that record details of lake level fluctuations and geologic events such as catastrophic floods drainage basin changes and isostatic rebound cores can potentially yield continuous records when sampled in stable parts of lake basins and can provide proxies for changes in lake level water temperature and chemistry and ecological conditions in the surrounding landscape however proxies such as stable isotopes may be influenced by several competing factors the relative effects of which may be difficult to assess and interpretations may be confounded by geologic events within the drainage basin that were unrecorded or not recognized in a core the best evidence for documenting absolute lake level changes lies within the shore nearshore and deltaic sediments that were deposited across piedmonts and at the mouths of streams as lake level rose and fell we review the different shorezone environments and resulting deposits used in such reconstructions and discuss potential estimation errors lake level studies based on deposits and landforms have provided paleohydrologic records ranging from general changes during the past million years to centennial scale details of fluctuations during the late pleistocene and holocene outcrop studies have documented the integration histories of several important drainage basins including the humboldt amargosa owens and mojave river systems that have evolved since the miocene within the active tectonic setting of the great basin these histories have influenced lake levels in terminal basins many pre late pleistocene lakes in the western great basin were significantly larger and record wetter conditions than the youngest lakes outcrop based lake level data provide important checks on core based proxy interpretations we discuss four such comparisons in some cases such as for lakes owens and manix outcrop and core data synthesis yields stronger and more complete records in other cases such as for bonneville and lahontan conflicts point toward reconsideration of confounding factors in interpretation of core based proxies

chronology of lake bonneville 30 000 to 10 000 yr b p abstract 368 radiocarbon ages between 30 and 10 cal ka for samples collected from outcrops and shorelines from the lake bonneville basin have been compiled for this paper samples include 1 organic materials from sources outside the lake such as wood charcoal plant fragments from emergent aquatics and dispersed wetland organics and 2 carbonate materials deposited in the lake such as mollusk shells tufa charophyte debris ooids and marl in general organic materials provide rather than precise ages for lacustrine events but their interpretation is less complicated than interpretation of results from carbonate samples the data set shows evidence of contamination of different carbonate samples with both younger carbon and older carbon for example a radiocarbon reservoir at certain places within the lake during the middle transgressive phase accounts for ages of mollusk shells that are older than basal wood ages at similar altitudes the large number of ages permits an accurate reconstruction of the lake chronology conflicts between ages can be detected and reliable ages can be meaningfully integrated in interpretations if fewer ages were available the chronology might look simpler but its accuracy would be unknown

on the effect of digital elevation model accuracy on hydrology and geomorphology this study compares published cartometric and photogrammetric digital elevation models dems of various grid spacings with a ground truth data set obtained by ground survey and studies the implications of these differences on key hydrologic statistics inferred catchment sizes and stream networks from published dems were found to be significantly different than those from the ground truth in most instances furthermore the width functions and cumulative area relationships determined from the published dems were found to fall consistently outside the 90 confidence limits determined from the ground truth for more than 60 of the relationship suggesting that these hydrologic properties are poorly estimated from published dems however the slope area relationships determined from published dems were found to be less sensitive to catchment shape size and stream network with the relationship falling outside the 90 confidence limits for less than 40 of the relationship for all catchments identified from the published dems a published relationship linking the horizontal resolution with the vertical accuracy of the dem was tested predicting a horizontal resolution of about 10 m for the published dems tested

holocene paleohydrology of central alberta testing the general circulation model climate simulations paleoecological investigations of 28 lakes and bogs provide the basis for a holocene paleohydrological record for central alberta shallow basins empty during the early holocene began flooding shortly after 8000 years ago most filled from 6500 to 4500 years ago and none filled later than 3000 years ago pollen of a hypersaline indicator genus ruppia was found in cores of lakes where the plant does not presently grow indicating lower lake levels evaporation stress and increased salinity from 8000 to 3000 years ago moore lake 54 30 n 100 30 w may have dropped 15 m from 9200 to 5900 years ago and lofty lake 54 44 n 112 29 w largely dried up between 8700 and 6300 years ago lower water levels are suggested for several lakes by shifts in fossil diatom populations from planktonic to benthic particularly epipelic species including pioneering taxa these shallow warm lakes were highly productive and sedimentary pigment levels reached maximum values between about 9000 and 4000 years ago the pale

reconstruction of past lake levels and their relation to groundwater hydrology in the parkers prairie sandplain west central minnesota abstract the parkers prairie sandplain lies along the prairie forest border in west central minnesota an area with ample palynological evidence of significant environmental change during the late pleistocene and holocene past lake levels provide independent evidence of environmental change over the same temporal and spatial scales as pollen analysis we used the method established by digerfeldt 1986 1988 to determine the lowest past levels of three closed basin lakes in the sandplain the method uses a combination of evidence from macrofossils coarse minerogenic matter and organic matter from a transect of cores taken in the near shore region to reconstruct past lake levels lower than present carbonate content magnetic susceptibility and pollen content in the cores provided supplementary evidence and aided cross correlation among the cores each of the three lakes had its lowest level between about 7200 and 6700 yr b p lake levels ranged from 2 8 m to 6 2 m lower than modern levels the greatest lake level lowering occurred for the lake farthest from the river that drains the sandplain the least lake level lowering occurred for the lakes close to the river this spatial pattern of lake level lowering implies that a lowering of the regional water table caused by a reduction in groundwater recharge was the ultimate cause of the lower lake levels the results demonstrate that surface water hydrology is not the sole determinant of lake level and that there are possible dangers in paleoclimatic interpretations of past lake levels based on the assumption of hydrologic isolation of basins

steering of westerly storms over western north america at the last glacial maximum the last glacial maximum hydroclimate over western north america differed from the modern climate a proxy model comparison suggests that the glacial storm track was squeezed and steered by atmospheric high pressure systems

pollen response surface estimates of late quaternary changes in the moisture balance of the northeastern united states abstract quantitative estimates of late quaternary climate in the northeastern united states are reconstructed from fossil pollen data to evaluate changes in the regional moisture balance inferred from water level fluctuations we use environmental response surfaces to calibrate modern pollen data for 17 different taxa to an index of effective soil moisture and mean annual precipitation we apply these response surfaces to fossil pollen data from 60 sites in the region to reconstruct changes in soil moisture and mean annual precipitation at 3000 yr intervals from 12 000 yr b p to present the mapped reconstructions of soil moisture and mean annual precipitation illustrate how the regional moisture balance of the northeast may have changed over the last 12 000 yr in response to changing climate reconstructions of annual precipitation show a gradual increase from 30 below modern values at 12 000 yr b p to near modern values by 6000 yr b p and then remain relatively constant thereafter reconstructed changes in the index of effective soil moisture however show a pattern of near modern values at 12 000 6000 and 3000 yr b p with significantly lower values estimated for 9000 yr b p the time of maximum pine pollen abundances in the northeast this pattern of change is similar to the change in regional moisture balance inferred from stratigraphic records of water level fluctuations these results confirm previous interpretations based on records of water level fluctuations that conditions in the northeast were significantly drier during the early to middle holocene than at other times during the last 12 000 yr

carving and adaptive drainage enforcement of grid digital elevation models 1 an effective and widely used method for removing spurious pits in digital elevation models consists of filling them until they overflow however this method sometimes creates large flat regions which in turn pose a problem for the determination of accurate flow directions in this study we propose to suppress each pit by creating a descending path from it to the nearest point having a lower elevation value this is achieved by carving i e lowering the terrain elevations along the detected path carving paths are identified through a flooding simulation starting from the river outlets the proposed approach allows for adaptive drainage enforcement whereby river networks coming from other data sources are imposed to the digital elevation model only in places where the automatic river network extraction deviates substantially from the known networks an improvement to methods for routing flow over flat regions is also introduced detailed results are presented over test areas of the danube basin

shear deformation and eclogite formation within granulite facies anorthosites of the bergen arcs western norway the anorthosites of the bergen arcs contain two distinct metamorphic mineral assemblages which are stable under different p t conditions a granulite facies assemblage which equilibrated at t 800 900 c p 10 kbar is defined by the minerals plagioclase al rich diopside cpx i garnet gnt i scapolite orthopyroxene dark brown hornblende this assemblage is linked by corona structures to a primary magmatic mineralogy where coexistence of olivine and plagioclase limits the maximum crystallization pressure of the anorthosites to 8 kbar at 1200 c the anorthosites are transected by a series of fine grained shear zones in which the granulite facies assemblage is replaced by an eclogite mineralogy defined by omphacite cpx ii jd40 50 garnet kyanite epidote clinozoisite phengite quartz amphibole with plagioclase sometimes remaining in the marginal parts mineral compositions indicate equilibration conditions for the eclogite paragenesis of t 700 800 c p 16 19 kbar a series of textural and mineralogical changes which can be observed as the strain increases into the shear zones demonstrates that the eclogite mineralogy is formed at the expense of the granulite facies minerals during the early stages of shear zone development the eclogite mineralogy defines a strong mylonitic fabric which again is truncated by centimetres thick secondary shear zones containing clasts of both granulite and eclogite minerals these latter fine grained ultramylonites may be connected to uplift of this deep seated complex the eclogite shear zones are highly oblique to the granulite structures indicating that the eclogite shear zones formed under a new stress regime rb sr dating shows that the granulite facies assemblage formed in sveconorwegian time or earlier while minerals from shear zones define ages around 400 ma indicating a caledonian formation of the eclogites this work demonstrates the critical role of deformation in triggering metamorphic reactions whether the shear stress also was important in stabilizing the eclogite assemblages is not clear from the available data

last glacial maximum and holocene lake levels of owens lake eastern california usa abstract stratigraphic investigations of fluvio deltaic and lacustrine sediments exposed in stream cuts quarry walls and deep trenches east of the sierra nevada in owens valley near lone pine california have enabled the reconstruction of pluvial owens lake level oscillations age control for these sediments is from 22 radiocarbon 14 c dates and the identification and stratigraphic correlation of a tephra which when plotted as a function of age versus altitude define numerous oscillations in the level of pluvial owens lake during the latest pleistocene and early holocene we have constructed a lake level altitude curve for the time interval 27 000 cal yr bp to present that is based on the integration of this new stratigraphic analysis with published surface stratigraphic data and subsurface core data pluvial owens lake regressed from its latest pleistocene highstands from 27 000 to 15 300 cal yr bp as recorded by 15 m of down cutting of the sill from the altitudes of 1160 to 1145 m by 11 600 cal yr bp the lake had dropped 45 m from the 1145 m sill this lowstand was followed by an early holocene transgression that attained a highstand near 1135 m before dropping to 1120 m at 7860 7650 cal yr bp that had not been recognized in earlier studies the lake then lowered another 30 m to shallow and near desiccation levels between 6850 and 4300 cal yr bp fluvial cut and fill relations north of lone pine and well preserved shoreline features at 1108 m indicate a minor lake level rise after 4300 cal yr bp followed by alkaline and shallow conditions during the latest holocene the new latest quaternary lake level record of pluvial owens lake offers insight to the hydrologic balance along the east side of the southern sierra nevada and will assist regional paleoclimatic models for the western basin and range

reconciling divergent trends and millennial variations in holocene temperatures cooling during most of the past two millennia has been widely recognized and has been inferred to be the dominant global temperature trend of the past 11 700 years the holocene epoch however long term cooling has been difficult to reconcile with global forcing and climate models consistently simulate long term warming the divergence between simulations and reconstructions emerges primarily for northern mid latitudes for which pronounced cooling has been inferred from marine and coastal records using multiple approaches here we show that temperatures reconstructed from sub fossil pollen from 642 sites across north america and europe closely match simulations and that long term warming not cooling defined the holocene until around 2 000 years ago the reconstructions indicate that evidence of long term cooling was limited to north atlantic records early holocene temperatures on the continents were more than two degrees celsius below those of the past two millennia consistent with the simulated effects of remnant ice sheets in the climate model community climate system model 3 ccsm3 ccsm3 simulates increases in growing degree days a measure of the accumulated warmth above five degrees celsius per year of more than 300 kelvin days over the holocene consistent with inferences from the pollen data it also simulates a decrease in mean summer temperatures of more than two degrees celsius which correlates with reconstructed marine trends and highlights the potential importance of the different subseasonal sensitivities of the records despite the differing trends pollen and marine based reconstructions are correlated at millennial to centennial scales probably in response to ice sheet and meltwater dynamics and to stochastic dynamics similar to the temperature variations produced by ccsm3 although our results depend on a single source of palaeoclimatic data pollen and a single climate model simulation they reinforce the notion that climate models can adequately simulate climates for periods other than the present day they also demonstrate that amplified warming in recent decades increased temperatures above the mean of any century during the past 11 000 years

ebsd measured lattice preferred orientations and seismic properties of eclogites we investigated the deformation mechanisms and the seismic properties of 10 eclogite samples from different localities alps norway mali and eastern china through the analysis of their microstructures and lattice preferred orientations lpo these samples are representative of various types and intensity of deformation under eclogitic metamorphic conditions omphacite and garnet lpo were determined from electron backscatter diffraction ebsd technique garnet appears to be almost randomly oriented whereas omphacite develops strong lpo characterized by the 001 axes concentrated sub parallel to the lineation and the 010 poles concentrated sub perpendicular to the foliation in order to analyze the deformation mechanisms that produced such omphacite lpo we compare our observations to lpo simulated by viscoplastic self consistent numerical models a good fit to the measured lpo is obtained for models in which the dominant slip systems are 1 2h110i 11 0 001 110 and 001 100 dominant activation of these slip systems is in agreement with tem studies of naturally deformed omphacite seismic properties of eclogite are calculated by combining the measured lpo and the single crystal elastic constants of omphacite and garnet although eclogite seismic anisotropies are very weak less than 3 for both p and s wave they are generally characterized by a maximum p wave velocity sub parallel to the lineation and by a minimum velocity approximately normal to foliation the mean p and s wave velocities are high respectively 8 6 and 4 9 km s the s wave anisotropy pattern displays complex relationships with the structural frame but the fast polarization plane generally tends to be parallel to the foliation calculated reflection coefficients show that an eclogite crust interface is generally a good reflector rc 0 1 whereas an eclogite body embedded in the upper mantle would be hardly detectable

estimating cell to cell land surface drainage paths from digital channel networks with an application to the amazon basin cell to cell surface flow paths are commonly derived from gridded digital elevation models dem by choosing the direction of steepest descent to one of the eight surrounding cells however adequate dems often are not available we developed a topography independent method for creating gridded land and stream drainage direction maps based on corrected vector river networks we applied it to the digital chart of the world river network in the amazon basin gridded at 0 0058 resolution in this basin low relief and poor topographic data have prevented the effective use of dem based methods we geo registered 224 hydrographic gages against the processed network and compared extracted vs published drainage areas drainage areas ranged from 227 to 4 620 000 km 2 median relative error was 4 5 increasing in smaller basins to 94 in basins 2000 km 2 the effective limit of reliability may differ from 2000 km 2 across the basin the drainage direction map and derivative datasets represent an improvement over existing datasets for regional research in the amazon basin methods exploiting vector networks complement terrain approaches and combined they may yield advances in the automated extraction of drainage maps and handling of topologically realistic river systems q 2005 elsevier b v all rights reserved

patterns of hydroclimatic change in the rocky mountains and surrounding regions since the last glacial maximum abstract the paleohydrologic record of western north america since the last glacial maximum reveals a wide range of hydroclimatic variability in time and space to improve the understanding of abrupt hydroclimatic shifts and millennial scale hydrologic changes in the central rocky mountains we reconstruct the lake level histories of two small lakes in the beartooth and bighorn mountains in northern wyoming over the past 17 ka to do so we use ground penetrating radar gpr and sediment cores to track the elevations of shoreline sediments within the lakes through time we compare the stratigraphies with those from four other lakes in wyoming and colorado and find widespread evidence for a terminal pleistocene drought from 15 to 11 ka an early holocene humid period from 11 to 8 ka and mid holocene aridity from 8 to 5 5 ka the northern wyoming lakes also provide evidence of high levels in the pleistocene possibly before ca 15 ka and rapid hydroclimatic changes that may have correlated with heinrich event 1 ca 16 8 ka we place the changes in a broad context by summarizing and mapping water level changes from 107 additional previously studied lakes important patterns include 1 extensive drying across the western u s after 15 ka 2 sub regional differences during the pleistocene holocene transition 3 a north south contrast from 9 to 6 ka consistent with a northward shift in storm tracks as the influence of the laurentide ice sheet diminished and 4 rapid increases in effective moisture across much of western north america from 6 to 4 ka

radiocarbon chronology of lake bonneville eastern great basin usa abstract lake bonneville occupied a series of connected topographically closed structural basins in the eastern great basin from about 30 ka to 12 ka the following synthesis of lake bonneville history is based on a critical evaluation of the stratigraphic and geomorphic contexts of 83 radiocarbon ages of a variety of samples including wood charcoal dispersed organic matter mollusk shells and tufa the lake began to rise from levels close to average holocene levels after about 28 ka by 22 ka it had transgressed approximately 100 m between 22 and 20 ka it regressed about 45 m in the stansbury oscillation and the stansbury shoreline was formed transgression after 20 ka proceeded in two phases a rapid phase from 20 to 18 ka and a slower phase from 18 to 15 ka the lake overflowed intermittently at its highest level the bonneville shoreline from about 15 to 14 5 ka then catastrophically dropped 100 m during the bonneville flood to the level of the provo shoreline which it occupied until about 14 ka subsequent closed basin regression was rapid and complete by 12 ka and was followed by a modest transgression to form the gilbert shoreline between 10 9 and 10 3 ka the lake bonneville record is an accurate proxy of the changing water balance in the bonneville basin during the late pleistocene although the nature of the climatic changes during this period are still uncertain

enhanced late holocene enso pdo expression along the margins of the eastern north pacific abstract pacific climate is known to have varied during the holocene but spatial patterns remain poorly defined this paper compiles terrestrial and marine proxy data from sites along the northeastern pacific margins and proposes that they indicate 1 suppressed enso conditions during the middle holocene between 8000 and 4000 cal bp with a north pacific that generally resembled a la nina like or more negative pdo phase and 2 a climate transition between 4200 and 3000 cal bp that appears to be the teleconnected expression to a more modern like enso pacific compared to modern day conditions the compiled data suggest that during the middle holocene the aleutian low was generally weaker during the winter and or located more to the west while the north pacific high was stronger during the summer and located more to the north coastal upwelling off california was more enhanced during the summer and fall but suppressed during the spring oregon and california sea surface temperatures ssts were cooler the santa barbara basin had an anomalous record suggesting warmer ssts late holocene records indicate a more variable el nino like and more positive pdo pacific the aleutian low became more intensified during the winter and or located more to the east the north pacific high became weaker and or displaced more to the south coastal upwelling off california intensified during the spring but decreased during the fall oregon and california ssts became warmer recording the shoreward migration of sub tropical gyre waters during the fall while spring upwelling cooler sst increased in the santa barbara basin the high resolution proxy records indicate enhanced enso and pdo variability after 4000 cal bp off southern california 3400 cal bp off northern california and by 2000 cal bp in southwestern yukon a progressively northward migration of the enso teleconnection during the late holocene is proposed

mid to late holocene climate change an overview the last 6000 years are of particular interest to the understanding of the earth system because the boundary conditions of the climate system did not change dramatically in comparison to larger glacial interglacial changes and because abundant detailed regional palaeoclimatic proxy records cover this period we use selected proxy based reconstructions of different climate variables together with state of the art time series of natural forcings orbital variations solar activity variations large tropical volcanic eruptions land cover and greenhouse gases underpinned by results from general circulation models gcms and earth system models of intermediate complexity emics to establish a comprehensive explanatory framework for climate changes from the mid holocene mh to pre industrial time the redistribution of solar energy due to orbital forcing on a millennial timescale was the cause of a progressive southward shift of the northern hemisphere nh summer position of the intertropical convergence zone itcz this was accompanied by a pronounced weakening of the monsoon systems in africa and asia and increasing dryness and desertification on both continents the associated summertime cooling of the nh combined with changing temperature gradients in the world oceans likely led to an increasing amplitude of the el nin o southern oscillation enso and possibly increasingly negative north atlantic oscillation nao indices up to the beginning of the last millennium on decadal to multi century timescales a worldwide coincidence between solar irradiance minima tropical volcanic eruptions and decadal to multi century scale cooling events was not found however reconstructions show that widespread decadal to multi century scale cooling events accompanied by advances of mountain glaciers occurred in the nh e g in scandinavia and the european alps this occurred namely during the little ice age lia between ad w1350 and 1850 when the lower summer insolation in the nh due to orbital forcing coincided with solar activity minima and several strong tropical volcanic eruptions the role of orbital forcing in the nh cooling the southward itcz shift and the desertifi cation of the sahara are supported by numerous model simulations other simulations have suggested that the fingerprint of solar activity variations should be strongest in the tropics but there is also evidence that changes in the ocean heat transport took place during the lia at high northern latitudes with possible additional implications for climates of the southern hemisphere sh

late glacial and holocene vegetation history and paleoclimate of the kaibab plateau arizona sediment cores spanning the last 13 500 calendar years cal yr were obtained from two lakes fracas lake 2518 m bear lake 2778 m on the kaibab plateau in northern arizona pollen and plant macrofossil records indicate that before 12 900 cal yr b p high elevation landscapes of the kaibab plateau near bear lake were covered by alpine tundra while lower elevations near fracas lake were occupied by picea woodland at12 900 cal yr b p picea engelmannii and abies lasiocarpa forest expanded upward to occupy the top of the plateau pinus ponderosa arrived near fracas lake 11 000 cal yr b p replacing picea forests since then pinus ponderosa forest has dominated the fracas lake area pinus ponderosa did not appear at bear lake until9730 cal yr b p mixed forests of picea mainly picea pungens abies lasiocarpa pinus ponderosa a ndpseudotsuga after8000 cal yr b p grew near bear lake for the remainder of the holocene picea engelmannii populations reexpanded near bear lake after 4000 cal yr b p charcoal records indicate that fire probably helpedpinus ponderosa to become established near bear lake climate changes on the kaibab plateau since the late glacial were inferred from lake levels and vegetation patterns the late glacial 11 000 cal yr b p was cold and probably wet the early holocene 11 000 to8000 cal yr b p was cooler than today and may have been the wettest period fracas lake and bear lake were probably deepest then during this period a strengthened summer monsoon brought in more moisture from the eastern pacific ocean and the gulf of mexico during the dry and warm mid holocene fracas lake and bear lake experienced water level declines the late holocene was relatively wet and cool again and aquatic plants were abundant in the two lakes increasing effective moisture in the late holocene was related to decreasing summer insolation 1999 elsevier science b v all rights reserved

drainage networks from grid digital elevation models current algorithms that deduce the drainage network from a digital elevation model dem represented by a regular array of surface elevations share a fault unless the terrain is rugged the derived water channels tend to flow in parallel lines along preferred directions engendered by the sampling grid orientation we present a simple solution to the problem a second difficulty is the presence of noise that creates artificial pits we briefly describe a method which deals with pits in what we believe to be a more efficient manner for virtual memory environments than previous efforts our system has treated dems of nearly 9 000 000 pixels we show how depth first search of the resulting drainage network permits segmentation of the dem into basins by various criteria analysis of stream sediment anomaly dilution profiles improved hydrological models and other applications

extraction and representation of nested catchment areas from digital elevation models in lake dominated topography this paper presents a new method for extracting flow directions contributing upslope areas and nested catchments from digital elevation models in lake dominated areas existing tools for acquiring descriptive variables of the topography such as surface flow directions and contributing areas were developed for moderate to steep topography these tools are typically difficult to apply in gentle topography owing to limitations in explicitly handling lakes and other flat areas this paper addresses the problem of accurately representing general topographic features by first identifying distinguishing features such as lakes in gentle topography areas and then using these features to guide the search for topographic flow directions and catchment marking lakes are explicitly represented in the topology of a watershed for use in water routing nonlake flat features help guide the search for topographic flow directions in areas of low signal to noise this combined feature based and grid based search for topographic features yields improved contributing areas and watershed boundaries where there are lakes and other flat areas lakes are easily classified from remotely sensed imagery which makes automated representation of lakes as subsystems within a watershed system tractable with widely available data sets

how does modifying a dem to reflect known hydrology affect subsequent terrain analysis summary many digital elevation models dems have difficulty replicating hydrological patterns in flat landscapes efforts to improve dem performance in replicating known hydrology have included a variety of soft i e algorithm based approaches and hard techniques such as stream burning or surface reconditioning e g agree or anudem using a representation of the known stream network these methods trench or mathematically warp the original dem to improve how accurately stream position stream length and catchment boundaries replicate known hydrological conditions however these techniques permanently alter the dem and may affect further analyses e g slope this paper explores the impact that commonly used hydrological correction methods stream burning agree aml and anudem v4 6 3 and anudem v5 1 have on the overall nature of a dem finding that different methods produce non convergent outcomes for catchment parameters such as catchment boundaries stream position and length and differentially compromise secondary terrain analysis all hydrological correction methods successfully improved calculation of catchment area stream position and length as compared to using the dem without any modification but they all increased catchment slope no single method performing best across all categories different hydrological correction methods changed elevation and slope in different spatial patterns and magnitudes compromising the ability to derive catchment parameters and conduct secondary terrain analysis from a single dem modification of a dem to better reflect known hydrology can be useful however knowledge of the magnitude and spatial pattern of the changes are required before using a dem for subsequent analyses

deformation induced clinopyroxene fabrics from eclogites abstract clinopyroxene omphacite microfabric studies in eclogites are useful in characterizing the deformation that has affected deeper parts of subduction zones we present the microfabrics and tem structures of four samples representative of the various deformation types known in eclogites the 001 axis tends to be parallel to the lineation while the 010 plane is close to the foliation but these general features undergo strong variations in relation to the shape fabric in strongly lineated samples 010 tends to be spread out in a great circle normal to the lineation whereas 001 is dispersed in the foliation in the case of strongly foliated samples the slip systems identified by tem are 001 100 001 110 and 1 2 110 110 operation of dislocation creep is testified by the presence of dislocations and dislocation walls but this process alone cannot explain all the observed features particularly some aspects of the crystallographic preferred orientation hence other deformation mechanisms must be involved these may be dynamic migration recrystallization by grain boundary migration as well as mass transfer processes this study illustrates that not every crystallographic mineral preferred orientation can be explained by dislocation glide and or twinning processes alone in addition it also demonstrates that eclogite microfabrics can be used to unravel deformation parameters in subduction and or collision zones

the assignment of drainage direction over flat surfaces in raster digital elevation models abstract drainage direction assignment over flat surfaces in raster digital elevation models dem has been a stubborn problem for dem processing methods a new approach that improves upon current methods of flat surface treatment is presented the approach is based on the recognition that in natural landscapes drainage is generally away from higher and towards lower terrain to produce such a drainage dem elevations of a flat surface are modified to impose two gradients one away from higher terrain and one towards lower terrain subsequent processing of the dem produces a drainage pattern over the flat surface that is topographically consistent and exhibits flow convergence properties the proposed approach is not restricted by the shape of the flat surface the number of outlets on its edge or the complexity of the surrounding topography a comparison with the drainage pattern of an established method that displays the parallel flow problem shows significant improvements in producing realistic drainage patterns the proposed approach extends automated dem processing to digital landscapes for which existing methods cannot provide adequate flow directions over flat surfaces to conduct a drainage analysis the necessary algorithmic details for implementation of the approach are provided 1997 elsevier science b v

the prediction of hillslope flow paths for distributed hydrological modelling using digital terrain models the accuracy of the predictions of distributed hydrological models must depend in part on the proper specification of flow pathways this paper examines some of the problems of deriving flow pathways from raster digital terrain data in the context of hydrological predictions using topmodel distributed moisture status is predicted in topmodel on the basis of spatial indices that depend on flow path definition the sensitivity of this index to flow path algorithm and grid size is examined for the case where the surface topography is a good indicator of local hydraulic gradients a strategy for the case where downslope subsurface flow pathways may deviate from those indicated by the surface topography is described with an example application

path based methods for the determination of nondispersive drainage directions in grid based digital elevation models 1 path based methods for the determination of nondispersive drainage directions in grid based digital elevation models are presented these methods extend the descriptive capabilities of the classical d8 method by cumulating the deviations between selected and theoretical drainage directions along the drainage paths it is shown that either angular or transversal deviations can be employed accordingly two classes of methods designated d8 lad eight drainage directions least angular deviation and d8 ltd eight drainage directions least transversal deviation are developed detailed tests on four synthetic drainage systems of known geometry and on the liro catchment central italian alps indicate that the proposed methods provide significant improvement over the d8 method for the determination of drainage directions and drainage areas