USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Journal”

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

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

At least 1,243 records · Page 69Linked to original sources
USGS2011

Modern thermokarst lake dynamics in the continuous permafrost zone, northern Seward Peninsula, Alaska

Quantifying changes in thermokarst lake extent is of importance for understanding the permafrost-related carbon budget, including the potential release of carbon via lake expansion or sequestration as peat in drained lake basins. We used high spatial resolution remotely sensed imagery from 1950/51, 1978, and 2006/07 to quantify changes in thermokarst lakes for a 700 km 2 area on the northern Seward Peninsula, Alaska. The number of water bodies larger than 0.1 ha increased over the entire observation period (666 to 737 or +10.7%); however, total surface area decreased (5,066 ha to 4,312 ha or -14.9%). This pattern can largely be explained by the formation of remnant ponds following partial drainage of larger water bodies. Thus, analysis of large lakes (>40 ha) shows a decrease of 24% and 26% in number and area, respectively, differing from lake changes reported from other continuous permafrost regions. Thermokarst lake expansion rates did not change substantially between 1950/51 and 1978 (0.35 m/yr) and 1978 and 2006/07 (0.39 m/yr). However, most lakes that drained did expand as a result of surface permafrost degradation before lateral drainage. Drainage rates over the observation period were stable (2.2 to 2.3 per year). Thus, analysis of decadal-scale, high spatial resolution imagery has shown that lake drainage in this region is triggered by lateral breaching and not subterranean infiltration. Future research should be directed toward better understanding thermokarst lake dynamics at high spatial and temporal resolution as these systems have implications for landscape-scale hydrology and carbon budgets in thermokarst lake-rich regions in the circum-Arctic.

Alaska↗
USGS2012

Linkages between denitrification and dissolved organicmatter quality, Boulder Creek watershed, Colorado

Dissolved organic matter (DOM) fuels the majority of in-stream microbial processes, including the removal of nitrate via denitrification. However, little is known about how the chemical composition of DOM influences denitrification rates. Water and sediment samples were collected across an ecosystem gradient, spanning the alpine to plains, in central Colorado to determine whether the chemical composition of DOM was related to denitrification rates. Laboratory bioassays measured denitrification potentials using the acetylene block technique and carbon mineralization via aerobic bioassays, while organic matter characteristics were evaluated using spectroscopic and fractionation methods. Denitrification potentials under ambient and elevated nitrate concentrations were strongly correlated with aerobic respiration rates and the percent mineralized carbon, suggesting that information about the aerobic metabolism of a system can provide valuable insight regarding the ability of the system to additionally reduce nitrate. Multiple linear regressions (MLR) revealed that under elevated nitrate concentrations denitrification potentials were positively related to the presence of protein-like fluorophores and negatively related to more aromatic and oxidized fractions of the DOM pool. Using MLR, the chemical composition of DOM, carbon, and nitrate concentrations explained 70% and 78% of the observed variability in denitrification potential under elevated and ambient nitrate conditions, respectively. Thus, it seems likely that DOM optical properties could help to improve predictions of nitrate removal in the environment. Finally, fluorescence measurements revealed that bacteria used both protein and humic-like organic molecules during denitrification providing further evidence that larger, more aromatic molecules are not necessarily recalcitrant in the environment.

Colorado↗
USGS2012

Peat accumulation in drained thermokarst lake basins in continuous, ice-rich permafrost, northern Seward Peninsula, Alaska

Thermokarst lakes and peat-accumulating drained lake basins cover a substantial portion of Arctic lowland landscapes, yet the role of thermokarst lake drainage and ensuing peat formation in landscape-scale carbon (C) budgets remains understudied. Here we use measurements of terrestrial peat thickness, bulk density, organic matter content, and basal radiocarbon age from permafrost cores, soil pits, and exposures in vegetated, drained lake basins to characterize regional lake drainage chronology, C accumulation rates, and the role of thermokarst-lake cycling in carbon dynamics throughout the Holocene on the northern Seward Peninsula, Alaska. Most detectable lake drainage events occurred within the last 4,000 years with the highest drainage frequency during the medieval climate anomaly. Peat accumulation rates were highest in young (50–500 years) drained lake basins (35.2 g C m −2 yr −1 ) and decreased exponentially with time since drainage to 9 g C m −2 yr −1 in the oldest basins. Spatial analyses of terrestrial peat depth, basal peat radiocarbon ages, basin geomorphology, and satellite-derived land surface properties (Normalized Difference Vegetation Index (NDVI); Minimum Noise Fraction (MNF)) from Landsat satellite data revealed significant relationships between peat thickness and mean basin NDVI or MNF. By upscaling observed relationships, we infer that drained thermokarst lake basins, covering 391 km 2 (76%) of the 515 km 2 study region, store 6.4–6.6 Tg organic C in drained lake basin terrestrial peat. Peat accumulation in drained lake basins likely serves to offset greenhouse gas release from thermokarst-impacted landscapes and should be incorporated in landscape-scale C budgets.

