Keywords: 1848, Iconoclasme, Italie, Mezzogiorno, Politisation populaire, Révolution 1848, Iconoclasia, Italia, Mezzogiorno, Politización popular, Revolución 1848, Iconoclasm, Italy, Mezzogiorno, Popular politicization, Revolution Carte ?

Les acteurs ordinaires face aux images de la monarchie. Autour de l'iconoclasme révolutionnaire dans les Calabres (1848)

Les périphéries rurales du royaume des Deux-Siciles ont vu s'affirmer, dès les premières décennies du XIX e siècle, une opposition politique dite « libérale » à la monarchie des Bourbons. Plus spontanée qu'idéologique, elle a largement exploité la conflictualité ordinaire et les coutumes comme autant de ressorts de l'action collective. En ce sens et par rapport aux modèles dominants de la « naissance de la politique moderne » qui marquent l'Europe méridionale au lendemain de l'époque révolutionnaire, elle constitue un espace alternatif du politique 1 . Dans une société marquée par une très forte religiosité populaire, bien décrite par l'anthropologie, la mobilisation du religieux a été un facteur de politisation décisif, qui s'est situé du côté des émotions protestataires. Chez les libéraux, elle a revêtu des formes diverses, qu'il s'agisse des prêches révolutionnaires, du culte des martyrs ou des appels à la « croisade » contre la tyrannie des Bourbons et contre la domination exercée par l'Autriche sur une grande partie de la péninsule italienne 2 . La révolution de 1848 en a vu les attestations les plus nombreuses, alors que les libéraux cherchaient à s'assurer le soutien des masses, à la suite d'un mouvement pris en charge depuis les années 1830 par les élites lettrées de l'opposition, qui ont érigé le peuple en objet de savoir afin de préciser l'orientation à donner aux mobilisations libérales 3 . L'iconoclasme politique relève de ces emprunts au religieux, même s'il procède de conditions différentes. Il est d'abord le fait d'acteurs ordinaires, et s'inscrit dans une généalogie complexe qui le relie à la fois aux formes chrétiennes de la violence et aux pratiques de la conflictualité sociale dont on trouve la trace dans les Calabres depuis l'époque moderne 4 . Lors des insurrections populaires de 1799 puis, en 1806, lors de la résistance locale à l'invasion des armées françaises, les destructions d'images ont été pratiquées chez les partisans de la contre-révolution. C'est plutôt son adaptation au cadre révolutionnaire qui est nouvelle : depuis la révolution de 1820-1821, l'iconoclasme fait partie du répertoire d'action courant des acteurs ordinaires, en particulier dans les franges les plus radicales de l'opposition. À plus large échelle, sa place dans les révolutions européennes du XIX e siècle suscite l'intérêt de l'historiographie récente 5 . Pour l'Italie, des travaux ont montré la récurrence de la pratique à la faveur des révolutions de 1848, en particulier dans des espaces politiques périphériques 6 . Ils en décrivent les mécanismes, la diversité des terrains d'application, l'articulation avec des mobilisations encore largement construites à l'échelle locale. Dans les Calabres, la destruction des images de la monarchie a répondu à des revendications précises autour de la souveraineté locale, qui ont suscité un large consensus populaire depuis la politique de centralisation engagée par la monarchie d'occupation du decennio francese (1806-1815) et maintenue par la suite par les Bourbons de Naples, au moment de la Restauration de 1815 7 . Dans une société au très fort fractionnement local, structurée autour des élites municipales ou ecclésiastiques traditionnelles, la diffusion de l'iconoclasme révolutionnaire est une manifestation de la diversité des scénarios de la politisation libérale, sur lesquels insiste l'historiographie actuelle de l'Italie 8 .

Les sources policières et surtout judiciaires, pléthoriques, renseignent sur la récurrence et la diversité des pratiques iconoclastes. Elles ont suscité les procès de 564 sujets dans la province de Cosenza et de 74 autres dans la province de Reggio, pour des actes perpétrés pour la plupart de façon collective et au cours de la seule année 1848. Ces procès politiques, inédits, intègrent donc la destruction des images au sein de pratiques contestataires plus larges. Ils ne présentent pas l'iconoclasme comme un chef d'accusation à part entière, le présentant comme une pratique nécessairement communautaire qui relève du trouble à l'ordre public, raison explicite invoquée pour justifier l'inculpation 9 . En effet, pour les acteurs méridionaux de l'époque, l'iconoclasme est qualifié selon la forme précise que prennent les destructions. On parle alors de bris de statues (infrangimento), de réduction en éclats (frantumazione), de lacération d'emblèmes (lacerazione), sans qu'il y ait dans les sources un terme générique pour désigner ce phénomène. À la différence de la France, la pratique est absente des quelques témoignages littéraires dont on dispose sur la révolution de 1848, qu'ils émanent de ses partisans ou de ses adversaires10 . Plus encore, l'acte n'est pas forcément perçu dans sa dimension politique : il relève le plus souvent de la police « générale », alors que les faits politiques et les atteintes à l'autorité de l'État dépendent de l'Alta Polizia, à la suite des réformes que le royaume connaît dans les premières décennies du XIX e siècle11 . C'est donc la politisation de l'iconoclasme en révolution qu'on voudrait interroger. Il faudra d'abord en saisir la nature, par-delà les dénominations multiples utilisées par les contemporains 7 Sur la Calabre, voir FERRARI, 2016. Sur la politique de centralisation du decennio francese et les réactions qu'elle a suscitées, voir DELPU, 2015.

8 DE SENSI SESTITO, PETRUSEWICZ (dir.), 2014. 9 Sur ces procès, voir DELPU, 2020. Les chiffres donnés s'appuient sur leur analyse au moyen d'une base de données relationnelles. En revanche, il s'agit d'estimations basses, la répression qui a suivi les révolutions de 1848 n'ayant touché qu'une partie des insurgés et n'étant documentée que de façon partielle (les provinces de Catanzaro et de Reggio, par exemple, ne disposent pas d'archives comparables).

et la diversité des cas, pour comprendre comment il s'inscrit dans des répertoires de protestation complexes qui visent la destitution symbolique du roi et l'imposition de communautés politiques locales autonomes.

Fracturer, pendre, fusiller : des gestes politiques Si l'iconoclasme n'est pas envisagé en tant que tel par les acteurs de l'époque, c'est surtout parce qu'il relève d'une attitude spontanée de la part de sujets a priori étrangers à l'espace politique. Le vocabulaire employé par les acteurs montre les formes de la violence dans leur diversité, du simple bris de statue ou de la lacération d'un drapeau à l'incendie d'objets signifiants, porteurs de représentations du pouvoir royal. Ces gestes destructeurs permettent de restituer le répertoire d'action des insurgés12 . Construits et contextualisés, ils sont porteurs d'un sens politique, qui se définit par rapport au régime en place. Ils visent des supports courants des représentations ordinaires de l'autorité de l'État. Parce qu'ils sont fortement liés au développement des camps de révoltés et parce qu'ils relèvent de démonstrations publiques, les actes d'iconoclasme s'inscrivent dans des stratégies de marquage de l'espace qui consacrent la mise en place de communautés locales autonomes. Les résistances qui leur ont été opposées ont été très minoritaires, y compris de la part des forces de l'ordre, qui ont largement soutenu l'insurrection. L'interposition d'un garde national, à Rogliano en juin 1848, lors de l'outrage d'un drapeau et d'une statue par plusieurs sujets en partance pour les camps de révoltés est demeurée une exception 53 . Il semble plutôt que les élites locales aient été des soutiens de ce processus, encourageant parfois les sujets de la province à détruire les emblèmes de la royauté. Outre les interventions répétées des curés de village, les autorités civiles ont parfois incité aux destructions des statues. C'est le cas du chancelier communal de Pellaro, dans la province de Reggio, Bruno Romeo, qui ordonne en avril 1848 à l'un des administrateurs des taxes de détruire l'inscription « Ferdinando Re » placée sur la façade de la mairie, « parce qu'il n'était plus roi »54 . Parce qu'il s'est inscrit dans les formes locales de la résistance à l'autorité du roi, l'iconoclasme a été un acte de défiance à l'égard de celle-ci, qu'il s'agissait de ne plus reconnaître. Les destructions se sont souvent accompagnées de proclamations de la mort du souverain, en conformité avec les revendications du par le mouvement libéral calabrais. Reggio, en juillet 1848, la destruction du télégraphe par quatre artisans de la commune a ainsi constitué un préalable à l'outrage aux statues du roi 61 . Au total, l'iconoclasme s'inscrit dans un processus insurrectionnel plus large, qui a cherché à entretenir le fonctionnement local des communautés en profitant à la fois de leur caractère auto-structuré et de leur faible ouverture sur l'extérieur. Mais il interroge la réalité des transferts de sacralité qu'induisent les destructions, entre élaboration de projets politiques propres à l'évolution des sociétés au XIX e siècle et maintien des logiques politiques traditionnelles de la révolte.

FIG. 1

 1 FIG. 1. -Typologie des supports des destructions d'emblèmes ou de symboles de l'autorité (provinces de Cosenza et de Reggio, 1848, élaboration propre) Pour les trois quarts d'entre elles, les destructions visent des blasons royaux, exposés dans l'espace public, ou des statuettes en grès du couple royal, postées dans les lieux de pouvoir. Dans des provinces

FIG. 2

 2 FIG. 2. -Répartition des lieux-cibles des attentats iconoclastes (provinces de Cosenza et de Reggio, 1848, élaboration propre) La répartition des lieux visés par les actes iconoclastes est moins le fait d'une volonté assumée de la part des acteurs que de la présence matérielle des emblèmes de la monarchie, systématiques dans les mairies et les postes de garde, plus rares ailleurs, où leur existence n'est imposée par aucun texte juridique. Mais les destructions sont ritualisées, les objets-cibles étant le plus souvent transportés sur la place communale, lieu essentiel de la vie sociale dans l'Italie du XIX e siècle 30 , pour y être données en spectacle, en particulier dans le cas des statues. Fin mars 1848 dans le village de Luzzi, l'avocat Gatti déplace en place publique la statue du roi qui se trouvait alors dans la mairie, afin de la fracturer et de la brûler publiquement. Lors de son procès en avril 1849, il justifie ce transfert par le fait que son village était « bien plus attaché au roi que tout autre », et ne pouvait être acquis aux libéraux que par des actes aisément compréhensibles de tous et à la très forte portée symbolique 31 . Dans des provinces rurales où l'implantation du courant libéral est récente et encore très parcellaire, une très importante conflictualité demeure entre partisans et adversaires de la monarchie, sans qu'il existe d'options idéologiques

