Соган · предресерч

Геофизика и рельеф

Кратко.

По открытым геофизическим и топографическим данным участок 62 га обеспечивается ровно двумя источниками с достаточным разрешением: Copernicus DEM GLO-30 (756 пикселей на участок; высоты 1597,5–1754,5 м, перепад 157 м, медианный уклон 18,9°, преобладание склонов северной и северо-восточной экспозиции) и архивом Sentinel-1 SLC (934 гранулы, три трека — один восходящий и два нисходящих, 12-суточный интервал, без зон перекрытия и тени на участке при типичном угле падения 39°). Всё остальное — глобальные магнитные (EMAG2v3 с ячейкой ~12,1 км², WDMAM с ячейкой ~27,2 км²) и гравитационные модели (WGM2012 и XGM2019e с шагом 2′, EGM2008 с 2,5′, GOCO06s с пределом 67 км), национальная аэромагнитная съёмка Ирана (галсы 7,5 км, отфильтрованы волны короче 15 км), модели Мохо (сетка 0,5°, СКО 4,9–6,2 км) и каталоги сейсмичности — работает исключительно как региональный контекст. Существенный региональный факт: участок лежит внутри выраженного положительного магнитного максимума (EMAG2v3 даёт +216,5 нТл на высоте 4 км при региональном фоне до −320 нТл), что согласуется с ожидаемым откликом серпентинизированных ультрамафитов комплекса Соган; и внутри зоны максимального утолщения коры Санандадж-Сирджан (Мохо 42–55 км). Район сейсмически активен: 671 событие в радиусе 100 км по ComCat, включая Mw 5,2 в 4,6 км от центра участка (1988) и Mw 6,0 в 26 км южнее (2006), 1964 события M ≥ 3 в радиусе 100 км по каталогу ISC. Открытая аэрогеофизика Ирана недоступна в виде скачиваемых данных: карты 1:250 000 существуют в бумажном/сканированном виде у коммерческого дистрибьютора, цифровые массивы держат GSI и NGDIR, чьи сайты из данной среды не отвечают. Ограничения 2026 г. на съёмку по Ирану касаются коммерческой оптики высокого разрешения; открытые Sentinel-1, Copernicus DEM и Landsat работают, что подтверждено фактическими загрузками. Главный методический риск, который надо не допустить в отчёте: выдавать значения глобальных моделей за характеристику участка — там один отсчёт сетки на весь объект.

0. Главная рамка: масштабный разрыв между глобальными геофизическими моделями и площадью 62 га

Ключевой и не подлежащий смягчению вывод предресерча: ни одна глобальная геофизическая модель (магнитная или гравитационная) не разрешает объекты размером с участок. Участок 1480 × 419 м = 0,62 км². Ячейка EMAG2v3 на широте 28,36° — 3,71 км по меридиану × 3,26 км по параллели ≈ 12,1 км²; весь участок целиком помещается внутрь ОДНОЙ ячейки и занимает около 5 % её площади. Ячейка WDMAM (3′) — примерно 5,56 × 4,89 км ≈ 27,2 км², участок — 2,3 % ячейки. Спутниковая гравиметрия GOCO06s (d/o 300) имеет предел разрешения ~67 км по полуволне — участок в 45 раз короче. Отсюда практический принцип: глобальные модели — только РЕГИОНАЛЬНЫЙ КОНТЕКСТ (положение офиолитового пояса, крупные структуры, региональные градиенты). Всё, что относится к самому участку, из открытых данных даёт по существу только рельеф (DEM 30 м и производные), радарная интерферометрия Sentinel-1 (пиксель порядка 5 × 20 м) и оптика (не тема этого раздела). Любое утверждение вида «по EMAG2 на участке выделяется аномалия» физически бессмысленно: там один отсчёт сетки.

Источники (2)

1. Глобальные магнитные модели: EMAG2 v3 и WDMAM v2.x — что реально даёт и чего не даёт

EMAG2v3 (Meyer, Saltus, Chulliat, 2017) — глобальная сетка аномалий полного вектора магнитного поля с шагом 2 угловые минуты. Два варианта: «Upward Continued» — пересчёт на постоянную высоту 4 км над эллипсоидом WGS84, и «Sea Level» — на уровне моря для океанских и прибрежных областей. Принципиальное отличие от EMAG2 (2009): версия 3 опирается ТОЛЬКО на измеренные данные и больше не «дорисовывает» аномалии по априорной геологии и возрасту океанской коры; в областях без съёмок она честно показывает отсутствие данных. Над Ираном источник коротковолновой части — та же аэромагнитная съёмка 1974–1977 гг. (см. раздел 2), из которой в национальной компиляции отфильтрованы длины волн короче 15 км; то есть даже «сырьё» EMAG2 над Ираном не содержит сигнала короче ~15 км.