Alaska↗
USGS2012

Persistence of soil organic matter in eroding versus depositional landform positions

Soil organic matter (SOM) processes in dynamic landscapes are strongly influenced by soil erosion and sedimentation. We determined the contribution of physical isolation of organic matter (OM) inside aggregates, chemical interaction of OM with soil minerals, and molecular structure of SOM in controlling storage and persistence of SOM in different types of eroding and depositional landform positions. By combining density fractionation with elemental and spectroscopic analyses, we showed that SOM in depositional settings is less transformed and better preserved than SOM in eroding landform positions. However, which environmental factors exert primary control on storage and persistence of SOM depended on the nature of the landform position considered. In an annual grassland watershed, protection of SOM by physical isolation inside aggregates and chemical association of organic matter (complexation) with soil minerals, as assessed by correlation with radiocarbon concentration, were more effective in the poorly drained, lowest-lying depositional landform positions, compared to well-drained landform positions in the upper parts of the watershed. Results of this study demonstrated that processes of soil erosion and deposition are important mechanisms of long-term OM stabilization.

Journal of Geophysical Research: Biogeosciences↗
USGS2012

Characterizing mercury concentrations and fluxes in a Coastal Plain watershed: Insights from dynamic modeling and data

Mercury (Hg) is one of the leading water quality concerns in surface waters of the United States. Although watershed-scale Hg cycling research has increased in the past two decades, advances in modeling watershed Hg processes in diverse physiographic regions, spatial scales, and land cover types are needed. The goal of this study was to assess Hg cycling in a Coastal Plain system using concentrations and fluxes estimated by multiple watershed-scale models with distinct mathematical frameworks reflecting different system dynamics. We simulated total mercury (Hg T , the sum of filtered and particulate forms) concentrations and fluxes from a Coastal Plain watershed (McTier Creek) using three watershed Hg models and an empirical load model. Model output was compared with observed in-stream Hg T . We found that shallow subsurface flow is a potentially important transport mechanism of particulate Hg T during periods when connectivity between the uplands and surface waters is maximized. Other processes (e.g., stream bank erosion, sediment re-suspension) may increase particulate Hg T in the water column. Simulations and data suggest that variable source area (VSA) flow and lack of rainfall interactions with surface soil horizons result in increased dissolved Hg T concentrations unrelated to DOC mobilization following precipitation events. Although flushing of DOC-Hg T complexes from surface soils can also occur during this period, DOC-complexed Hg T becomes more important during base flow conditions. TOPLOAD simulations highlight saturated subsurface flow as a primary driver of daily Hg T loadings, but shallow subsurface flow is important for Hg T loads during high-flow events. Results suggest limited seasonal trends in Hg T dynamics.

South Carolina↗
USGS2012

Landscape controls on total and methyl Hg in the Upper Hudson River basin, New York, USA

Approaches are needed to better predict spatial variation in riverine Hg concentrations across heterogeneous landscapes that include mountains, wetlands, and open waters. We applied multivariate linear regression to determine the landscape factors and chemical variables that best account for the spatial variation of total Hg (THg) and methyl Hg (MeHg) concentrations in 27 sub-basins across the 493 km 2 upper Hudson River basin in the Adirondack Mountains of New York. THg concentrations varied by sixfold, and those of MeHg by 40-fold in synoptic samples collected at low-to-moderate flow, during spring and summer of 2006 and 2008. Bivariate linear regression relations of THg and MeHg concentrations with either percent wetland area or DOC concentrations were significant but could account for only about 1/3 of the variation in these Hg forms in summer. In contrast, multivariate linear regression relations that included metrics of (1) hydrogeomorphology, (2) riparian/wetland area, and (3) open water, explained about 66% to >90% of spatial variation in each Hg form in spring and summer samples. These metrics reflect the influence of basin morphometry and riparian soils on Hg source and transport, and the role of open water as a Hg sink. Multivariate models based solely on these landscape metrics generally accounted for as much or more of the variation in Hg concentrations than models based on chemical and physical metrics, and show great promise for identifying waters with expected high Hg concentrations in the Adirondack region and similar glaciated riverine ecosystems.

