SEARCH · Search USGS

Results for “Journal of Geophysical Research - Atmospheres”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25Linked to original sources

USGS2005

Interferometric synthetic aperture radar study of Okmok volcano, Alaska, 1992-2003: Magma supply dynamics and postemplacement lava flow deformation

Okmok volcano, located in the central Aleutian arc, Alaska, is a dominantly basaltic complex topped with a 10-km-wide caldera that formed circa 2.05 ka. Okmok erupted several times during the 20th century, most recently in 1997; eruptions in 1945, 1958, and 1997 produced lava flows within the caldera. We used 80 interferometric synthetic aperture radar (InSAR) images (interferograms) to study transient deformation of the volcano before, during, and after the 1997 eruption. Point source models suggest that a magma reservoir at a depth of 3.2 km below sea level, located beneath the center of the caldera and about 5 km northeast of the 1997 vent, is responsible for observed volcano-wide deformation. The preeruption uplift rate decreased from about 10 cm yr −1 during 1992–1993 to 2 ∼ 3 cm yr −1 during 1993–1995 and then to about −1 ∼ −2 cm yr −1 during 1995–1996. The posteruption inflation rate generally decreased with time during 1997–2001, but increased significantly during 2001–2003. By the summer of 2003, 30 ∼ 60% of the magma volume lost from the reservoir in the 1997 eruption had been replenished. Interferograms for periods before the 1997 eruption indicate consistent subsidence of the surface of the 1958 lava flows, most likely due to thermal contraction. Interferograms for periods after the eruption suggest at least four distinct deformation processes: (1) volcano-wide inflation due to replenishment of the shallow magma reservoir, (2) subsidence of the 1997 lava flows, most likely due to thermal contraction, (3) deformation of the 1958 lava flows due to loading by the 1997 flows, and (4) continuing subsidence of 1958 lava flows buried beneath 1997 flows. Our results provide insights into the postemplacement behavior of lava flows and have cautionary implications for the interpretation of inflation patterns at active volcanoes.

Alaska↗

USGS2004

Quantifying probabilities of volcanic events: The example of volcanic hazard at Mount Vesuvius

We describe an event tree scheme to quantitatively estimate both long- and short-term volcanic hazard. The procedure is based on a Bayesian approach that produces a probability estimation of any possible event in which we are interested and can make use of all available information including theoretical models, historical and geological data, and monitoring observations. The main steps in the procedure are (1) to estimate an a priori probability distribution based upon theoretical knowledge, (2) to modify that using past data, and (3) to modify it further using current monitoring data. The scheme allows epistemic and aleatoric uncertainties to be dealt with in a formal way, through estimation of probability distributions at each node of the event tree. We then describe an application of the method to the case of Mount Vesuvius. Although the primary intent of the example is to illustrate the methodology, one result of this application merits special mention. The present emergency response plan for Mount Vesuvius is referenced to a maximum expected event (MEE), the largest out of all the possible eruptions within the next few decades. Our calculation suggest that there is a nonnegligible (1-20%) chance that the next eruption could be larger than that stipulated in the present MEE. The methodology allows all assumptions and thresholds to be clearly identified and provides a rational means for their revision if new data or information are obtained. Copyright 2004 by the American Geophysical Union.

Mount Vesuvius↗

USGS2005

Significance of stress transfer in time-dependent earthquake probability calculations

A sudden change in stress is seen to modify earthquake rates, but should it also revise earthquake probability? Data used to derive input parameters permits an array of forecasts; so how large a static stress change is require to cause a statistically significant earthquake probability change? To answer that question, effects of parameter and philosophical choices are examined through all phases of sample calculations, Drawing at random from distributions of recurrence-aperiodicity pairs identifies many that recreate long paleoseismic and historic earthquake catalogs. Probability density funtions built from the recurrence-aperiodicity pairs give the range of possible earthquake forecasts under a point process renewal model. Consequences of choices made in stress transfer calculations, such as different slip models, fault rake, dip, and friction are, tracked. For interactions among large faults, calculated peak stress changes may be localized, with most of the receiving fault area changed less than the mean. Thus, to avoid overstating probability change on segments, stress change values should be drawn from a distribution reflecting the spatial pattern rather than using the segment mean. Disparity resulting from interaction probability methodology is also examined. For a fault with a well-understood earthquake history, a minimum stress change to stressing rate ratio of 10:1 to 20:1 is required to significantly skew probabilities with >80-85% confidence. That ratio must be closer to 50:1 to exceed 90-95% confidence levels. Thus revision to earthquake probability is achievable when a perturbing event is very close to the fault in question or the tectonic stressing rate is low.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Rupture dynamics with energy loss outside the slip zone