Мной проверено фактически: GeoTIFF EMAG2v3 UpCont читается напрямую по HTTP (через /vsicurl), сетка 10800 × 5399, шаг 0,03333° (3,71 км по широте). В точке центра участка значение аномалии ΔT = +216,5 нТл (на высоте 4 км). В окне 1° × 1° вокруг участка значения от +0,5 до +421,2 нТл — то есть участок сидит внутри отчётливого регионального ПОЛОЖИТЕЛЬНОГО магнитного максимума, что согласуется с ожиданием для серпентинизированных ультрамафитов офиолитового комплекса. В окне 4° × 4° диапазон от −320,9 до +421,2 нТл, СКО 98,3 нТл. Это полезный региональный факт, но он относится к телу размером в единицы–десятки километров, а не к участку.

WDMAM (World Digital Magnetic Anomaly Map) — проект IAGA/CGMW; текущая версия 2.2 (Choi et al., 2025), сетка 3 угловые минуты, приведение к высоте 5 км над уровнем моря. Распространяется в том числе как готовые сетки в GMT (earth_wdmam, разрешения от 1° до 3′). WDMAM ещё грубее EMAG2 и на 62 га неприменима вовсе; её ценность — общая рамка «внешний офиолитовый пояс Загроса как магнитная граница».

Честный вывод: обе модели годятся, чтобы (а) подтвердить, что участок лежит в региональной положительной магнитной аномалии, характерной для ультрамафитового комплекса, и (б) очертить границы этого комплекса с точностью в километры. Для оконтуривания рудных тел хромититов внутри 62 га они бесполезны — там нужна наземная (или БПЛА) магнитная съёмка с шагом 10–25 м.

Источники (7)

2. Аэрогеофизика Ирана: что снято, с каким шагом и доступно ли это за рубежом

Национальная аэромагнитная съёмка Ирана выполнена компанией Aeroservice (Хьюстон) по заказу Геологической службы Ирана (GSI) в 1974–1977 гг. Параметры: двухмоторный самолёт, цезиевый магнитометр чувствительностью 0,02 нТл, СРЕДНЕЕ РАССТОЯНИЕ МЕЖДУ ГАЛСАМИ 7,5 км, 62 съёмочных блока, барометрическая высота полёта 1000–3600 м, что соответствует 500–1000 м над поверхностью земли. Saleh (2006) выполнил повторное уравнивание всех 62 блоков и собрал единую композитную сетку Ирана 1 × 1 км, при этом ИЗ СЕТКИ ОТФИЛЬТРОВАНЫ ДЛИНЫ ВОЛН КОРОЧЕ 15 км. Это принципиально: сетка «1 км» — только шаг интерполяции, реальное содержание — от 15 км и длиннее. При шаге галсов 7,5 км и высоте 500–1000 м никакой объект масштаба 1,5 × 0,4 км в этих данных не проявится.

Эта съёмка была использована для двух известных региональных работ: Teknik & Ghods (2017, GJI) — глубина магнитного фундамента Ирана по фрактальному спектральному анализу; Teknik et al. (2020, Tectonophysics) — новая тектоническая карта Иранского плато с выделением магматических дуг и ОФИОЛИТОВЫХ ПОЯСОВ именно по аэромагнитным данным. Последняя работа релевантна напрямую: внешний офиолитовый пояс Загроса (куда входит Эсфандаге–Фарьяб) прослеживается как магнитный элемент регионального масштаба.

Публикационная форма: карты аэромагнитной интенсивности масштаба 1:250 000 покрывают практически всю страну (цветные надпечатки на серой топооснове; покрытие завершено к 1989 г.). Отдельно есть современная вертолётная магнитно-радиометрическая съёмка Центрального Ирана (заказчик Eagle Mining and Exploration, 2010 г.): около 140 000 км², шаг галсов 500 м (в северо-западной части, район Анарак — 250 м), от Исфахана на западе до района южнее Мешхеда, через Дехшир и Шахр-э-Бабак, до юга Йезда и Кермана. Это высококачественные данные, НО их описанный контур — Урумие-Дохтарский магматический пояс и Центральный Иран; район Урзуийе / Эсфандаге лежит юго-западнее, в зоне Санандадж-Сирджан, и покрытие участка этой съёмкой НЕ подтверждено.

Доступность: открытого международного портала загрузки иранской аэрогеофизики не найдено. Формально данные держат GSI (gsi.ir) и Национальная база геонаук Ирана NGDIR (ngdir.ir, создана в 1999 г. при Министерстве промышленности, рудников и торговли). Из данной рабочей среды оба домена не ответили за 40 с (HTTP 000) — это может быть и сетевое ограничение среды, и недоступность сайтов; проверить нужно отдельно. Бумажные и сканированные листы 1:250 000 (аэромагнитные) и геологические карты 1:25 000 / 1:100 000 / 1:250 000 продаёт коммерческий дистрибьютор East View Geospatial. Публикации иностранных авторов по иранской аэромагнитике почти всегда указывают, что данные получены «by courtesy of GSI» — то есть по индивидуальному запросу, а не из открытого каталога.

Источники (6)

3. Гравитационные модели: WGM2012, GOCE/GOCO, XGM2019e, EGM2008, GGMplus — разрешение и подводные камни

