Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Polar Research”

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 307 records · Page 17Linked to original sources

A revised surface age for the North Polar Layered Deposits of Mars

The North Polar Layered Deposits (NPLD) of Mars contain a complex stratigraphy that has been suggested to retain a record of past eccentricity- and obliquity-forced climate changes. The surface accumulation rate in the current climate can be constrained by the crater retention age. We scale NPLD crater diameters to account for icy target strength and compare surface age using a new production function for recent small impacts on Mars to the previously used model of Hartmann (2005). Our results indicate that ice is accumulating in these craters several times faster than previously thought, with a 100 m diameter crater being completely infilled within centuries. Craters appear to have a diameter-dependent lifetime, but the data also permit a complete resurfacing of the NPLD at ~1.5 ka.

Geophysical Research Letters↗

Improving gas-derived parameterization of groundwater using free phase gas measurements

Dissolved atmogenic gasses in groundwater provide significant information about recharge conditions, flowpath, and age. Free phase gas in aquifers is largely ignored in these analyses and there is a lack of quantitative analysis for gas flux mechanisms. Many related fields encountering multiphase flow acknowledge that the presence of bubbles allows for the rapid exsolution of dissolved gasses and volatile compounds through diffusive and polar forces. By measuring the mass flow of the exsolved gas at a spring, coupled with compositional analysis in the free and dissolved phases, we show that not incorporating the effects of the free gas phase of bubbling springs introduces error in the estimation of total gas quantities, particularly light noble gasses. This can significantly affect the corresponding estimation of noble gas temperature (NGT) and apparent age. We examine the transport of free and dissolved gas from the recharge zone, using water level variation data, to the discharge location where the gases are measured. This technique of using the free gas phase for assessing aquifer dynamics will improve groundwater conceptual models, particularly in karstic aquifers where rapid fluctuations in the water table facilitate the development of excess air, generating multiphase spring discharge.

Environmental Science: Water Research & Technology↗

Paleomagnetic record determined in cores from deep research wells in the Quaternary Santa Clara basin, California

Paleomagnetic study of cores from six deep wells provides an independent temporal framework for much of the alluvial stratigraphy of the Quaternary basin beneath the Santa Clara Valley. This stratigraphy consists of 8 upward-fining cycles in the upper 300 m of section and an underlying 150 m or more of largely fine-grained sediment. The eight cycles have been correlated with the marine oxygen isotope record, thus providing one means of dating the section. The section has also proved to contain a rich paleomagnetic record despite the intermittent sedimentation characteristic of alluvial environments. Each well was designed to reach a depth of ~300 m, although 2 were terminated at shallower depth where bedrock was encountered and one (GUAD) was deepened to bedrock at 407.2 m. Cores were taken at intermittent intervals in most of the wells, composing ~20%–25% of their depths. In GUAD an attempt was made to core the entire upper 300 m, with core recovery of 201.8 m (67%). The paleomagnetic framework ranges from the 32 ka Mono Lake excursion near the top of the second sedimentary cycle to below the 780 ka Brunhes-Matuyama geomagnetic reversal beneath the eighth cycle. These ages nicely fit those assigned to the section based on correlation with the marine oxygen isotope record. Several episodes of anomalous magnetic inclinations were also found within the cyclic section in some of the wells. Some of the episodes of anomalous magnetic inclinations are only separated by short normal intervals in a pattern similar to that described for some welldocumented excursions. We consider that a geomagnetic excursion was likely only if the anomalous inclinations were found at approximately the same stratigraphic position in more than one drill hole. A deeper time constraint is provided by the upper boundary (990 ka) of the Jaramillo Normal Polarity Subchron recognized at a depth of 302 m in one deeply penetrating well (GUAD). Approximately 100 m of normal Jaramillo section is evident below that in wells GUAD and EVGR. The reversal that we identify as the 780 ka Brunhes-Matuyama boundary, found at depths of 291–303 m in three wells, indicates an average rate of deposition in this upper section of ~37 cm/k.y. In GUAD, the top of the underlying normally polarized section, which we assign to the upper part of the Jaramillo Normal Polarity Subchron, was found between 301.8 and 304.5 m. The resultant 10 m of reversed polarity section above the Jaramillo seems anomalously short for this 210 k.y. part of the Matuyama Chron, during which several times that thickness of section probably should have accumulated. This observation indicates that a significant unconformity should be present in that short section between the Jaramillo Subchron and the Brunhes-Matuyama boundary. Deeper cores in two wells (GUAD and EVGR) all have normal polarity and seem to represent much of the Jaramillo Subchron, although no base for that subchron was found. The resultant minimum rate of sedimentation for this lower section beneath the unconformity is 170 cm/k.y. The Mono Lake (ca. 32 ka), Pringle Falls (ca. 210 ka), and Big Lost (ca. 565 ka) geomagnetic excursions all seem to be represented in the Santa Clara Valley wells. Possible correlations to the Laschamp (ca. 40 ka) and Blake (ca. 110 ka) excursions are also noted. Three additional excursions that have apparently not been previously reported from western North America occur within cycle 6 (between 536 and 433 ka), near the base of cycle 5 (after 433 ka), and near the middle of cycle 2 (before ca. 75 ka).