51

  Archivio di Stato (Napoli), Archivio Borbone, b. 1044, f. 38. 52 Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 88 bis, f. 517. 53 Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 35 bis, f. 220. On trouverait un autre cas à Montalto, où deux gardes nationaux se sont interposés lors de la destruction des statues du couple royal. Voir Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 88 bis, f. 517.

  au fusil sur les restes d'un blason royal qu'ils avaient préalablement jeté à terre, fracturé puis piétiné, à l'été 184819 . Les attestations de ce type de destructions sont nombreuses, et s'inscrivent parfois dans des scénarios plus complexes qui articulent plusieurs cibles et plusieurs modes de destruction.Lorsque quatre sujets d'Aprigliano, en périphérie sud-est de Cosenza, cherchent à épurer l'espace du village de tous les symboles royaux, à l'été 1848, ils exécutent d'abord au fusil les deux bustes du couple royal situés sur le poste de garde, puis lacèrent, à l'aide d'un couteau de boucher, l'emblème de la monarchie placé sur la façade d'une boutique privée 20 . La rareté des emblèmes de l'autorité centrale facilite ce processus d'élimination, qui prend la forme d'un rituel collectif impliquant le peuple en armes.Les sources renseignent en effet sur les formes de la prise d'armes. L'interdiction de la circulation des armes à feu, entérinée par une importante législation prise au début des années 1830 afin de prémunir le royaume contre une révolution potentielle à l'image de celles qui ont traversé l'Europe quelques années prêtre Ferdinando Vigna, à Aprigliano, qui exécute la statue de la reine de plusieurs coups de pistolet, le 2 juillet 1848 25 . Cet usage des armes interdites n'est pas propre aux Calabres et connaît des équivalents dans la plupart des révolutions européennes, où il contribue à former l'espace public d'opposition26 . Par son aspect collectif, il joue un rôle de socialisation politique indéniable, et suppose parfois l'appui d'insurgés rompus au métier des armes, en particulier des gardes nationaux. À Majera dans la province de Cosenza, les frères Francesco et Tommaso Brani, anciens gardes nationaux pendant la révolution de 1820-1821 qui ont été réinvestis dans leurs fonctions en 1848, ont fourni des armes à des sujets de la commune en vue de détruire les emblèmes de la royauté 27 . Ces soutiens ont été facilités par l'appui des Gardes Nationales calabraises à l'insurrection 28 . Mais leur investissement dans les pratiques iconoclastes n'a été qu'occasionnel, celui-ci relevant la plupart du temps de l'expérience des acteurs ordinaires. Il montre comment les insurgés ont investi l'espace public à travers la mise en scène des destructions, qui représentent des tyrannicides symboliques et ritualisés.

	périphériques, ce sont souvent les seules représentations visibles de l'autorité. Les visites qu'y a
	effectuées le roi Ferdinand II, en 1833, 1844 et 1848, n'ont pas suffi à y diffuser durablement et
	massivement l'image du pouvoir en place. Dans la mairie de Paternò, par exemple, les deux statues du
	roi et de la reine détruites en juin 1848 sont signalées comme les deux seules du village 13 . La rareté des
	représentations de la monarchie explique la difficile appropriation des signes politiques par les acteurs
	ordinaires, habitués aux autonomies locales qui ont persisté jusqu'à l'abolition de la féodalité par la
	monarchie du decennio francese 14 . C'est donc le constat d'un vide politique et administratif qui explique
	l'iconoclasme des insurgés, perçu à la fois comme la destruction des images du pouvoir et comme
	l'entérinement d'un état de fait. Cette double fonction explique que la chronologie des destructions
	recoupe celle des crispations de la monarchie constitutionnelle. Elle se concentre en effet sur la séquence
	avril-juillet 1848. Les actes d'iconoclasme connaissent une augmentation notable au lendemain des
	massacres du 15 mai 1848, alors que les pratiques répressives de la monarchie ont été portées au centre
	du débat national par les libéraux et ont constitué l'un des moteurs principaux de la résistance populaire
	dans les Calabres 15 .
	Les formes prises par la destruction des images sont renseignées par les sources judiciaires, au moins
	dans leurs grandes lignes, et s'inscrivent dans le répertoire de la violence ordinaire qui a prévalu dans
	les révoltes locales pendant l'époque moderne 16 . Le fait de jeter les statues à terre ou de les réduire en
	pièces en constitue la modalité la plus fréquente, toujours attestée comme une pratique collective et
	dirigée vers les deux statues du couple royal. Elle vise donc autant l'institution monarchique que la
	personne du roi. Les lacérations des drapeaux de la monarchie et des décrets royaux placardés dans
	l'espace public le confirment. D'autres pratiques interviennent de façon plus occasionnelle. À Villa San
	20 Archivio di Stato (Cosenza), Intendenza, Processi politici, I, b. 35 bis, f. 220.
	21 CADET, 2015, p. 206. 25 Archivio di Stato (Naples), Archivio Borbone, b. 1044, f. 38.
	22 Plusieurs mentions chez des propriétaires ou des prêtres, dans la province de Cosenza, à Rogliano (Archivio di Stato [Cosenza], Intendenza, Processi politici, I, b. 35 bis, f. 220), à Montalto (ibid., b. 88 bis, f. 517) ou à Cassano 26 HINCKER, 2008 ; FOURNIER, 2016.

Giovanni près de Reggio, l'archiprêtre Giovanni Zagarella pend à l'aide d'une écharpe tricolore le portrait du roi qui se trouvait alors dans l'église, en juin 1848

17

. Il arrive parfois que les formes de l'outrage se combinent, montrant l'acharnement des insurgés à l'encontre des représentations visuelles de la monarchie. À Luzzi dans la province de Cosenza, le 28 mars 1848, le juriste Alessandro Gatti brûle ainsi la statue du roi, qu'il a préalablement fracturée d'un coup de bâton 18 . Un autre procédé courant consiste à achever au moyen de coups de feu l'exécution de statues ou d'emblèmes déjà détruits, simulant ainsi l'exécution du roi. À Palmi dans la province de Reggio, quatre artisans de la commune

13

Archivio di Stato (Cosenza),

Intendenza, Processi politici, I, b. 32, f. 195. 14 FERRARI, 2016.

15 BUTTIGLIONE, inédite.

16 CECERE, 2013. 17

Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 212. tirent plus tôt, explique qu'elles aient circulé clandestinement et que les mentions dans les sources de police en soient rares. Dans la plupart des cas, l'armement est traditionnel, souvent très sommaire, et correspond, dans ses grandes lignes, à celui observé par Nicolas Cadet chez les insurgés qui défendaient les provinces calabraises contre les troupes napoléoniennes en 1806 21 . Le couteau est l'arme la plus fréquemment utilisée. Il est un ustensile d'usage quotidien dans les sociétés rurales calabraises, y compris chez des sujets qui appartiennent aux élites locales 22 ; il est par ailleurs associé à la figure du brigand calabrais, topos littéraire du romantisme italien 23 . L'usage d'un bâton est parfois mentionné, comme dans le cas d'Alessandro Gatti, avocat italo-albanais qui fracture les statues du couple royal à la fin mars 1848 24 . Parce qu'il s'inscrit dans les schémas traditionnels de la révolte, l'armement accrédite l'aspect spontané des actes iconoclastes, qui apparaissent le plus souvent improvisés et conjoncturels. L'usage des armes à feu est plus occasionnel, qu'il s'agisse des fusils à baïonnette, relevés dans 93 cas pour les deux provinces, ou des pistolets, utilisés pour exécuter des statues. Le recours à ces armes relève de logiques sociales distinctes : la très grande majorité des détenteurs de fusils, très courants dans les provinces calabraises, sont des paysans, alors que l'usage des pistolets est plus courant parmi les élites locales. C'est le cas des frères Antonio Maria et Domenico Marchiano, propriétaires du village italoalbanais de San Demetrio Corone où ils ont été parmi les principaux acteurs de l'insurrection, ou encore 19 Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 344. (Archivio di Stato [Naples], Archivio Borbone, b. 1044, f. 38), tous en juin 1848. du Le spectacle de la mise à mort du roi Parce qu'il s'effectue dans l'espace public, l'acte iconoclaste cherche à signifier auprès des masses la destitution symbolique du pouvoir en place et, indirectement, le transfert au peuple de la souveraineté dont il est détenteur. Il constitue en cela un procédé de démocratisation, qui passe par le meurtre symbolique du roi. En ce sens, des travaux récents l'ont montré pour d'autres espaces, il est producteur d'un langage politique, par le choix des cibles et par l'investissement de l'espace public, et ne saurait être réduit à un simple acte de vandalisme, malgré son impréparation matérielle apparente 29 . La répartition des lieux-cibles en est révélatrice, montrant la volonté de purger la totalité de l'espace public local des références à la monarchie. 27 Archivio di Stato (Naples), Archivio Borbone, b. 1044, f. 38. 28 BUTTIGLIONE, inédite, p. 64.

  formalisées. La destruction des images de la royauté, dès lors, réduit le programme politique des opposants à leur expression la plus simple, qui est précisément le meurtre du tyran.La profération de cris séditieux tyrannicides, parmi lesquels Morte al Tiranno et Abbasso al Borbone sont les plus couramment évoqués par la documentation, contribuent à expliciter et à justifier ces actes, présentés comme l'expression de la volonté du peuple. Punis par la loi du 5 mai 1827 qui porte sur les injures sacrilèges en public et en privé, ils ont donc une place centrale dans le processus de mise à mort symbolique du roi. Si l'accusation de tyrannie est la plus courante, certaines désignations font référence à la trahison du peuple napolitain, surtout après la répression des émeutes du 15 mai 1848 32 . Ailleurs, au mois de juillet, le roi est qualifié de « cornu, voleur et assassin » par référence à sa politique fiscale centralisée et au caractère répressif de sa pratique du pouvoir, qui en font l'une des figures emblématiques du malgoverno, notion utilisée dans la culture politique italienne depuis le Moyen Âge pour désigner un régime tyrannique associant la corruption à l'inefficacité et à la cruauté33 . À la a la charge, incite ses fidèles à tuer le roi et le général Nunziante et à manger leur chair 45 . À Stilo, le gendarme Carmine Anania déclare, en décapitant les deux statues du roi et de la reine, vouloir les hisser sur un piquet sur la place du Mercato à Naples, puis boire leur sang 46 . On perçoit la ritualisation du massacre, qui combine les pratiques de la violence ordinaire calabraise au répertoire des révoltes nationales, à travers l'investissement du principal théâtre des révoltes urbaines depuis l'époque moderne, qui a notamment vu éclater la révolte de Masaniello en 1647. Anania réutilise ici le procédé par lequel les soldats du roi ont ensuite exécuté Masaniello, incluant ce dernier, au même titre que les victimes calabraises des Bourbons, dans une communauté de la souffrance et du martyre, envisagée à l'échelle de la nation napolitaine 47 . La mise en scène du tyrannicide articule donc des pratiques hétérogènes, locales et nationales, qui posent la question du projet politique porté par les insurgés. Elles s'inscrivent en effet dans la tradition méridionale de revendication des autonomies locales, qui atteint son paroxysme avec les révolutions survenues à l'été 1848 dans les périphéries du royaume. 'expliquent à la fois par la volonté consciente d'épurer l'espace public des symboles de l'autorité centrale et par des réflexes spontanés qui témoignent de la haine ordinaire à l'égard des Bourbons.Dès lors, les premiers jours de juin 1848, qui ont constitué l'un des temps forts de la résistance calabraise, ont vu se multiplier les destructions publiques des statues par des acteurs ordinaires ayant pris les armes et se dirigeant vers les camps. À San Benedetto Ullano, Raffaele Conforti outrage les statues du couple royal avec quatre autres membres de la bande armée qu'il dirige, lors du départ pour le camp de Paola 51 .À la fin juillet 1848, le retour des camps s'est accompagné de destructions comparables, quoique de façon plus ponctuelle. À Montalto dans la province de Cosenza, l'emblème royal situé sur la façade d'une boutique privée est décroché et lacéré par des insurgés de retour du camp de Paola. L'acte s'accompagne de chants nocturnes en procession publique qui dénoncent les autorités royale et pontificale, dont les révoltés dénoncent la traîtrise, les promesses non-tenues envers les libéraux et le tournant conservateur 52 .

	Résistance locale et transferts de souveraineté
	32 Par exemple Luigi Panetta, propriétaire à Ardore (province de Reggio), le 4 juin 1848. Voir Archivio di Stato
	(Reggio Calabria), Atti di Polizia, I, b. 1, f. 141.
	33 Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 32, f. 194. 47 Sur la révolte de Masaniello en 1647, voir HUGON, 2011.
	34 Par exemple dans la province de Reggio, où les cris séditieux faisant référence à l'armée bourbonienne et à Nunziante sont nombreux lors des actes d'iconoclasme. Voir Archivio di Stato (Reggio Calabria), Atti di Polizia, 48 DELLA PERUTA, 1997.
	I, b. 1, f. 19. 49 Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 12, f. 64.
	35 Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 20.