WGM2012 (Bureau Gravimétrique International, CGMW, 2012) — первые глобальные сетки аномалий силы тяжести в сферической геометрии: полная сферическая аномалия Буге, изостатическая (Эри–Хейсканен, Tc = 30 км), поверхностная аномалия в свободном воздухе (Молоденский), плюс гравитационные возмущения поверхностных масс из ETOPO1. Шаг сетки — 2′ × 2′. Исходники — EGM2008 и DTU10; топографические поправки с разрешением 1′ по ETOPO1 (1 угловая минута ≈ 1,8 км). Доступ свободный через BGI по DOI.

GOCO06s (2021) — «чисто спутниковая» модель, d/o 300, что соответствует ~67 км по полуволне; собрана из 1,16·10⁹ наблюдений с 19 спутников (GOCE TIM6, GRACE ITSG-Grace2018s, кинематические орбиты Swarm/TerraSAR-X/TanDEM-X/CHAMP, лазерная локация LAGEOS и др.). Средние–короткие волны от степени 120 определяются градиентометром GOCE; на степенях 280–300 отношение сигнал/шум приближается к 1. Для участка это заведомо только «плита в масштабе провинции».

XGM2019e — комбинированная модель, сфероидальные гармоники до d/o 5399, номинально 2′ (~4 км). Состав: GOCO06s на длинных волнах до d/o 300 + наземные аномалии NGA с разрешением 15′ (~28 км) + топографически смоделированная гравитация EARTH2014 на суше + альтиметрия DTU13 на океане. Важнейшая оговорка самих авторов: на высоких частотах над сушей модель опирается на ТОПОГРАФИЧЕСКОЕ ПРЯМОЕ МОДЕЛИРОВАНИЕ и в хорошо изученных районах уступает моделям с реальными гравиметрическими измерениями. Лицензия CC BY 4.0, доступ через ICGEM.

EGM2008 — d/o 2190 (2159 полный), исходник — 5′ × 5′ аномалии; BGI отдаёт производные сетки аномалий в свободном воздухе и Буге с шагом 2,5′. Оценка по Ирану: СКО расхождений аномалий в свободном воздухе EGM2008 против наземного банка BGI — 17,6 мГал (для сравнения EGM96 — 60,1 мГал). Это ошибка регионального уровня и она больше, чем весь искомый эффект от небольшого рудного тела.

GGMplus (Hirt et al., 2013; обновление GGMplus2023 у BGI) — «сверхвысокое разрешение» ~200 м (шаг 7,2″), суша и прибрежье в полосе ±60° широты. Формально это единственная гравитационная модель, чьё разрешение сопоставимо с участком. Но её короткие волны (от ~10 км до 200 м) получены ПРЯМЫМ МОДЕЛИРОВАНИЕМ ОТ ТОПОГРАФИИ SRTM при постоянной плотности; корреляция коротковолновой гравитации с топографией в GGMplus близка к 100 %. Иначе говоря, на масштабе участка GGMplus не несёт независимой информации о плотности недр — это фактически перекрашенный рельеф. Использовать его как «гравиразведку по участку» нельзя; это была бы методическая ошибка.

Отдельная практическая заметка: если когда-нибудь делать наземную гравиразведку на участке, рельефная поправка будет большой — перепад высот 157 м в пределах прямоугольника и 676 м в окне 10 × 10 км (см. раздел 6), а это в горной местности легко даёт поправки в единицы–десятки мГал.

Источники (8)

4. Сейсмичность как признак активной тектоники: фактические цифры по участку (USGS ComCat, ISC, IRSC)

Здесь получены самые конкретные результаты предресерча — прямыми запросами к открытым каталогам.

USGS ANSS ComCat через FDSN-сервис (запрос с явным starttime=1900-01-01; ВНИМАНИЕ: по умолчанию сервис отдаёт только последние 30 суток — это ловушка, из-за которой первый прогон дал «4 события»): — в радиусе 30 км от центра участка: 125 событий, из них 17 с M ≥ 5 и 1 с M ≥ 6; — в радиусе 50 км: 267 событий, 42 с M ≥ 5; — в радиусе 100 км: 671 событие, 588 с M ≥ 4, 111 с M ≥ 5, 4 с M ≥ 6; — в радиусе 200 км: 217 событий с M ≥ 5, 15 с M ≥ 6.

Ближайшие сильные события: — Mw 5,2, 09.06.1988, 28,386° с.ш. 56,868° в.д., глубина 26,8 км — примерно 4,6 км от центра участка (практически под ним); — Mw 6,0, 28.02.2006, 28,120° с.ш. 56,865° в.д., глубина 18 км — около 26 км южнее; — Mw 5,9, 22.12.1964, 28,114° с.ш. 56,906° в.д., глубина 24,4 км; — Mw 5,9, 18.12.1987, 28,191° с.ш. 56,677° в.д., глубина 19,8 км; — серия 04.03.1999 (Mw 5,2 / 5,1 / 5,0) на удалении 10–20 км. Самое сильное в радиусе 200 км — Ms 6,9, 21.03.1977, 27,609° с.ш. 56,393° в.д. (район Бендер-Аббаса, Хургу).