California↗

Two million years of climate-driven cave-sediment aggradation and valley incision in the southern Ozark Plateau from Fitton Cave, northern Arkansas, USA

Landscape evolution in karst terrains affects both subterranean and surface settings. For better understanding of controlling processes and connections between the two, multiple geochronometers were used to date sediments and speleothems in upper-level passages of Fitton Cave adjacent to the Buffalo River, northern Arkansas, within the southern Ozark Plateau. Burial cosmogenic-nuclide dating of coarse sediments indicates that gravel pulses washed into upper passages at 2.2 Ma and 1.25 Ma. These represent the oldest epigenetic cave deposits documented in this region. Associated sands and clay-rich sediments mostly have reversed magnetic polarity and thermally transferred optically stimulated luminescence dates of 1.2 to 1.0 Ma. Abandonment of these upper passages began before 0.72 Ma, when coarse sediment was deposited in a passage incised below older sediment. Maximum U-series dates of 0.7–0.4 Ma for flowstones capping clastic deposits mark the stabilization of older sediments and a change to vadose conditions that allowed post–0.4 Ma stalagmite growth. Resulting valley incision rates since 0.85 Ma are estimated at 27 m/Ma. Coarse cave-sediment pulses correlate to Laurentide glacial tills about 300 km to the north, suggesting climate influence on periglacial sediment production. Dated cave sediments also may correlate with undated older strath terraces preserved at similar heights above the Buffalo River.

Arkansas↗

The Kiowa core, a continuous drill core through the Denver Basin bedrock aquifers at Kiowa, Elbert County, Colorado

The Kiowa core was obtained as a component of the Denver Basin Project, a cooperative research effort to study the evolution of the Denver Basin, Colorado. The Kiowa core provides a virtually continuous stratigraphic record of the Upper Cretaceous and lower Tertiary strata of the Denver Basin. The upper portion of the core recovered strata conventionally referred to as the Arapahoe and Denver Formations and the Dawson Arkose. A prominent unconformity marked by a mature paleosol breaks these strata into two unconformity-bounded sequences; the lower sequence is termed Dl and the upper sequence, D2. Beneath these units and also penetrated by the core occur the Laramie Formation, Fox Hills Sandstone, and Pierre Shale. The site for coring was selected in order to obtain fine-grained strata suitable for both palynological and paleomagnetic analyses. The coring effort recovered 93 percent of the 2,256 ft of rock penetrated, resulting in a nearly continuous record of the sedimentary rocks recording the retreat of the Cretaceous Interior Seaway and the subsequent uplift of the Front Range portion of the Rocky Mountains. Palynological data constrain the Cretaceous-Tertiary boundary to a depth between 878 and 880 ft in the core. The palynological data also serve to bracket the age of the paleosol marking the unconformity between the Dl and D2 sequences to between middle Paleocene and earliest Eocene. The paleomagnetic data are interpreted to represent polarity intervals ranging from polarity subchrons 31r to 28n and polarity subchron 24r. Hydrologic analyses indicate variable aquifer characteristics across the State-defined bedrock aquifers. Individual aquifer units exhibit generally lower water-yield potential than was identified to the west in a core drilled by the U.S. Geological Survey (USGS) in 1987 at Castle Pines, Colorado. Downhole temperature measurements indicate a normal geothermal gradient of 30°C/ km. Perturbations of the gradient may represent active fluid flow through the aquifers penetrated by the core. Petrographic examination of the cored sandstone and mudstone units document both the clay-rich character of the paleosol series marking the boundary between the Dl and D2 sequences, and variation in sandstone composition with depth. The lower sequence (Dl) is characterized by litharenites with a significant volcaniclastic component, while the upper sequence (D2) is more arkosic. Extensive lignite beds occur in Dl in the cored interval and these appear as strong reflectors on the seismic line that passes near the core hole. A set of electric logs, core descriptions, and derived data sets accompany this report.

