Forced waves in the martian atmosphere from MGS TES nadir data

We have analyzed the temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) nadir spectra to yield latitude-height resolved maps of various atmospheric forced wave modes as a function of season for a full Mars year. Among the isolated wave modes is the zonal mean, time mean temperature, which we used to derive zonal mean zonal winds and stationary wave quasi-geostrophic indices of refraction, diagnostic of their propagation. The diurnal Kelvin wave was isolated in the data, with results roughly consistent with models (Wilson and Hamilton, 1996, J. Atmos. Sci. 33, 1290-1326). The s = 1 and s = 2 stationary waves were found to have significant amplitude in ducts extending up the winter polar jets, while the s = 3 stationary wave was found to be confined to near the surface. The s = 1 stationary wave was found to have little phase tilt with height during northern winter, but significant westward phase tilt with height in the southern winter. This indicates that the wave carries heat poleward, slightly more than that found in Barnes et al. (1996; J. Geophys. Res. 101, 12,753-12,776). The s = 1 stationary wave is likely the dominant mechanism for eddy meridional heat transport for the southern winter. We noted that the phase of the s = 2 stationary wave is nearly constant with time, but that the s = 1 stationary wave moved 90° of longitude from fall to winter and back in spring in the North. While interannual variability is not yet addressed, overall, these results provide the first comprehensive benchmark for forced waves in Mars’s atmosphere against which future atmospheric models of Mars can be compared.     

Surficial properties in Melas Chasma, Mars, from Mars Odyssey THEMIS data

We examine the nature of the surface layer in a small area of the Melas Chasma region as determined from high-resolution thermal and visible Mars Odyssey Thermal Emission Imaging System (THEMIS) data as well as how our conclusions compare to past analyses. At THEMIS resolution, the thermal structure is dominated by local control and all significant thermal variations can be linked to morphology. Thus, THEMIS provides us with detailed images that contain thermophysical information as well, allowing us to create a surficial geologic map intended to reflect the surface structure of the region. Eight units have been defined: (i) blanketed plateaus with thermally distinct craters and fractures, (ii) blanketed canyon walls with rocky edges, (iii) indurated and/or rocky canyon wall slide material partially covered by aeolian material, (iv) an anomalous wall region with fluvial-like depressions partially filled with particulate material, (v) indurated and/or rocky ridged and non-ridged canyon floor landslide material mingled with aeolian material, (vi) sand sheets, (vii) indurated and/or rocky rounded blocks intermingled with small areas of aeolian material, and (viii) transverse dunes. The THEMIS thermal data support conclusions from previous studies but also reveal much more structure than was seen in the past. We have found that all significant thermal variations in this region can be linked to morphology but all morphological variations cannot be linked to significant thermal variations. THEMIS visible images provide an intermediate resolution that bridges the gap between MOC and Viking and allow for a more meaningful interpretation of the geologic context of a region. Surfaces indicate that landslides were an important geologic process long ago, shaping the canyon walls and floor, while aeolian processes have subsequently altered the surface layer in many locations and may still be active.     

Decoupling the host and nuclear spectra of type I AGNs using integral field spectroscopy: A test on 3C 120

We have developed a new technique to decouple the spectra of the host and the nucleus of type I AGNs using integral field spectroscopy data. The technique is a simple extension of methods widely tested in 2D imaging. We present here the results from applying the technique to data taken with INTEGRAL at the 4.2 m William Herschel Telescope on the Seyfert 1 radio-galaxy 3C 120. We obtained, for the first time, a clean spectrum of the host galaxy, without contamination from the nuclear source.     

Pore-scale characteristics of multiphase flow in porous media: A comparison of air–water and oil–water experiments

Studies of NAPL dissolution in porous media have demonstrated that measurement of saturation alone is insufficient to describe the rate of dissolution. Quantification of the NAPL–water interfacial area provides a measure of the expected area available for mass transfer and will likely be a primary determinant of NAPL removal efficiency. To measure the interfacial area, we have used a synchrotron-based CMT technique to obtain high-resolution 3D images of flow in a Soltrol–water–glass bead system. The interfacial area is found to increase as the wetting phase saturation decreases, reach a maximum, and then decrease as the wetting phase saturation goes to zero. These results are compared to previous findings for an air–water–glass bead study; The Soltrol–water interfacial areas were found to peak at similar saturations as those measured for the air–water system (20–35% saturation range), however, the peak values were in some cases almost twice as high for the oil-water system. We believe that the observed differences between the air–water and oil–water systems to a large degree can be explained by the differences in interfacial tensions for the two systems.     

Quantification and regionalization of groundwater–surface water interaction along an alluvial stream