ISC Bulletin (Международный сейсмологический центр), круговой поиск радиусом 100 км, 1960–2026, M ≥ 3: найдено 1964 события. Разница с ComCat закономерна: ISC собирает данные всех агентств, включая иранские сети, и полнота по слабым событиям заметно выше. Оговорка ISC: данные с 01.07.2024 ещё не отревизованы.

Иранские каталоги: IRSC (Iranian Seismological Center, Тегеранский университет) работает с 2006 г., каталог содержит более 187 000 событий за 2006–2024 гг.; с 2015 г. налажен обмен данными с IIEES, что уменьшило ошибки локации. Известные проблемы каталога IRSC, зафиксированные в рецензируемых работах: неоднородность магнитуды полноты Mc в пространстве и во времени (с 2006 г. Mc улучшилась в разных районах на 0–0,5 единицы магнитуды), и наличие взрывов — статистический алгоритм отнёс к карьерным и рудничным взрывам около 24 % всех событий каталога. При работе с иранскими данными это надо чистить.

Интерпретация для участка: район бесспорно сейсмически активен (внутренняя часть коллизии Аравия–Евразия, юго-восточное окончание Загроса и переход к Макрану). Глубины преимущественно 10–33 км, то есть коровые; часть глубин в каталоге фиксирована (10, 15, 33 км) и не должна восприниматься как измеренная. Для 62 га это означает: сейсмичность — фактор проектирования (сейсмостойкость сооружений, устойчивость откосов и отвалов), а НЕ поисковый признак. Разрешение эпицентров (горизонтальная ошибка в приведённом примере 2026 г. — 15,1 км) на два порядка грубее размеров участка.

Источники (6)

5. Мощность коры и глубина Мохо: только региональный уровень

Прямых измерений Мохо под участком нет и быть не может — ближайшие сейсмические станции и профили отстоят на десятки–сотни километров. Что известно по региону:

Профильные функции приёмника через Центральный Загрос (Paul et al., 2006, GJI; 45 временных станций на профиле длиной 620 км, среднее расстояние между станциями 14 км): от побережья Персидского залива до 25 км юго-западнее Главного Загросского надвига (MZT) Мохо почти горизонтальна на глубине около 45 км; затем мощность коры резко нарастает до максимума ~70 км под метаморфической зоной САНАНДАДЖ-СИРДЖАН (в 50–90 км северо-восточнее выхода MZT на поверхность); далее к северо-востоку Мохо поднимается до ~42 км под Урумие-Дохтарским магматическим поясом и югом Центрально-Иранского микроконтинента. Гравитационное моделирование в той же работе согласуется с поддвигом коры Загроса под кору Центрального Ирана вдоль MZT как структуры кораового масштаба.

Комбинированные гравиметрически-сейсмические модели Мохо по Ирану (коллокация): сетка 0,5° × 0,5° (то есть ячейка ~55 × 49 км), исходники — спутниковая модель GOCO03S до d/o 180, рельеф SRTM30_PLUS, кора CRUST1.0, 213 сейсмических точек Мохо (140 в инверсию, 73 на контроль). В зоне Санандадж-Сирджан решение даёт до 55 км, контрольные сейсмические точки 42–55 км. Точность: чисто гравиметрический вариант — СКО 6,2 км, комбинированный — 4,9 км. Отдельная работа по быстрой коллокации по данным GOCE использует средний уровень Мохо 40 км как опорный для устранения неоднозначности инверсии.

Вывод: цифра «Мохо под участком» из открытых моделей получается интерполяцией по сетке 0,5° с погрешностью около 5 км — это допустимо только как фраза «кора здесь мощная, порядка 45–55 км, район находится в области максимального утолщения коры зоны Санандадж-Сирджан». Никакой поисковой ценности для 62 га это не имеет; ценность — контекстная (тектонический режим, обдукция офиолитов, поздняя коллизия).

Источники (4)

6. Рельеф: какой DEM брать и что он показывает НА САМОМ УЧАСТКЕ (измерено)

Рекомендуемая основа — Copernicus DEM GLO-30 (шаг 1″). Это цифровая модель ПОВЕРХНОСТИ (DSM) на базе WorldDEM/TanDEM-X, съёмка 2011–2015 гг., с заполнением пропусков из ASTER, SRTM, GMTED2010, TerraSAR-X и ALOS World 3D-30m. Заявленная точность: абсолютная по высоте <4 м (LE90), относительная <2 м на склонах до 20 % и <4 м на склонах круче 20 % (в пределах ячеек 1° × 1°), плановая <6 м (CE90). Выпущено девять редакций с декабря 2019 г. по июль 2024 г. (последняя 2024_1). Лицензия свободная, требуется атрибуция «© DLR e.V. 2010-2014 and © Airbus Defence and Space GmbH 2014-2018 provided under COPERNICUS by the European Union and ESA».