Colorado↗

Pollen in polar ice implies eastern Canadian forest dynamics diverged from climate after European settlement

Rapid warming and human exploitation threaten boreal forests. Understanding links among vegetation, climate, and people in this vast biome requires highly resolved long-term records that integrate regional inputs. We developed an 850-year pollen-based record of supraregional vegetation change using a southern Greenland ice core and atmospheric modeling that identified the boreal and mixed-conifer forests of eastern Canada as the dominant pollen source regions. Conifer pollen increased ∼1400 CE at the onset of the cooler and drier Little Ice Age. A subsequent decline began ∼1650 CE and a statistically significant pollen change after 1760 CE suggests ecological consequences of the Little Ice Age cooling and initial human exploitation that persisted until recent decades. These supraregional changes are broadly consistent with local records and demonstrate intensification of human impacts on northern forests, suggesting a shift from a climate-modulated to an increasingly human-controlled system during recent centuries.

Geophysical Research Letters↗

Stratigraphy of the north polar layered deposits of Mars from high-resolution topography

The stratigraphy of the layered deposits of the polar regions of Mars is theorized to contain a record of recent climate change linked to insolation changes driven by variations in the planet's orbital and rotational parameters. In order to confidently link stratigraphic signals to insolation periodicities, a description of the stratigraphy is required based on quantities that directly relate to intrinsic properties of the layers. We use stereo Digital Terrain Models (DTMs) from the High Resolution Imaging Science Experiment (HiRISE) to derive a characteristic of North Polar Layered Deposits (NPLD) strata that can be correlated over large distances: the topographic protrusion of layers exposed in troughs, which is a proxy for the layers’ resistance to erosion. Using a combination of image analysis and a signal-matching algorithm to correlate continuous depth-protrusion signals taken from DTMs at different locations, we construct a stratigraphic column that describes the upper ~500 m of at least 7% of the area of the NPLD, and find accumulation rates that vary by factors of up to two. We find that, when coupled with observations of exposed layers in orbital images, the topographic expression of the strata is consistently continuous through large distances in the top 300 – 500 m of the NPLD, suggesting it is better related to intrinsic layer properties than brightness alone.

Journal of Geophysical Research E: Planets↗

Non-double-couple mechanisms of microearthquakes induced by hydraulic fracturing

We have inverted polarity and amplitude information of representative microearthquakes to investigate source mechanisms of seismicity induced by hydraulic fracturing in the Carthage Cotton Valley, east Texas, gas field. With vertical arrays of four and eight three-component geophones in two monitoring wells, respectively, we were able to reliably determine source mechanisms of the strongest events with the best signal-to-noise ratio. Our analysis indicates predominantly non-double-couple source mechanisms with positive volumetric component consistent with opening cracks oriented close to expected hydraulic fracture orientation. Our observations suggest the induced events are directly the result of opening cracks by fluid injection, in contrast to many previous studies where the seismicity is interpreted to be primarily shearing caused by pore pressure diffusion into the surrounding rock or associated with shear stresses created at the hydraulic fracture tip.

Journal of Geophysical Research B: Solid Earth↗

Mapping invasive Phragmites australis in the coastal Great Lakes with ALOS PALSAR satellite imagery for decision support