Energy loss in a fault damage zone, outside the slip zone, contributes to the fracture energy that determines rupture velocity of an earthquake. A nonelastic two-dimensional dynamic calculation is done in which the slip zone is modeled as a fault plane and material off the fault is subject to a Coulomb yield condition. In a mode 2 crack-like solution in which an abrupt uniform drop of shear traction on the fault spreads from a point, Coulomb yielding occurs on the extensional side of the fault. Plastic strain is distributed with uniform magnitude along the fault, and it has a thickness normal to the fault proportional to propagation distance. Energy loss off the fault is also proportional to propagation distance, and it can become much larger than energy loss on the fault specified by the fault constitutive relation. The slip velocity function could be produced in an equivalent elastic problem by a slip-weakening friction law with breakdown slip Dc increasing with distance. Fracture energy G and equivalent Dc will be different in ruptures with different initiation points and stress drops, so they are not constitutive properties; they are determined by the dynamic solution that arrives at a particular point. Peak slip velocity is, however, a property of a fault location. Nonelastic response can be mimicked by imposing a limit on slip velocity on a fault in an elastic medium.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

New constraints on mechanisms of remotely triggered seismicity at Long Valley Caldera

Regional-scale triggering of local earthquakes in the crust by seismic waves from distant main shocks has now been robustly documented for over a decade. Some of the most thoroughly recorded examples of repeated triggering of a single site from multiple, large earthquakes are measured in geothermal fields of the western United States like Long Valley Caldera. As one of the few natural cases where the causality of an earthquake sequence is apparent, triggering provides fundamental constraints on the failure processes in earthquakes. We show here that the observed triggering by seismic waves is inconsistent with any mechanism that depends on cumulative shaking as measured by integrated energy density. We also present evidence for a frequency-dependent triggering threshold. On the basis of the seismic records of 12 regional and teleseismic events recorded at Long Valley Caldera, long-period waves (>30 s) are more effective at generating local seismicity than short-period waves of comparable amplitude. If the properties of the system are stationary over time, the failure threshold for long-period waves is ~0.05 cm/s vertical shaking. Assuming a phase velocity of 3.5 km/s and an elastic modulus of 3.5 x 10 10 Pa, the threshold in terms of stress is 5 kPa. The frequency dependence is due in part to the attenuation of the surface waves with depth. Fluid flow through a porous medium can produce the rest of the observed frequency dependence of the threshold. If the threshold is not stationary with time, pore pressures that are >99.5% of lithostatic and vary over time by a factor of 4 could explain the observations with no frequency dependence of the triggering threshold.

California↗

USGS2005

Evolution of melt-vapor surface tension in silicic volcanic systems: Experiments with hydrous melts

We evaluate the melt‐vapor surface tension (σ) of natural, water‐saturated dacite melt at 200 MPa, 950–1055°C, and 4.8–5.7 wt % H 2 O. We experimentally determine the critical supersaturation pressure for bubble nucleation as a function of dissolved water and then solve for σ at those conditions using classical nucleation theory. The solutions obtained give dacite melt‐vapor surface tensions that vary inversely with dissolved water from 0.042 (±0.003) J m −2 at 5.7 wt % H 2 O to 0.060 (±0.007) J m −2 at 5.2 wt % H 2 O to 0.073 (±0.003) J m −2 at 4.8 wt % H 2 O. Combining our dacite results with data from published hydrous haplogranite and high‐silica rhyolite experiments reveals that melt‐vapor surface tension also varies inversely with the concentration of mafic melt components (e.g., CaO, FeO total , MgO). We develop a thermodynamic context for these observations in which melt‐vapor surface tension is represented by a balance of work terms controlled by melt structure. Overall, our results suggest that cooling, crystallization, and vapor exsolution cause systematic changes in σ that should be considered in dynamic modeling of magmatic processes.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

Stress orientations at intermediate angles to the San Andreas Fault, California