Проверено фактически: тайл N28 E056 доступен в открытом бакете AWS (copernicus-dem-30m), запрос диапазона отдаёт HTTP 206 — то есть НИКАКИХ ограничений доступа по Ирану у GLO-30 нет, вопреки распространённому опасению (ограничения когда-то касались Армении, Азербайджана и Молдовы, они сняты в редакции 2023_1 от июля 2024 г.).

ИЗМЕРЕНИЯ ПО САМОМУ УЧАСТКУ (мои прямые чтения GLO-30, окно по четырём заданным углам UTM 40N, пересчитанным в WGS84: 56,829484–56,844592° в.д., 28,353297–28,357096° с.ш.): — размер окна 14 × 54 пикселя ≈ 756 ячеек; шаг пикселя 27,2 м по долготе и 30,9 м по широте; — высоты: минимум 1597,5 м, максимум 1754,5 м, среднее 1664,5 м, ПЕРЕПАД 157,0 м на прямоугольник 1480 × 419 м; — уклоны: медиана 18,9°, 75-й процентиль 26,0°, 90-й процентиль 30,9°, максимум 38,8°; — экспозиция склонов: С 19,4 %, СВ 18,3 %, В 12,2 %, ЮВ 8,7 %, Ю 7,4 %, ЮЗ 8,1 %, З 11,2 %, СЗ 14,7 % — то есть преобладают склоны северной и северо-восточной экспозиции. В окне 10 × 10 км вокруг центра: высоты 1349,5–2025,9 м, перепад 676,5 м, среднее 1535,1 м.

Это крутой горный рельеф. Практические следствия: (1) любые наземные работы затруднены, шаг сети придётся привязывать к склонам; (2) для гравиразведки понадобятся серьёзные рельефные поправки; (3) для магнитной съёмки с БПЛА нужен режим следования рельефу; (4) оптические и радарные данные будут иметь заметные ракурсные искажения и тени в зависимости от экспозиции.

Альтернативные DEM: NASADEM и SRTM v3 (1″, вертикальная точность порядка ±9 м), ALOS AW3D30 (1″, ±5 м; в независимых сравнениях по геодезическим точкам в горах Рифа RMSE 7,77 м против 9,09 м у Copernicus — то есть в отдельных горных тестах AW3D30 выигрывает), ASTER GDEM3 (заметно шумнее), FABDEM (Copernicus с удалёнными лесом и застройкой — для голых аридных гор преимущества минимальны), TanDEM-X 90 м (свободно от DLR) и TanDEM-X 12 м (только по научной заявке в DLR). Разумная стратегия: базовый анализ на GLO-30, кросс-проверка морфометрии на AW3D30 и NASADEM, чтобы отделить реальные формы рельефа от артефактов конкретной модели.

Источники (6)

7. Производные DEM для структурного анализа: лайнаменты, дренаж, морфометрия — алгоритмы и на чём считать

Это самый содержательный из открытых методов на масштабе участка и его ближайшего обрамления. Схема работ, подтверждённая литературой:

ЛАЙНАМЕНТЫ. Классика — автоматическое выделение по многоазимутальным теневым отмывкам (hillshade) с DEM: набор отмывок при разных азимутах и углах освещения (обычно 8 азимутов через 45° при высоте солнца 30–45°), детектирование краёв (Sobel, Canny) и связывание в линии. Реализации: модуль LINE в PCI Geomatica (комбинация детектирования краёв и извлечения кривых, параметры RADI/GTHR/LTHR/FTHR/ATHR/DTHR), преобразование Хафа в MATLAB, метод MHHC (Multi-Hillshade Hierarchic Clustering) — иерархическая кластеризация «протолайнаментов» с разных отмывок и классификация на положительные и отрицательные лайнаменты. Обязателен контроль: отмывка при одном азимуте систематически подавляет линии, параллельные направлению освещения — отсюда требование многоазимутальности. Дальше — розы-диаграммы направлений, плотность лайнаментов, пересечения (узлы) как потенциальные структурные ловушки.

ДРЕНАЖ. Стандартный конвейер: заполнение локальных понижений → направление стока (D8 или D-infinity) → аккумуляция стока → пороговое выделение водотоков → порядок по Стралеру → водосборы. Инструменты открытого стека: WhiteboxTools, TauDEM, SAGA GIS, GRASS GIS (r.watershed, r.stream.*), PySheds/RichDEM в Python. Анализируются: плотность дренажной сети, тип рисунка (дендритовый / решётчатый / параллельный — последние два указывают на структурный контроль), коленчатые изгибы и обезглавленные долины как признаки разрывов, асимметрия водосборов.

МОРФОТЕКТОНИЧЕСКИЕ ИНДЕКСЫ (для района, не для участка): индекс градиента длины водотока SL и его нормированная форма SLk (чувствителен к перегибам профиля от разломов и дифференциальных поднятий), гипсометрический интеграл HI, отношение ширины днища долины к её высоте Vf, извилистость горного фронта Smf, коэффициент асимметрии водосбора AF, индексы крутизны ksn и χ-анализ. Есть готовые наборы функций MATLAB для извлечения гипсометрии, SL, ksn, χ и swath-профилей из DEM.