The invasive variety of Phragmites australis (common reed) forms dense stands that can cause negative impacts on coastal Great Lakes wetlands including habitat degradation and reduced biological diversity. Early treatment is key to controlling Phragmites, therefore a map of the current distribution is needed. ALOS PALSAR imagery was used to produce the first basin-wide distribution map showing the extent of large, dense invasive Phragmites-dominated habitats in wetlands and other coastal ecosystems along the U.S. shore of the Great Lakes. PALSAR is a satellite imaging radar sensor that is sensitive to differences in plant biomass and inundation patterns, allowing for the detection and delineation of these tall (up to 5 m), high density, high biomass invasive Phragmites stands. Classification was based on multi-season ALOS PALSAR L-band (23 cm wavelength) HH and HV polarization data. Seasonal (spring, summer, and fall) datasets were used to improve discrimination of Phragmites by taking advantage of phenological changes in vegetation and inundation patterns over the seasons. Extensive field collections of training and randomly selected validation data were conducted in 2010–2011 to aid in mapping and for accuracy assessments. Overall basin-wide map accuracy was 87%, with 86% producer's accuracy and 43% user's accuracy for invasive Phragmites. The invasive Phragmites maps are being used to identify major environmental drivers of this invader's distribution, to assess areas vulnerable to new invasion, and to provide information to regional stakeholders through a decision support tool.

Lake Erie;Lake Huron;Lake Michigan;Lake Ontario;La↗

Wavefield properties of a shallow long-period event and tremor at Kilauea Volcano, Hawaii

The wavefields of tremor and a long-period (LP) event associated with the ongoing eruptive activity at Kilauea Volcano, Hawaii, are investigated using a combination of dense small-aperture (300 m) and sparse large-aperture (5 km) arrays deployed in the vicinity of the summit caldera. Measurements of azimuth and slowness for tremor recorded on the small-aperture array indicate a bimodal nature of the observed wavefield. At frequencies below 2 Hz, the wavefield is dominated by body waves impinging the array with steep incidence. These arrivals are attributed to the oceanic microseismic noise. In the 2-6 Hz band, the wavefield is dominated by waves propagating from sources located at shallow depths (<1 km) beneath the eastern edge of the Halemaumau pit crater. The hypocenter of the LP event, determined from frequency-slowness analyses combined with phase picks, appears to be located close to the source of tremor but at a shallower depth (<0.1 km). The wavefields of tremor and LP event are characterized by a complex composition of body and surface waves, whose propagation and polarization properties are strongly affected by topographic and structural features in the summit caldera region. Analyses of the directional properties of the wavefield in the 2-6 Hz band point to the directions of main scattering sources, which are consistent with pronounced velocity contrasts imaged in a high-resolution three-dimensional velocity model of the caldera region. The frequency and Q of the dominant peak observed in the spectra of the LP event may be explained as the dominant oscillation mode of a crack with scale length 20-100 m and aperture of a few centimeters filled with bubbly water. The mechanism driving the shallow tremor appears to be consistent with a sustained excitation originating in the oscillations of a bubbly cloud resulting from vesiculation and degassing in the magma. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

The widespread influence of Great Lakes microseisms across the United States revealed by the 2014 polar vortex

During the winter of 2014, a weak polar vortex brought record cold temperatures to the north‐central (“Midwest”) United States, and the Great Lakes reached the highest extent of ice coverage (92.5%) since 1979. This event shut down the generation of seismic signals caused by wind‐driven wave action within the lakes (termed “lake microseisms”), giving an unprecedented opportunity to isolate and characterize these novel signals through comparison with nonfrozen time periods. Using seismic records at 72 broadband stations, we observe Great Lakes microseism signals at distances >300 km from the lakes. In contrast to conventional oceanic microseisms, there is no clear relationship between the frequency content of the seismic signals (observed from ~0.5–5‐s period) and the dominant swell period or resonance periods of the lakes based on their bathymetric profiles. Thus, the exact generation mechanism is not readily explained by conventional microseism theory and warrants further investigation.

Geophysical Research Letters↗

Seismic velocity anisotropy and heterogeneity beneath the Mantle Electromagnetic and Tomography Experiment (MELT) region of the East Pacific Rise from analysis of P and S body waves