New York↗
USGS2012

Carbon and sediment accumulation in the Everglades (USA) during the past 4000 years: rates, drivers, and sources of error

Tropical and sub-tropical wetlands are considered to be globally important sources for greenhouse gases but their capacity to store carbon is presumably limited by warm soil temperatures and high rates of decomposition. Unfortunately, these assumptions can be difficult to test across long timescales because the chronology, cumulative mass, and completeness of a sedimentary profile are often difficult to establish. We therefore made a detailed analysis of a core from the principal drainage outlet of the Everglades of South Florida, to assess these problems and determine the factors that could govern carbon accumulation in this large sub-tropical wetland. Accelerator mass spectroscopy dating provided direct evidence for both hard-water and open-system sources of dating errors, whereas cumulative mass varied depending upon the type of method used. Radiocarbon dates of gastropod shells, nevertheless, seemed to provide a reliable chronology for this core once the hard-water error was quantified and subtracted. Long-term accumulation rates were then calculated to be 12.1 g m -2 yr -1 for carbon, which is less than half the average rate reported for northern and tropical peatlands. Moreover, accumulation rates remained slow and relatively steady for both organic and inorganic strata, and the slow rate of sediment accretion ( 0.2 mm yr -1 ) tracked the correspondingly slow rise in sea level (0.35 mm yr -1 ) reported for South Florida over the past 4000 years. These results suggest that sea level and the local geologic setting may impose long-term constraints on rates of sediment and carbon accumulation in the Everglades and other wetlands.

Florida↗
USGS2012

Spatiotemporal analysis of black spruce forest soils and implications for the fate of C

Post-fire storage of carbon (C) in organic-soil horizons was measured in one Canadian and three Alaskan chronosequences in black spruce forests, together spanning stand ages of nearly 200 yrs. We used a simple mass balance model to derive estimates of inputs, losses, and accumulation rates of C on timescales of years to centuries. The model performed well for the surface and total organic soil layers and presented questions for resolving the dynamics of deeper organic soils. C accumulation in all study areas is on the order of 20–40 gC/m 2 / yr for stand ages up to ∼200 yrs. Much larger fluxes, both positive and negative, are detected using incremental changes in soil C stocks and by other studies using eddy covariance methods for CO 2 . This difference suggests that over the course of stand replacement, about 80% of all net primary production (NPP) is returned to the atmosphere within a fire cycle, while about 20% of NPP enters the organic soil layers and becomes available for stabilization or loss via decomposition, leaching, or combustion. Shifts toward more frequent and more severe burning and degradation of deep organic horizons would likely result in an acceleration of the carbon cycle, with greater CO 2 emissions from these systems overall.

Journal of Geophysical Research↗
USGS2012

Reoccupation of floodplains by rivers and its relation to the age structure of floodplain vegetation

River channel dynamics over many decades provide a physical control on the age structure of floodplain vegetation as a river occupies and abandons locations. Floodplain reoccupation by a river, in particular, determines the interval of time during which vegetation can establish and mature. A general framework for analyzing floodplain reoccupation and a time series model are developed and applied to five alluvial rivers in the United States. Channel dynamics in these rivers demonstrate time-scale dependence with short-term oscillation in active channel area in response to floods and subsequent vegetation growth and progressive lateral movement that accounts for much of the cumulative area occupied by the rivers over decades. Rivers preferentially reoccupy locations recently abandoned causing a decreasing probability of reoccupation with time since abandonment. For a typical case, a river is 10 times more likely to reoccupy an area it abandoned in the past decade than it is to reoccupy an area it abandoned 30 yrs ago. The decreasing probability of reoccupation over time is consistent with observations of persistent stands of late seral stage floodplain forest. A power function provides a robust approach for estimating the cumulative area occupied by a river and the age structure of riparian forests resulting from a specific historical sequence of streamflow in comparison to either linear or exponential alternatives.

Journal of Geophysical Research↗
USGS2012