УКЛОНЫ И УСТОЙЧИВОСТЬ. Уклон, кривизна (плановая и профильная), TPI, TRI/шероховатость, геоморфоны — карты классов форм рельефа. По моим измерениям уклонов (медиана 18,9°, p90 30,9°, максимум 38,8°) участок относится к категории крутых склонов: это и повышенный риск осыпей/оползней, и жёсткое ограничение на технику.

ВАЖНАЯ ОГОВОРКА О МАСШТАБЕ. На 30-метровом DEM минимальный надёжно выделяемый линейный элемент — порядка 3–5 пикселей, то есть ~90–150 м. Значит внутри участка (419 м в поперечнике = 14 пикселей) можно распознать лишь единичные крупные линеаменты; статистика направлений и плотности набирается только на площади в десятки–сотни км². Практический вывод: строить структурную схему по DEM надо для окна 10–30 км, а на участок переносить уже готовые направления и проверять их наземно. Если нужен структурный анализ внутри 62 га — требуется съёмка с БПЛА (фотограмметрия SfM, DEM 5–20 см), это уже не открытые данные.

Источники (5)

8. InSAR Sentinel-1: реальное покрытие участка, геометрия съёмки и что можно измерить

ПРОВЕРЕНО ПО АРХИВУ. Запрос к поисковому API Alaska Satellite Facility по точке центра участка (POINT 56.83704 28.35520), уровень SLC: — всего в архиве 934 гранулы Sentinel-1 SLC (на 27.07.2026); — 234 гранулы с 01.01.2024; — 142 гранулы с 01.01.2025, распределённые по трём трекам: • восходящий трек 57, кадр 88 — 46 сцен (Sentinel-1A 43, Sentinel-1D 3); • нисходящий трек 93, кадр 499 — 48 сцен (S1A 45, S1D 3); • нисходящий трек 166, кадр 499 — 48 сцен (S1A 46, S1D 2). То есть участок покрыт одним восходящим и двумя нисходящими треками, съёмка регулярная, интервал по каждому треку около 12 суток. Отчётливо видно, что с 2026 г. в работу вошёл Sentinel-1D (S1D-гранулы с апреля–июля 2026), Sentinel-1C введён в эксплуатацию 26.03.2025 после запуска 05.12.2024.

ГЕОМЕТРИЯ. Sentinel-1 в режиме IW: углы падения 29–46°, разрешение SLC порядка 5 × 20 м (наклонная дальность × азимут), длина волны C-диапазона 5,6 см. Я посчитал геометрические искажения по GLO-30 при типичном угле падения 39°: — восходящая орбита (взгляд на восток): медианный локальный угол падения 37,0°, минимум 4,4°; площадь в перекрытии (layover) 0,0 %, в тени 0,0 %; — нисходящая орбита (взгляд на запад): медианный локальный угол падения 44,8°, минимум 9,4°; перекрытия и тени 0,0 %. То есть, несмотря на крутизну (максимум 38,8°), из-за преимущественно северной/северо-восточной экспозиции склонов критических искажений на участке нет. Нисходящая геометрия выглядит предпочтительнее (больший локальный угол падения → лучше сохраняется когерентность и меньше сжатие в дальности). Минимум 4,4° на восходящей орбите означает отдельные пиксели на грани перекрытия — их надо маскировать.

ЧТО РЕАЛЬНО ИЗМЕРИМО. Аридная, почти безрастительная местность — благоприятный случай для C-диапазона: когерентность сохраняется долго. Достижимая точность средних скоростей смещения по линии визирования при многолетних рядах — единицы мм/год. Прикладные задачи для 62 га: медленные оползни и крип склонов, осадки/подвижки бортов и отвалов при горных работах, просадки над подземными выработками, деформации от откачки подземных вод. Ограничения: измеряется только проекция на линию визирования (по одному восходящему и одному нисходящему треку раскладывается на восточную и вертикальную компоненты; северная компонента принципиально недоступна); быстрые подвижки >~14 мм за 12 суток (четверть длины волны на пиксель) приводят к срыву развёртки фазы; тропосферная задержка в горах требует коррекции (GACOS).

ГОТОВЫЕ ПРОДУКТЫ. Портал COMET-LiCS (Университет Лидса) автоматически считает интерферограммы и карты когерентности по кадрам (LiCSAR на базе GAMMA), продукты выкладываются в открытый доступ; поверх них считаются временные ряды пакетом LiCSBAS. По Ирану уже сделан национальный продукт: Watson et al. (2024, GRL) — поле скоростей InSAR+GNSS для площади 1,7 млн км², 7 лет наблюдений Sentinel-1, около 85 000 интерферограмм из ~18 500 снимков по 9 восходящим и 10 нисходящим трекам, обработка COMET-LiCSAR + GACOS + LiCSBAS, привязка к системе отсчёта с неподвижной Евразией, разложение на восточную и вертикальную компоненты с ШАГОМ 500 м. Это готовая региональная канва: по ней видно, есть ли на площади аномальные скорости, а участок в ней занимает 3 × 1 ячейки — то есть служит только индикатором «есть/нет крупного сигнала», детальную картину надо считать самостоятельно по SLC.