31

Archivio di Stato (Cosenza), Intendenza, Processi politici,

I, b. 32, f. 196.

véritablement dénonciation du roi s'ajoute parfois celle de ses soutiens, comme le général Nunziante, principal acteur de la répression dans les Calabres, régulièrement présenté, à l'été 1848, comme l'un des instruments de la tyrannie bourbonienne et associé aux destructions des images de la royauté 34 . Ces appels au tyrannicide expriment donc des réactions spontanées face à la monarchie et témoignent de la haine ordinaire des Bourbons, qui justifient des appels publics à l'iconoclasme. C'est le cas d'un propriétaire de Reggio, Carmine Cuzzocrea, proche des chefs révolutionnaires Antonino et Agostino Plutino, qui aurait modifié le nom de son frère Ferdinando en Emilio « par haine du nom du roi » et aurait exhorté de jeunes gens de la province à détruire les symboles de la royauté

35

. Ce sont parfois les prêtres qui ont assumé ce rôle, intermédiaires décisifs de la révolution tout au long du premier XIX e siècle. Les sources en donnent plusieurs exemples, à l'image de Filippo Stilo, à Fiumara, ou de l'archiprêtre Antonio Giampaolo, à San Luca, auteurs de prêches incitant leurs paroissiens à détruire les statues du roi à l'été 1848

36

.

Ainsi construits et justifiés, les actes d'iconoclasme tyrannicide se présentent sous la forme d'exécutions simulées sur le modèle du jugement, qui recourent à des pratiques informelles de justice populaire 37 . Là 36 Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 81 et f. 252. Voir, pour la province de Cosenza, VALENTE, 1956. 37 Sur ces pratiques, courantes dans les sociétés sud-européennes de l'époque, voir DELIVRÉ, BERGER (dir.), 2014. il Les pratiques d'iconoclasme révolutionnaire s'inscrivent dans des mobilisations protestataires collectives et ritualisées, qui s'effectuent pour la plupart d'entre elles en bandes armées. Cette forme militaire, qui s'est développée dans les rangs des démocrates mazziniens, a trouvé un écho dans les structures villageoises calabraises 48 . Elle a rencontré des pratiques collectives du combat, celles des « troupes en masse » (truppe a massa) déjà mobilisées en 1799 et en 1806 avant d'être récupérées par les carbonari dans les années 1810. La province de Cosenza en offre les exemples les plus nombreux, comme à Figline, le 3 juillet 1848, où seize sujets du village détruisent publiquement les statues du roi et de la reine situées dans la mairie 49 . L'épisode montre comment l'acte iconoclaste s'inscrit dans un rituel plus large : il s'effectue au moment du départ des insurgés pour le camp de révoltés (campo rivoltoso) de Castrovillari 50 . Le départ pour ces structures dissidentes, où les insurgés ont cherché à organiser la résistance des provinces calabraises face à l'État centralisé, a été largement mis en scène auprès des populations locales et a recouru à des pratiques d'iconoclasme : les destructions intervenues 45 Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 102. 46 Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 187. 50 DELPU, 2019, p. 311-314. Les camps de révoltés (campi rivoltosi), apparus en Calabre à partir de juin 1848, ont permis d'organiser la résistance locale face aux troupes bourboniennes commandées par le général Nunziante, venues de Naples pour défendre l'ordre public. Cette orientation défensive explique leur mise en place aux entrées stratégiques de la province de Cosenza, par voie terrestre et maritime, afin de l'isoler du reste du royaume et en particulier de la capitale.

s

  De la tyrannie des Bourbons aux « rois du bas peuple », des transferts de sacralité incertains Par le répertoire d'action sur lequel il s'appuie et par son intégration dans le mouvement calabrais des revindiche, l'iconoclasme révolutionnaire est le support d'un projet politique qui consiste à transférer l'autorité du roi au peuple et celle de la nation à la communauté locale. Les cris séditieux, les actes des procès intentés aux insurgés après 1848 permettent d'en préciser le contenu. Les actes iconoclastes se sont en effet accompagnés d'appels fréquents à proclamer la république. Cette réclamation, qui avait perdu de l'importance chez les libéraux du royaume après l'écrasement traumatique de la République napolitaine, en 1799, par les armées des Bourbons, a été de plus en plus courante chez les opposants à la monarchie à partir des années 1830, sans être dotée d'un sens politique précis. À Cassano, le propriétaire Luigi Nigro dit devoir proclamer la république lorsqu'il fracture la statue du roi en 1848 et poignarde son effigie sur une monnaie 62 . Dans les trois provinces calabraises, ces revendications s'appuient sur des motivations diverses, qui vont de la déploration de la tyrannie bourbonienne à la revendication de l'abolition des taxes. Dans la province de Reggio, c'est le prétexte fiscal qui pousse Giuseppe Chianese, sujet de Sant'Ilario, à proclamer la république lorsqu'il fracture les statues du roi : il voit dans ce régime le seul capable d'abolir la fondiaria 63 . Plus que l'intérêt pour un gouvernement républicain donc, les appels à fonder ce régime se fondent d'abord sur des considérations fiscales, liées à la situation économique et sociale des sujets calabrais. Ces intentions, portées pour l'essentiel par des notables locaux, sont en effet dépourvues d'un contenu idéologique précis et correspondent plutôt à une posture d'opposition, destinée à porter au pouvoir des acteurs choisis par le peuple. Elle se situe donc à distance des républiques qui ont vu le jour dans les révolutions européennes de 1848, en France et à Rome, définies par une idéologie et des institutions. Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 205.initiatives n'aient été coordonnées ou que l'échelle à laquelle il fallait fonder la république soit précisée.Dans la très grande majorité des cas, elles sont le fait d'élites villageoises, occupant souvent des charges administratives ou judiciaires qui contribuent à leur légitimité auprès des populations locales. C'est le cas du substitut du juge Francesco Greco, auteur de la destruction des statues du roi à Cassano, suivi par dix hommes du village avec lesquels il proclame la république. Celle mise en place présente la particularité de confier le pouvoir à des notables locaux, comme les démocrates Carlo Mileti et Domenico Mauro ou encore Greco lui-même, intronisés rois par le peuple. Cette fonction, qui se retrouverait pour d'autres espaces périphériques du royaume des Deux-Siciles et notamment dans les provinces des Pouilles, correspond à une institution informelle, qui repose à la fois sur la notabilité, le charisme et le parcours politique64 . La diversité des désignations le confirme, des « petits rois » (piccoli re) aux « rois du bas peuple » (re del basso popolo), si ce n'est qu'elles décrivent l'origine démocratique de leur autorité. Dans la province de Cosenza, Francesco Greco a ainsi participé, aux côtés des insurgés calabrais, à la résistance locale contre l'occupation française en 1806, puis a participé à la révolution de 1820-1821 comme carbonaro 65 . Les nominations de piccoli re ont fréquemment accompagné la destruction des emblèmes royaux lors du départ pour les camps de révoltés. À Fuscaldo en juin 1848, un propriétaire de la commune, Giovanni Battista Carnevale, obtient ainsi le titre de roi dans sa commune, lors de la destruction d'un buste du roi Ferdinand II lors du départ pour le camp de Paola 66 . Saracena en juin 1848, sont détruits successivement les statues du couple royal placées dans la mairie, les blasons bourboniens affichés sur le poste de garde, dans la chancellerie communale et dans deux boutiques privées. Les insurgés appellent à tuer le tyran et à proclamer la république 67 . L'épuration des images de la monarchie passe donc par le contrôle des lieux de pouvoir communaux et par la diffusion des symboles de la révolution. Elle produit donc des contre-symboles, qui sont un marqueur de la transition politique 68 . Les arbres de la liberté, régulièrement évoqués dans les sources de police, en sont les plus visibles, sans qu'il soit possible de les quantifier précisément, mais l'espace insurgé se dote parfois d'images nouvelles de la souveraineté, vidées de leur contenu . C'est ce qui explique, le 11 avril 1848, la substitution d'un portrait du roi Charles-Albert de Piémont-Sardaigne à celui de Ferdinand II dans l'église de San Lorenzo del Vallo, au nord de la province de Cosenza 69 . L'acte, qui demeure isolé, n'en est pas moins chargé d'une valeur symbolique importante, parce qu'il propose de remplacer le souverain du royaume par un autre dirigeant italien. Il témoigne des sympathies italophiles qui ont alors circulé dans les Calabres, sans pour autant s'inscrire véritablement dans les mobilisations pour l'indépendance italienne 70 . Plus qu'un ralliement aux options majoritaires du Risorgimento, qui se construisent au même moment, l'acte exprime la volonté d'en finir avec « l'esclavage imposé par les tyrans » : Charles-Albert apparaît davantage comme un symbole d'émancipation que comme un souverain potentiel des Calabres intégrées à l'Italie.Les figures alternatives de la royauté sont donc mobilisées pour leur valeur symbolique, extraites de leur réalité idéologique et de leur contexte. L'horizon envisagé par les insurgés iconoclastes a été celui de la souveraineté locale. Les exhortations à fonder des républiques, qui ont accompagné la destruction des emblèmes de la royauté, ont davantage été perçues comme des appels à vaincre la tyrannie qu'à fonder un régime républicain. Leur éclatement géographique, alors qu'il n'y a pas eu de projet républicain global, y contribue. Par ailleurs, les actes des procès montrent que la notion de république n'a pas toujours été précisément comprise par ses promoteurs locaux : à Saracena dans la province de Cosenza, l'un des acteurs de la destruction des statues, Antonio Tedeschi, cordonnier de 36 ans, dit avoir contribué à la protestation collective sans savoir exactement ce que désignait le terme repubblica 71 . Cette méconnaissance est fréquente chez les insurgés des classes populaires, en particulier les petits commerçants et les artisans. Elle rappelle que l'iconoclasme s'est situé du côté des économies morales propres aux sociétés traditionnelles plutôt que dans les cadres de la politique moderne, alors inégalement consolidés dans la société du royaume.Pour la majorité des acteurs en effet, la destruction des statues du couple royal permet de répondre à un état de fait, la mort du roi et la vacance du pouvoir, connues par la circulation des rumeurs 72 . La structuration locale de la société calabraise a largement contribué à leur diffusion et à leur succès. Parce que les fausses nouvelles ont généralement annoncé la mort du roi ou sa destitution, l'iconoclasme a souvent permis, dès les derniers mois de l'année 1847, d'acter une transition politique en détruisant les effigies d'un pouvoir disparu. C'est la raison invoquée, parmi des cas très nombreux dans la province de Reggio, par le chancelier communal de Pellaro, puis par l'archiprêtre Zagarella de Villa San 69 Archivio di Stato (Napoli), Archivio Borbone, b. 1044, f. 35.70 Les chroniqueurs locaux du XIX e siècle ont rappelé la forte diffusion des sympathies italiennes dans les trois provinces calabraises, qui ont conduit à la mobilisation de volontaires calabrais dans le corps expéditionnaire napolitain parti en Lombardie-Vénétie à partir d'avril 1848. Voir par exempleANDREOTTI, 1874, vol. 3, p. 295.

	Les proclamations républicaines effectuées dans les Calabres ont relevé d'acteurs isolés, sans que leurs 64 LUCARELLI, 1953.
	65 Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 12, f. 66.
	61 Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 344.