Dissolved organic carbon and chromophoric dissolved organic matter properties of rivers in the USA

Dissolved organic carbon (DOC) concentration and chromophoric dissolved organic matter (CDOM) parameters were measured over a range of discharge in 30 U.S. rivers, covering a diverse assortment of fluvial ecosystems in terms of watershed size and landscape drained. Relationships between CDOM absorption at a range of wavelengths (a 254 , a 350 , a 440 ) and DOC in the 30 watersheds were found to correlate strongly and positively for the majority of U.S. rivers. However, four rivers (Colorado, Colombia, Rio Grande and St. Lawrence) exhibited statistically weak relationships between CDOM absorption and DOC. These four rivers are atypical, as they either drain from the Great Lakes or experience significant impoundment of water within their watersheds, and they exhibited values for dissolved organic matter (DOM) parameters indicative of autochthonous or anthropogenic sources or photochemically degraded allochthonous DOM and thus a decoupling between CDOM and DOC. CDOM quality parameters in the 30 rivers were found to be strongly correlated to DOM compositional metrics derived via XAD fractionation, highlighting the potential for examining DOM biochemical quality from CDOM measurements. This study establishes the ability to derive DOC concentration from CDOM absorption for the majority of U.S. rivers, describes characteristics of riverine systems where such an approach is not valid, and emphasizes the possibility of examining DOM composition and thus biogeochemical function via CDOM parameters. Therefore, the usefulness of CDOM measurements, both laboratory-based analyses and in situ instrumentation, for improving spatial and temporal resolution of DOC fluxes and DOM dynamics in future studies is considerable in a range of biogeochemical studies.

Journal of Geophysical Research G: Biogeosciences↗
USGS2012

Uncovering the nonadiabatic response of geosynchronous electrons to geomagnetic disturbance

We describe an energy spectrum method for scaling electron integral flux, which is measured at a constant energy, to phase space density at a constant value of the first adiabatic invariant which removes much of the variation due to reversible adiabatic effects. Applying this method to nearly a solar cycle (1995 - 2006) of geosynchronous electron integral flux (E>2.0MeV) from the GOES satellites, we see that much of the diurnal variation in electron phase space density at constant energy can be removed by the transformation to phase space density at constant μ (4000 MeV/G). This allows us a clearer picture of underlying non-adiabatic electron population changes due to geomagnetic activity. Using scaled phase space density, we calculate the percentage of geomagnetic storms resulting in an increase, decrease or no change in geosynchronous electrons as 38%, 7%, and 55%, respectively. We also show examples of changes in the electron population that may be different than the unscaled fluxes alone suggest. These examples include sudden electron enhancements during storms which appear during the peak of negative Dst for μ-scaled phase space density, contrary to the slow increase seen during the recovery phase for unscaled phase space density for the same event.

Journal of Geophysical Research↗
USGS2012

Problem of the Love‐Gannon relation between the asymmetric disturbance field and Dst

Love and Gannon (2009) discovered that statistically, over a fifty year period the difference in the dawn and dusk disturbance‐field H component at low latitudes (hourly averaged) is linearly proportional to Dst. If the difference is designated by δ DD in units of nT/R E , then the Love‐Gannon (L‐G) relation is δ DD = − 0.2 Dst. At any time departures from the relation can be large. Nonetheless, the relation is evident for all values of Dst and persists throughout magnetic storms, both the main phase and the recovery phase. The Love‐Gannon discovery presents a problem to current understanding of the relation between the causes of δ DD and Dst because the dawn dusk asymmetry in the disturbance field is presumably governed by a long‐established magnetosphere‐ionosphere coupling theory which predicts a characteristic time scale (the shielding time) of less than an hour whereas the characteristic time scale for Dst (the ring current decay time) is more like ten hours. Thus, without forcing both time scales toward each other to the limits of their ranges, a linear proportionality between δ DD and Dst cannot be derived from the current understanding of the causes of the asymmetry and the ring current. This conclusion is the paper's main contribution. In addition, we attempt to get around the conflict of time scales by looking at other possibilities for generating δ DD that depend directly on the ring current. The most promising of these is the possibility that the ring current decay mechanism creates a quasi‐permanent, local‐time modification of the ring current compared to what it would be in the absence of the decay mechanism and that this modification causes a field‐aligned current that closes through the ionosphere and generates the asymmetry δ DD . This idea has the virtue of coupling the asymmetry directly to the ring current and of accounting for the persistence of the L‐G proportionality through the recovery phase of magnetic storms.