Источники (8)

9. Геологический контекст участка, релевантный геофизике (что уже опубликовано по комплексу Соган)

Комплекс Соган (Soghan mafic–ultramafic complex) — одно из крупнейших хромитовых полей юго-востока Ирана, часть комплексов Эсфандаге и Фарьяб (провинция Керман). Породы: ритмичное переслаивание дунитов, хромититов, гарцбургитов, с подчинёнными пироксенитами, верлитами и габбро; примерно 70 % хромитовых залежей — прерывистые слои, остальное — линзовидные тела. Годовая добыча по комплексам Эсфандаге и Фарьяб — порядка 50 000 т хромитовой руды; Эсфандаге считается вторым производителем хромита в Иране после Фарьяба. По геохимии ЭПГ и золота (Sci. Direct, 2014) в хромититах Согана концентрации IPGE выше 100 ppb, что указывает на бонинитовый расплав и альпинотипный тип оруденения.

По дистанционному зондированию комплекс изучен: работы по литологическому картированию Согана и Эсфандаге спектральными преобразованиями ASTER и Landsat TM (коэффициент корреляции CC, индекс оптимальности OIF, отношения каналов, метод главных компонент). Это подтверждает, что открытые оптические данные над районом работают, но задача этого раздела — геофизика и рельеф.

Геофизическое следствие: ультрамафиты с серпентинизацией дают сильный магнитный контраст, а хромитит сам по себе слабомагнитен и ПЛОТЕН (плотность хромитита ~4,0–4,5 г/см³ против ~2,7–3,3 г/см³ у вмещающих серпентинизированных перидотитов). Это классическое сочетание для локальной наземной магниторазведки (картирование серпентинизированных зон и структурного каркаса) и микрогравиразведки (прямой поиск плотных линз). Никакая глобальная модель этого масштаба не увидит; речь именно о наземных работах на 62 га.

Источники (4)

10. Доступ к данным в 2026 году: что работает, что ограничено

Что проверено и РАБОТАЕТ из этой среды на 27.07.2026: — USGS FDSN Event Web Service (каталог ANSS ComCat) — отвечает, выдаёт CSV/GeoJSON; — ISC Bulletin (веб-форма поиска по каталогу) — отвечает, выдаёт CSV; — ASF Search API (архив Sentinel-1) — отвечает, 934 SLC над участком; — открытый бакет AWS Copernicus DEM GLO-30 и GLO-90 — тайл N28 E056 отдаётся (HTTP 206); — NOAA NCEI EMAG2v3 GeoTIFF — читается по HTTP диапазонами.

Что НЕ ОТВЕТИЛО: gsi.ir и www.ngdir.ir — таймаут 40 с, код HTTP 000. Причина не установлена: это может быть блокировка исходящего трафика в данной среде, геоблокировка со стороны Ирана, или временная недоступность. Требуется отдельная проверка из другой сети.

Отдельный сюжет 2026 года: ограничения на КОММЕРЧЕСКУЮ спутниковую съёмку высокого разрешения по Ирану. По сообщениям (SpaceNews, Bloomberg, CBS, GIJN, апрель 2026) правительство США попросило операторов добровольно ввести бессрочную задержку публикации снимков широкой полосы Ближнего Востока задним числом с 9 марта 2026 г.; Planet Labs объявила о 14-суточной задержке публикации новых снимков Ирана, Персидского залива и зон конфликта; Vantor заявила о «расширенном контроле доступа» в периоды геополитических конфликтов. ВАЖНО: эти ограничения касаются коммерческих операторов оптики высокого разрешения. На открытые программы (Copernicus Sentinel, Landsat, SRTM, Copernicus DEM) ограничения не распространяются, и это подтверждено эмпирически: Sentinel-1 SLC над участком доступны вплоть до 24.07.2026, GLO-30 скачивается. Тем не менее для планирования работ это риск, который надо держать в виду: политика может меняться.

Лицензии: Copernicus DEM GLO-30/GLO-90 — свободная всемирная лицензия с обязательной атрибуцией DLR/Airbus/Copernicus/ESA; Copernicus Sentinel — свободный полный открытый доступ (с оговоркой в законодательстве ЕС о возможных ограничениях в редких случаях угроз безопасности); XGM2019e — CC BY 4.0; EMAG2v3 и WGM2012 — свободно, требуется цитирование.

Источники (4)

11. Итоговая матрица: что осмысленно на 62 га, а что только на район в сотни км²

ПРИМЕНИМО К САМОМУ УЧАСТКУ (0,62 км²): 1. Copernicus DEM GLO-30 — 756 пикселей на участок. Даёт высоты, уклоны, экспозицию, кривизну, микроводоразделы. Уже посчитано: 1597,5–1754,5 м, перепад 157 м, медианный уклон 18,9°. 2. Sentinel-1 InSAR — 934 SLC в архиве, три трека, 12-суточный интервал; при пикселе 5 × 20 м на участок приходятся сотни независимых отсчётов. Реально: скорости смещений с точностью единицы мм/год, мониторинг склонов и бортов. 3. Отдельные крупные лайнаменты (>100–150 м) по DEM — но только как элементы, продолжающиеся за пределы участка.