62

Archivio di Stato (Napoli), Archivio

Borbone, b. 1044, f. 38. 63

Les destructions d'emblèmes du roi ont donc permis, symboliquement, la réhabilitation de l'autorité de figures notabiliaires locales, contre la tyrannie centralisatrice des Bourbons. Cette réhabilitation suppose l'investissement des lieux de pouvoir à travers des processions ritualisées, qui consacrent le rôle d'acteurs politiques locaux. Lors d'un cortège investissant le chef révolutionnaire Domenico Mauro dans ses fonctions de roi, à 66 Archivio di Stato (Napoli), Archivio Borbone, b. 1044, f. 38. 67 Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 37, f. 227. 68 Sur ce point, voir sur un autre contexte FUREIX, 2015.

idéologique71 Archivio di Stato (Cosenza), Intendenza, Processi politici, b. 37, f. 227. 72 CAGLIOTI, 2001.

On entend par politique moderne un ensemble de codes et de pratiques qui permettent l'affirmation d'un cadre politique unifié fondé sur la souveraineté nationale. Le modèle, défini par plusieurs historiens autour de Jean-Philippe Luis pour le cas espagnol (LUIS,

2005), s'applique notamment à la péninsule italienne et à la France méridionale, entre les révolutions de la fin du XVIII e siècle et les années 1850. 2 DELPU, 2019. 3 Sur celui-ci, voir PETRUSEWICZ, 2005. 4 Sur les formes chrétiennes de la violence, voir BUC, 2017. Sur les pratiques de la conflictualité dans les Calabres, voir CECERE, 2013, CADET, 2015 et ADDANTE, 2021. 5 FUREIX (dir.), 2014 et, sur l'importance et la fonction de cette pratique dans la France du XIX e siècle, FUREIX, 2019.Pour une approche comparative, voir FUREIX, ARISI ROTA (dir.), 2018. 6 ARISI ROTA, 2018 et, pour le cas des provinces des Pouilles, CALEFATI, inédite.

Chez les libéraux, voir ANDREOTTI, 1874 ; chez les contre-révolutionnaires, voir DE' SIVO, 1868. En français, le terme est attesté dès 1834 (FUREIX (dir.), 2014, p. 14). En italien, il conserve tout au long du XIX e siècle son sens traditionnel, qui le relie à l'histoire du Moyen Âge byzantin et aux formes de la dévotion chrétienne.

Celles-ci sont bien connues de l'historiographie spécialisée. Voir par exempleFIORE, inédite, et DI FIORE, 2013.

Voir, sur la méthodologie, MUCHEMBLED, 1987.

Archivio di Stato (Cosenza), Intendenza, Processi politici, I, b. 32, f.

Par exemple et parmi une riche production, l'Errico de Domenico Mauro, paru en 1834.

Archivio di Stato (Cosenza), Intendenza, Processi politici, I, b. 32, f. 196.

Voir surtout les travaux d'Emmanuel Fureix sur la France de la même époque[START_REF] Fureix | L'iconoclasme : une pratique politique ? (1814-1848)[END_REF].

Archivio di Stato (Reggio Calabria), Atti di Polizia, I, b. 1, f. 103.


Le Teuff

A Lucas
email: lucas@ipgp.fr

Possibly seismically triggered avalanches after the S1222a Marsquake and S1000a impact event

On May 4, 2022, a major martian seismic event was recorded

• We identified possibly seismically induced dust avalanches in the area of the estimated epicenter

• We discuss avalanche triggering conditions and derive a possible epicenter location based on avalanche spatial density

Introduction

On May 4, 2022, a major seismic event named S1222a (Kawamura et al., 2023) was recorded by the SEIS instrument [START_REF] Lognonné | SEIS: Insight's Seismic Experiment for Internal Structure of Mars[END_REF] of the InSight mission [START_REF] Banerdt | Initial results from the insight mission on mars[END_REF]. It was an unprecedented marsquake in the SEIS recording period with an estimated moment magnitude of M 𝑀 𝑎 𝑊 4.7 (InSight Marsquake Service, 2022). In comparison, 95% of events recorded by SEIS since landing in November 2018 had a magnitude below 3.5 [START_REF] Clinton | The marsquake catalogue from insight, sols 0-478[END_REF][START_REF] Böse | Magnitude Scales for Marsquakes Calibrated from InSight Data[END_REF][START_REF] Ceylan | The marsquake catalogue from insight, sols 0-1011[END_REF][START_REF] Knapmeyer | The global seismic moment rate of mars after event s1222a[END_REF]. As for some of the InSight events, a location was estimated with a back-Azimuth (bearing from the event toward InSight) of 101 • (96 • -112 • ) and an epicentral distance Δ = 37 • (±1.6 • ) which places the event epicenter at the location of 3.0 • S,171.9 • E (Kawamura et al., 2023) (green star on Fig. 1). Other nearby locations for the epicenter have also been proposed (Panning et al., 2023;[START_REF] Kim | Structure along the martian dichotomy constrained by rayleigh and love waves and their overtones[END_REF] (Fig. 1). No new impact crater has been reported that could be the source of this event (Fernando et al., 2023). The region shows many topographic features including a few tectonic structures expressed as north-south wrinkle ridges [START_REF] Knapmeyer | Working models for spatial distribution and level of Mars' seismicity[END_REF] and impact craters (Fig. 1). To the east of this region, the only major structures are Appollinaris Patera, a Noachian volcano [START_REF] Tanaka | Geologic map of mars[END_REF] about 200 km in diameter, and a large alluvial fan spanning southwards from the volcano's rim.

From orbital images, we identified dust avalanches (also known as slope streaks) in this region (orange symbols on Fig. 1). These are known active mass wasting processes occurring on Mars in several contexts [START_REF] Ferguson | Dark streaks on talus slopes, Mars[END_REF]Sullivan et al., 2001;Aharonson et al., 2003;[START_REF] Schorghofer | Slope streaks on Mars: Correlations with surface properties and the potential role of water: Slope streaks on Mars[END_REF][START_REF] Schorghofer | Three decades of slope streak activity on mars[END_REF][START_REF] Schorghofer | Sporadic formation of slope streaks on Mars[END_REF][START_REF] Gerstell | A distinct class of avalanche scars on Mars[END_REF][START_REF] Baratoux | The role of the wind-transported dust in slope streaks activity: Evidence from the HRSC data[END_REF][START_REF] Chuang | HiRISE observations of slope streaks on mars[END_REF]Bergonio et al., 2013;[START_REF] Heyer | Seasonal formation rates of martian slope streaks[END_REF][START_REF] Heyer | Dust devil triggering of slope streaks on mars[END_REF][START_REF] Valantinas | CaSSIS color and multi-angular observations of Martian slope streaks[END_REF]. They appear as relatively dark or bright streaks on steep dust-covered slopes and occur in regions with a high albedo and low to very low thermal inertia (Sullivan et al., 2001;Aharonson et al., 2003). Dust avalanches on Mars typically appear darker than the surrounding terrain. This is likely due to the removal of lighter-colored surface dust by the avalanches. When a slope streak is formed, loose dust and sand on the surface are mobilized and cascade down the slope, exposing the darker, underlying material [START_REF] Malin | Context Camera Investigation on board the Mars Reconnaissance Orbiter[END_REF][START_REF] Dundas | Geomorphological evidence for a dry dust avalanche origin of slope streaks on mars[END_REF]. This material may be darker due to several factors, such as the presence of iron-rich minerals or alteration by weathering processes (Christensen et al., 2001). In addition, the removal of surface dust by the avalanches may expose a rougher, more textured surface, which can scatter and absorb more light, making the streak appear even darker. Many studies discuss possible triggering conditions and emplacement mechanisms. Purely dry avalanches of fine dust have been explored from the perspective of both observations [START_REF] Schorghofer | Three decades of slope streak activity on mars[END_REF][START_REF] Phillips | Mass movement within a slope streak on mars[END_REF][START_REF] Dundas | Geomorphological evidence for a dry dust avalanche origin of slope streaks on mars[END_REF], and numerical simulations (Lucas, 2010). Spring discharge involving salty groundwater and/or brines in the shallow subsurface has been proposed [START_REF] Ferris | Dark slope streaks on mars: Are aqueous processes involved[END_REF][START_REF] Miyamoto | Fluid dynamical implications of anastomosing slope streaks on Mars[END_REF][START_REF] Head | Slope Streaks in the Antarctic Dry Valleys: Characteristics, Candidate Formation Mechanisms, and Implications for Slope Streak Formation in the Martian Environment[END_REF][START_REF] Kreslavsky | Slope streaks on Mars: A new "wet" mechanism[END_REF][START_REF] Bhardwaj | Martian slope streaks as plausible indicators of transient water activity[END_REF][START_REF] Bhardwaj | Are slope streaks indicative of global-scale aqueous processes on contemporary mars?[END_REF]. Other possible triggers include wind [START_REF] Baratoux | The role of the wind-transported dust in slope streaks activity: Evidence from the HRSC data[END_REF][START_REF] Heyer | Seasonal formation rates of martian slope streaks[END_REF], seismic activity from impacts or internal forces, or boulder track [START_REF] Chuang | HiRISE observations of slope streaks on mars[END_REF] have been proposed.

While previous studies looked at boulder falls and associated tracks triggered by possible paleo-seismic activity [START_REF] Roberts | Possible evidence of paleomarsquakes from fallen boulder populations, Cerberus Fossae, Mars[END_REF][START_REF] Brown | Possible evidence for variation in magnitude for marsquakes from fallen boulder populations, Grjota Valles, Mars[END_REF], no previous work could have directly tested the possibility of seismically induced mass wasting on Mars due to a lack of seismic event records before the InSight mission. In the framework of the recent seismic events S1000a and S1222a, we investigate the effects of the induced ground acceleration aftermaths as a potential triggering mechanism for dust avalanches in the vicinity of the located epicenter. To do so, we conduct regional mapping of the avalanches from pre-event and post-event imagery in order to estimate the effect of the marsquake and impact crater on the rate of avalanches. We take into account possible biases due to the limited number of images, the time span between images, the sub-surface properties through thermal behavior, and the various sensitivities of each camera sensor.