In the present study, groundwater seepage to an alluvial stream and two tributary streams was examined at nine field sites using hydrological, geophysical, and geomorphological observations. The data indicate that seepage enters the streams in the following ways: (i) directly through the streambed; (ii) as nearly superficial flow from diffuse discharge areas on the flood plains or; (iii) as a combination of (i) and (ii). At about 40% of the sites more than 50% of seepage flows through the streambed. Moreover, it was found that the ratio C, defined as the width of the wet zone of the flood plain divided by the effective width of the stream, can be used as an indicator of the percentage of water entering the stream directly through the streambed. When C is small streambed seepage is large, while when C is large streambed seepage is small and ground water enters the stream mainly as nearly superficial or over-bank flow from the wet zone.     

The rupture process and tectonic implications of the great 1964 Prince William Sound earthquake

We have determined the rupture history of the March 28, 1964, Prince Williams Sound earthquake (M _w=9.2) from long-period WWSSNP-wave seismograms. Source time functions determined from the long-periodP waves indicate two major pulses of moment release. The first and largest moment pulse has a duration of approximately 100 seconds with a relatively smooth onset which reaches a peak moment release rate at about 75 seconds into the rupture. The second smaller pulse of moment release starts at approximately 160 seconds after the origin time and has a duration of roughly 40 seconds. Because of the large size of this event and thus a deficiency of on-scale, digitizableP-wave seismograms, it is impossible to uniquely invert for the location of moment release. However, if we assume a rupture direction based on the aftershock distribution and the results of surface wave directivity studies we are able to locate the spatial distribution of moment along the length of the fault. The first moment pulse most likely initiated near the epicenter at the northeastern down-dip edge of the aftershock area and then spread over the fault surface in a semi-circular fashion until the full width of the fault was activated. The rupture then extended toward the southwest approximately 300 km (Ruff andKanamori, 1983). The second moment pulse was located in the vicinity of Kodiak Island, starting at 500 km southwest of the epicenter and extending to about 600 km. Although the aftershock area extends southwest past the second moment pulse by at least 100 km, the moment release remained low. We interpret the 1964 Prince William Sound earthquake as a multiple asperity rupture with a very large dominant asperity in the epicentral region and a second major, but smaller, asperity in the Kodiak Island region.The zone that ruptured in the 1964 earthquake is segmented into two regions corresponding to the two regions of concentrated moment release. Historical earthquake data suggest that these segments behaved independently during previous events. The Kodiak Island region appears to rupture more frequently with previous events occurring in 1900, 1854, 1844, and 1792. In contrast, the Prince William Sound region has much longer recurrence intervals on the order of 400–1000 years.     

Field measurements of shear stress and friction in the surf zone

Measurements of near-bed shear stress were undertaken in the shallow subtidal zone at Durras Beach, NSW, Australia using a sideways-looking acoustic velocity meter installed within the wave boundary layer. The wave climate was swell-dominated and wave conditions comprised shoaling and breaking waves as well as surf bores. The sediment at the field site was medium-grained sand, and observations of bedform geometry were conducted using a pencilbeam-sonar system. Using frequency-filtering techniques, the measured stresses were partitioned into terms representing turbulent (Reynolds) stress, stresses due to gravity and infragravity-scale oscillatory motions, and wave-turbulence-mean current cross-terms. Gravity wave-orbital scale motions contributed the largest fraction of the stresses, comprising 24% on average, followed by long-wave advection of vertical orbital motion (16%). The presence of wave orbital-scale motions near or at the water/sediment interface was likely due to the porous nature of the seabed, facilitating interfacial flow. Shear stresses did not scale with bed roughness but exhibited a linear relationship with the relative wave height. This indicates that for the experimental conditions, surf zone processes overwhelmed bed roughness effects on shear stress and friction. Calculations of the wave friction factor, f_w, showed that in a natural surf zone, this was a factor 3–4 larger than conventional predictions. © 2020 John Wiley & Sons, Ltd.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.     

Evaluation of MODFLOW-LGR in connection with a synthetic regional-scale model

This work studies costs and benefits of utilizing local-grid refinement (LGR) as implemented in MODFLOW-LGR to simulate groundwater flow in a buried tunnel valley interacting with a regional aquifer. Two alternative LGR methods were used: the shared-node (SN) method and the ghost-node (GN) method. To conserve flows the SN method requires correction of sources and sinks in cells at the refined/coarse-grid interface. We found that the optimal correction method is case dependent and difficult to identify in practice. However, the results showed little difference and suggest that identifying the optimal method was of minor importance in our case. The GN method does not require corrections at the models' interface, and it uses a simpler head interpolation scheme than the SN method. The simpler scheme is faster but less accurate so that more iterations may be necessary. However, the GN method solved our flow problem more efficiently than the SN method. The MODFLOW-LGR results were compared with the results obtained using a globally coarse (GC) grid. The LGR simulations required one to two orders of magnitude longer run times than the GC model. However, the improvements of the numerical resolution around the buried valley substantially increased the accuracy of simulated heads and flows compared with the GC simulation. Accuracy further increased locally around the valley flanks when improving the geological resolution using the refined grid. Finally, comparing MODFLOW-LGR simulation with a globally refined (GR) grid showed that the refinement proportion of the model should not exceed 10% to 15% in order to secure method efficiency.