ПРИМЕНИМО ТОЛЬКО К РАЙОНУ (десятки–сотни км² и больше): 4. EMAG2v3 (12,1 км² на ячейку) и WDMAM (27,2 км²) — положение и контур магнитной аномалии офиолитового комплекса. Установлено: участок внутри регионального положительного максимума, ΔT ≈ +216 нТл на высоте 4 км. 5. Национальная аэромагнитная съёмка Ирана (галсы 7,5 км, фильтр <15 км) — региональный структурный каркас, границы поясов. 6. WGM2012 (2′), XGM2019e (2′), EGM2008 (2,5′), GOCO06s (67 км) — региональное гравитационное поле, изостазия, крупные плотностные неоднородности. 7. GGMplus (200 м) — формально разрешает участок, но коротковолновая часть является производной топографии, а не измерений; как гравиразведка непригодна. 8. Мохо и мощность коры (сетки 0,5°, СКО 4,9–6,2 км) — только тектонический контекст. 9. Сейсмичность и активная тектоника (горизонтальные ошибки эпицентров порядка 10–15 км) — оценка сейсмической опасности, но не поиск. 10. Морфотектонические индексы SL, HI, Vf, Smf, AF — требуют бассейнов и горных фронтов, то есть площади в десятки–сотни км².

ЧЕГО ОТКРЫТЫЕ ДАННЫЕ НЕ ДАЮТ СОВСЕМ: — локальных магнитных и гравитационных аномалий от рудных тел внутри 62 га; — электроразведочных данных (ЭМ, ВЭЗ, ВП) — открытых глобальных аналогов не существует; — сейсморазведки любых масштабов; — скважинной информации; — субметрового рельефа для микроструктурного анализа (нужна съёмка с БПЛА).

Источники (2)

Что дальше по этому направлению

  1. Скачать тайл Copernicus DEM GLO-30 N28 E056 и построить полный набор производных на окно 20 × 20 км вокруг участка: уклон, экспозиция, плановая и профильная кривизна, TPI, TRI, геоморфоны; отдельно — карту зон уклона >30° как ограничение для техники и признак осыпей
  2. Сделать автоматическое выделение лайнаментов по многоазимутальным отмывкам (8 азимутов через 45°, высота солнца 35°) с кросс-проверкой на AW3D30 и NASADEM, чтобы отсеять артефакты одной модели; построить розу-диаграмму направлений и карту плотности лайнаментов, сопоставить с известным простиранием офиолитового пояса Эсфандаге–Фарьяб
  3. Извлечь дренажную сеть (WhiteboxTools или GRASS r.watershed: заполнение понижений → D8 → аккумуляция → порог) на окно 20 × 20 км, оценить рисунок сети и коленчатые изгибы как индикаторы структурного контроля; посчитать SL/SLk и гипсометрический интеграл по бассейнам
  4. Собрать временной ряд InSAR по нисходящему треку 93 или 166 (у них выше локальный угол падения — 44,8° против 37,0° у восходящего): выгрузить готовые интерферограммы с портала COMET-LiCS, обработать LiCSBAS с коррекцией тропосферы GACOS, получить карту средних скоростей по линии визирования и оценить, есть ли на участке сигнал выше уровня шума
  5. Скачать сетки WGM2012 (аномалия Буге и изостатическая) и EMAG2v3, вырезать окно 2° × 2° и подготовить региональные карты-подложки — но с явной подписью разрешения ячейки и обязательной оговоркой, что участок занимает 5 % одной ячейки
  6. Выгрузить полный каталог ISC по кругу радиусом 100 км (1964 события M ≥ 3) и построить карту эпицентров с эллипсами ошибок, разрез глубин, кумулятивный график Гутенберга–Рихтера; отделить фиксированные глубины (10/15/33 км) от измеренных
  7. Проверить доступность gsi.ir и ngdir.ir из другой сети; если сайты живы — установить, есть ли в открытом виде листы аэромагнитных карт 1:250 000 на район Урзуийе/Эсфандаге и цифровые данные вертолётной съёмки 2010 г.; параллельно уточнить у East View Geospatial наличие соответствующего листа
  8. Установить точное название и номер листа геологической карты GSI 1:100 000 на район участка (кандидаты — Соган, Эсфандаге, Урзуийе) и листа аэромагнитной карты 1:250 000; это нужно для любых дальнейших запросов
  9. Запланировать наземные работы, которые действительно решают задачу на 62 га: магнитная съёмка (наземная или с БПЛА, шаг 10–25 м по профилям вкрест простирания) для картирования серпентинизированных зон, и микрогравиразведка для прямого поиска плотных хромититовых линз; при съёмке с БПЛА обязателен режим следования рельефу при уклонах до 39°
  10. Заказать или снять фотограмметрию с БПЛА (SfM, DEM 5–20 см) на весь участок — только это даст структурный анализ внутри 62 га; открытые данные этого не заменят