Methods

Orbital data and mapping

As soon as the S1222a event was detected by SEIS and an estimate of the epicenter location was provided, we investigated orbital observations provided by the Context (CTX)

and High Resolution Imaging Science Experiment (HiRISE) cameras [START_REF] Malin | Context Camera Investigation on board the Mars Reconnaissance Orbiter[END_REF][START_REF] Mcewen | Mars Reconnaissance Orbiter's High Resolution Imaging Science Experiment (HiRISE)[END_REF], both on board the Mars Reconnaissance Orbiter (MRO). Along with MRO imagery, we examine images from the Mars Global Surveyor (MGS)/Mars Orbiter Camera (MOC) [START_REF] Malin | Mars observer camera[END_REF], and THEMIS-Vis/Odyssey [START_REF] Fergason | High-resolution thermal inertia derived from the Thermal Emission Imaging System (THEMIS): Thermal model and applications[END_REF].

This led to a set of hundreds of images acquired before the seismic event. In addition, we requested new MRO observation over areas where we mapped avalanches inside the uncertainty area (Kawamura et al., 2023) (Fig. 1, Supp. Info text S1). At the time of writing this paper, a dozen HiRISE images and thirty new CTX observations were obtained, all acquired after the S1222a seismic event. In addition to imagery, we used Digital Terrain Models (DTMs) from both Mars Orbiter Laser Altimeter (MOLA, Smith et al. (2001)) and High

Resolution Stereo Camera (HRSC, Neukum and Jaumann (2004)), the geological map from [START_REF] Tanaka | Geologic map of mars[END_REF] and the thermal inertia map [START_REF] Christensen | The thermal emission imaging system (themis) for the mars 2001 odyssey mission[END_REF]) (See Supp. Info.

Text S2), which all provide contextual information. Note that Figure 1 specifically depicts the area covered by post-S1222a event CTX imagery. All the data have been combined into a Geographical Information System (GIS) in order to manually map all avalanches in the region of interest (Fig. 1), by two independent people (see Supp. Info. Text S1 for details on the imagery processing and mapping). The older observations, provided by both MOC and THEMIS-Vis, were only used for confirming the very low fading rate (Sullivan et al., 2001), being in good agreement with the dust activity reported in this region (Battalio & Wang, 2021).

Estimates of avalanche rate and statistics

Avalanche rate 𝑞 is obtained from equation provided in Aharonson et al. (2003):

𝑞 = 100 × Δ𝑛 𝑛Δt , ( 1
)
where 𝑛 is the total number of avalanches observed in both the two overlapping images, Δ𝑛 being the newly observed avalanches on the recent image and not in the older image, and Δt being the time span between the two observations in Martian years. This rate 𝑞 is expressed in % of new events/Martian year (Aharonson et al., 2003). This method has also been used by recent work [START_REF] Heyer | Seasonal formation rates of martian slope streaks[END_REF]. The time periods between overlapping images in our database range from ∼0.3 to almost 7 martian years.

Finally we agglomerate avalanches in the same location (i.e. crater) and hence to compute the avalanche rate in each area where new events can be observed between two overlapping images. As opposed to a squared binning, hexagons are more similar to circles, hence they better translate data aggregation around the bin center. As most areas covered by avalanches in this region are impact craters, this provides a more valuable way to decipher the avalanche coverage.

Results and discussion

Avalanches triggered after S1000a impact event

Before discussing S1222a event, we investigated S1000a impact event which occurred on September 18 2021, and left a crater over 150 m in diameter at 38.1 • N;-79.87 • E (Fig. 2).

This event was recorded by SEIS and then orbital imagery revealed its actual location. Its magnitude was estimated to be around M 𝑀 𝑎 𝑤 4.1, hence about 25 times smaller that S1222a in energy [START_REF] Ceylan | The marsquake catalogue from insight, sols 0-1011[END_REF][START_REF] Posiolova | Largest recent impact craters on mars: Orbital imaging and surface seismic co-investigation[END_REF]. It should be noted that the location estimate from the seismic signal analysis was felt at approximately 130 km away from the actual location, as discussed in [START_REF] Posiolova | Largest recent impact craters on mars: Orbital imaging and surface seismic co-investigation[END_REF]. By analysing all pre-event images including CTX, and HRSC and post-event HiRISE images, we could map a very large number of avalanches not seen in pre-event imagery. By looking back in time using all available images, including MOC/MGS, we observed that these areas were not covered by dust avalanches prior to the impact event (Fig. 2.).

We looked at the the density distribution of the new avalanches (as seen on the postevent images and having the same brightness, hence the same age) as a function of their respective distance to the impact crater (histogram inset in Fig. 2). This distribution follows a bell-shaped curve. As seen on Earth, seismically triggered mass-wasting is absent very close to the epicenter, and increases at farther distances until it decreases again at the farthest distances (e.g., Tatard, 2010;Livio & Ferrario, 2020). Nonetheless, the mechanism here is different. It is very likely that the avalanches are triggered by secondary impacts, and not seismic waves. As an example of a typical scenario, ejecta leaving the primary impact at a velocity v = 200 m.s -1 , with a launch angle of 𝜃 = 45 shows that an impact with a seismic magnitude M 𝑀 𝑎 𝑤 4.1 can trigger a very large number of avalanches on Mars. As such, this is likely be discussed more thoroughly in a following work which would evaluate the ballistic recomposition in order to evaluate potential effects of secondary impacts on the dust avalanche triggering. However, the ground accelerations caused by a surface impact and a deep earthquake are not the same. So, in view of our results for the S1000a event, we discuss our results for S1222a in the following sections.

Avalanche rate increase in post-marsquake S1222a images

We analyzed all image pairs over the whole area of interest near the S1222a estimated epicenter. We identified 4532 avalanches (orange symbols in Figure 1). More than 200 avalanches were identified on pre-event images (over the 2005-2021 period), and 122 were identified on the post-event CTX images with respect to their 2005-2021 period counterparts respectively.

An example is given in Figure 3-a. Note that, while doubtful avalanches may have been detected (e.g., yellow symbols in Fig. 3-a), we only took into account the robust observations of new avalanches (e.g., red symbols in Fig. 3-a). For the statistical robustness, we then derived avalanche rates 𝑞 for each CTX/CTX pair only. When times series were available, we derived avalanche rate chronicles (Fig. 3-b). As exemplified on Fig. 3-b, a strong increase of 𝑞 is observed after the S1222a event. Indeed, over the whole area of interest (Fig. 1), the pre-event rates (circles in fig. 3-c) lie around 2.6%.MYear -1 with a maximum value of 6%.MYear -1 , accounting for uncertainties following Aharonson et al. (2003). These values are in agreement with in previous estimated by Aharonson et al. (2003), and avalanche rates do not differ substantially across the region covered by our study. In contrast, post-event values of 𝑞 show a significantly different distribution both spatially and in amplitude (Figure 3c,d). While most rates still fall below 10%, we observe that in 9 places, the rates are >10%, as high as 40% (excluding outlier, Figure 3-d). If we keep only the sub-10% values, the average is the same as that before the seismic event (2.6%.MYear -1 ), and there is also no dependence on the epicentral distance. Interestingly, the highest post-event 𝑞 (>20%) are found at the smallest distances from the epicenter of the S1222a event proposed by Kawamura et al. (2023). When relating the derived avalanche rate 𝑞 to the epicentral distance Δ with respect to the estimated location from Kawamura et al. (2023), we obtained a slight decreasing trend of 𝑞 with Δ. Finally we also verified that temporal sampling of the orbital images (Δt) does not bias the avalanche rate estimates (Fig. 3-e).

To address the limited number of observations, we employed a permutation test, also known as bootstrapping [START_REF] Efron | An introduction to the bootstrap[END_REF][START_REF] Davison | Bootstrap methods and their application[END_REF]. This nonparametric approach does not rely on specific distribution assumptions about the data. We began by calculating the avalanche rate for each CTX/CTX pair for both pre-event and postevent observations, determining the mean difference as our observed statistic. Then, we merged the pre-event and post-event rates, treating them as a combined dataset without distinction of their original times. This pooled data was randomly shuffled to create new groups, preserving the original group sizes. We calculated the permuted test statistic by assessing the avalanche rate in this permuted data. This permutation process was iterated a million times, generating a distribution of test statistics under the hypothesis of no marsquake influence.

Comparing our observed statistic to the 95% confidence interval derived from bootstrapping, we found that the post-event avalanche rates in all CTX observations exceeded 95% of the bootstrap statistic distribution. This indicates a significant increase in avalanche activity following the seismic event. However, it is important to note that the area studied includes locations possibly too distant from the epicenter to be affected by the marsquake. Focusing on rates exceeding 6%.MYear -1 , the post-marsquake rates surpassed the 99.98% confidence level. These findings, along with the detailed bootstrapping procedure, are outlined in Algorithm 1 and illustrated in Figure 4.

Effect of the relative thermal inertia

Subsurface properties at shallow depths can be analyzed through thermal inertia, which et al. (2004). Due to significant variations between orbits, we calculate the ratio of the apparent value of thermal inertia at avalanche scar location with respect to the median value on the surrounding plains, and named hereafter Γ * = Γ 𝑎𝑣𝑎𝑙𝑎𝑛𝑐ℎ𝑒 /Γ 𝑝𝑙𝑎𝑖𝑛 (see Supp. Info.

Text S2). By examining Γ * , we found that areas with the lowest values experience the most significant increases in avalanche rates (see Fig. 3-c). Specifically, when Γ * >> 1.5, postevent avalanche rates do not exceed pre-event rates. Conversely, an increase in 𝑞 is observed when Γ * < 1.5 This leads us to conclude that post-event avalanche susceptibility on Mars is primarily influenced by scarp locations with steep slopes and the lowest apparent thermal inertia. Such conditions correspond to the most unconsolidated terrains or areas with fine granular material.

Epicentral distance and possible sources of the quake

Although the epicentral distance is far from being the only parameter that controls the avalanche rates, it remains an important control factor (Tatard, 2010;Livio & Ferrario, 2020) (Supp. Info. Text S3 and S5). The reason is that the transition between a static state and a flowing state is modelled by introducing a threshold allowing the material to flow. This has been shown to quantitatively capture debris and rock avalanche morphodynamics on Mars (Lucas, 2010;Lucas & Mangeney, 2007;Lucas et al., 2011Lucas et al., , 2014) (see Supp. Info. Text S3).