We use teleseismic P and S delay times and shear wave splitting measurements to constrain isotropic and anisotropic heterogeneity in the mantle beneath the southern East Pacific Rise (SEPR). The data comprise 462 P and S delay times and 18 shear wave splitting observations recorded during the Mantle Electromagnetic and Tomography (MELT) Experiment. We estimate the mantle melt content (F) and temperature (T) variation from the isotropic velocity variation. Our results indicate that the maximum variation in F beneath our array is between zero and ???1.2%, and maximum variation in T is between zero and ???100 K. We favor an explanation having partial contributions from both T and F. We approximate the seismic anisotropy of the upper mantle with hexagonal symmetry, consistent with the assumption of two dimensionality of mantle flow. Our new tomographic technique uses a nonlinear inversion of P and slow S polarization delay times to simultaneously solve for coupled VP and VS heterogeneity throughout the model and for the magnitude of anisotropy within discrete domains. The domain dimensions and the dip of the anisotropy are fixed for each inversion but are varied in a grid search, obtaining the misfit of the models to the body wave delay data and to split times of vertically propagating S waves. The data misfit and the isotropic heterogeneity are sensitive to domain dimensions and dip of anisotropy. In a region centered beneath the SEPR the best average dip of the hexagonal symmetry axis is horizontal or dipping shallowly (<30??) west. Given the resolution of our data, a subaxial region characterized by vertically aligned symmetry axes may exist but is limited to be <80 km deep. We infer that the mantle flow beneath the SEPR is consistent with shallow asthenospheric return flow from the direction of the South Pacific superswell.

Journal of Geophysical Research B: Solid Earth↗

Characterizing physical properties of streambed interface sediments using in situ complex electrical conductivity measurements

Streambed sediment physical properties such as surface area, are difficult to quantify in situ but exert a high‐level control on a wide range of biogeochemical processes and sorption of contaminants. We introduce the use of complex electrical conductivity (CC) methods (also known as spectral induced polarization (SIP)) that measure both real and imaginary conductivity to non‐invasively and efficiently characterize shallow streambed sediments. We explore the method through synthetic modeling, laboratory, and field measurements to demonstrate the sensitivity of imaginary conductivity to sediment surface area, controlled in part by fine‐grained iron oxides produced by anoxic groundwater discharge. Laboratory measurements verify expected relationships between CC parameters and sediment properties. Synthetic modeling using a 1D analytical model illustrates the influence of water layer depth and conductivity on the field CC measurements made at the streambed‐stream water interface. Specifically, the inverted sediment imaginary conductivity is less impacted by uncertainty in the water layer depth and conductivity relative to the real conductivity and phase shift. Field CC measurements along a landfill‐impacted river reveal discrete streambed zones with enhanced bulk surface area generally corresponding to anoxic groundwater discharges zones with high concentrations of fine‐grained iron oxide precipitates.

Water Resources Research↗

The grass is not always greener on the other side: Seasonal reversal of vegetation greenness in aspect-driven semiarid ecosystems

Our current understanding of semiarid ecosystems is that they tend to display higher vegetation greenness on polar-facing slopes (PFS) than on equatorial-facing slopes (EFS). However, recent studies have argued that higher vegetation greenness can occur on EFS during part of the year. To assess whether this seasonal reversal of aspect-driven vegetation is a common occurrence, we conducted a global scale analysis of vegetation greenness on a monthly time scale over an 18-year period (2000-2017). We examined the influence of climate seasonality on the normalised difference vegetation index (NDVI) values of PFS and EFS at 60 different catchments with aspect-controlled vegetation located across all continents except Antarctica. Our results show that an overwhelming majority of sites (70%) display seasonal reversal, associated with transitions from water-limited to energy-limited conditions during wet winters. These findings highlight the need to consider seasonal variations of aspect-driven vegetation patterns in ecohydrology, geomorphology, and earth system models.

Geophysical Research Letters↗

Late Quaternary variations in relative sea level due to glacial cycle polar wander

Growth and decay of continental ice sheets can excite significant motion of the Earth's rotation pole and cause a complex spatio-temporal pattern of changes in relative sea level. These two effects have generally been considered separately, but may interact in important ways. In particular, a simple model of the melting of the Laurentide ice sheet causes a uniform eustatic sea level rise of 55 m, and also induces a motion of the rotation pole by 0.1 to 1 degree, depending on viscosity structure in the mantle. This motion produces a secular pole tide, which is a spherical harmonic degree 2, order 1 component of the relative sea level pattern, with peak-to-peak amplitude of 20 to 40 m. The maximum effect is along the great circle passing through the path of the pole and at latitudes of ±45°. This secular pole tide has been ignored in most previous attempts to estimate ice sheet loading history and mantle viscosity from global patterns of relative sea level change. It has a large influence along the East coast of North America and the West coast of South America, and significantly contributes to present day rates of relative sea level change.