There are currently two competing models for the frictional strength of the San Andreas Fault in California: the strong-fault model and the weak-fault model. The strong-fault model predicts the maximum horizontal compressive stress axis to be at low angles to the fault, while the relatively weak fault model predicts it to be at high angles. Previous studies have disagreed as to which model is supported by observed stress orientations. We review and compare these studies and present results from several new focal mechanism stress inversions. We find that the observed stress orientations of different studies are generally consistent, implying that the disagreement is one of interpretation. The majority of studies find compressive stress orientations at intermediate angles to the fault, not strictly consistent with either current model. The strong-fault model is acceptable if the San Andreas is assumed to be a nonoptimally orientated fault that fails because optimally oriented, preexisting planes are not present. The relatively weak fault model is not consistent with the stress orientations. We propose two alternative models to better explain the observed intermediate stress orientations: an intermediate-strength San Andreas model and a model in which all major active faults are weak.

California↗

USGS2005

Forearc structure beneath southwestern British Columbia: A three-dimensional tomographic velocity model

This paper presents a three-dimensional compressional wave velocity model of the forearc crust and upper mantle and the subducting Juan de Fuca plate beneath southwestern British Columbia and the adjoining straits of Georgia and Juan de Fuca. The velocity model was constructed through joint tomographic inversion of 50,000 first-arrival times from earthquakes and active seismic sources. Wrangellia rocks of the accreted Paleozoic and Mesozoic island arc assemblage underlying southern Vancouver Island in the Cascadia forearc are imaged at some locations with higher than average lower crustal velocities of 6.5-7.2 km/s, similar to observations at other island arc terranes. The mafic Eocene Crescent terrane, thrust landward beneath southern Vancouver Island, exhibits crustal velocities in the range of 6.0-6.7 km/s and is inferred to extend to a depth of more than 20 km. The Cenozoic Olympic Subduction Complex, an accretionary prism thrust beneath the Crescent terrane in the Olympic Peninsula, is imaged as a low-velocity wedge to depths of at least 20 km. Three zones with velocities of 7.0-7.5 km/s, inferred to be mafic and/or ultramafic units, lie above the subducting Juan de Fuca plate at depths of 25-35 km. The forearc upper mantle wedge beneath southeastern Vancouver Island and the Strait of Georgia exhibits low velocities of 7.2-7.5 km/s, inferred to correspond to ???20% serpentinization of mantle peridotites, and consistent with similar observations in other warm subduction zones. Estimated dip of the Juan de Fuca plate beneath southern Vancouver Island is ???11??, 16??, and 27?? at depths of 30, 40, and 50 km, respectively. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

The calcite → aragonite transformation in low-Mg marble: Equilibrium relations, transformations mechanisms, and rates

Experimental transformation of a rather pure natural calcite marble to aragonite marble did not proceed via the expected straightforward polymorphic replacement. Instead, the small amount of Mg in the starting material (0.36 wt %) was excluded from the growing aragonite and diffused preferentially into the remaining calcite grains, producing Mg-rich calcite rods that persisted as relicts. Nucleation of aragonite occurred exclusively on grain boundaries, with aragonite [001] oriented subparallel to calcite [0001]. The aragonite crystals preferentially consumed the calcite crystal on which they nucleated, and the reaction fronts developed preferentially along the {010} and {110} planes of aragonite. Each aragonite neoblast that grew was nearly free of Mg (typically <0.1 wt %). The excess Mg was taken up by the calcite grains in between, stabilizing them and causing a few volume percent rodlike relicts of Mg-enriched calcite (up to 10 wt % MgO) to be left behind by the advancing reaction front. The aragonite growth rates are approximately linear and range from &sim;3 × 10 −11 m s −1 at 600°C to &sim;9 × 10 −9 m s −1 at 850°C, with an apparent activation enthalpy of 166 ± 91 kJ mol −1 . This reaction mechanism and the resultant texture are akin to cellular precipitation reactions in metals. Similar transformation textures have been reported from high-Mg marbles in Japan and China that disproportionated to low-Mg calcite and dolomite.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Magnitude and location of historical earthquakes in Japan and implications for the 1855 Ansei Edo earthquake