Nonetheless, local geology, fractures, aftershocks and historical events will have a significant effect on the aftermaths of an earthquake by leading the slopes close to failure (Tatard, 2010;Livio & Ferrario, 2020;[START_REF] Chen | Distribution pattern of coseismic landslides triggered by the 2017 jiuzhaigou ms 7.0 earthquake of china: Control of seismic landslide susceptibility[END_REF]Rosser et al., 2021;[START_REF] Lombardo | From scenario-based seismic hazard to scenario-based landslide hazard: fast-forwarding to the future via statistical simulations[END_REF]. Taking into account all these considerations, the rate would not be expected to be controlled only by epicentral distance. However, our constraints on the characteristics of the marsquake are weak, especially in terms of depth, focal mechanisms, and therefore the resulting ground acceleration. Our knowledge on the geological heterogeneity is also poorly constrained. Also, compared to terrestrial standards, this marsquake remains a small event.

Nonetheless, small seismic events have shown to significantly increase the rate of landslides on Earth [START_REF] Martino | Increase in landslide activity after a low-magnitude earthquake as inferred from dinsar interferometry[END_REF]. Indeed, recent studies show that even very small amplitude seismicity may trigger instabilities on metastable slopes [START_REF] Bontemps | Rain and small earthquakes maintain a slow-moving landslide in a persistent critical state[END_REF]Durand et al., Minor revision).

Nonetheless, under the hypothesis that event S1222a did trigger avalanches, we considered the empirical model proposed by Livio and Ferrario (2020) which relates the distribution of triggered avalanches N 𝑎𝑣𝑎 with the epicentral distance Δ:

𝐺 (m) = 𝑁 𝑎𝑣𝑎 = 𝑎 exp - Δ -𝑏 𝑐 2 , ( 2
)
where a is the amplitude of the distribution, b the distance of the peak amplitude and c the width of the distribution. While we do not have images just before and just after the event,

we derived an estimation of the number of triggered avalanches from this relationship:

𝑁 𝑎𝑣𝑎 = Δ𝑛 -q × 𝑛Δ𝑡/100, ( 3
)
where q is the long-term avalanche rate (i.e., we conservatively considered 6%.MYear -1 ).