Geophysical Research Letters↗

Conceptual model for the removal of cold-trapped H2O ice on the Mars northern seasonal springtime polar cap

The transport of H 2 O ice along the retreating north polar seasonal CO 2 ice cap has previously been modeled and observed. Spectral observations show that H 2 O ice forms on the interior of the seasonal cap, while thermal observations show these regions to be consistent with CO 2 ice. Prior to the sublimation of the seasonal CO 2 , the observed H 2 O ice deposits are diminished—and because H 2 O ice sublimation rates are extremely slow while in direct thermal contact with CO 2 ice, an alternate removal process must be operating. We propose a model where the process of removing these H 2 O deposits starts with insolation‐induced basal sublimation of the underlying CO 2 ice. This sublimed gas would “seep” upward and into the interface between the two ices, increasing pressure until the gas pressure fractures the cold‐trapped H 2 O ice. Small fragments would be suspended while larger fragments would be pushed aside, exposing the underlying CO 2 ice.

Geophysical Research Letters↗

Mapping of the Moon by Clementine

The "faster, cheaper, better" Clementine spacecraft mission mapped the Moon from February 19 to May 3, 1994. Global coverage was acquired in 11 spectral bandpasses from 415 to 2792 nm and at resolutions of 80-330 m/pixel; a thermal-infrared camera sampled ???20% of the surface; a high-resolution camera sampled selected areas (especially the polar regions); and a lidar altimeter mapped the large-scale topography up to latitudes of ??75??. The spacecraft was in a polar, elliptical orbit, 400-450 km periselene altitude. Periselene latitude was -28.5?? for the first month of mapping, then moved to +28.5??. NASA is supporting the archiving, systematic processing, and analysis of the ???1.8 million lunar images and other datasets. A new global positional network has been constructed from 43,000 images and ???0.5 million match points; new digital maps will facilitate future lunar exploration. In-flight calibrations now enable photometry to a high level of precision for the uv-visible CCD camera. Early science results include: (1) global models of topography, gravity, and crustal thicknesses; (2) new information on the topography and structure of multiring impact basins; (3) evidence suggestive of water ice in large permanent shadows near the south pole; (4) global mapping of iron abundances; and (5) new constraints on the Phanerozoic cratering rate of the Earth. Many additional results are expected following completion of calibration and systematic processing efforts. ?? 1997 COSPAR. Published by Elsevier Science Ltd.

Advances in Space Research↗

Microseismicity at the North Anatolian Fault in the Sea of Marmara offshore Istanbul, NW Turkey

The North Anatolian Fault Zone (NAFZ) below the Sea of Marmara forms a “seismic gap” where a major earthquake is expected to occur in the near future. This segment of the fault lies between the 1912 Ganos and 1999 İzmit ruptures and is the only NAFZ segment that has not ruptured since 1766. To monitor the microseismic activity at the main fault branch offshore of Istanbul below the Çınarcık Basin, a permanent seismic array (PIRES) was installed on the two outermost Prince Islands, Yassiada and Sivriada, at a few kilometers distance to the fault. In addition, a temporary network of ocean bottom seismometers was deployed throughout the Çınarcık Basin. Slowness vectors are determined combining waveform cross correlation and P wave polarization. We jointly invert azimuth and traveltime observations for hypocenter determination and apply a bootstrap resampling technique to quantify the location precision. We observe seismicity rates of 20 events per month for M < 2.5 along the basin. The spatial distribution of hypocenters suggests that the two major fault branches bounding the depocenter below the Çınarcık Basin merge to one single master fault below ∼17 km depth. On the basis of a cross-correlation technique we group closely spaced earthquakes and determine composite focal mechanisms implementing recordings of surrounding permanent land stations. Fault plane solutions have a predominant right-lateral strike-slip mechanism, indicating that normal faulting along this part of the NAFZ plays a minor role. Toward the west we observe increasing components of thrust faulting. This supports the model of NW trending, dextral strike-slip motion along the northern and main branch of the NAFZ below the eastern Sea of Marmara.

Istanbul↗