Japan Meteorological Agency (JMA) intensity assignments IJMA are used to derive intensity attenuation models suitable for estimating the location and an intensity magnitude Mjma for historical earthquakes in Japan. The intensity for shallow crustal earthquakes on Honshu is equal to -1.89 + 1.42MJMA - 0.00887?? h - 1.66log??h, where MJMA is the JMA magnitude, ??h = (??2 + h2)1/2, and ?? and h are epicentral distance and focal depth (km), respectively. Four earthquakes located near the Japan Trench were used to develop a subducting plate intensity attenuation model where intensity is equal to -8.33 + 2.19MJMA -0.00550??h - 1.14 log ?? h. The IJMA assignments for the MJMA7.9 great 1923 Kanto earthquake on the Philippine Sea-Eurasian plate interface are consistent with the subducting plate model; Using the subducting plate model and 226 IJMA IV-VI assignments, the location of the intensity center is 25 km north of the epicenter, Mjma is 7.7, and MJMA is 7.3-8.0 at the 1?? confidence level. Intensity assignments and reported aftershock activity for the enigmatic 11 November 1855 Ansei Edo earthquake are consistent with an MJMA 7.2 Philippine Sea-Eurasian interplate source or Philippine Sea intraslab source at about 30 km depth. If the 1855 earthquake was a Philippine Sea-Eurasian interplate event, the intensity center was adjacent to and downdip of the rupture area of the great 1923 Kanto earthquake, suggesting that the 1855 and 1923 events ruptured adjoining sections of the Philippine Sea-Eurasian plate interface.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Response of seismicity to Coulomb stress triggers and shadows of the 1999 Mw=7.6 Chi-Chi, Taiwan, earthquake

The correlation between static Coulomb stress increases and aftershocks has thus far provided the strongest evidence that stress changes promote seismicity, a correlation that the Chi-Chi earthquake well exhibits. Several studies have deepened the argument by resolving stress changes on aftershock focal mechanisms, which removes the assumption that the aftershocks are optimally oriented for failure. Here one compares the percentage of planes on which failure is promoted after the main shock relative to the percentage beforehand. For Chi-Chi we find a 28% increase for thrust and an 18% increase for strike-slip mechanisms, commensurate with increases reported for other large main shocks. However, perhaps the chief criticism of static stress triggering is the difficulty in observing predicted seismicity rate decreases in the stress shadows, or sites of Coulomb stress decrease. Detection of sustained drops in seismicity rate demands a long catalog with a low magnitude of completeness and a high seismicity rate, conditions that are met at Chi-Chi. We find four lobes with statistically significant seismicity rate declines of 40-90% for 50 months, and they coincide with the stress shadows calculated for strike-slip faults, the dominant faulting mechanism. The rate drops are evident in uniform cell calculations, 100-month time series, and by visual inspection of the M ??? 3 seismicity. An additional reason why detection of such declines has proven so rare emerges from this study: there is a widespread increase in seismicity rate during the first 3 months after Chi-Chi, and perhaps many other main shocks, that might be associated with a different mechanism. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Distribution of stress drop, stiffness, and fracture energy over earthquake rupture zones