Because the avalanche susceptibility is not evenly distributed (i.e., steep slopes only located inside impact craters, non-homogeneous surface/sub-surface properties), we only consider observations that meet the following criteria: Γ * < 1.93, and Δt < 1.5 MYear, to only account for the lowest thermal inertia (see Fig. 3) and the smallest time span between images to reduce biases. Then, we used a Monte Carlo method to invert the most probable epicenter location using a maximum likelihood function with a Laplacian distribution of errors (Mosegaard & Tarantola, 1995) (See Supp. Info. Text S4). The resulting probability distribution of the epicenter under all of these considerations is given in Figure 5. It is situated in between the locations obtained from both body and surface waves analysis respectively (Kawamura et al., 2023;[START_REF] Panning | Locating the largest event observed on mars with multi-orbit surface waves[END_REF][START_REF] Kim | Structure along the martian dichotomy constrained by rayleigh and love waves and their overtones[END_REF], then included in the uncertainty ellipses of epicentral locations (green contours in Fig. 5).

This distribution can lead us to two different interpretations regarding the source mechanism of the quake, mainly related to internal tectonic activity. A first hypothesis would be based on the fact that our distribution is slightly shifted toward the East from the wrinkle ridges, on the flanks of Apollinaris Patera. It is now well supported that Mars still hosts remnant volcano-tectonic activity, especially along Cerberus Fossae [START_REF] Giardini | The Seismicity of Mars[END_REF][START_REF] Horvath | Evidence for geologically recent explosive volcanism in Elysium Planitia[END_REF][START_REF] Perrin | Geometry and segmentation of cerberus fossae, mars: Implications for marsquake properties[END_REF][START_REF] Stähler | Tectonics of Cerberus Fossae unveiled by marsquakes[END_REF], possibly due to the presence of a plume [START_REF] Broquet | Geophysical evidence for an active mantle plume underneath Elysium Planitia on Mars[END_REF], and associated with normal slip motion (Brinkman et al., 2021;[START_REF] Jacob | Seismic sources of insight marsquakes and seismotectonic context of elysium planitia, mars[END_REF]. While the moment tensor analysis of the S1222a event can give very different slip motions, NNW-SSE normal faulting is a possible solution [START_REF] Maguire | Moment Tensor Estimation of Event S1222a and Implications for Tectonics Near the Dichotomy Boundary in Southern ELysium Planitia Mars[END_REF], highlighting a possible activity of Apollinaris Patera at depth. However, unlike Cerberus Fossae, Apollinaris Patera is an old Noachian volcano, thus it seems unlikely that remnant volcanic activity would be present at shallow depth. A second hypothesis would be related to the 450 km long wrinkle ridge, trending NNE-SSW, and cross-cutting the Hesperian terrains between the two epicentral locations (black lines in figure 5). The probability distribution of the epicenter inferred from the avalanche rate is about 30 to 60 km East of this major structure. The shape of the topographic profile across the ridge is an asymmetric arch-ridge, with a steep slope facing West and a shallow slope facing East (Fig. 1), which would imply a main East-dipping thrust at depth (Andrews-Hanna, 2020). Assuming a fault dip of 34 • to 42 • for arch-ridges (Andrews-Hanna, 2020), a probability distribution situated about 30 to 60 km East of the wrinkle ridge would lead to a hypocentral depth ranging from 20 to 54 km. This range of depth is in agreement with the best solutions found by [START_REF] Maguire | Moment Tensor Estimation of Event S1222a and Implications for Tectonics Near the Dichotomy Boundary in Southern ELysium Planitia Mars[END_REF]. Along with another study by [START_REF] Brinkman | Applications of time-frequency domain polarization filtering to insight seismic data[END_REF], their preferred solutions present mainly reverse slip motions striking E-W to NW-SE, which is not optimally oriented with the overall wrinkle ridge observed from orbital imagery. However, local large variations in fault strikes are possible along a wrinkle ridge. Note that the wrinkle ridges are cross-cutting a large E-W bulge situated at about -5 • latitude, connecting the flanks of Apollinaris Patera and a large crater in the west (Fig. 1). This bulge presents hundreds of meters of difference in elevation and slight apparent thermal inertia anomalies that could indicate a bedrock affected by an old tectonic structure. Interestingly, the bulge's azimuth is aligned with our probability distribution of the epicenter. More work would be needed to understand the origin of this structure and a possible link with the source of the marsquake.

It should also be noted that source locations obtained from other methods such as surface waves or coda characteristics give different locations [START_REF] Kim | Structure along the martian dichotomy constrained by rayleigh and love waves and their overtones[END_REF][START_REF] Panning | Locating the largest event observed on mars with multi-orbit surface waves[END_REF][START_REF] Menina | Stratification of heterogeneity in the lithosphere of mars from envelope modeling of event s1222a and near impacts: Interpretation and implications for very-highfrequency events[END_REF]. Both studies using surface waves predict source locations more towards the south as shown in 5. This is due to different back azimuth they obtained for surface waves compared from that described in (Kawamura et al., 2023) using body waves.

Panning et al. ( 2023) also discusses the possibility that the source location could be in the southern hemisphere. Interestingly, [START_REF] Menina | Stratification of heterogeneity in the lithosphere of mars from envelope modeling of event s1222a and near impacts: Interpretation and implications for very-highfrequency events[END_REF] conclude that they need a thick ( 60km) diffusive layer to explain the coda shape of S1222a. This could imply that either the source location could be in the highlands of the southern hemisphere [START_REF] Wieczorek | Insight constraints on the global character of the martian crust[END_REF], or that thermal anomalies at depth are present in the Appollinaris area.

Our work leads us to propose that the source of the quake is likely due to thermal contraction due to Mars' cooling through time. The peak of thermal contraction and wrinkle ridge formation occurred during the early Hesperian and decreased progressively until now [START_REF] Watters | Compressional tectonism on Mars[END_REF]. Even if the wrinkle ridge in figure 5 is well expressed in morphology, its surface trace ends to the north, near the transition between Hesperian and Amazonian terrains [START_REF] Tanaka | Geologic map of mars[END_REF]. This indicates that the ridge has not been active in recent times.

However, thermal contraction is still ongoing on Mars and might re-activate local mechanical weaknesses in the martian crust, such as wrinkle ridges, over larger recurrence time periods.

If such activity is real, microseismicity should be associated with it.

Conclusions

In our comprehensive study of surface features surrounding the S1000a and S1222a seismic events on Mars, we utilized MRO orbital data to assess the associated avalanche rates. Our findings reveal a substantial increase in avalanches following the S1000a impact event, suggesting its indirect aftermaths, likely via secondary impacts. The S1222a event presented a more complex scenario, necessitating thorough investigation. We established pre-event avalanche rates in line with global estimates from Aharonson et al. (2003) and those near Olympus Mons obtained by [START_REF] Heyer | Seasonal formation rates of martian slope streaks[END_REF], ranging between 1 and 6%.MYears -1 .

These rates, when compared to post-event rates of up to 40%.MYear -1 near the estimated epicenter (Kawamura et al., 2023;[START_REF] Panning | Locating the largest event observed on mars with multi-orbit surface waves[END_REF][START_REF] Kim | Structure along the martian dichotomy constrained by rayleigh and love waves and their overtones[END_REF], underscore a significant increase in areas of lower apparent thermal inertia. This leads us to propose that the S1222a marsquake could be the driving factor behind the observed increase in avalanches.

This analysis also enabled us to estimate a probable epicenter for the marsquake, considering the Γ * threshold of 1.5 and a radial ground acceleration. This inferred location is intriguingly situated near a volcanic edifice and a North-South wrinkle ridge, highlighting the geological complexity of the region. Our study not only might suggest that current seismic activity on Mars can initiate mass wasting processes like dust avalanches but also opens avenues for exploring regions with observed avalanches and other seismic events detected by the InSight mission. The increased rates of avalanches in areas with historical seismic sources suggest that ground deformation plays a pivotal role in these phenomena. This methodology can be invaluable in future seismic event analyses, where visible aftermaths such as avalanches can offer significant insights into epicenter locations. Overall, our findings demonstrate that avalanches on Mars serve as a crucial tool for documenting rapid processes, from discrete surface perturbations like impacts to more continuous events like quakes. This understanding significantly enhances our ability to study and interpret the dynamic surface and subsurface processes of Mars.

Acknowledgments

Authors Figure 1. Regional map of dust avalanches near the S1222a event estimated location (green star with associated green contours, (Kawamura et al., 2023). The white star is the location estimated by multi-orbit surface waves [START_REF] Panning | Locating the largest event observed on mars with multi-orbit surface waves[END_REF]. The yellow star shows the estimated location according to surface waves (Kim et al., 2022). Orange symbols are all avalanches mapped. Red symbols show where avalanches are observed on post-event images. Basemap is the MOLA elevation map (Smith et al., 2001). Symbols show number of avalanches observed on post-event images. Color scales with the rate 𝑞. The green star (upper center) is the maximum peak of the estimated epicenter and its uncertainty ellipses (green contours) obtained from body waves (Kawamura et al., 2023), the white star, the location estimated from multi-orbit surface waves [START_REF] Panning | Locating the largest event observed on mars with multi-orbit surface waves[END_REF], and the yellow star is the estimated epicenter derived from the surface waves [START_REF] Kim | Structure along the martian dichotomy constrained by rayleigh and love waves and their overtones[END_REF]. Black lines are detailed surface traces of the main wrinkle ridge in the vicinity of the epicentral area. Background map is the thermal inertia Γ. (Inset) Expected avalanche density distribution with confidence interval from Monte Carlo inversion using equation 2 with respect to the number of avalanche 𝑁 𝑎𝑣𝑎 derived from equation 3. Rosser, N., Kincey, M., Oven, K., Densmore, A., Robinson, T., Pujara, D. S., . . . Dhital, M. R. ( 2021). Changing significance of landslide hazard and risk after the 2015 mw 7.8 gorkha, nepal earthquake. Progress in Disaster Science, 10, 100159. Schorghofer, N., [START_REF] Schorghofer | Three decades of slope streak activity on mars[END_REF]. Three decades of slope As well known, the brightness of these streaks is usually inversely proportional to their age. The darker they are, the more recently they formed (Sullivan et al.;2001). Note that we observe a very slow fading rate in this area, as some image pairs spanning 18 Earth years still show the same avalanches. Such a low fading rate allows us to estimate the long-term avalanche formation rate without concern that ongoing fading would significantly affect our measured rates (Aharonson et al.;2003). Additionally, this slow fading rate is in good agreement with the dust activity reported in this region (Battalio and Wang;[START_REF] Banerdt | First focal mechanisms of marsquakes[END_REF].

For the sake of comparison, the fading rate of dust devil tracks at the InSight landing site was also measured. There, dust devil tracks fade in a few months to almost two terrestrial years (Perrin et al.;2020). Notwithstanding, these are not formed by the same processes and thus might not fade at the same rates, in particular since avalanches remove a thicker layer of dust cover compared to dust devils. Other types of albedo features fade at different rates on Mars. For example, blast zones around new impacts have a median fading lifetime of 8 Mars years [START_REF] Daubar | Changes in blast zone albedo patterns around new martian impact craters[END_REF][START_REF] Daubar | Changes in blast zone albedo patterns around new martian impact craters[END_REF], while dust devil tracks and rover tracks fade much more quickly as they remove dust more superficially [START_REF] Balme | Mars: Dust devil track survey in argyre planitia and hellas basin[END_REF]2003;[START_REF] Verba | Observations from the high resolution imaging science experiment (hirise): Martian dust devils in gusev and russell craters[END_REF]2010;[START_REF] Geissler | Gone with the wind: Eolian erasure of the mars rover tracks[END_REF]2010). To compare the same type of feature as studied here, other slope streaks have shown gradual fading over approximately 20 Mars years (e.g., Schorghofer et al.;2007;Bergonio et al.;2013). Avalanche motion treated as a single-phase dry granular flow with Coulomb-type behaviour can be described by a two-dimensional Saint-Venant system. Hence, the transition between a static state and a flowing state is modelled by introducing a threshold allowing the material to flow. The friction forces f (x, y, t) must satisfy ||f (x, y, t)|| ≤ µρgh,

where µ is the friction coefficient, ρ the density, g acceleration due to the gravity, and h the thickness. This has been shown to quantitatively capture debris and rock avalanche morphodynamics on Mars (Lucas;2010;Lucas and Mangeney;2007;Lucas et al.;2011, 2014). Hence, this transition is sensitive to ground acceleration or dust deposition. To evaluate the stability of a slope composed of Martian dust, taking into account the potential cohesion of dust due to electrostatic forces, one can use the Coulomb failure criterion

τ = c + σ n tan(ϕ), (2)
where: τ is the shear stress, c is the cohesion of the material, σ is the normal stress, and ϕ is the angle of internal friction (where µ = tan(ϕ)). The normal stress on a slope is calculated as:

σ n = ρgh cos 2 (θ), (3)
where ρ is the density of the material, g is the acceleration due to the gravity, h is the depth of the material, and θ is the slope angle.

• Triggering by dust deposition

The additional shear stress δτ due to an additional height δh of deposition is estimated by: δτ = ρgδh sin(θ) cos(θ) (4)

By rearranging the aforementioned equations, the additional deposit height necessary for slope failure, denoted as δh, can be calculated. One can perform a parametric exploration accounting for the following parameter ranges estimated from soil mechanics experiments recently performed with HP 3 on InSight Spohn et al. (2022a,b): initial dust thickness h = {10 cm, 2m}, slope angle θ = 30 • , friction angle ϕ ={12 • ,25 • ,35 • } and cohesion c = {0,12} (kPa). The additional dust height δh required to trigger failure ranges from millimetric to metric scales as illustrated in the Fig. S4. These ranges can be compared to the typical dust deposition rates observed on the ground from both rovers and landers. As a matter of fact, dust deposition of a few particles radius per 100 sols are observed at the MER rover locations [START_REF] Kinch | Dust deposition on the mars exploration rover panoramic camera (pancam) calibration targets[END_REF]2007). By considering the 140 sols between the closest image pair (Ls = {80;220}) and the typical grain size radii discussed by [START_REF] Kinch | Dust deposition on the mars exploration rover panoramic camera (pancam) calibration targets[END_REF], would lead to 140 sols x 3 radii/100 sols x 1.5µm = 6.30µm, which is 3 folds smaller than the lower bound discussed previously. Accounting for larger particle radii (r=60µm) as it was assessed around the InSight landing site [START_REF] Chen-Chen | Characterisation of deposited dust particles on mars insight lander instrument context camera (icc) lens[END_REF][START_REF] Banerdt | A Tectonic Origin for the Largest Marsquake Observed by InSight[END_REF], the dust accumulation would be around 250 µm over the 140 sols period, which is comparable to the dust settling rate which was recorded to be 1µm/sols at Phoenix landing site [START_REF] Drube | Magnetic and optical properties of airborne dust and settling rates of dust at the phoenix landing site[END_REF]2010). Still, such extreme values of dust rate would be insufficient.

5

The dust rate at the S1222a location can be estimated from GCMs results, by dividing the total dust deposit over the 140 sols period by the density of the material. Considering a of 1600 kg.m -3 and a typical mean dust deposition of 0.5 × 10. -9 kg.m -2 .s -1 (Figure S5), the deposit accumulated is of about 7 µm, which in agreement to MER's location observations as shown in [START_REF] Kinch | Dust deposition on the mars exploration rover panoramic camera (pancam) calibration targets[END_REF]. Additionally, pre-event images display wind streaks oriented northward, associated with small impact craters. However, post-event images do not show any changes, casting doubt on the wind's ability to transport dust within the observed time span (Fig S6). Predictions from GCM indicate surface wind speed below 3 m.s -1 if we except slope winds over Apollinaris Mons (Figure S7). Finally, the orientation of the slopes on which avalanches are observed, whether in images before or after event S1222a, show the same distribution as that of steep slopes (>10 • ), which allows us to rule out, once again, any connection with the action of the wind (Fig. S8).

9

• Triggering by atmospheric pressure variation In order to evaluate the atmospheric pressure variation required to initiate slope failure, we need to focus on the change in shear stress. The pressure variation, in this context, is essentially the increase in shear stress needed to reach the failure point. To find the pressure variation required for failure, one needs to calculate the difference between the current shear stress due to the weight of the material and the shear strength (given by the Coulomb failure criterion). This difference represents the additional shear stress required to reach the failure condition. For the specified parameters (cohesion of 1000 Pa, friction angle of 30 degrees, density of 1600 kg/m 3 , height of 1 meter, slope angle of 30 degrees, and Martian gravity), the calculated pressure variation (additional shear stress) required to initiate slope failure is approximately 1000 Pascals. In comparison, the typical pressure variations due to turbulence observed at the InSight landing site are of about a few 10's of Pa, which is two folds smaller than the required pressure variations.

• Triggering by ground acceleration

To evaluate how ground acceleration from seismic activity, could contribute to reaching the slope failure condition, we need to consider the additional shear stress induced by this acceleration. One can approach this by considering the pseudostatic analysis, where the seismic force is simplified as a static force acting on the slope. The additional shear stress (δτ ) due to ground acceleration (a) can be estimated using the formula:

δτ = ρah sin(θ), ( 5
)
where a is the horizontal ground acceleration. The total shear stress acting on a slope is the cumulative effect of the inherent shear stress, which is a consequence of gravitational forces, and the additional shear stress induced by ground acceleration. A slope is prone to failure when this combined shear stress surpasses the material's shear strength, as defined by the Coulomb failure criterion. Continuing with the same parameters as before -cohesion (c) at 1000 Pa, internal friction angle at 30 degrees, density at 1600 kg.m -3 , and a slope angle of 30 degrees -we consider a hypothetical ground acceleration of 0.5 m.s -2 , which is 15% of Martian gravity (A typical M w 5 on Earth leads to 10% of terrestrial gravity, (e.g., [START_REF] Sandron | Seismic parameter design assessment for the kribi deep seaport in cameroon[END_REF][START_REF] Andrews-Hanna | The tectonic architecture of wrinkle ridges on Mars[END_REF]). Under these conditions, an initial material thickness of 3 meters results in the total shear stress beginning to surpass the shear strength. This indicates a condition likely to lead to slope failure.

A key factor influencing the outcome in terms of both total shear stress and shear strength is the initial height of the material on the slope. Within our model's framework, the gravitational shear stress and the extra shear stress from ground acceleration are both directly proportional to the material's height. Consequently, an increase in height leads to a corresponding rise in both types of stress. However, the rates at which they increase may vary. This variance plays a crucial role in determining how changes in height influence the probability of slope failure. This aspect can be related to the fact the thermal properties, as shown in the main text, seems to control the increase of the avalanche rate.

The ground acceleration associated with an earthquake is generally evaluated through the PGA (Peak Ground Acceleration), the maximum amplitude recorded during an earthquake. On the zeroth-order, this value will depend on the epicentral distance. Following [START_REF] Gomberg | Earthquake dynamic triggering and ground motion scaling[END_REF]; [START_REF] Gomberg | A model of earthquake triggering probabilities and application to dynamic deformations constrained by ground motion observations[END_REF], we can relate the ground motion G with epicentral distance ∆ through this simple relation:

P GA(∆) = K (α + ∆) 2 , (6) 10
where K and α are empirical parameters (Gomberg and Felzer;[START_REF] Gomberg | A model of earthquake triggering probabilities and application to dynamic deformations constrained by ground motion observations[END_REF]. Nonetheless, as on Earth, avalanches or rockfalls occur in response to a combination and cumulative forcings acting at different time scales including gravity, weathering and chemical effects, thermal, meteorological, seismic or volcanic activities (e.g., Tatard;2010).

Text S4. Note on the inversion method and uncertainties

Following Mosegaard and Tarantola (1995), we can use inverse problem method to estimate what is the epicentral location that best explain the avalanche rate derived from orbital imagery.

As the "direct problem" consisting in predicting the d-values that we should observe when we make observations on a given system, such system can be described by a vector m of parameters, such as m = {m 1 ,m 2 ,. . . ,m i }. We can define a probability density σ(m) that combines the a priori ρ(m) information with the information provided by the experimental measurements (the data vector) and the information provided by the physics of the problem (i.e., the model). Hence, σ(m) is called the posterior information. It is the combination of all the information we have and reads:

σ(m) = kρ(m)L(m), (7)
with k being a normalization constant (=1) and L(m) the likelihood function accounting for a Laplacian experimental uncertainties, which reads:

L(m) = exp - i ∥G i (m) -d i obs ∥ σ i , (8)
where G i (m) is the prediction by the model G(m) of the i th observation, d i obs is the i th observation with its associated uncertainties σ i . In this study, the model G(m) is provided after Livio and Ferrario (2020) which relates the avalanche (or landslides in their paper) distribution with the epicentral distance ∆:

G(m) = N ava = a exp - ∆ -b c 2 , ( 9
)
where N ava is the avalanches number triggered by the seismic event, where m = {a, b, c}, with a the peak amplitude of the distribution, b the distribution peak distance to the epicenter, and c the distribution width. Hence, we randomly explore the longitude-latitude space (Lon ∈ [166,175] and Lat ∈ [-10,-1]) over 1 million runs. For each run, compute the epicentral distance for each avalanche's area that meet both conditions (∆t¡1.5 MYear and Γ * < 450). Then, we fit the equation 9 on the observed N ava vs ∆ (expressed in euclidean distances and accounting for Mars flattening) using a non-linear least squares method. Finally, we compute the maximum likelihood function from equation 8 which accounts for a Laplacian distribution of uncertainties on the avalanche number.

The uncertainties on the avalanche rate q is computed after Aharonson et al. (2003), which reads:

σ = ⟨q -⟨q⟩⟩ 1/2 2 = ⟨q⟩ 1 + i ∆n i , (10)
11 with ⟨q⟩ being the expectation value of q, which reads:

⟨q⟩ = 1 + ∆n i i n i ∆t i . (11)

Figure 2 .Figure 3 .Figure 4 .Figure 5 .

 2345 Figure 2. (top) Pre-event CTX mosaic around the impact location of S1000a (dated from 2018-09-12) showing absence of any dust avalanches. (bottom) Post-event HiRISE mosaic on top of CTX images around the impact location of S1000a event (red star) with associated triggered avalanches (orange areas). Insets show close-up on the crater, the avalanches areas and slopes without new avalanches (from top-left, to bottomright, respectively). The density distribution of avalanches with respect to the epicentral distance is shown in the bottom-left inset.

Figure S1 :

 S1 Figure S1: Location of the pre-event avalanche mapping once an estimate of the S1222a location was provided. (left panel) close up on CTX G10 021978 1748 XI 05S192W (2011-04-05).

Figure S2 :

 S2 Figure S2: Slope map derived from MOLA 128ppd (Smith et al.; 2001). Orange symbols show location of avalanches over the 2005-2021 period.

Figure S4 :Figure S5 :

 S4S5 Figure S4: Estimation of the required additional dust heights for reaching slope failure under various conditions. The bounds for each computation are corresponding to an initial dust thickness ranging from 10 cm (lower bounds) up to 2 meters (upper bounds) and a slope angle of 30 • .

Figure S7 :Figure S8 :

 S7S8 Figure S7: Surface wind speed predicted from GCM at different Ls. Plots generated from the Mars Climate Database (c) LMD/LATMOS/OU/IAA/ESA/CNES.

  • , will have a ballistic flight time t 𝑓 of 76 sec (i.e., t 𝑓 = 2 × 𝑣 sin 𝜃/𝑔), and will land at a distance d 𝑙 = 10.78 km (neglecting the air friction, d 𝑙 = 𝑣 cos 𝜃 ×𝑡 𝑓 ). Hence, the histogram in the inset of figure 2 is similar to the statis-

	tical distribution of secondary ejecta impacting the ground. This correlation indicates those
	secondary impacts are a likely source for the avalanches. Of course, the S1000a event is an
	ideal case. First of all, we know the position of the epicenter perfectly well, thanks to the or-
	bital imagery revealing the source crater. What's more, the presence of northeast-southwest
	trending ripples implies the presence of uniformly distributed topographic slopes as mov-
	ing away from the impact crater, hence the avalanche susceptibility. Note that Burleigh et
	al. (2012) demonstrated that impact blast can trigger slope streaks. The S1000a event also

  Define the observed test statistic: Δ𝑛 𝑜𝑏𝑠 (change in the number of avalanches), 𝑛 𝑜𝑏𝑠 (number of avalanches), and Δ𝑡 𝑜𝑏𝑠 (time difference between images);Combine the before and after data into a single pool, disregarding their original labels;Perform resampling with replacement: Randomly sample, with replacement, from the pooled data to create a bootstrap sample of the same size as the original data set. Repeat this process to generate a large number of bootstrap samples.Calculate the test statistic for each bootstrap sample: Compute the rate of avalanches for each bootstrap sample, given by Δ𝑛 𝑏𝑜𝑜𝑡 𝑛 𝑏𝑜𝑜𝑡 Δ𝑡 𝑏𝑜𝑜𝑡 .; Calculate the bootstrap statistic distribution: Collect the calculated test statistics from step 4 to form the bootstrap distribution of the test statistic;Calculate the confidence interval: Determine the desired confidence level (e.g., 95%).Compute the lower and upper percentiles of the bootstrap distribution corresponding to the chosen confidence level;

indicates how solar energy absorption and subsequent subsurface heat propagation relate to material properties. Thermal inertia is represented as Γ ≡ √ 𝜅 𝑒 (1 -𝑝) 𝜌𝐶 (𝑇), where 𝜅 𝑒 denotes effective thermal conductivity, 𝑝 is porosity, 𝜌 represents density, and 𝐶 (𝑇) is the specific heat capacity. Therefore, low thermal inertia can indicate high porosity, low den-Algorithm 1: Assessing Influence of marsquake on avalanche rate Input: Pre-event image pairs (𝑛, Δ𝑛, Δ𝑡) Input: Post-event image pairs (𝑛, Δ𝑛, Δ𝑡) Output: Assessment of marsquake influence sity, small grain size, or a combination of these factors. As Sullivan et al. (2001); Aharonson et al. (2003) previously demonstrated, dust avalanches on Mars typically occur on steep slopes and are found in areas with low absolute thermal inertia. It's important to note that thermal inertia values are derived from models and assumptions, as detailed by Christensen

  thank Menina S., Margerin L., Kim D., Malystskyy D., Stähler S., Wieczorek M., Panning M. for InSightful discussions on the S1222a event. Authors thank Dundas C.M. for his help in accessing the post-event HiRISE images for both events. The authors express their gratitude to Ferrari C. for providing valuable insights into thermal inertia, to A. Spiga for his expertise on the output of General Circulation Models (GCMs), and to P. Souloumiac

	3°S
		for her expert knowledge on slope failure conditions. Authors thank anonymous review-
		ers and the associate editor for their constructive feedback that substantially improved the
		manuscript. All authors declare no conflict of interest. French co-authors acknowledge the
		French Space Agency CNES and ANR (ANR-19-CE31-0008). AL, TK, PL, SR, GS, AM
		acknowledge Idex Paris Cité (ANR-18-IDEX-0001). IJD was funded by NASA InSight PSP
		grant 80NSSC20K0971.
	5°S	6 Open Research The orbital data are available online: HRSC are available at ESA's Planetary Science
		Archive (https://www.cosmos.esa.int/web/psa/mars-express). THEMIS
		data are available at Arizona State University's repository (https://themis.asu.edu).
	7°S
	9°S
		168°E	170°E	172°E	174°E

MOC images are available at the PDS Imaging Node (https://pds-imaging.jpl .nasa.gov/data/mgs-m-moc-na_wa-2-sdp-l0-v1.0/). MOLA data are available at the PDS Geosciences Node (https://pds-geosciences.wustl.edu/missions/ mgs/mola.html). Apparent thermal inertia map is provided by the USGS (https:// astrogeology.usgs.gov/maps/mars-themis-derived-global-thermal -inertia-mosaic). HiRISE data, including the post-event images, are available at the University of Arizona's dedicated website (https://www.uahirise.org). CTX image are available at the Imaging PDS Node (https://pds-imaging.jpl.nasa.gov/ data/mro/mars_reconnaissance_orbiter/ctx/). The post-event CTX images will be posted on the NASA PDS by MSSS by the time of publication. Meanwhile, referee's can have access to the mosaic at https://www.dropbox.com/sh/u1cykaotwxvi7ga/ AAAsDcqw4FrkGDqjb4HTFmjka?dl=0. The avalanche catalogue is available on Zenodo (doi:10.5281/zenodo.7679315). The InSight seismic event catalogue version 9 (InSight Marsquake Service, 2022) and waveform data

(InSight Mars SEIS Data Service, 2019a,b)

are available from the IPGP Datacenter and IRIS-DMC, as are previous catalogue versions.

Seismic waveforms are also available from NASA PDS. The crustal thickness grid is available on Zenodo (doi:10.5281/zenodo.6477509).

Supporting Information for "Discussion of possible seismically triggered avalanches after the S1222a Marsquake and S1000a impact event" Contents of this file

Text S1. Note on post-marsquake image request and mapping

Once the S1222a event were recorded and a estimate of the location was provided, we investigate the orbital image archive for checking if either this location were falling into a region showing specific surface features (tectonics or avalanche signatures). It appears that this fairly flat region showing a few impact craters, into which, dust avalanches were detectable (Fig. S1). We therefore initiate the mapping using CTX and HiRISE images in order to established the list of location of interest. Then we request to MRO's team to target those specific places. As we focus on maximising the temporal coverage, we utilized all released images from MOC, THEMIS-Vis, CTX and HiRISE sensors. This leads to various conditions of observation (i.e., emission angle, local time etc.) and hence large parallax effects, especially on steep slopes areas where avalanches occur (i.e., typically from 30 to 8 degrees).

Consequently, the accuracy of the ortho-rectification performed with ISIS using either MOLA (Smith et al.;2001), and/or HRSC DTMs (Neukum and Jaumann;[START_REF] Miyamoto | Fluid dynamical implications of anastomosing slope streaks on Mars[END_REF], is not sufficient for accurate sub-pixel co-registration between all the images, even with bundle adjustment. This is an important aspect because as a result no automatic method could be used successfully (i.e., all of our principal component analysis (PCA) and Convolutional Neural Network (CNN) attempts failed to reach the required accuracy that could be achieved with human inspection).

We manually map all avalanches in the region shown in Fig. 1 in a Geographical Information System (GIS). Then each pair of overlapping images (mostly CTX/CTX pairs, but we also considered CTX/HiRISE and HiRISE/HiRISE overlapping observations) leads to a detection of new (or absence of new) avalanches. This mapping work has been conducted by two independent people to compare and validate the results. Although non-MRO images offer a greater historical depth (∼15 Earth years), the resolutions are very different between these sensors. Hence, over the 2008-2022 period, only CTX and HiRISE are considered. We down-sample the HiRISE data to 6 m/pixel, a pixel scale close to CTX's. The thermal inertia (Γ) is a complex data set that is not directly obtained from observations (see https: //astrogeology.usgs.gov/maps/mars-themis-derived-global-thermal-inertia-mosai Numerous factors contribute to thermal inertia. In our study, however, the emphasis is on the consistency of observations. We have conducted thorough verifications to ensure homogeneity in the slopes within the areas where avalanches were mapped. This implies that any topographic effects would uniformly impact these areas. By using flat terrain as a reference for comparison, the validity of our comparative analysis of the Γ * = Γ avalanche /Γ plain across the avalanche-affected areas is maintained. The map of the Γ is given in Figure S3.

Figure S3: Apparent thermal inertial map derived from THEMIS-IR over the region of interest provided by the USGS. Green star is the epicentral location estimated from Kawamura et al. (2023). Inset shows a close-up on a crater presenting dust avalanches with thermal inertial and its visible counterpart from CTX. 4