Journal of Geophysical Research A: Space Physics↗
USGS2012

Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity

A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.

British Columbia;California;Oregon;Washington↗
USGS2012

Laboratory triggering of stick-slip events by oscillatory loading in the presence of pore fluid with implications for physics of tectonic tremor

The physical mechanism by which the low-frequency earthquakes (LFEs) that make up portions of tectonic (also called non-volcanic) tremor are created is poorly understood. In many areas of the world, tectonic tremor and LFEs appear to be strongly tidally modulated, whereas ordinary earthquakes are not. Anomalous seismic wave speeds, interpreted as high pore fluid pressure, have been observed in regions that generate tremor. Here we build upon previous laboratory studies that investigated the response of stick-slip on artificial faults to oscillatory, tide-like loading. These previous experiments were carried out using room-dry samples of Westerly granite, at one effective stress. Here we augment these results with new experiments on Westerly granite, with the addition of varying effective stress using pore fluid at two pressures. We find that raising pore pressure, thereby lowering effective stress can significantly increase the degree of correlation of stick-slip to oscillatory loading. We also find other pore fluid effects that become important at higher frequencies, when the period of oscillation is comparable to the diffusion time of pore fluid into the fault. These results help constrain the conditions at depth that give rise to tidally modulated LFEs, providing confirmation of the effective pressure law for triggering and insights into why tremor is tidally modulated while earthquakes are at best only weakly modulated.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Clustering of GPS velocities in the Mojave Block, southeastern California

We find subdivisions within the Mojave Block using cluster analysis to identify groupings in the velocities observed at GPS stations there. The clusters are represented on a fault map by symbols located at the positions of the GPS stations, each symbol representing the cluster to which the velocity of that GPS station belongs. Fault systems that separate the clusters are readily identified on such a map. The most significant representation as judged by the gap test involves 4 clusters within the Mojave Block. The fault systems bounding the clusters from east to west are 1) the faults defining the eastern boundary of the Northeast Mojave Domain extended southward to connect to the Hector Mine rupture, 2) the Calico-Paradise fault system, 3) the Landers-Blackwater fault system, and 4) the Helendale-Lockhart fault system. This division of the Mojave Block is very similar to that proposed by Meade and Hager. However, no cluster boundary coincides with the Garlock Fault, the northern boundary of the Mojave Block. Rather, the clusters appear to continue without interruption from the Mojave Block north into the southern Walker Lane Belt, similar to the continuity across the Garlock Fault of the shear zone along the Blackwater-Little Lake fault system observed by Peltzer et al. Mapped traces of individual faults in the Mojave Block terminate within the block and do not continue across the Garlock Fault [Dokka and Travis, ].

California↗
USGS2012

Coupled atmosphere-ocean-wave simulations of a storm event over the Gulf of Lion and Balearic Sea

The coastal areas of the North-Western Mediterranean Sea are one of the most challenging places for ocean forecasting. This region is exposed to severe storms events that are of short duration. During these events, significant air-sea interactions, strong winds and large sea-state can have catastrophic consequences in the coastal areas. To investigate these air-sea interactions and the oceanic response to such events, we implemented the Coupled Ocean-Atmosphere-Wave-Sediment Transport Modeling System simulating a severe storm in the Mediterranean Sea that occurred in May 2010. During this event, wind speed reached up to 25 m.s-1 inducing significant sea surface cooling (up to 2°C) over the Gulf of Lion (GoL) and along the storm track, and generating surface waves with a significant height of 6 m. It is shown that the event, associated with a cyclogenesis between the Balearic Islands and the GoL, is relatively well reproduced by the coupled system. A surface heat budget analysis showed that ocean vertical mixing was a major contributor to the cooling tendency along the storm track and in the GoL where turbulent heat fluxes also played an important role. Sensitivity experiments on the ocean-atmosphere coupling suggested that the coupled system is sensitive to the momentum flux parameterization as well as air-sea and air-wave coupling. Comparisons with available atmospheric and oceanic observations showed that the use of the fully coupled system provides the most skillful simulation, illustrating the benefit of using a fully coupled ocean-atmosphere-wave model for the assessment of these storm events.

Journal of Geophysical Research↗