Using information provided by slip models and the methodology of McGarr and Fletcher (2002), we map static stress drop, stiffness (k = ????/u, where ???? is static stress drop and u is slip), and fracture energy over the slip surface to investigate the earthquake rupture process and energy budget. For the 1994 M6.7 Northridge, 1992 M7.3 Landers, and 1995 M6.9 Kobe earthquakes, the distributions of static stress drop show strong heterogeneity, emphasizing the importance of asperities in the rupture process. Average values of static stress drop are 17, 11, and 4 Mpa for Northridge, Landers, and Kobe, respectively. These values are substantially higher than estimates based on simple crack models, suggesting that the failure process involves the rupture of asperities within the larger fault zone. Stress drop as a function of depth for the Northridge and Landers earthquakes suggests that stress drops are limited by crustal strength. For these two earthquakes, regions of high slip are surrounded by high values of stiffness. Particularly for the Northridge earthquake, the prominent patch of high slip in the central part of the fault is bordered by a ring of high stiffness and is consistent with expectations based on the failure of an asperity loaded at its edge due to exterior slip. Stiffness within an asperity is inversely related to its dimensions. Estimates of fracture energy, based on static stress drop, slip, and rupture speed, were used to investigate the nature of slip weakening at four locations near the hypocenter of the Kobe earthquake for comparison with independent results based on a dynamic model of this earthquake. One subfault updip and to the NE of the hypocenter has a fracture energy of 1.1 MJ/m2 and a slip-weakening distance, Dc, of 0.66 m. Right triangles, whose base and height are Dc and the dynamic stress drop, respectively, approximately overlie the slip-dependent stress given by Ide and Takeo (1997) for the same locations near the hypocenter. The total fracture energy for the Kobe earthquake, 3.7 ?? 1014 J, is about the same as the seismic energy (Ea = 3.2 ?? 1014 J.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

A frictional population model of seismicity rate change

We study models of seismicity rate changes caused by the application of a static stress perturbation to a population of faults and discuss our results with respect to the model proposed by Dieterich (1994). These models assume distribution of nucleation sites (e.g., faults) obeying rate-state frictional relations that fail at constant rate under tectonic loading alone, and predicts a positive static stress step at time to will cause an immediate increased seismicity rate that decays according to Omori's law. We show one way in which the Dieterich model may be constructed from simple general idead, illustratted using numerically computed synthetic seismicity and mathematical formulation. We show that seismicity rate change predicted by these models (1) depend on the particular relationship between the clock-advanced failure and fault maturity, (2) are largest for the faults closest to failure at to, (3) depend strongly on which state evolution law faults obey, and (4) are insensitive to some types of population hetrogeneity. We also find that if individual faults fail repeatedly and populations are finite, at timescales much longer than typical aftershock durations, quiescence follows at seismicity rate increase regardless of the specific frictional relations. For the examined models the quiescence duration is comparable to the ratio of stress change to stressing rate ????/??,which occurs after a time comparable to the average recurrence interval of the individual faults in the population and repeats in the absence of any new load may pertubations; this simple model may partly explain observations of repeated clustering of earthquakes. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Vertical motions of the Puerto Rico Trench and Puerto Rico and their cause

The Puerto Rico trench exhibits great water depth, an extremely low gravity anomaly, and a tilted carbonate platform between (reconstructed) elevations of +1300 m and -4000 m. I argue that these features are manifestations of large vertical movements of a segment of the Puerto Rico trench, its forearc, and the island of Puerto Rico that took place 3.3 m.y. ago over a time period as short as 14-40 kyr. I explain these vertical movements by a sudden increase in the slab's descent angle that caused the trench to subside and the island to rise. The increased dip could have been caused by shearing or even by a complete tear of the descending North American slab, although the exact nature of this deformation is unknown. The rapid (14-40 kyr) and uniform tilt along a 250 km long section of the trench is compatible with scales of mantle flow and plate bending. The proposed shear zone or tear is inferred from seismic, morphological, and gravity observations to start at the trench at 64.5??W and trend southwestwardly toward eastern Puerto Rico. The tensile stresses necessary to deform or tear the slab could have been generated by increased curvature of the trench following a counterclockwise rotation of the upper plate and by the subduction of a large seamount.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

An efficient algorithm for double-difference tomography and location in heterogeneous media, with an application to the Kilauea volcano

Improving our understanding of crustal processes requires a better knowledge of the geometry and the position of geological bodies. In this study we have designed a method based upon double-difference relocation and tomography to image, as accurately as possible, a heterogeneous medium containing seismogenic objects. Our approach consisted not only of incorporating double difference in tomography but also partly in revisiting tomographic schemes for choosing accurate and stable numerical strategies, adapted to the use of cross-spectral time delays. We used a finite difference solution to the eikonal equation for travel time computation and a Tarantola-Valette approach for both the classical and double-difference three-dimensional tomographic inversion to find accurate earthquake locations and seismic velocity estimates. We estimated efficiently the square root of the inverse model's covariance matrix in the case of a Gaussian correlation function. It allows the use of correlation length and a priori model variance criteria to determine the optimal solution. Double-difference relocation of similar earthquakes is performed in the optimal velocity model, making absolute and relative locations less biased by the velocity model. Double-difference tomography is achieved by using high-accuracy time delay measurements. These algorithms have been applied to earthquake data recorded in the vicinity of Kilauea and Mauna Loa volcanoes for imaging the volcanic structures. Stable and detailed velocity models are obtained: the regional tomography unambiguously highlights the structure of the island of Hawaii and the double-difference tomography shows a detailed image of the southern Kilauea caldera-upper east rift zone magmatic complex. Copyright 2005 by the American Geophysical Union.

Hawaii↗

USGS2005

Slicing up the San Francisco Bay Area: Block kinematics and fault slip rates from GPS-derived surface velocities

Observations of surface deformation allow us to determine the kinematics of faults in the San Francisco Bay Area. We present the Bay Area velocity unification (BA??VU??, "bay view"), a compilation of over 200 horizontal surface velocities computed from campaign-style and continuous Global Positioning System (GPS) observations from 1993 to 2003. We interpret this interseismic velocity field using a three-dimensional block model to determine the relative contributions of block motion, elastic strain accumulation, and shallow aseismic creep. The total relative motion between the Pacific plate and the rigid Sierra Nevada/Great Valley (SNGV) microplate is 37.9 ?? 0.6 mm yr-1 directed toward N30.4??W ?? 0.8?? at San Francisco (??2??). Fault slip rates from our preferred model are typically within the error bounds of geologic estimates but provide a better fit to geodetic data (notable right-lateral slip rates in mm yr-1: San Gregorio fault, 2.4 ?? 1.0; West Napa fault, 4.0 ?? 3.0; zone of faulting along the eastern margin of the Coast Range, 5.4 ?? 1.0; and Mount Diablo thrust, 3.9 ?? 1.0 of reverse slip and 4.0 ?? 0.2 of right-lateral strike slip). Slip on the northern Calaveras is partitioned between both the West Napa and Concord/ Green Valley fault systems. The total convergence across the Bay Area is negligible. Poles of rotation for Bay Area blocks progress systematically from the North America-Pacific to North America-SNGV poles. The resulting present-day relative motion cannot explain the strike of most Bay Area faults, but fault strike does loosely correlate with inferred plate motions at the time each fault initiated. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Wave spectral energy variability in the northeast Pacific

The dominant characteristics of wave energy variability in the eastern North Pacific are described from NOAA National Data Buoy Center (NDBC) buoy data collected from 1981 to 2003. Ten buoys at distributed locations were selected for comparison based on record duration and data continuity. Long‐period (LP) [ T > 12] s, intermediate‐period [6 ≤ T ≤ 12] s, and short‐period [ T < 6] s wave spectral energy components are considered separately. Empirical orthogonal function (EOF) analyses of monthly wave energy anomalies reveal that all three wave energy components exhibit similar patterns of spatial variability. The dominant mode represents coherent heightened (or diminished) wave energy along the West Coast from Alaska to southern California, as indicated by composites of the 700 hPa height field. The second EOF mode reveals a distinct El Niño‐Southern Oscillation (ENSO)‐associated spatial distribution of wave energy, which occurs when the North Pacific storm track is extended unusually far south or has receded to the north. Monthly means and principal components (PCs) of wave energy levels indicate that the 1997–1998 El Niño winter had the highest basin‐wide wave energy within this record, substantially higher than the 1982–1983 El Niño. An increasing trend in the dominant PC of LP wave energy suggests that storminess has increased in the northeast Pacific since 1980. This trend is emphasized at central eastern North Pacific locations. Patterns of storminess variability are consistent with increasing activity in the central North Pacific as well as the tendency for more extreme waves in the south during El Niño episodes and in the north during La Niña.

Journal of Geophysical Research C: Oceans↗

USGS2005

Numerical modeling of an estuary: A comprehensive skill assessment

Numerical simulations of the Hudson River estuary using a terrain-following, three-dimensional model (Regional Ocean Modeling System (ROMS)) are compared with an extensive set of time series and spatially resolved measurements over a 43 day period with large variations in tidal forcing and river discharge. The model is particularly effective at reproducing the observed temporal variations in both the salinity and current structure, including tidal, spring neap, and river discharge-induced variability. Large observed variations in stratification between neap and spring tides are captured qualitatively and quantitatively by the model. The observed structure and variations of the longitudinal salinity gradient are also well reproduced. The most notable discrepancy between the model and the data is in the vertical salinity structure. While the surface-to-bottom salinity difference is well reproduced, the stratification in the model tends to extend all the way to the water surface, whereas the observations indicate a distinct pycnocline and a surface mixed layer. Because the southern boundary coindition is located near the mouth the estuary, the salinity within the domain is particularly sensitive to the specification of salinity at the boundary. A boundary condition for the horizontal salinity gradient, based on the local value of salinity, is developed to incorporate physical processes beyond the open boundary not resolved by the model. Model results are sensitive to the specification of the bottom roughness length and vertical stability functions, insofar as they influence the intensity of vertical mixing. The results only varied slightly between different turbulence closure methods of k-??, k-??, and k-kl. Copyright 2005 by the American Geophysical Union.

New York↗