Cellulose nanocrystals (CNCs) are one of nature's most abundant structural material building blocks and possess outstanding mechanical properties such as a tensile modulus comparable to Kevlar. It remains challenging to upscale these properties in CNC neat films and nanocomposites due to the difficulty of characterizing interfacial bonding between CNCs that govern stress transfer under deformation. Here we present new analyses based on atomistic simulations of shear and tensile failure of the interfaces between Iβ CNCs in elementary fibrils, providing new insight into factors governing the mechanical behavior of hierarchical nanocellulose materials. We compare the two most relevant crystal surfaces and find that hydrogen bonded surfaces have greater tensile strength compared to the surfaces governed by weaker interactions. On the contrary, shearing simulations show that friction between the atomic interfaces depends not only on surface energy but also the energy landscape along the shear direction. While being a weaker interface, the intersheet plane exhibits greater energy barriers to shear. The molecular roughness of this interface, characterized by a greater energy barrier, exhibits a stick-slip deformation behavior as opposed to a continuous sliding mechanism observed for the interfaces with hydrogen bonds. Analytical models to describe the energy landscapes are developed using energy scaling relations for van der Waals surfaces in combination with a modification of the Prandtl-Tomlinson model for atomic friction. Our simulations pave the way for tailoring hierarchical CNC materials by taking a similar approach to techniques employed for describing metals, where mechanical properties can be tuned through a deeper understanding of grain boundary physics and nanoscale interfaces. 3 Graphical Abstract: 4